I tried to use Implicit domain meshing MMG3D to get an adaptive mesh based on 0 value of the level-set function
After reading the mesh into FreeFem++, I need to identify the boundary label in order to apply the boundary condition. Could you please give me some suggestion about this?
Thanks a lot!
Here is one way to do it.
mesh Th = square(10, 10);
int[int] lab = labels(Th);
varf onG(u, v) = on(lab, u = label);
fespace Vh(Th, P1);
u = onG(0, Vh, tgv = -1);
Thanks so much for you kind help.
It does work!
I tried your method for 2D and 3D problem, it works fine!