I have been using mmgls as part of a FEM. Each time step I use the level set function provided by mmgls (5.5.2 - tried the develop branch too) but I find that sometimes I get the following error when calling mmgls using the API:
mmg3d/tools_3d.c:552: MMG5_BezierNom: Assertion `( p0->tag & MG_NOM ) && ( p1->tag & MG_NOM )’ failed.
Unexpected error: *** potential lack of memory.
I have provided two examples in attached zip file along with a .F90 script which tests the levelset function using the same parameters as set in the model. I often avoid the error when testing it like this by changing the Hausdorff distance but since I cannot do this when using the FEM, I was wondering if there was a more robust solution. Additionally, is it possible to handle such errors in mmg with an error output so an exception can be written to save the model falling over?
mmgexamples.zip (282.7 KB)
Thanks in advance for any help!