Assertion fail:(ie==nbe) in msh2.cpp

General Discussion
1

Dear all,

I have a 3D mesh generated by GID, which has an internal surface as shown in the following figure.
Using one processor there is no problem of reading the mesh and executing my parallel code. However, the code cannot deal with the mesh correctly if I use more that one processor, and the error is:

image
image1122×295 11.8 KB

image
image873×684 74.7 KB

image
image862×591 59.6 KB

Does this error occur when executing buildDmesh(Th) ?
Does any one have an ideas what’s going wrong?

Best,
Yongxing.

2

Most probably your mesh is not a closed surface and there is a problem with the boundary elements.
If you just try to load the mesh and do Th = trunc(Th, 1, split = 2); in a sequential script, do you have an error?

3

Thank you very much! I did find a problem the surface mesh.
Yongxing.