dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

Thanks, this should be easy to fix. I'll fix it after breakfast...

/Anders


On Thu, Mar 30, 2006 at 01:08:59PM +0000, Alexander Jarosch wrote:
> Hi everybody,
> 
> sorry to bring up the subject again but I still have problems with gmsh 
> meshes and dolfin-convert.
> 
> a really simple geo file is
> 
> ----------------------
> Point(1) = {0,0,0,1};
> Point(2) = {1,0,0,1};
> Point(3) = {1,1,0,1};
> Point(4) = {0,1,0,1};
> Point(5) = {0,0,1,1};
> Point(6) = {1,0,1,1};
> Point(7) = {1,1,1,1};
> Point(8) = {0,1,1,1};
> Line(1) = {1,2};
> Line(2) = {2,3};
> Line(3) = {4,3};
> Line(4) = {4,1};
> Line(5) = {5,6};
> Line(6) = {6,7};
> Line(7) = {7,8};
> Line(8) = {8,5};
> Line(9) = {5,1};
> Line(10) = {4,8};
> Line(11) = {6,2};
> Line(12) = {3,7};
> Line Loop(13) = {1,-11,-5,9};
> Plane Surface(14) = {13};
> Line Loop(15) = {11,2,12,-6};
> Plane Surface(16) = {15};
> Line Loop(17) = {2,-3,4,1};
> Plane Surface(18) = {17};
> Line Loop(19) = {5,6,7,8};
> Plane Surface(20) = {19};
> Line Loop(21) = {7,-10,3,12};
> Plane Surface(22) = {21};
> Line Loop(23) = {10,8,9,-4};
> Plane Surface(24) = {23};
> Surface Loop(25) = {22,-20,-14,18,-16,24};
> Volume(26) = {25};
> Physical Volume(27) = {26};
> -----------
> 
> if I mesh it with
> gmsh test.geo -3 -o test.msh
> 
> i get this mesh file
> 
> -----------------------
> $NOD
> 9
> 1 0 0 0
> 2 1 0 0
> 3 1 1 0
> 4 0 1 0
> 5 0 0 1
> 6 1 0 1
> 7 1 1 1
> 8 0 1 1
> 17 0.5 0.5 0.5
> $ENDNOD
> $ELM
> 12
> 1 4 27 26 4 1 17 2 4
> 2 4 27 26 4 17 1 2 5
> 3 4 27 26 4 1 17 4 5
> 4 4 27 26 4 2 17 3 4
> 5 4 27 26 4 17 2 3 6
> 6 4 27 26 4 2 17 5 6
> 7 4 27 26 4 17 3 4 8
> 8 4 27 26 4 3 17 6 7
> 9 4 27 26 4 3 17 7 8
> 10 4 27 26 4 17 4 5 8
> 11 4 27 26 4 17 5 6 8
> 12 4 27 26 4 17 6 7 8
> $ENDELM
> -------------------
> 
> and if I use dolfin-convert to convert it I get
> 
> Converting from Gmsh format (.msh, .gmsh) to DOLFIN XML format
> Expecting 9 vertices
> Found all vertices
> Expecting 12 cells
> *** Vertex 17 of tetrahedron 1 not previously defined.
> *** Make sure to only export the physical mesh from Gmsh.
> 
> I am pretty sure I only exported the physical mesh and if i look into 
> the mesh vertex 17 was defined...
> 
> so whats going wrong in that conversion? Am I missing something in the 
> use of gmsh?
> 
> thanks for the help
> 
> /Alex
> 

-- 
Anders Logg
Research Assistant Professor
Toyota Technological Institute at Chicago
http://www.tti-c.org/logg/