<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jan 10, 2014 at 3:44 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
        I have a question and feature request regarding DMPlexCreateGmsh.<br>
<br>
Gmsh, mesh format outputs, physical regions in the mesh as<br>
For example in the a particular 2 d mesh with eight physical regions<br>
it may look like the one below<br>
$PhysicalNames<br>
8<br>
1 1 "RegionName1"<br>
1 2 "RegionName2"<br>
1 3 "RegionName3"<br>
2 4 "RegionName4"<br>
2 5 "RegionName5"<br>
2 6 "RegionName6"<br>
2 7 "RegionName7"<br>
2 8 "RegionName8"<br>
$EndPhysicalNames<br>
<br>
It would be nice to have the DMPLexCreateGmsh process this information.<br>
<br>
The first number in each row of name is the topological dimension (tdim),<br>
<br>
If the code can group all the regions with dimension < mesh dimension<br>
(meshDim) as "Boundary"  and all regions with dimension =<br>
meshDimension as "Region"<br>
<br>
Gmsh outputs (if requested) the lower dimension cells in regions with<br>
tdim  < meshDim. If interpolate mesh is enabled, i would like to have<br>
cells of boundary region, which in the above example will be lines and<br>
points, included in the boundary strata with global numbering for the<br>
facets.<br></blockquote><div><br></div><div>I really hate GMsh right now. The format is a complete cluster of idiocy. They allow</div><div>outputting facets as separate lower-dimensional cells, unconnected to the original</div>
<div>cell, so we would have to initiate a search for every insertion. Not to mention that sizes</div><div>are not declared up front so that we have to run through the file in multiple passes. I</div><div>really really hate these guys. Is there an overwhelming reason to waste my time messing</div>
<div>around with a format designed by the mentally infirm?</div><div><br></div><div>I am willing to listen to reason, but reading Gmsh files would make anyone unreasonably angry.</div><div><br></div><div>    Matt</div><div>
 </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
This ways the total number of Cells in the mesh will be : NumNodes +<br>
NumLines+ NumTriangles or NumNodes+ NumTriangles (if interpolate is<br>
off)<br>
<br>
<br>
In the Current implementation, correct me if i am wrong, I see the the<br>
total number of cells = numNodes + (some Lines which are elements of<br>
the lower dimensional regions) + numTriangles<br>
<br>
<br>
Thanks<br>
Reddy<br>
<br>
<br>
--<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>