That's the thing; you WANT a MFFD matrix in the snes->jacobian slot but not the snes->jacobian_pre slot.
Somewhere the PC's getting the wrong one and putting it on the smoother. - Peter On Wed, Mar 6, 2013 at 3:44 PM, Mark F. Adams <mark.adams at columbia.edu>wrote: > > > > 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. > > > > #0 0x0000000100fca059 in MatCreate_MFFD () > #1 0x0000000100eedc34 in MatSetType () > #2 0x0000000100fc7bc3 in MatCreateMFFD () > #3 0x00000001012acc90 in MatCreateSNESMF () > #4 0x00000001012c092e in SNESSetUp () > #5 0x00000001012c277b in SNESSolve () > > I will keep digging. This is the first one and it is not in the PC ... -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130306/f6871cd9/attachment.html>
