[petsc-users] (Fast?) assembly of upper triangle and MatZeroRows
Barry Smith
bsmith at mcs.anl.gov
Mon Aug 8 17:16:24 CDT 2011
On Aug 8, 2011, at 2:54 PM, Jed Brown wrote:
> On Mon, Aug 8, 2011 at 14:44, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > How would you use this interface with SNES?
>
> The user calls this as part the FormJacobian() (in the same way as if they explicitly remove those rows and columns (what Matt likes) in the matrix assembly.
>
> But they need access to the right hand side at the same time. It seems too fragile to reach into the SNES and modify the residual vector during Jacobian evaluation.
Add a KSP/SNES./TSSetDirichlet(ksp,snes,ts,IS listofnodes,Vec valuesfornodes)?
Barry
More information about the petsc-users
mailing list