How to identify the boundary label for an import mesh

Dear developers,

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);
Vh u;
u[] = onG(0, Vh, tgv = -1);
plot(u);
3 Likes

Thanks so much for you kind help.
It does work!
I tried your method for 2D and 3D problem, it works fine!

Thanks agian.

Sincerely yours,
Hao