Coming back to this after a while. It would still be good if I
could provide a callback for user-defined adjacency. Does
https://bitbucket.org/petsc/petsc/pull-requests/690/plex-support-user-defined-adjacencies-via/diff
On Fri, Aug 4, 2017 at 9:37 AM, Lawrence Mitchell <
lawrence.mitch...@imperial.ac.uk> wrote:
>
> > On 6 Jun 2017, at 09:01, Lawrence Mitchell ac.uk> wrote:
> >
> >> If you do really want to prune them, then I guess overriding the
> DMPlexGetAdjacency() as you
> On 6 Jun 2017, at 09:01, Lawrence Mitchell
> wrote:
>
>> If you do really want to prune them, then I guess overriding the
>> DMPlexGetAdjacency() as you propose is probably the best way. I would
>> be willing to put it in. Please send me a reminder email
> On 5 Jun 2017, at 23:01, Matthew Knepley wrote:
>
> To do a FEM integral on a facet f I need:
>
> i) to evaluate coefficients at quadrature points (on the facet)
> ii) to evaluate basis functions at quadrature points (on the facet)
>
> for (i), I need all the dofs in
On Mon, Jun 5, 2017 at 12:18 PM, Lawrence Mitchell <
lawrence.mitch...@imperial.ac.uk> wrote:
>
> > On 2 Jun 2017, at 05:09, Matthew Knepley wrote:
> >
> >
> > Coming back to this, I think I understand the problem a little better.
> >
> > Consider this mesh:
> >
> > ++
> >
> On 2 Jun 2017, at 05:09, Matthew Knepley wrote:
>
>
> Coming back to this, I think I understand the problem a little better.
>
> Consider this mesh:
>
> ++
> |\ 3 |
> | \ |
> |2 \ |
> | \|
> ++
> |\ 1 |
> | \ |
> |0 \ |
> | \|
> ++
>
> Let's say I run