I have imported a 3D mesh from gmsh with different physical volumes inside it in freefem, all are imported in one mesh.
I need to define different epsilon on each region. I can get my physical volumes tages by command
fespace Ph(Th, P1); Ph reg = region; plot(reg); real[int] tags; Unique(reg, tags); cout << "tags"<<tags << endl;
May you suggest me a command so to define epsilons for each region based on physical volumes tags?
Thanks in advance.