Hello everyone. I'm encountering this basic problem with GridIn::read_msh().
For the sake of simplicity I brought a very simple mesh, created with Gmsh and exported as .msh ASCII 4 file. - Rectangle 2x1 - Divided in two square regions 1x1, both meshed in a structured way, uniform spacing. When I try to import it in deal.II with: const std::string filename = "mesh.msh"; std::ifstream input_file(filename); Triangulation<2> triangulation; GridIn<2> grid_in; grid_in.attach_triangulation(triangulation); grid_in.read_msh(input_file); I get this error: terminate called after throwing an instance of 'dealii::StandardExceptions::ExcMessage' what(): -------------------------------------------------------- An error occurred in line <2594> of file </usr/src/dealii-v9.5.0/source/grid/tria.cc> in function static void dealii::internal::TriangulationImplementation::Implementation::process_subcelldata(const dealii::internal::TriangulationImplementation::CRS<T>&, dealii::internal::TriangulationImplementation::TriaObjects&, const std::vector<dealii::CellData<dim> >&, const std::vector<dealii::Point<spacedim> >&) [with int structdim = 1; int spacedim = 2; T = unsigned int] The violated condition was: boundary_id != numbers::internal_face_boundary_id Additional information: The input arguments for creating a triangulation specified a boundary id for an internal face. This is not allowed. The object in question has vertex indices 1,40, which are located at positions (1 0)(1 0.25). Of course I only defined as Physical Curves only two segments at the boundary. Can anybody help me to fix the issue? I attach a basic script to reproduce the error, together with the .geo and .msh files and the CMakeLists to compile the snippet of code. Thanks in advance. -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/1fb44be2-808e-4eee-a43f-4a22ab7c9e3fn%40googlegroups.com.
