<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, May 20, 2016 at 9:30 AM, Juha Jaykka <span dir="ltr"><<a href="mailto:juhaj@iki.fi" target="_blank">juhaj@iki.fi</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">> > It does not. And I should have mentioned earlier, that I tried -snes_mf, -<br>
> > snes_mf_operator, -snes_fd and -snes_fd_color already and none of those<br>
> > converges. Your suggested options result in<br>
><br>
> For solver questions, I always need to see the result of<br>
><br>
>   -snes_view -snes_monitor -ksp_converged_reason -ksp_monitor_true_residual<br>
> -snes_converged_reason<br>
><br>
> Turn on LU for the linear solver so it plays no role.<br>
<br>
</span>I'm not sure what you mean by the last sentence, but I ran with -pc_type lu<br>
and the other options you suggested. Output is in<br>
<a href="http://paste.debian.net/688191/" rel="noreferrer" target="_blank">http://paste.debian.net/688191/</a><br>
<br>
I don't see anything wrong in it except it does not converge. :) And it<br>
settles in a weird oscillatory solution (looks almost like a sin()) whereas<br>
the real solution is 4*arctan(exp(x)).<br></blockquote><div><br></div><div>Okay, so the linear solves are accurate and we think the Jacobian is correct because</div><div>you ran with -snes_test, or ran using -snes_fd and it made no difference.</div><div><br></div><div>Nonlinear Gauss-Siedel (coordinate descent) converges but slowly, and Newton does</div><div>not. This usually indicates a bad guess for Newton. I would now advocate trying at least</div><div>two things:</div><div><br></div><div>  1) Grid sequencing: This is easy if you use a DMDA. You just use -snes_grid_sequence</div><div>      and its automatic. Since you report that smaller grids converge, this is usually enough.</div><div><br></div><div>  2) NPC: If the above fails, I would try preconditioning Newton with GS. You can do this</div><div>      just using multiplicative composition, or left preconditioning.</div><div><br></div><div>If 1 works, then its likely that FAS would work even better.</div><div><br></div><div>  Thanks,</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">
Thanks for your help.<br>
<br>
Cheers,<br>
Juha<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>