Re: [petsc-users] Mat indices for DMPlex jacobian

2024-09-23 Thread Matthew Knepley
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

[petsc-users] Mat indices for DMPlex jacobian

2024-09-23 Thread Matteo Semplice via petsc-users
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