[petsc-dev] sor smoothers

[Jed Brown]

"Mark F. Adams" <mfadams at lbl.gov> writes:

> So I should add a residual method to Mat, get rid of PCMGResidualDefault, make a MatResidualDefault, and have SetUp_MatAIJ set this default.
>
> I guess PCMGSetResidual would just set the residual for Mat to keep
> the same interface or should we just get rid of it and make a
> MatSetResidual?

PCMGSetResidual (if we continue to support it) should not modify the Mat
because we can't let its effects escape the solve.  That is, we should
be able to create two solvers and not have the presence of the first
affect the algorithm used when we pass the matrix on to the second.

Note that there is MatShellSetOperation.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130908/9cd5e9f0/attachment.sig>