[petsc-dev] PetscSection

[Jed Brown]

On Sat, Nov 10, 2012 at 2:28 PM, Matthew Knepley <knepley at gmail.com> wrote:

> > point that is a member. What about if a given level happens to be empty?
>
> DMComplexGetLabelIdIS(). Its not problem if a level is empty.
>

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?

How does one get an IS containing _only_ the labeled points?


> > Can you answer how you can distinguish a quad from a tet in a
> > non-interpolated mesh?
>
> You can't, but that is the point. You are not supposed to distinguish them.


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?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121110/21e5e09f/attachment.html>