Chen,
I want to import the mesh file exported by gmsh.
[...]
*Could anyone let me known how to import mesh ?*
First, our apologies for not replying to this issue earlier. These
interactions with other programs are really thorny :-(
I've taken a look at your input VTK file. The first issue is that our reader
in GridIn expects (maybe unreasonably so) that cells are listed before lines
in the CELLS section. This is relatively easy to fix by swapping the lines and
cells in the CELLS and CELL_TYPES section of your input file.
But the more relevant part is that one then runs into the following error:
...........................................
An error occurred in line <2690> of file
</home/bangerth/p/deal.II/1/dealii/source/grid/tria.cc> in function
static void
dealii::internal::TriangulationImplementation::Implementation::create_triangulation(const
std::vector<dealii::Point<dim> >&, const std::vector<dealii::CellData<2> >&,
const dealii::SubCellData&, dealii::Triangulation<2, spacedim>&) [with int
spacedim = 2]
The violated condition was:
line->boundary_id() != numbers::internal_face_boundary_id
Additional information:
The input data for creating a triangulation contained information about a
line with indices 1 and 49 that is described to have boundary indicator 0.
However, this is an internal line not located on the boundary. You cannot
assign a boundary indicator to it.
If this happened at a place where you call
Triangulation::create_triangulation() yourself, you need to check the
SubCellData object you pass to this function.
If this happened in a place where you are reading a mesh from a file, then you
need to investigate why such a line ended up in the input file. A typical case
is a geometry that consisted of multiple parts and for which the mesh
generator program assumes that the interface between two parts is a boundary
when that isn't supposed to be the case, or where the mesh generator simply
assigns 'geometry indicators' to lines at the perimeter of a part that are not
supposed to be interpreted as 'boundary indicators'.
...........................................
The description of the error message already suggests the problem. You need to
somehow convince gmsh to not output information about line segments at all.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/2b99a9fc-aaac-13fd-6da0-01dbfd212816%40colostate.edu.
