[petsc-users] Regarding ksp ex42 - Citations

[Lawrence Mitchell]

> On 21 Jul 2016, at 10:55, domenico lahaye <domenico_lahaye at yahoo.com> wrote:
> 
> Apologies for being not sufficient clear in my previous message.
> 
> I would like to be able to Galerkin coarsen A^h to obtain A^H
> and to separately Galerkin coarsen M^h to obtain M^H.
> 
> So, yes, the way in which I currently (partially) understand your
> description of the new DMCreateMatrices would do the job.

If you want to separately coarsen A and M via Galerkin, I think it will be easier to just change the code in PCSetUp_MG to handle the case where A and M are different on the coarse levels.  Effectively you just need to replicate the code that computes the coarse grid "B" matrix to separately compute coarse grid A and B matrices and pass them in to KSPSetOperators.

Cheers,

Lawrence

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 455 bytes
Desc: Message signed with OpenPGP using GPGMail
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160721/93008bf6/attachment.pgp>