<div class="gmail_extra">On Sat, Nov 10, 2012 at 2:28 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">> point that is a member. What about if a given level happens to be empty?<br>
<br>
</div>DMComplexGetLabelIdIS(). Its not problem if a level is empty.<br></blockquote><div><br></div><div>Ah, I misinterpreted what that routine did. Aren't these trivially short arrays (max length of 4 for meshes in 3D) that are expected to be identical on every process?</div>
<div><br></div><div>How does one get an IS containing _only_ the labeled points?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"></div>
<div class="im">
> Can you answer how you can distinguish a quad from a tet in a<br>
> non-interpolated mesh?<br>
<br>
</div>You can't, but that is the point. You are not supposed to distinguish them.</blockquote></div><br></div><div class="gmail_extra">Hmm, so these are the same at the category-theory level, but definitely not at the mesh topology level. It doesn't seem right to not distinguish them at all. How do you write DMComplexInterpolate for mixed-dim strata?</div>