[petsc-users] Solving multiple linear systems with similar matrix sequentially

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:
>    Freeze the hierarchy and coarse grid interpolations of GAMG but
>    compute the new coarse grid operators RAP for each new linear
>    system (this is a much cheaper operation). Use
>    PCGAMGSetReuseInterpolation() to freeze/unfreeze the hierarchy.

The RAP is often more than 50% of PCSetUp, so this might not save much.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150223/719456e1/attachment.pgp>