Segmentation error with fine mesh

General Discussion
1

Plane_5.zip (5.7 KB)

Hello all,
I have the following mesh. When I try to open it like this:

mesh Th1=gmshload("Plane_5.msh");

I get the following errror:

Assertion fail : (area > 0) line :303, in file ../seq/include/fem.hpp Assertion fail : (area > 0) line :303, in file ../seq/include/fem.hpp err code 6 , mpirank 0
This error happens only when the mesh is refined. Indeed, if I consider coarser elements, this problem does not occur.

Could someone help to understand what is the problem ?

Thank you in advance,

Best regards,

Loïc,

2

I attach also the same mesh with coarser elements.
Plane_5.zip (3.6 KB)

I have no problem to open this one.

Best regards,

Loïc,

3

Are you sure about your mesh?

load "PETSc"
DM dm("Plane_5.msh");
mesh Th = dm;
fespace Vh(Th, P1);
Vh u = x + y;
plot(u);

Screenshot 2022-04-13 at 5.09.40 PM
Screenshot 2022-04-13 at 5.09.40 PM3638×2202 820 KB

4

The second mesh looks better, no weird corner on the bottom left.

5

Hello @prj ,

Thank you for this vizualization.
I see my error now, the weird corner on the bottom left should not be there.

Best regards,

Loïc,