8 Aug
2011
8 Aug
'11
5:16 p.m.
On Aug 8, 2011, at 2:54 PM, Jed Brown wrote:
On Mon, Aug 8, 2011 at 14:44, Barry Smith <[email protected]> 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