[CIG-SHORT] Fw: cells connected to a fault vertex

[Matthew Knepley]

On Mon, Nov 16, 2015 at 10:05 AM, Birendra jha <bjha7333 at yahoo.com> wrote:

> Hi Matt
>
> I was wondering if you could help me with this. I tried to follow Sieve's
> documentation to get the indices of cells connected to a fault vertex using
> orientedStar but it is not working (please see below). Do you have any
> suggestions?
>

Okay, you want the regular, not cohesive, cells connected to a fault
vertex, which is not clamped, and you want both sides.

1) As long as a vertex is not clamped, it will only give you cells on one
side of the fault, so you can just use

  DMPlexGetTransitiveClosure(dm, v, PETSC_FALSE, ...)

http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexGetTransitiveClosure.html

and then filter out anything not in [0, cMax) where you get cMax from

  DMPlexGetHybridBounds(dm, &cMax, ...)

http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexGetHybridBounds.html

  2) You will need to get the shadow vertex of this vertex on the other
side of the fault in order to get the cells over there.
      This is stored inside PyLith, but Brad will remember better than I do
which structure holds this info.

  Thanks,

     Matt

Thank you very much
> Birendra
>
> --- On Sun, 11/15/15, Birendra jha <bjha7333 at yahoo.com> wrote:
>
> > From: Birendra jha <bjha7333 at yahoo.com>
> > Subject: cells connected to a fault vertex
> > To: cig-short at geodynamics.org
> > Date: Sunday, November 15, 2015, 2:19 PM
> > Hi
> >
> > How can I get cells connected to a given fault vertex, both
> > positive and negative side cells? I used sieve's
> > orientedStar with v_positive and v_negative but this gives
> > all faces at the vertex plus the vertex number of that
> > vertex. Do I need to use sieve's support on individual
> > faces, as the second step, to get the cells on both sides of
> > each face? This will give duplicate cell numbers.
> >
> > Thanks
> > Birendra
> >
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20151116/1c2a47ae/attachment.html>