<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, May 20, 2016 at 10:11 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="">> not. This usually indicates a bad guess for Newton. I would now advocate<br>
<br>
</span>Oh, did I forget to mention, that other initial conditions DO converge. The<br>
funny thing is the initial condition which does NOT converge is actually<br>
CLOSER to the solution than ones which DO. Had it been the other way around, I<br>
probably would not have thought twice of it.<br>
<span class=""><br>
>   1) Grid sequencing: This is easy if you use a DMDA. You just use<br>
> -snes_grid_sequence<br>
>       and its automatic. Since you report that smaller grids converge, this<br>
> is usually enough.<br>
<br>
</span>Unfortunately, no.</blockquote><div><br></div><div>Isn't this a 1D problem with no geometry? You should use DMDA. It would make</div><div>it easier.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
>   2) NPC: If the above fails, I would try preconditioning Newton with GS.<br>
> You can do this<br>
>       just using multiplicative composition, or left preconditioning.<br>
<br>
</span>I'm not sure what you mean here. Do you mean I should run normal newton line<br>
search snes with SNESSetNPC(snes_of_type_gs)?<br></blockquote><div><br></div><div>Yes.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
If so, does the inner snes_of_type_gs and the outer snes_of_type_newtonls get<br>
all the same routines? I guess so, but just checking.</blockquote><div><br></div><div>The outer SNES get the same routines, and the inner one is setup automatically.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
> If 1 works, then its likely that FAS would work even better.<br>
<br>
</span>I wish it was so easy: I did try all the snes types before posting the first<br>
post. Also all KSP types. Except those that need things I don't have, of<br>
course.<br></blockquote><div><br></div><div>Not sure what you mean here. To use FAS, you can either use DMDA, or provide</div><div>interpolation operators between grids.</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">
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>