<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof"><span style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">Hi Matt,</span></div>
<div class="elementToProof"><span style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">After investigation, it fails because, at some point, the boolean needH
is set to PETSC_FALSE when initializing the BNK method with TAOBNKInitialize (line 103 of $PETSC_DIR/src/tao/bound/impls/bnk/bntr.c). The Hessian and the precondtitioner are thus not updated throughout the TAO iterations.</span><span style="color: rgb(0, 0, 0); font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;"> It
has something to do with the option BNK_INIT_INTERPOLATION set by default. It works when I choose BNK_INIT_CONSTANT<span style="background-color:rgb(255, 255, 255);display:inline !important" class="ContentPasted0">. In my case, in all the successful calls
of TAOSolve,<span style="background-color:rgb(255, 255, 255);display:inline !important" class="ContentPasted1"> the computed trial objective value is better than the current value which implies needH = PETSC_TRUE
<span style="background-color:rgb(255, 255, 255);display:inline !important" class="ContentPasted2">
within </span><span class="ContentPasted1 ContentPasted2" style="margin:0px;background-color:rgb(255, 255, 255);display:inline !important">TAOBNKInitialize</span>. At some point, the trial value becomes equal to the current objective value up to machine precision
and then, needH = PETSC_FALSE. <span style="background-color:rgb(255, 255, 255);display:inline !important" class="ContentPasted3">I have to admit I am struggling understanding how that boolean needH is computed when BNK is initialized with<span class="Apple-converted-space ContentPasted3"> </span></span><span class="ContentPasted0 ContentPasted3" style="margin:0px;background-color:rgb(255, 255, 255);display:inline !important">BNK_INIT_INTERPOLATION.
Can you help me with that?</span></span></span></span></div>
<div id="appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)" class="elementToProof">
Thanks a lot.</div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)" class="elementToProof">
Alexis</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>De :</b> Alexis Marboeuf <alexis.marboeuf@hotmail.fr><br>
<b>Envoyé :</b> samedi 14 janvier 2023 05:24<br>
<b>À :</b> Matthew Knepley <knepley@gmail.com><br>
<b>Cc :</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Objet :</b> RE: [petsc-users] Nonconforming object sizes using TAO (TAOBNTR)</font>
<div> </div>
</div>
<div dir="ltr">
<div class="x_elementToProof"><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 class="x_elementToProof"><span class="x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">Indeed, it fails on 1 process with the same error. The source code is available
here: <a href="https://github.com/bourdin/mef90" data-auth="NotApplicable" id="LPlnkOWALinkPreview">https://github.com/bourdin/mef90</a> (branch marboeuf/vdef-tao-test)</span></div>
<div class="x__Entity x__EType_OWALinkPreview x__EId_OWALinkPreview x__EReadonly_1">
<div id="LPBorder_GTaHR0cHM6Ly9naXRodWIuY29tL2JvdXJkaW4vbWVmOTA." class="x_LPBorder354162" style="width:100%; margin-top:16px; margin-bottom:16px; max-width:800px; min-width:424px">
<table id="LPContainer354162" role="presentation" style="padding:12px 36px 12px 12px; width:100%; border:1px solid rgb(200,200,200); border-radius:2px">
<tbody>
<tr valign="top" style="border-spacing:0px">
<td>
<div id="LPImageContainer354162" style="margin-right:12px; height:120px; overflow:hidden; width:240px">
<a href="https://github.com/bourdin/mef90" data-auth="NotApplicable" id="LPImageAnchor354162"><img id="LPThumbnailImageId354162" alt="" height="120" width="240" style="display:block" src="https://opengraph.githubassets.com/8f51eb183957c4e2f2dd59e2733f43a7bc667a50d4aaad934ebb3ac8f25a17ab/bourdin/mef90"></a></div>
</td>
<td style="width:100%">
<div id="LPTitle354162" style="font-size:21px; font-weight:300; margin-right:8px; font-family:wf_segoe-ui_light,"Segoe UI Light","Segoe WP Light","Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; margin-bottom:12px">
<a href="https://github.com/bourdin/mef90" data-auth="NotApplicable" id="LPUrlAnchor354162" style="text-decoration:none">GitHub - bourdin/mef90: Official repository for mef90/vDef</a></div>
<div id="LPDescription354162" style="font-size:14px; max-height:100px; color:rgb(102,102,102); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif; margin-bottom:12px; margin-right:8px; overflow:hidden">
mef90 / vDef: A reference implementation of the variational approach to fracture, as described in: Francfort, G. and Marigo, J.-J. (1998). Revisiting brittle fracture as an energy minimization problem.</div>
<div id="LPMetadata354162" style="font-size:14px; font-weight:400; color:rgb(166,166,166); font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif">
github.com</div>
</td>
</tr>
</tbody>
</table>
</div>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0); background-color:rgb(255,255,255)">
I can share the details (installation + command line) for running it. But the ideal would be to reproduce this error with a Petsc example so it's easier for you to investigate. I looked for a bound constraint minimization problem with TAO and TS but I didn't
find it. What example could I use? </div>
<div id="x_appendonsend"></div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Thanks!</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Alexis</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>De :</b> Matthew Knepley <knepley@gmail.com><br>
<b>Envoyé :</b> samedi 14 janvier 2023 02:38<br>
<b>À :</b> Alexis Marboeuf <alexis.marboeuf@hotmail.fr><br>
<b>Cc :</b> petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><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 3:21 PM Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr" data-auth="NotApplicable">alexis.marboeuf@hotmail.fr</a>> wrote:<br>
</div>
<div class="x_x_gmail_quote">
<blockquote class="x_x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_x_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="x_x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_x_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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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-top: 0px; margin-bottom: 0px;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="x_x_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="x_x_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" data-auth="NotApplicable">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" data-auth="NotApplicable">alexis.marboeuf@hotmail.fr</a>><br>
<b>Cc :</b> <a href="mailto:petsc-users@mcs.anl.gov" data-auth="NotApplicable">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" data-auth="NotApplicable">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" data-auth="NotApplicable">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/" data-auth="NotApplicable">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/" data-auth="NotApplicable">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/" data-auth="NotApplicable">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="x_x_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/" data-auth="NotApplicable">https://www.cse.buffalo.edu/~knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>