Hello mmg community,

First of all, thank you so much for all your efforts rendering the task of creating complex geometries soooo easy!!

I have a question regarding implicit meshing, that hopefully you can help me with. I am trying to implement 3D shape optimization for inverse homogenization using FreeFem++, By nature of the application, I have to use 3d periodic FE spaces in a cube with FreeFem++, and FreeFem’s constraint is that for us to be able to define periodic boundary conditions ont two given faces, these faces must be absolutely identical. My first attempt was to use a fixed cubic mesh, and do shape optimization with a level set to define my interface (inclusion of two crossing cylinders inside the cubic cell), but the results are not so accurate.

I would like to try and mesh my interface better using mesh evolution, and it seems the only obstacle is to make sure mmg will only remesh the interior of my domain, and leave the cube surfaces unchanged, or ideally, if we were to split two identical cube surfaces with respect to two identical inclusions (a disk), we get two identical remeshed surfaces.

Thank you in advance for any help or explanation you can point me to.

Sincerly your,

Salah