<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:10.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
p.m6963147566296146480p1, li.m6963147566296146480p1, div.m6963147566296146480p1
{mso-style-name:m_6963147566296146480p1;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
p.m6963147566296146480m477101460960849173p1, li.m6963147566296146480m477101460960849173p1, div.m6963147566296146480m477101460960849173p1
{mso-style-name:m_6963147566296146480m477101460960849173p1;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.m6963147566296146480m477101460960849173apple-converted-space
{mso-style-name:m_6963147566296146480m477101460960849173apple-converted-space;}
p.m6963147566296146480m477101460960849173p2, li.m6963147566296146480m477101460960849173p2, div.m6963147566296146480m477101460960849173p2
{mso-style-name:m_6963147566296146480m477101460960849173p2;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.m6963147566296146480gmailsignatureprefix
{mso-style-name:m_6963147566296146480gmailsignatureprefix;}
span.gmailsignatureprefix
{mso-style-name:gmail_signature_prefix;}
span.EmailStyle24
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-GB" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US">Thank you for your response.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US">I will test the Jacobian.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US">The linear solve in the second step convergences. Is it possible that Jacobian had been zeroed out?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US">Best,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US">Karthik.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div id="mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Matthew Knepley <knepley@gmail.com><br>
<b>Date: </b>Thursday, 7 September 2</span><span style="font-size:12.0pt"> ot<span style="color:black">023 at 14:19<br>
<b>To: </b>Chockalingam, Karthikeyan (STFC,DL,HC) <karthikeyan.chockalingam@stfc.ac.uk><br>
<b>Cc: </b>Stefano Zampini <stefano.zampini@gmail.com>, petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Subject: </b>Re: [petsc-users] Non-linear solve: DIVERGED_LINE_SEARCH<o:p></o:p></span></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">On Thu, Sep 7, 2023 at 8:49 AM Karthikeyan Chockalingam - STFC UKRI via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<o:p></o:p></span></p>
</div>
<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>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">I just noticed that the line search diverges even for the linear Poisson equation (where the Jacobian is a constant). I am not sure, what I am doing
wrong.<o:p></o:p></span></p>
</div>
</div>
</div>
</blockquote>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">You need to test the Jacobian. We have facilities for this. Run using<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"> --snes_test_jacobian<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">This will tell you at each step whether the Jacobian matches the finite difference approximant. It looks like your Jacobian at the second step is wrong. Is it possible that you are not zeroing out the matrix
before adding in updates? That would make the second step wrong.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"> Thanks,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"> Matt<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"> <o:p></o:p></span></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>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Best,<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Karthik.<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div id="m_6963147566296146480mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Chockalingam, Karthikeyan (STFC,DL,HC) <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>><br>
<b>Date: </b>Thursday, 7 September 2023 at 12:19<br>
<b>To: </b>Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>><br>
<b>Cc: </b><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject: </b>Re: [petsc-users] Non-linear solve: DIVERGED_LINE_SEARCH</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">I am trying to solve the non-linear scalar problem<o:p></o:p></span></p>
<p class="m6963147566296146480p1">-\nabla \cdot \big( (1+u^2)\nabla u\big)=f<o:p></o:p></p>
<p class="m6963147566296146480p1">Using finite element formulation.<o:p></o:p></p>
<p class="m6963147566296146480p1"> <o:p></o:p></p>
<p class="m6963147566296146480p1">The linear iteration converged – does that mean my Jacobian is correct??<o:p></o:p></p>
<p class="m6963147566296146480p1"> <o:p></o:p></p>
<p class="m6963147566296146480p1">I am hoping the non-linear residual goes down to 1-e11.<o:p></o:p></p>
<p class="m6963147566296146480p1"> <o:p></o:p></p>
<p class="m6963147566296146480p1">Kind regards,<o:p></o:p></p>
<p class="m6963147566296146480p1">Karthik.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div id="m_6963147566296146480mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>><br>
<b>Date: </b>Thursday, 7 September 2023 at 11:49<br>
<b>To: </b>Chockalingam, Karthikeyan (STFC,DL,HC) <<a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" target="_blank">karthikeyan.chockalingam@stfc.ac.uk</a>><br>
<b>Cc: </b><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject: </b>Re: [petsc-users] Non-linear solve: DIVERGED_LINE_SEARCH</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">The solver did not diverge.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">It was the line search that was not able to make further progress in minimizing the 2-norm of the residual.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">This is common in nonlinear solvers. It would help if you tell us what you are trying to solve.<o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Note that at the first step, your residual norm is already 1.e-6. What kind of accuracy do you want?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
</div>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Il giorno gio 7 set 2023 alle ore 13:08 Karthikeyan Chockalingam - STFC UKRI via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>>
ha scritto:<o:p></o:p></span></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Hello,<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">The non-linear solution diverged. The final solution seems right and I believe the Jacobian is correct (not 100% certain).
<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">I am not sure if I doing something wrong in the solver setting.
<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black"> </span></span><span style="color:black">0 SNES Function norm 3.890991210938e-03<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">0 KSP Residual norm 9.037762538598e+00<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">1 KSP Residual norm 2.120375403775e-01<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">2 KSP Residual norm 5.155439334511e-03<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">3 KSP Residual norm 1.394364169369e-04<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">4 KSP Residual norm 9.233543407204e-06<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Linear solve converged due to CONVERGED_RTOL iterations 4</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p2"><span style="color:black"> </span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Using full step: fnorm 3.890991210938e-03 gnorm 7.701511565083e-06</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">1 SNES Function norm 7.701511565083e-06<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">0 KSP Residual norm 5.630229687829e-03<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">1 KSP Residual norm 1.030475601271e-04<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">2 KSP Residual norm 2.576454714319e-06<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">3 KSP Residual norm 6.669316846898e-08<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">4 KSP Residual norm 3.215810984829e-09<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Linear solve converged due to CONVERGED_RTOL iterations 4</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p2"><span style="color:black"> </span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: gnorm after quadratic fit 1.805500533481e-05</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 2.563759884284e-05 lambda=3.0804668685096816e-02</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 3.332721829751e-05 lambda=3.0804668685096817e-03</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 4.102754045833e-05 lambda=3.0804668685096822e-04</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 4.872893294880e-05 lambda=3.0804668685096822e-05</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 5.643043250787e-05 lambda=3.0804668685096822e-06</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 6.413194279696e-05 lambda=3.0804668685096827e-07</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 7.183345417492e-05 lambda=3.0804668685096828e-08</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 7.953496567312e-05 lambda=3.0804668685096829e-09</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 8.723647719173e-05 lambda=3.0804668685096831e-10</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 9.493798871875e-05 lambda=3.0804668685096832e-11</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 1.026395002516e-04 lambda=3.0804668685096835e-12</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: Cubic step no good, shrinking lambda, current gnorm 1.103410117889e-04 lambda=3.0804668685096835e-13</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: unable to find good step length! After 12 tries<span class="m6963147566296146480m477101460960849173apple-converted-space"> </span></span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span class="m6963147566296146480m477101460960849173apple-converted-space"><span style="color:black">
</span></span><span style="color:black">Line search: fnorm=7.7015115650831560e-06, gnorm=1.1034101178892401e-04, ynorm=6.2052357872955976e-03, minlambda=9.9999999999999998e-13, lambda=3.0804668685096835e-13, initial slope=-5.9313318983096354e-11</span><o:p></o:p></p>
<p class="m6963147566296146480m477101460960849173p1"><span style="color:black">Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Thank you for your help.<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Kind regards,<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Karthik.<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">-- <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#2D2D62">Dr. Karthik Chockalingam</span></b><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black">Senior Research Software Engineer</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black">High Performance Systems Engineering Group</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black">Hartree Centre | Science and Technology Facilities Council</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"><a href="mailto:karthikeyan.chockalingam@stfc.ac.uk" target="_blank"><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#0563C1">karthikeyan.chockalingam@stfc.ac.uk</span></a><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:black"> </span><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#414141"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#414141"> </span></b><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#1F497D"><b>Error!
Filename not specified.</b></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"><br clear="all">
<br>
<span class="m6963147566296146480gmailsignatureprefix">-- </span><o:p></o:p></span></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Stefano<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><br clear="all">
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span class="gmailsignatureprefix"><span style="font-size:11.0pt">--
</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">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<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>