[petsc-users] Preconditioning in matrix free SNES

[Jed Brown]

I applaud your attention to detail if you worked out that Jacobian without
AD. There has been some success preconditioning high order methods with
nonlinear reconstruction using a matrix assembled with first order upwinding
(sparser and less messy to work out). I works be interested to hear if that
works well for you.

On Feb 4, 2011 8:17 PM, "Travis C. Fisher" <travis.fisher at nasa.gov> wrote:

Jed,

Thanks for the response. I think I have resolved this particular issue. My
jacobian was incorrect. I use generalized coordinates and
I missed a scaling factor. I found this by generating the matrix finite
difference coloring and calculating
the approximate jacobian, which was much easier than I originally
anticipated. The jacobian I am trying to create is exact for convective
terms based on the numerical methods (ESWENO finite difference schemes).

Travis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110205/38856f3e/attachment.htm>