[petsc-users] DMPlexCreateGmsh

[Dharmendar Reddy]

Hello,
        I have a question and feature request regarding DMPlexCreateGmsh.

Gmsh, mesh format outputs, physical regions in the mesh as
For example in the a particular 2 d mesh with eight physical regions
it may look like the one below
$PhysicalNames
8
1 1 "RegionName1"
1 2 "RegionName2"
1 3 "RegionName3"
2 4 "RegionName4"
2 5 "RegionName5"
2 6 "RegionName6"
2 7 "RegionName7"
2 8 "RegionName8"
$EndPhysicalNames

It would be nice to have the DMPLexCreateGmsh process this information.

The first number in each row of name is the topological dimension (tdim),

If the code can group all the regions with dimension < mesh dimension
(meshDim) as "Boundary"  and all regions with dimension =
meshDimension as "Region"

Gmsh outputs (if requested) the lower dimension cells in regions with
tdim  < meshDim. If interpolate mesh is enabled, i would like to have
cells of boundary region, which in the above example will be lines and
points, included in the boundary strata with global numbering for the
facets.

This ways the total number of Cells in the mesh will be : NumNodes +
NumLines+ NumTriangles or NumNodes+ NumTriangles (if interpolate is
off)


In the Current implementation, correct me if i am wrong, I see the the
total number of cells = numNodes + (some Lines which are elements of
the lower dimensional regions) + numTriangles


Thanks
Reddy


--