[petsc-users] Using Pmat different from Amat

Chung-Kan Huang ckhuangf at gmail.com
Fri Feb 21 14:41:03 CST 2014


In my application I like to use
PetscErrorCode  KSPSetOperators(KSP ksp,Mat Amat,Mat Pmat,MatStructure flag)
and have Pmat different from Amat

if Amat = L + D + U
then Pmat  = Amat  - L* - U* + rowsum(L* + U*)
where L* and U* is a portion of L and U and rowsum(L* + U*) will add into D

I know I can explicitly construct Pmat step-by-step but wonder what will be
the most easy way to do this in PETSC?


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140221/653a89aa/attachment-0001.html>

More information about the petsc-users mailing list