It is because you have not set ids for all boundary faces. Add this and try
Physical Curve("others", 3) = {2,4,5,6};
best
praveen
> On 15 Oct 2024, at 2:43 PM, Matteo Malvestiti <[email protected]> wrote:
>
> [Errata corrige]
> Hello everyone.
> I'm getting an error when importing msh files in deal.II with
> GridIn::read_msh().
>
> For the sake of simplicity I brought a very basic grid:
> - Rectangle 2x1
> - Divided in two squares 1x1, both uniformly meshed, transfinite, with
> different granularity
>
> When I 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 experience 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 defined as physical curves only two edges at the boundary.
>
> Can anybody help me to fix the issue please?
> I attach the simple script to reproduce the error, together with the .geo and
> .msh files, as well as the CMakeLists to compile the 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]
> <mailto:[email protected]>.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/69854a0c-5be6-40db-8919-48c736d64baen%40googlegroups.com
>
> <https://groups.google.com/d/msgid/dealii/69854a0c-5be6-40db-8919-48c736d64baen%40googlegroups.com?utm_medium=email&utm_source=footer>.
> <Ticket-error-read_msh.zip>
--
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/4E845272-A1F0-4C48-AF80-7A185BF6B228%40gmail.com.
