<div dir="ltr"><div><div><div><div><div><div><div><div><div><div>Hello,<br></div>I have merged two .geo geometries, see attachment.<br></div>I have meshed the result, I have got the mesh in attachment.<br></div>I compute in deal.ii fem library, it threw error <br><br><span><br></span><div style="margin-left:40px"><div style="margin-left:40px"><span>
n error occurred in line <1447> of file</span><br><span>
</home/marek/SECONDARY/Librari<wbr>es/numerics/candi/candi_reinst<wbr>_install/tmp/unpack/deal.II-<wbr>v8.4.1/source/grid/grid_in.cc></span><br><span>
in function</span><br><span>
void dealii::GridIn<dim, spacedim>::read_msh(std::istre<wbr>am&) [with int dim</span><br><span>
= 3; int spacedim = 3; std::istream = std::basic_istream<char>]</span><br><span>
The violated condition was:</span><br><span>
false</span><br><span>
The name and call sequence of the exception was:</span><br><span>
ExcInvalidVertexIndex(cell, subcelldata.boundary_lines.bac<wbr>k().vertices[i])</span><br><span>
Additional Information:</span><br><span>
While creating cell 55, you are referencing a vertex with index 0 but no</span><br><span>
vertex with this index has been described in the input file.</span><br></div><span></span><br></div><div style="margin-left:40px"><span></span></div><span><br></span></div><span>meaning, that the numbering of vertices starts with 1<br><br></span>$MeshFormat<br>
2.2 0 8<br>
$EndMeshFormat<br>
$Nodes<br>
9708<br>
1 0 0 0<br>
2 1 0 0<br>
3 0 1 0<br>
...<br>
<br>
That is, it starts numbering vertices at 1 (first column). But then, later on, it describes cells as follows:<br>
<br>
$Elements<br>
11703<br>
1 15 2 0 1 1<br>
2 15 2 0 2 2<br>
3 15 2 0 3 3<br>
...<br><b>
56 1 2 0 5 0 52</b><br><br></div>which means, that it references vertex with number 0<br></div>Maybe I have merged the ,geo files incorrectly.<br></div>Could You advice me, how to merge the files in the right way.<br></div>Or should I modify them?<br><br></div>Thank You<br><br></div>Marek<br></div>