<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Jan 27, 2014 at 1:58 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">On Mon, Jan 27, 2014 at 11:18 AM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>

> On Thu, Jan 23, 2014 at 11:05 PM, Dharmendar Reddy <<a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a>><br>
> wrote:<br>
>><br>
>> Hello,<br>
>>          I am seeing inconsistent results using get stratum call.  For a<br>
>> 2d mesh<br>
><br>
><br>
> Sorry, I was swamped last week. You had not interpolated your mesh. I have<br>
> included modified source:<br>
<br>
Agreed, But if the mesh is not interpolated, depth = topological<br>
dimension  (dim) should refer to dim-cells, is not that true ?<br></blockquote><div><br></div><div>No, depth is depth in the DAG, not topological dimension.</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">

><br>
> Successfully created DM From Cell List<br>
>  Successfully interpolated mesh<br>
> DM_0x84000000_0 in 2 dimensions:<br>
>   0-cells: 9<br>
>   1-cells: 12<br>
>   2-cells: 4<br>
> Labels:<br>
>   depth: 3 strata of sizes (9, 12, 4)<br>
>  stratum size for value           1 is          12<br>
> Number of indices in set 12<br>
> 0 13<br>
> 1 14<br>
> 2 15<br>
> 3 16<br>
> 4 17<br>
> 5 18<br>
> 6 19<br>
> 7 20<br>
> 8 21<br>
> 9 22<br>
> 10 23<br>
> 11 24<br>
>  Node id start=           4 node id end=          12<br>
>  edge id start=          13 edge id end=          24<br>
>  element id start=           0 element id end=           3<br>
>  printing the cell ids based on height stratum<br>
>  element id start=           0 element id end=           3<br>
><br>
><br>
>    Matt<br>
><br>
>><br>
>> 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>
><br>
><br>
><br>
<span class="HOEnZb"><font color="#888888">><br>
> --<br>
> What most experimenters take for granted before they begin their experiments<br>
> is infinitely more interesting than any results to which their experiments<br>
> lead.<br>
> -- Norbert Wiener<br>
</font></span></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>