dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #04019
Re: dolfin-convert: bug with conversion from gmsh?
Looks ok now! Thanks.
/Anders
On Thu, Dec 07, 2006 at 03:03:34PM +0100, Angelo Simone wrote:
> I have fixed the gmsh2xml conversion.
>
> angelo
>
>
> >I've fixed the xml_old2xml conversion also, so only the gmsh2xml
> >conversion remains. I'll leave it as it is for now and wait for
> >another patch (bundle) from you.
> >
> >Make sure you make the bundle against the current dolfin tree (which
> >include my updates to mesh2xml and xml_old2xml).
> >
> >Thanks
> >/Anders
> >
> >
> >On Thu, Dec 07, 2006 at 12:03:48PM +0100, Anders Logg wrote:
> >
> >>On Thu, Dec 07, 2006 at 01:18:49AM +0100, angelo simone wrote:
> >>
> >>>Here's the bundle. dolfin-convert is now in pure python (i.e. there are
> >>>no calls to grep) and there are no temporary files. Is it ok now?
> >>>The script has been tested with the attached mesh files:
> >>>gmsh.msh (2d)
> >>>old-xml.xml.gz (2d)
> >>>tetgen.mesh.gz (3d)
> >>>
> >>Thanks, but I don't really like the write_output stuff. I have added
> >>your changes and modified mesh2xml so it first scans the file for the
> >>cell type and then does the conversion. Could you supply a patch that
> >>does the same thing for gmsh2xml and xml_old2xml?
> >>
> >>
> >>>I have tried the output of the first two conversions with a Dolfin
> >>>example and they worked out fine. The tetgen mesh, which comes from the
> >>>example.poly file distributed with tetgen 1.4.1, has not been tested
> >>>against Dolfin. The output does look ok though.
> >>>
> >>Thanks. I have added the test meshes and your updated conversion
> >>script, but I don't like the write_output stuff
> >>
> >>
> >>>As far as I can see from the tetgen web page
> >>>(http://tetgen.berlios.de/), tetgen cannot generate 2d meshes. Is there
> >>>a reason why in the script there is a check on num_dims == 2?
> >>>
> >>I haven't tested, but presumably the conversion works for meshes
> >>generated by Triangle. I think Triangle and TetGen share the same format.
> >>
> >>
> >>>I will soon send a bundle to update section 8.3 of the
> >>>dolfin-user-manual. A bundle on the procedure to create a bundle will be
> >>>ready soon as well.
> >>>
> >>Very nice.
> >>
> >>/Anders
> >>
> >>
> >>
> >>>angelo
> >>>
> >>>
> >>>
> >>>
> >>>Anders Logg wrote:
> >>>
> >>>>You could try generating some with tetgen.
> >>>>
> >>>>When you submit the patch, please also include some (very small) test
> >>>>meshes that you have run the code on. We can add these for future
> >>>>reference.
> >>>>
> >>>>Thanks
> >>>>/Anders
> >>>>
> >>>>
> >>>>On Tue, Dec 05, 2006 at 01:56:00AM +0100, angelo simone wrote:
> >>>>
> >>>>
> >>>>>Exactly. Anyway, I got rid of the ugly hack... ;-)
> >>>>>
> >>>>>Before sending a bundle, I would like to test the new version on 2d
> >>>>>and 3d Medit meshes, which I do not have. Does anyone have them?
> >>>>>Thanks!
> >>>>>
> >>>>>angelo
> >>>>>
> >>>>>
> >>>>>Anders Logg wrote:
> >>>>>
> >>>>>
> >>>>>>I think because the new mesh format requires celltype in the header of
> >>>>>>mesh, so the mesh needs to be scanned before writing to find out the
> >>>>>>celltype.
> >>>>>>
> >>>>>>Instead of doing this, the patch first generates an illegal xml file
> >>>>>>missing the celltype attribute and then scans the generated file to
> >>>>>>grep for the cell type and then fixes the header...
> >>>>>>
> >>>>>>It works but it's not a nice solution. Angelo, please submit an
> >>>>>>updated patch that first scans the file for the cell type (preferably
> >>>>>>not using grep).
> >>>>>>
> >>>>>>/Anders
> >>>>>>
> >>>>>>PS: Yes, I know you looked at xml_old2xml() that I wrote and I do use
> >>>>>>grep, but it's still an ugly hack. :-)
> >>>>>>
> >>>>>>
> >>>>>>On Mon, Dec 04, 2006 at 11:06:21AM +0100, Garth N. Wells wrote:
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>>Why are temporary files being created? Isn't updating to the new
> >>>>>>>mesh xml format just a case of changing the output formatting?
> >>>>>>>
> >>>>>>>Garth
> >>>>>>>
> >>>>>>>angelo simone wrote:
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>>hi,
> >>>>>>>>
> >>>>>>>>I have fixed dolfin-convert (patch attached). Now the output is in
> >>>>>>>>the current format for both Gmsh and Medit. I have tested
> >>>>>>>>dolfin-convert with triangular Gmsh meshes but is should work fine
> >>>>>>>>with all the remaining combinations as well.
> >>>>>>>>
> >>>>>>>>Regarding the addition to the user manual, maybe we should just add
> >>>>>>>>that on top of the steps described in the Gmsh mesh generation
> >>>>>>>>tutorial
> >>>>>>>>(http://www.geuz.org/gmsh/doc/gui_tutorial/gui_tutorial.pdf) one
> >>>>>>>>needs to create physical surfaces/volumes as well to generate a
> >>>>>>>>dolfin-friendly mesh -- yesterday I have used Gmsh for the first
> >>>>>>>>time... that's why I did that mistake. Anyway, dolfin-convert
> >>>>>>>>complains if the mesh is not in the correct format.
> >>>>>>>>
> >>>>>>>>angelo
> >>>>>>>>
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
Follow ups
References