[petsc-dev] PetscSection

Jed Brown jedbrown at mcs.anl.gov
Sat Nov 10 14:29:09 CST 2012

On Sat, Nov 10, 2012 at 2:09 PM, Blaise A Bourdin <bourdin at lsu.edu> wrote:

> Nope... and please, don't break that...
>  Here is an example of a mesh with 4 element blocks. In an analysis code,
> for each block , I would specify the element type, type of physics,
> coefficients etc in a petscBag. That said, I may be doing it wrong...
>  block 10 is made of tri
> block 22 is made of quad
> blocks 100 and 102 are made of bars

I'm not suggesting any demotion of this ability. I'm trying to understand
what is the advantage of having a stratum with unsorted mixed dimension. If
the blocks were sorted by dimension, you would still access them using a
label (which gives you an index set, perhaps contiguous). You would still
be able to do all the normal operations with those points. If you want to
have one loop over all elements of the various types, you would union or
concatenate the index sets of each type. Am I missing something?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20121110/4bc9dffd/attachment.html>

More information about the petsc-dev mailing list