<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Nov 30, 2017 at 7:41 AM, Buesing, Henrik <span dir="ltr"><<a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div lang="DE" link="blue" vlink="purple">
<div class="m_-2496744380777257325WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">Dear Matt,
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">I wanted to try out NGMRES. I am using options [1] from your tutorial slides. This works fine with SNES ex19, but fails with error [2] when
 I use my code. When I use my code with SNES options “-snes_type newtonls -snes_linesearch_type l2” it runs fine.
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">Is this me doing something wrong in FormJacobian or what could this be? Is there a SNES example with MatSetValuesStencil I could try with NGMRES?</span></p></div></div></blockquote><div><br></div><div>1) There should be more to the stack frame. I assume you call MatSetValuesStencil(). Do you use PetscFunctionBegin/Return() in that function? It would add to the stack</div><div><br></div><div>2) I am guessing that you are using a DMDA. There should be an ISL2G there automatically. It does not make sense to me that its gone,</div><div>and ex19 works that way. The first thing to do to track it down is to get the full stack. Then I think we probably need an example small</div><div>enough to run here to track through everything, unless you are good with teh debugger.</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"><div lang="DE" link="blue" vlink="purple"><div class="m_-2496744380777257325WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">Thank you!<br>
Henrik<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[1]<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">-snes_type ngmres -npc_snes_max_it 1 -snes_converged_reason -npc_snes_type fas -npc_fas_coarse_snes_<wbr>converged_reason -npc_fas_levels_snes_type
 newtonls -npc_fas_levels_snes_max_it 6 -npc_fas_levels_snes_<wbr>linesearch_type basic -npc_fas_levels_snes_max_<wbr>linear_solve_fail 30 -npc_fas_levels_ksp_max_it 20 -npc_fas_levels_snes_<wbr>converged_reason -npc_fas_coarse_snes_<wbr>linesearch_type basic<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[2]<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: Null argument, when expecting valid pointer<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: Null Object: Parameter # 1<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: Petsc Release Version 3.8.2, unknown
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: /rwthfs/rz/cluster/work/<wbr>hb111949/Descramble/1D_Model/<wbr>steam/JUBE/ngmres/debug/shem_<wbr>fw64gnu_steam.x on a arch-linux2-c-debug
 named <a href="http://linuxbmc0004.rz.RWTH-Aachen.DE" target="_blank">linuxbmc0004.rz.RWTH-Aachen.DE</a> by hb111949 Thu Nov 30 14:35:56 2017<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: Configure options --download-fblaslapack --with-cc=mpicc --with-fc=mpif90 --with-cxx=mpicxx --download-hypre --download-superlu_dist
 --download-suitesparse --download-scalapack --download-blacs --download-hdf5 --download-parmetis --download-metis --with-debugging=1 --download-mumps<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: #1 ISLocalToGlobalMappingApply() line 639 in /rwthfs/rz/cluster/home/<wbr>hb111949/Code/petsc/src/vec/<wbr>is/utils/isltog.c<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: #2 MatSetValuesLocal() line 2139 in /rwthfs/rz/cluster/home/<wbr>hb111949/Code/petsc/src/mat/<wbr>interface/matrix.c<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">[0]PETSC ERROR: #3 MatSetValuesStencil() line 1550 in /rwthfs/rz/cluster/home/<wbr>hb111949/Code/petsc/src/mat/<wbr>interface/matrix.c<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">--
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">Dipl.-Math. Henrik Büsing<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">Institute for Applied Geophysics and Geothermal Energy<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">E.ON Energy Research Center<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Courier New"">RWTH Aachen University<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">------------------------------<wbr>------------------------<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""><a href="https://maps.google.com/?q=Mathieustr.+10&entry=gmail&source=g">Mathieustr. 10</a>            |    Tel <a href="tel:+49%20241%208049907" value="+492418049907" target="_blank">+49 (0)241 80 49907</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">52074 Aachen, Germany     |    Fax <a href="tel:+49%20241%208049889" value="+492418049889" target="_blank">+49 (0)241 80 49889</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">------------------------------<wbr>------------------------<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""><a href="http://www.eonerc.rwth-aachen.de/GGE" target="_blank">http://www.eonerc.rwth-aachen.<wbr>de/GGE</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""><a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.de</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New"">------------------------------<wbr>------------------------<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Courier New""><u></u> <u></u></span></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #e1e1e1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Von:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Matthew Knepley [mailto:<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>]
<br>
<b>Gesendet:</b> 07 November 2017 12:54<br>
<b>An:</b> Buesing, Henrik <<a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.<wbr>de</a>><br>
<b>Cc:</b> petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Betreff:</b> Re: [petsc-users] Newton methods that converge all the time<u></u><u></u></span></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<p class="MsoNormal">On Tue, Nov 7, 2017 at 4:19 AM, Buesing, Henrik <<a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.<wbr>de</a>> wrote:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">Dear all,
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New""> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">I am solving a system of nonlinear, transient PDEs. I am using Newton’s method in every time step to solve the nonlinear algebraic
 equations. Of course, Newton’s method only converges if the initial guess is sufficiently close to the solution.
<br>
<br>
This is often not the case and Newton’s method diverges. Then, I reduce the time step and try again. This can become prohibitively costly, if the time steps get very small. I am thus looking for variants of Newton’s method, which have a bigger convergence radius
 or ideally converge all the time. </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New""> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">I tried out the pseudo-timestepping described in
<a href="http://www.mcs.anl.gov/petsc/petsc-current/src/ts/examples/tutorials/ex1f.F.html" target="_blank">
http://www.mcs.anl.gov/petsc/<wbr>petsc-current/src/ts/examples/<wbr>tutorials/ex1f.F.html</a>.
</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New""> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">However, this does converge even worse. I am seeing breakdown when I have phase changes (e.g. liquid to two-phase).
</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New""> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">I was under the impression that pseudo-timestepping should converge better. Thus, my question:
<br>
<br>
Am I doing something wrong or is it possible that Newton’s method converges and pseudo-timestepping does not?<br>
<br>
Thank you for any insight on this. </span><u></u><u></u></p>
</div>
</div>
</blockquote>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Hi Hendrik,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I would try using NGMRES as a nonlinear preconditioner. I have an example in my tutorial slides for using it with SNES ex19.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">I hope this will work because I suspect that around the phase boundary Newton directions are noisy, since sometimes you<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">step into the other phase. NGMRES takes a few directions (you set the m) and then picks the best one.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">  Hopefully this helps,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">     Matt<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">Henrik</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span lang="EN-US" style="font-family:"Courier New""><br>
<br>
</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New""> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">--
</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">Dipl.-Math. Henrik Büsing</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">Institute for Applied Geophysics and Geothermal Energy</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">E.ON Energy Research Center</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US" style="font-family:"Courier New"">RWTH Aachen University</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New"">------------------------------<wbr>------------------------</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New""><a href="https://maps.google.com/?q=Mathieustr.+10&entry=gmail&source=g" target="_blank">Mathieustr. 10</a>            |    Tel
<a href="tel:+49%20241%208049907" target="_blank">+49 (0)241 80 49907</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New"">52074 Aachen, Germany     |    Fax
<a href="tel:+49%20241%208049889" target="_blank">+49 (0)241 80 49889</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New"">------------------------------<wbr>------------------------</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New""><a href="http://www.eonerc.rwth-aachen.de/GGE" target="_blank">http://www.eonerc.rwth-aachen.<wbr>de/GGE</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New""><a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.de</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-family:"Courier New"">------------------------------<wbr>------------------------</span><u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><br>
<br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">-- <u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">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<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~<wbr>knepley/</a><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>