[Dune] GMSH reader fails in reading the same file for ALUSimplexGrid<2, 3> and not ALUSimplexGrid<3, 3>

[Sacconi, Andrea]

Hi everyone,

I have a problem with the GMSH reader.

I'm trying to build a new grid by reading a .msh file (attached); ideally, the computational mesh should be made only of triangles located in R^3, so it can be interpreted either as a 3D mesh with only boundary elements or as a 2D mesh embedded in R^3.

When I try the first solution, that is

#define ALUGRID_SIMPLEX
#define GRIDDIM 3
#define GRIDWORLD 3

everything works properly.

Instead, if I modify

#define ALUGRID_SIMPLEX
#define GRIDDIM 2
#define GRIDWORLD 3

I get this error:

Reading 2d Gmsh grid...
version 2.2 Gmsh file detected
file contains 304 nodes
file contains 693 elements
number of real vertices = 303
number of boundary elements = 84
number of elements = 602
terminate called after throwing an instance of 'Dune::GridError'
Abort trap

Does anyone know if there are problems when passing from < 3 , 3 > to < 2 , 3> ?
I've attached the file.

Thanks in advance!

Andrea
__________________________________________________________

Andrea Sacconi
PhD student, Applied Mathematics
AMMP Section, Department of Mathematics, Imperial College London,
London SW7 2AZ, UK
a.sacconi11 at imperial.ac.uk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20120606/224a8d7b/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sphere_surface.msh
Type: application/octet-stream
Size: 33156 bytes
Desc: sphere_surface.msh
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20120606/224a8d7b/attachment.obj>