<div dir="ltr"><div dir="ltr">On Fri, Jan 13, 2023 at 3:21 PM Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr">alexis.marboeuf@hotmail.fr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg5004949584998710548">
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Hi Matt,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Here is the output from Petsc when I view the TAO solver:</span></div></div></div></blockquote><div><br></div><div>Okay, there is no sophisticated caching going on. So, first I would get it to fail on1 process. It should if it just depends on</div><div>the convergence (I hope). Then send the source so we can run it. It should be simple for us to find where system size</div><div>changes and the KSP is not reset (if that indeed is what happens).</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg5004949584998710548"><div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Tao Object: (Damage_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: bntr</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Tao Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: bncg</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>CG Type: ssml_bfgs</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Skipped Stepdirection Updates: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scaled gradient steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Pure gradient steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Not a descent direction: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Line search fails: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Matrix has not been preallocated yet</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>TaoLineSearch Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: more-thuente</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>maximum function evaluations=30</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of function evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of gradient evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of function/gradient evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Termination reason: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Active Set subset type: subvec</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>convergence tolerances: gatol=1e-08, <span>
</span>steptol=0., <span>
</span>gttol=0.</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Residual in Function/Gradient:=0.</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Objective value=0.</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of iterations=0,<span>
</span>(max: 2000)</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Solver terminated: 0 <span> </span></span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Rejected BFGS updates: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>CG steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Newton steps: 11</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>BFGS steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scaled gradient steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Gradient steps: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>KSP termination reasons:</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>atol: 4</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>rtol: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>ctol: 7</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>negc: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>dtol: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>iter: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>othr: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>TaoLineSearch Object: (Damage_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: more-thuente</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>maximum function evaluations=30</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of function evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of gradient evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of function/gradient evaluations=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>using variable bounds</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Termination reason: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>KSP Object: (Damage_tao_bnk_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: stcg</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>maximum iterations=10000, nonzero initial guess</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>tolerances:<span> </span>relative=1e-08, absolute=1e-08, divergence=1e+10</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>left preconditioning</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>using UNPRECONDITIONED norm type for convergence test</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>PC Object: (Damage_tao_bnk_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: lmvm</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Mat Object: (Damage_tao_bnk_pc_lmvm_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: lmvmbfgs</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>rows=30634, cols=30634</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scale type: DIAGONAL</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scale history: 1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scale params: alpha=1., beta=0.5, rho=1.</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Convex factors: phi=0., theta=0.125</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Max. storage: 5</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Used storage: 5</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of updates: 11</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of rejects: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of resets: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Mat Object: (Damage_tao_bnk_pc_lmvm_J0_) 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: lmvmdiagbroyden</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>rows=30634, cols=30634</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scale history: 1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Scale params: alpha=1., beta=0.5, rho=1.</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Convex factor: theta=0.125</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Max. storage: 1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Used storage: 1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of updates: 11</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of rejects: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Number of resets: 0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>linear system matrix = precond matrix:</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Mat Object: 4 MPI processes</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>type: mpiaij</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>rows=468, cols=468</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total: nonzeros=2932, allocated nonzeros=2932</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of mallocs used during MatSetValues calls=0</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>not using I-node (on process 0) routines</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total KSP iterations: 103</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Active Set subset type: subvec</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>convergence tolerances: gatol=0.0001, <span>
</span>steptol=0., <span> </span>
gttol=1e-05</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Residual in Function/Gradient:=9.11153e-05</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Objective value=0.00665458</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of iterations=11,<span>
</span>(max: 50)</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of function evaluations=17,<span>
</span>max: -1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of gradient evaluations=13,<span>
</span>max: -1</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>total number of Hessian evaluations=12</span></p>
<p style="margin:0px;font:11px Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span>
</span>Solution converged:<span> </span>
||g(X)|| <= gatol</span></p>
<br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Thanks again for your help!</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Alexis</span></div>
<div id="m_-7317894617606137971appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_-7317894617606137971divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>De :</b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
<b>Envoyé :</b> samedi 14 janvier 2023 01:38<br>
<b>À :</b> Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr" target="_blank">alexis.marboeuf@hotmail.fr</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>Objet :</b> Re: [petsc-users] Nonconforming object sizes using TAO (TAOBNTR)</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">On Fri, Jan 13, 2023 at 11:22 AM Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr" target="_blank">alexis.marboeuf@hotmail.fr</a>> wrote:<br>
</div>
<div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Hi all,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span style="background-color:rgb(255,255,255);display:inline">In a variational approach of brittle fracture setting, </span>I
try to solve a bound constraint minimization problem using TAO. I checkout on the main branch of Petsc. Minimization with respect to the bounded variable (damage) is achieved through the Bounded Newton Trust Region (TAOBNTR). All other TAO parameters are set
by default. On a Linux machine, I get the following error with a 4 processors run:</span></div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Can you view the solver?</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<div style="font-weight:normal;color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Menlo,Monaco,"Courier New",monospace;font-size:12px;line-height:18px">
<span><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: ---------------------
</span><span style="color:rgb(206,145,120);font-weight:bold">Error</span><span> Message --------------------------------------------------------------</span></span>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Nonconforming object sizes</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Preconditioner number of local rows
</span><span style="color:rgb(86,156,214)">1122</span><span> does not equal input vector size
</span><span style="color:rgb(86,156,214)">1161</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: See </span>
<span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/" target="_blank">https://petsc.org/release/faq/</a></span><span> for trouble shooting.</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:
</span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span> GIT Date:
</span><span style="color:rgb(106,153,85)">2023-01-04</span><span> </span><span style="color:rgb(106,153,85)">13:37:04</span><span> +</span><span style="color:rgb(86,156,214)">0000</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan
</span><span style="color:rgb(86,156,214)">12</span><span> </span><span style="color:rgb(106,153,85)">16:55:18</span><span> </span><span style="color:rgb(86,156,214)">2023</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: ---------------------
</span><span style="color:rgb(206,145,120);font-weight:bold">Error</span><span> Message --------------------------------------------------------------</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --CXXOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --FOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --download-fblaslapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-mumps=</span><span style="color:rgb(86,156,214)">1</span><span> --download-chaco=</span><span style="color:rgb(86,156,214)">1</span><span> --download-exodusii=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hypre=</span><span style="color:rgb(86,156,214)">1</span><span> --download-ml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-triangle
--download-scalapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-superlu=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing-cc=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/gcc
--download-sowing-cxx=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp --download-sowing-cxxcpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp
--download-yaml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-bison=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hdf5=</span><span style="color:rgb(86,156,214)">1</span><span> --download-metis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-parmetis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-netcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-pnetcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-zlib=</span><span style="color:rgb(86,156,214)">1</span><span> --with-cmake=</span><span style="color:rgb(86,156,214)">1</span><span> --with-debugging=</span><span style="color:rgb(86,156,214)">0</span><span> --with-mpi-dir=/opt/HPC/mvapich2/</span><span style="color:rgb(86,156,214)">2</span><span>.</span><span style="color:rgb(86,156,214)">3</span><span>.</span><span style="color:rgb(86,156,214)">7</span><span>-</span><span style="color:rgb(86,156,214)">gcc11.2.1</span><span> --with-ranlib=ranlib
--with-shared-libraries=</span><span style="color:rgb(86,156,214)">1</span><span> --with-sieve=</span><span style="color:rgb(86,156,214)">1</span><span> --download-p4est=</span><span style="color:rgb(86,156,214)">1</span><span> --with-pic --with-mpiexec=srun
--with-x11=</span><span style="color:rgb(86,156,214)">0</span><span> PETSC_ARCH=arch-darwin-c</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">1</span><span> PCApply() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:</span><span style="color:rgb(86,156,214)">434</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">2</span><span> KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:</span><span style="color:rgb(86,156,214)">380</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">3</span><span> KSPCGSolve_STCG() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:</span><span style="color:rgb(86,156,214)">76</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">4</span><span> KSPSolve_Private() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">898</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">5</span><span> KSPSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">1070</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">6</span><span> TaoBNKComputeStep() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:</span><span style="color:rgb(86,156,214)">459</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">7</span><span> TaoSolve_BNTR() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:</span><span style="color:rgb(86,156,214)">138</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">8</span><span> TaoSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:</span><span style="color:rgb(86,156,214)">177</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Nonconforming object sizes</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Preconditioner number of local rows
</span><span style="color:rgb(86,156,214)">1229</span><span> does not equal input vector size
</span><span style="color:rgb(86,156,214)">1254</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: See </span>
<span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/" target="_blank">https://petsc.org/release/faq/</a></span><span> for trouble shooting.</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:
</span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span> GIT Date:
</span><span style="color:rgb(106,153,85)">2023-01-04</span><span> </span><span style="color:rgb(106,153,85)">13:37:04</span><span> +</span><span style="color:rgb(86,156,214)">0000</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan
</span><span style="color:rgb(86,156,214)">12</span><span> </span><span style="color:rgb(106,153,85)">16:55:18</span><span> </span><span style="color:rgb(86,156,214)">2023</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --CXXOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --FOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --download-fblaslapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-mumps=</span><span style="color:rgb(86,156,214)">1</span><span> --download-chaco=</span><span style="color:rgb(86,156,214)">1</span><span> --download-exodusii=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hypre=</span><span style="color:rgb(86,156,214)">1</span><span> --download-ml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-triangle
--download-scalapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-superlu=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing-cc=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/gcc
--download-sowing-cxx=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp --download-sowing-cxxcpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp
--download-yaml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-bison=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hdf5=</span><span style="color:rgb(86,156,214)">1</span><span> --download-metis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-parmetis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-netcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-pnetcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-zlib=</span><span style="color:rgb(86,156,214)">1</span><span> --with-cmake=</span><span style="color:rgb(86,156,214)">1</span><span> --with-debugging=</span><span style="color:rgb(86,156,214)">0</span><span> --with-mpi-dir=/opt/HPC/mvapich2/</span><span style="color:rgb(86,156,214)">2</span><span>.</span><span style="color:rgb(86,156,214)">3</span><span>.</span><span style="color:rgb(86,156,214)">7</span><span>-</span><span style="color:rgb(86,156,214)">gcc11.2.1</span><span> --with-ranlib=ranlib
--with-shared-libraries=</span><span style="color:rgb(86,156,214)">1</span><span> --with-sieve=</span><span style="color:rgb(86,156,214)">1</span><span> --download-p4est=</span><span style="color:rgb(86,156,214)">1</span><span> --with-pic --with-mpiexec=srun
--with-x11=</span><span style="color:rgb(86,156,214)">0</span><span> PETSC_ARCH=arch-darwin-c</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">1</span><span> PCApply() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:</span><span style="color:rgb(86,156,214)">434</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">2</span><span> KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:</span><span style="color:rgb(86,156,214)">380</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">3</span><span> KSPCGSolve_STCG() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:</span><span style="color:rgb(86,156,214)">76</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">4</span><span> KSPSolve_Private() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">898</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">5</span><span> KSPSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">1070</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">6</span><span> TaoBNKComputeStep() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:</span><span style="color:rgb(86,156,214)">459</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">7</span><span> TaoSolve_BNTR() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:</span><span style="color:rgb(86,156,214)">138</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">8</span><span> TaoSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:</span><span style="color:rgb(86,156,214)">177</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">9</span><span> /home/marboeua/Developpement/mef90/vDef/vDefTAO.</span><span style="color:rgb(86,156,214)">F90:370</span></div>
<div><span>application called MPI_Abort(MPI_COMM_SELF, </span><span style="color:rgb(86,156,214)">60</span><span>) - process
</span><span style="color:rgb(86,156,214)">0</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC </span>
<span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">9</span><span> /home/marboeua/Developpement/mef90/vDef/vDefTAO.</span><span style="color:rgb(86,156,214)">F90:370</span></div>
<div><span>application called MPI_Abort(MPI_COMM_SELF, </span><span style="color:rgb(86,156,214)">60</span><span>) - process
</span><span style="color:rgb(86,156,214)">0</span></div>
<div><span>slurmstepd: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> *** STEP
</span><span style="color:rgb(86,156,214)">5034</span><span>.</span><span style="color:rgb(86,156,214)">0</span><span> ON bb01 CANCELLED AT
</span><span style="color:rgb(106,153,85)">2023-01-12T17:21:07</span><span> ***</span></div>
<div><span>srun: Job step aborted: Waiting up to </span><span style="color:rgb(86,156,214)">32</span><span> seconds for job step to finish.</span></div>
<div><span>srun: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> bb01: tasks
</span><span style="color:rgb(86,156,214)">0</span><span>-</span><span style="color:rgb(86,156,214)">1</span><span>: Killed</span></div>
<span><span>srun: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> bb01: tasks
</span><span style="color:rgb(86,156,214)">2</span><span>-</span><span style="color:rgb(86,156,214)">3</span><span>: Exited with exit code
</span><span style="color:rgb(86,156,214)">1</span></span></div>
<br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
The error is raised in the middle of the computation after many successful calls of TAOSolve and TAO iterations. My guess is that TAO computes the preconditioner during its first iteration with all variables in the active set. But the preconditioner is never
updated when some variables are moved to the inactive set during the next TAO iterations. Am I right? Can you help me with that?</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Thanks a lot for your help and your time.</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Regards,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Alexis</div>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr">
<div dir="ltr">
<div>
<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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>