<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi everyone,<br>
<br>
I'm experiencing a problem, again, with the GMSH reader.<br>
<br>
You will find attached two .msh files.<br>
<br>
I managed to import a .msh file for a boundaryless surface (representing a sphere, an ellipsoid, ...) into ALUSimplexGrid<2, 3>, eliminating from the .msh all the elements that are not triangles (otherwise the segments would be interpreted as boundary objects,
but here there is no boundary at all). Thanks to Andreas Dedner for the useful suggestion.<br>
<br>
Anyway, if I try to use the same file with AlbertaGrid<2, 3>, this error message appears on my screen:<br>
<br>
Loading macro grid: ../data/ellipsoid.msh<br>
Reading 2d Gmsh grid...<br>
version 2.2 Gmsh file detected<br>
file contains 92 nodes<br>
file contains 178 elements<br>
number of real vertices = 91<br>
number of boundary elements = 0<br>
number of elements = 178<br>
Error: GridError [setOrientation:macrodata.cc:259]: Surface grid cannot be oriented.<br>
<br>
It seems that the same file gives problems with the orientation of the elements, but only in the AlbertaGrid case.<br>
Here my questions:<br>
<br>
1) are there any specific requirements about the orientation I should pay attention to? Check also the (long, sorry!) post scriptum.<br>
2) why doesn't ALUGrid detect the problem?<br>
<br>
Cheers,<br>
<br>
Andrea<br>
<div>________________________<br>
<div style="font-family:Tahoma; font-size:13px"><br>
PS, if I try with another file, the problem doesn't appear. For example:<br>
<br>
Loading macro grid: ../data/mesh_GMSH_coarse_new.msh<br>
Reading 2d Gmsh grid...<br>
version 2.2 Gmsh file detected<br>
file contains 19 nodes<br>
file contains 32 elements<br>
number of real vertices = 18<br>
number of boundary elements = 0<br>
number of elements = 32<br>
macro_test_2d: WARNING in ./../2d/macro_2d.c, line 451<br>
There is a cycle beginning in macro element 14.<br>
Correcting refinement edges....<br>
reorder_2d: Current elem: 0, best_edges: 0 3 3<br>
Current elem: 1, best_edges: 2 3 3<br>
Current elem: 2, best_edges: 1 3 3<br>
Current elem: 3, best_edges: 2 3 3<br>
Current elem: 4, best_edges: 0 3 3<br>
Current elem: 5, best_edges: 2 3 3<br>
Current elem: 6, best_edges: 0 1 2<br>
Looking at best_edges[6][0]...<br>
Found compatibly divisible neighbour 7!<br>
Current elem: 8, best_edges: 0 3 3<br>
Current elem: 9, best_edges: 1 3 3<br>
Current elem: 10, best_edges: 1 3 3<br>
Current elem: 11, best_edges: 2 3 3<br>
Current elem: 12, best_edges: 0 3 3<br>
Current elem: 13, best_edges: 1 3 3<br>
Current elem: 14, best_edges: 1 3 3<br>
Current elem: 15, best_edges: 2 3 3<br>
Current elem: 16, best_edges: 0 3 3<br>
Current elem: 17, best_edges: 0 1 2<br>
Looking at best_edges[17][0]...<br>
Found compatibly divisible neighbour 19!<br>
Current elem: 18, best_edges: 1 3 3<br>
Current elem: 20, best_edges: 0 3 3<br>
Current elem: 21, best_edges: 1 3 3<br>
Current elem: 22, best_edges: 1 3 3<br>
Current elem: 23, best_edges: 2 3 3<br>
Current elem: 24, best_edges: 0 3 3<br>
Current elem: 25, best_edges: 2 3 3<br>
Current elem: 26, best_edges: 2 3 3<br>
Current elem: 27, best_edges: 1 3 3<br>
Current elem: 28, best_edges: 0 3 3<br>
Current elem: 29, best_edges: 1 3 3<br>
Current elem: 30, best_edges: 1 3 3<br>
Current elem: 31, best_edges: 2 3 3<br>
check_mesh: checking mesh 'DUNE AlbertaGrid'<br>
checking done; no error detected<br>
Grid file read correctly.<br>
<br>
__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial College London,<br>
London SW7 2AZ, UK<br>
a.sacconi11@imperial.ac.uk<br>
</div>
</div>
</div>
</div>
</body>
</html>