On Mon, Sep 23, 2024 at 5:33 AM Matteo Semplice via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Dear petsc,
>
> I need to hand-code a jacobian and I can't figure out how to
> translate the DMplex points/fields to matrix indices.
>
> My DMPlex has a section with m fields per cell (which ma
Dear petsc,
I need to hand-code a jacobian and I can't figure out how to
translate the DMplex points/fields to matrix indices.
My DMPlex has a section with m fields per cell (which makes for n>m dof
per cell since some fields are vector). Say I want to insert an nxn
block for the row cor