I’m having some troubles in the identifications of different regions into a mesh imported from Gmsh with the format .mesh with the line:
mesh3 Th = readmesh3(“Trial1.mesh”);
The loading of the mesh is perfect since I can plot and visualize it and is exactly what I’m looking for, but my mesh is formed by two volumes that are identified as two separate “physical groups” in gmsh, however if I try to access the label of the regions in FreeFEM with the line:
int r1 = Th(xc,yc,zc).region;
int r2 = Th(xs,ys,zs).region;
(where the (x,y,z) tripletused before are representative of two points belonging to different volumes).
what I actually get is that r1 = r2, differently from what I was expecting and from what I need unfortunately!
Is some of you able to tell me some ways in which I can solve this problem? Should I do something particoular during the mesh generation procedure?
If it can help the two volumes are in contact since they are representative of a cable and its jacket.
Thanks a lot to everyone will help me!