<div dir="ltr"><div>That's the thing; you WANT a MFFD matrix in the snes->jacobian slot but not the snes->jacobian_pre slot.<br><br></div><div>Somewhere the PC's getting the wrong one and putting it on the smoother.<br>
</div><div><br></div>- Peter<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 6, 2013 at 3:44 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">><br>
>  Yikes, the logic of SNES and MG (with dm etc) is getting a bit too convoluted. I would run in the debugger with a break point for MatCreate_MFFD() this will give a hint why it is being used.<br>
><br>
<br>
</div>#0  0x0000000100fca059 in MatCreate_MFFD ()<br>
#1  0x0000000100eedc34 in MatSetType ()<br>
#2  0x0000000100fc7bc3 in MatCreateMFFD ()<br>
#3  0x00000001012acc90 in MatCreateSNESMF ()<br>
#4  0x00000001012c092e in SNESSetUp ()<br>
#5  0x00000001012c277b in SNESSolve ()<br>
<br>
I will keep digging.  This is the first one and it is not in the PC ...</blockquote></div><br></div>