I’m quite new, i’d like to use freefem with gmsh file mesh.
LoadMesh function make easy to import a 2D gmsh mesh file and labels(Th) function give all physical boundary tags (dim=1).
1 - I’d like to know how to get all physical region tags (dim=2), is there a function return that list ?
2 - I want to solve the heat equation with different kappa for each region, how to define a piecewise constant value for each subdomain identified by their physical tag ? All examples tackles the case where subdomains are easily defined by the indicator domain function, but in my case each domain are complex and only the physical region tag can help to indicate the label region. Is there an instruction to define kappa as the sum of indicator region function lika kappa = 4on(region1) + 2on(region2) ?