[petsc-users] Fwd: Implementing a homotopy solver

[zakaryah]

OK, thanks, I understand now.

FormJacobian() has two matrices as arguments, the Jacobian Amat and the
preconditioner Pmat, which may be the same.  My FormJacobian() routine sets
the values of the Jacobian in Pmat.  So, to calculate the effective
nullspace of the submatrix, I want to solve Pmat u = b.  To do this within
FormJacobian(), can I assemble Pmat, call SNESGetKSP(snes,&ksp), and then
KSPSolve(ksp,b,u)?  Is it safe to assume that the KSP matrix is Pmat?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180722/adc827e8/attachment.html>