Hi, thanks for your time.
I implemented the suggested fix, but unfortunately it didn't change
anything.
Still get the 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).
Il giorno martedì 15 ottobre 2024 alle 16:27:00 UTC+2 Praveen C ha scritto:
> 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].
> 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/394e4d09-5d50-45ed-9022-07f9fd3546ben%40googlegroups.com.
