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

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:

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


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?

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