<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Feb 24, 2015 at 12:57 PM, Lawrence Mitchell <span dir="ltr"><<a href="mailto:lawrence.mitchell@imperial.ac.uk" target="_blank">lawrence.mitchell@imperial.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">-----BEGIN PGP SIGNED MESSAGE-----<br>
Hash: SHA1<br>
<br>
Hi all,<br>
<br>
I'd like to be able to provide a matrix that is different to the<br>
actual operator for preconditioning purposes on coarse levels.  I<br>
think I see how to do this for PCMG, I do KSPSetComputeOperators which<br>
migrates to the DM and will be called to setup the operators on first<br>
KSPSolve calls.<br>
<br>
On the other hand, for FAS, it looks like there is no way to do this:<br>
<br>
The SNESen on each level will call SNESSetupMatrices which, if the<br>
snes doesn't have a jacobian matrix set, will call DMCreateMatrix once<br>
to create snes->jacobian and use that for both the jacobian and the<br>
preconditioner for the jacobian.  Is this right, or do I have the<br>
wrong end of the stick?<br></blockquote><div><br></div><div>That is what happens. It does check for snes->jacobian, but that would not be</div><div>set on coarse levels anyway. I think we should have an option that makes</div><div>SNESSetUpMatrices() create two matrices.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Cheers,<br>
<br>
Lawrence<br>
-----BEGIN PGP SIGNATURE-----<br>
Version: GnuPG v1<br>
<br>
iQEcBAEBAgAGBQJU7MmXAAoJECOc1kQ8PEYvgdMH/ip3m0z2g6wvURkiP680QZ63<br>
peAKK9R1F5TJZt1OjZR9d99YcplHwsTi7EemEAJ6kaDlo1QVJ+aGTh00EcaR9yiT<br>
lT0zlJQ8vsSa7Wwn12OGa02kZ6jh9wdRJJOhpSJoqdSih2zMhlcBYjgR1VupxpbT<br>
ALF9TMrxi1086y7237ukbcuQh4MCidzY8t5pwglNGjGVzGMRprWGhW+0nv9Wie7f<br>
gY6r0NAfyGc3Yi7FQvPupi07WMo1hDwnRaWVah3w11JCqVRPZXs3XxKD/KaZO/cN<br>
0HJ70WDNJWCR1gnQ2fdFIR5grm5Dr6xIJ2rtUbBCRwldi3ZhdcmHVfjMD3o75r0=<br>
=AB63<br>
-----END PGP SIGNATURE-----<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>