Marek,
I have prepared slightly better mesh - without "doubling" in the intersection,
see attachment.
However I got error on reading in the mesh file.
n error occurred in line <1447> of file
</home/marek/SECONDARY/Libraries/numerics/candi/candi_reinst_install/tmp/unpack/deal.II-v8.4.1/source/grid/grid_in.cc>
in function
void dealii::GridIn<dim, spacedim>::read_msh(std::istream&) [with int dim
= 3; int spacedim = 3; std::istream = std::basic_istream<char>]
The violated condition was:
false
The name and call sequence of the exception was:
ExcInvalidVertexIndex(cell, subcelldata.boundary_lines.back().vertices[i])
Additional Information:
While creating cell 55, you are referencing a vertex with index 0 but no
vertex with this index has been described in the input file.
Do You have any advice, how to cope with this problem.
The error message is correct. If you look at the file, it starts with
$MeshFormat
2.2 0 8
$EndMeshFormat
$Nodes
9708
1 0 0 0
2 1 0 0
3 0 1 0
...
That is, it starts numbering vertices at 1 (first column). But then, later on,
it describes cells as follows:
$Elements
11703
1 15 2 0 1 1
2 15 2 0 2 2
3 15 2 0 3 3
...
56 1 2 0 5 0 52
...
Here, each line corresponds to an point, edge, face, or cell with properties.
Column 1 has the number of the object, column 2 the type of the object
(line=1, quadrangle=2, hexahedron=5, point=15). After the type is a number
that indicates how many properties there are for this object (here, 2
properties for points and lines), followed by this many property values,
followed by the vertex indices. Indeed, object 56 (a line, with two properties
with values 0 and 5) consists of vertices 0 and 52. But we don't have a vertex
with number 0, since they were numbered starting at 1.
I don't know how this mesh was created, but the content of the file indeed
does not seem correct to me. You probably need to talk to the people who wrote
gmsh what this is supposed to mean.
Best
Wolfgang
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
