<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jan 23, 2014 at 11:05 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:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hello,<br>
         I am seeing inconsistent results using get stratum call.  For a 2d mesh<br></blockquote><div><br></div><div>Sorry, I was swamped last week. You had not interpolated your mesh. I have included modified source:</div>
<div><br></div><div><div>Successfully created DM From Cell List</div><div> Successfully interpolated mesh</div><div>DM_0x84000000_0 in 2 dimensions:</div><div>  0-cells: 9</div><div>  1-cells: 12</div><div>  2-cells: 4</div>
<div>Labels:</div><div>  depth: 3 strata of sizes (9, 12, 4)</div><div> stratum size for value           1 is          12</div><div>Number of indices in set 12</div><div>0 13</div><div>1 14</div><div>2 15</div><div>3 16</div>
<div>4 17</div><div>5 18</div><div>6 19</div><div>7 20</div><div>8 21</div><div>9 22</div><div>10 23</div><div>11 24</div><div> Node id start=           4 node id end=          12</div><div> edge id start=          13 edge id end=          24</div>
<div> element id start=           0 element id end=           3</div><div> printing the cell ids based on height stratum</div><div> element id start=           0 element id end=           3</div></div><div><br></div><div>
<br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
call DMPlexGetDepthStratum(dmMesh,2,globalidStart,globalidEnd,ierr)<br>
<br>
prints: globalIdStart = 0 and globalIdEnd = 0<br>
<br>
<br>
where as<br>
<br>
call DMPlexGetHeightStratum(dmMesh,0,globalidStart,globalidEnd,ierr)<br>
<br>
prints globalIdStart = 0 and globalIdEnd = 4<br>
<br>
For a mesh with 4 cells. I have attached the test case. I am on petsc<br>
version v3.4.3-2568-gb45093d<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>