> As Matt notes you absolutely need to run with -pc_type lu
> -ksp_monitor_true_residual -ksp_converged_reason to make sure that the
Why lu? I was running with bjacobi for whatever reason, which I no longer
recall. Probably gave the fastest convergence before introducing the horribly
nonlinear terms.
I just tried lu again: it converges, too, but snes still diverges. Running
with -pc_type none requires too many KSP iterations, I had a limit of 1e6 and
ran out. I will revisit this while I wait for any further replies, but as it
takes so long, I will finish this reply first.
> linear system solve is not the issue. It is also possible that a bug in
> the function evaluation (using uninitialized values or something) is
> causing the difficulty: run with valgrind
Valgrind complains a LOT about openmpi, which guess is not the issue here. It
also complains about PetscStrcasecmp (sic!) and some other Petsc* functions,
but my code is clean: the last complaint from valgrind (apart from the summary
at the end) comes from SNESSetFromOptions, which is called before I construct
the initial guess.
Cheers,
-Juha
--
-----------------------------------------------
| Juha J?ykk?, juhaj at iki.fi |
| http://www.maths.leeds.ac.uk/~juhaj |
-----------------------------------------------
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 836 bytes
Desc: This is a digitally signed message part.
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110623/417a20ff/attachment.pgp>
