I’m relatively new to MMG, but already managed to get some very good results re-meshing sub-domains with Level set functions. However, for my work the set up seems a bit complicated and I’m hoping to improve it.
My current set-up takes an input mesh, then (1) re-meshes it in MMG using a fixed size (hsiz) which is then sampled with an implicit function. (2) That mesh is then cut at an iso-value (I do this in VTK as I need to preserve boundaries and labels (text strings) which works better within the current setup). (3) The result from (2) is then optimized in MMG to obtain the final result with all sub-domains preserved correctly. I repeat the process for multiple sub-domains and the results are really good.
However, in my field (geology), area extents can be 100’s of meter or even kilometres, whereas the domain to model can be 10’s of centimetres up to a few metres. This means step (1) needs to sample at a very fine resolution across the whole space to be able to construct the isosurface correctly.
My question is if in step (1) I can locally insert only new vertices into the existing mesh close to the zero crossing of the level set function instead of using the -hsiz option? This could save a lot or memory, and potentially processing time.
BTW, I’m using the library so advice that shows how to do it using the API are much appreciated.
Many thanks in advance,