<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div>Hi Barry,</div>
<div><br>
</div>
That is not ideal of course, but at least the error message will clarify the fact that Tao is not designed for repeated solves.
<div><br>
</div>
<div>Thanks,</div>
<div>Blaise</div>
<div><br>
<div><br>
<blockquote type="cite">
<div>On Jan 17, 2023, at 7:07 PM, Barry Smith <bsmith@petsc.dev> wrote:</div>
<br class="Apple-interchange-newline">
<div>
<div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div><br>
</div>
It appears that Tao is not written to allow multiple TaoSolve() on the same Tao object; this is different from KSP, SNES, and TS.
<div><br>
</div>
<div> If you look at the converged reason at the beginning of the second TaoSolve() you will see it is the reason that occurred in the first solve and all the Tao data structures are in the state they were in when the previous TaoSolve ended. Thus it is using
incorrect flags and previous matrices incorrectly.</div>
<div><br>
</div>
<div> Fixing this would be a largish process I think. </div>
<div><br>
</div>
<div> I added an error check for TaoSolve that checks if converged reason is not iterating (meaning the Tao object was previously used and left in a bad state) so that this same problem won't come up for other users. <a href="https://gitlab.com/petsc/petsc/-/merge_requests/5986">https://gitlab.com/petsc/petsc/-/merge_requests/5986</a></div>
<div><br>
</div>
<div> Barry</div>
<div><br>
<div><br>
<blockquote type="cite">
<div>On Jan 16, 2023, at 5:07 PM, Blaise Bourdin <bourdin@mcmaster.ca> wrote:</div>
<br class="Apple-interchange-newline">
<div>
<div>
<div style="line-break:after-white-space">Hi,
<div><br>
</div>
<div>I am attaching a small modification of the eptorsion1.c example that replicates the issue. It looks like this bug is triggered when the upper and lower bounds are equal on enough (10) degrees of freedom.</div>
<div> </div>
<div></div>
</div>
<div style="line-break:after-white-space">
<div></div>
<div><br>
</div>
<div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">---- Elastic-Plastic Torsion Problem -----</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">mx: 10 my: 10 </span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code"; min-height: 16px;">
<span style="font-variant-ligatures:no-common-ligatures"></span><br>
</div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 0</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 1</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 2</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 3</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 4</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 5</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 6</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 7</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 8</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">i: 9</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code"; color: rgb(180, 36, 25);">
<span style="font-variant-ligatures:no-common-ligatures"><b>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</b></span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: Nonconforming object sizes</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: Preconditioner number of local rows 44 does not equal input vector size 54</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: Petsc Development GIT revision: v3.18.2-242-g4615508c7fc GIT Date: 2022-11-28 10:21:46 -0600</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: ./eptorsion1 on a ventura-gcc12.2-arm64-g64 named bblaptop.math.mcmaster.ca by blaise Mon Jan 16 17:06:57 2023</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: Configure options --CFLAGS="-Wimplicit-function-declaration -Wunused" --FFLAGS="-ffree-line-length-none -fallow-argument-mismatch -Wunused" --download-ctetgen=1 --download-exodusii=1 --download-hdf5=1
--download-hypre=1 --download-metis=1 --download-netcdf=1 --download-mumps=1 --download-parmetis=1 --download-pnetcdf=1 --download-scalapack --download-triangle=1 --download-zlib=1 --with-64-bit-indices=1 --with-debugging=1 --with-exodusii-fortran-bindings
--with-shared-libraries=1 --with-x11=0</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #1 PCApply() at /opt/HPC/petsc-main/src/ksp/pc/interface/precon.c:434</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #2 KSP_PCApply() at /opt/HPC/petsc-main/include/petsc/private/kspimpl.h:380</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #3 KSPCGSolve_STCG() at /opt/HPC/petsc-main/src/ksp/ksp/impls/cg/stcg/stcg.c:76</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #4 KSPSolve_Private() at /opt/HPC/petsc-main/src/ksp/ksp/interface/itfunc.c:898</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #5 KSPSolve() at /opt/HPC/petsc-main/src/ksp/ksp/interface/itfunc.c:1070</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #6 TaoBNKComputeStep() at /opt/HPC/petsc-main/src/tao/bound/impls/bnk/bnk.c:459</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #7 TaoSolve_BNTR() at /opt/HPC/petsc-main/src/tao/bound/impls/bnk/bntr.c:138</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #8 TaoSolve() at /opt/HPC/petsc-main/src/tao/interface/taosolver.c:177</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: #9 main() at eptorsion1.c:166</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">[0]PETSC ERROR: No PETSc Option Table entries</span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code"; color: rgb(180, 36, 25);">
<span style="font-variant-ligatures:no-common-ligatures"><b>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</b></span></div>
<div style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: "Fira Code";">
<span style="font-variant-ligatures:no-common-ligatures">Abort(60) on node 0 (rank 0 in comm 16): application called MPI_Abort(MPI_COMM_SELF, 60) - process 0</span></div>
<div><span style="font-variant-ligatures:no-common-ligatures"><br>
</span></div>
<div><span style="font-variant-ligatures:no-common-ligatures">I hope that this helps.</span></div>
<div><span style="font-variant-ligatures:no-common-ligatures"><br>
</span></div>
<div><span style="font-variant-ligatures:no-common-ligatures">Blaise</span></div>
<div><span style="font-variant-ligatures:no-common-ligatures"><br>
</span></div>
<div><br>
<blockquote type="cite">
<div>On Jan 16, 2023, at 3:14 PM, Alexis Marboeuf <alexis.marboeuf@hotmail.fr> wrote:</div>
<br class="x_Apple-interchange-newline">
<div>
<div class="x_elementToProof" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<span class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">Hi Matt,</span></div>
<div class="x_elementToProof" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<span class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">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="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 class="x_ContentPasted0" style="background-color:rgb(255,255,255); display:inline!important">. In my case, in all the successful calls of TAOSolve,<span class="x_ContentPasted1" style="background-color:rgb(255,255,255); display:inline!important"> the
computed trial objective value is better than the current value which implies needH = PETSC_TRUE<span class="x_Apple-converted-space"> </span><span class="x_ContentPasted2" style="background-color:rgb(255,255,255); display:inline!important">within </span><span class="x_ContentPasted1 x_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 class="x_ContentPasted3" style="background-color:rgb(255,255,255); display:inline!important">I have to admit I am struggling
understanding how that boolean needH is computed when BNK is initialized with<span class="x_Apple-converted-space x_ContentPasted3"> </span></span><span class="x_ContentPasted0 x_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="x_appendonsend" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
</div>
<div class="x_elementToProof" style="font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
Thanks a lot.</div>
<div class="x_elementToProof" style="font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
Alexis</div>
<hr tabindex="-1" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; display:inline-block; width:1418.046875px">
<span style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; float:none; display:inline!important"></span>
<div id="x_divRplyFwdMsg" dir="ltr" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>De :</b><span class="x_Apple-converted-space"> </span>Alexis Marboeuf <alexis.marboeuf@hotmail.fr><br>
<b>Envoyé :</b><span class="x_Apple-converted-space"> </span>samedi 14 janvier 2023 05:24<br>
<b>À :</b><span class="x_Apple-converted-space"> </span>Matthew Knepley <knepley@gmail.com><br>
<b>Cc :</b><span class="x_Apple-converted-space"> </span>petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Objet :</b><span class="x_Apple-converted-space"> </span>RE: [petsc-users] Nonconforming object sizes using TAO (TAOBNTR)</font>
<div> </div>
</div>
<div dir="ltr" style="font-family:Helvetica; font-size:12px; font-style:normal; font-variant-caps:normal; font-weight:400; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<div class="x_x_elementToProof"><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">Hi Matt,</span></div>
<div class="x_x_elementToProof"><span class="x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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" id="LPlnkOWALinkPreview">https://github.com/bourdin/mef90</a> (branch
marboeuf/vdef-tao-test)</span></div>
<div class="x_x__Entity x_x__EType_OWALinkPreview x_x__EId_OWALinkPreview x_x__EReadonly_1">
<div id="LPBorder_GTaHR0cHM6Ly9naXRodWIuY29tL2JvdXJkaW4vbWVmOTA." class="x_x_LPBorder354162" style="width:800px; margin-top:16px; margin-bottom:16px; max-width:800px; min-width:424px">
<table id="LPContainer354162" role="presentation" style="padding:12px 36px 12px 12px; width:800px; 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" id="LPImageAnchor354162"><img id="LPThumbnailImageId354162" alt="" height="120" width="240" src="https://opengraph.githubassets.com/8f51eb183957c4e2f2dd59e2733f43a7bc667a50d4aaad934ebb3ac8f25a17ab/bourdin/mef90" style="display:block"></a></div>
</td>
<td style="width:488px">
<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" 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_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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_x_appendonsend"></div>
<div class="x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
Thanks!</div>
<div class="x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt">
Alexis</div>
<hr tabindex="-1" style="display:inline-block; width:1418.046875px">
<div id="x_x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt"><b>De :</b><span class="x_Apple-converted-space"> </span>Matthew Knepley <knepley@gmail.com><br>
<b>Envoyé :</b><span class="x_Apple-converted-space"> </span>samedi 14 janvier 2023 02:38<br>
<b>À :</b><span class="x_Apple-converted-space"> </span>Alexis Marboeuf <alexis.marboeuf@hotmail.fr><br>
<b>Cc :</b><span class="x_Apple-converted-space"> </span>petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov><br>
<b>Objet :</b><span class="x_Apple-converted-space"> </span>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">alexis.marboeuf@hotmail.fr</a>> wrote:<br>
</div>
<div class="x_x_x_gmail_quote">
<blockquote class="x_x_x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div class="x_x_x_msg5004949584998710548">
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">Hi Matt,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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_x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div class="x_x_x_msg5004949584998710548">
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures">Tao Object: (Damage_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: bntr</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Tao Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: bncg</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>CG Type: ssml_bfgs</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Skipped Stepdirection Updates: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scaled gradient steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Pure gradient steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Not a descent direction: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Line search fails: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Matrix has not been preallocated yet</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>TaoLineSearch Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: more-thuente</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>maximum function evaluations=30</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of function evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of gradient evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of function/gradient evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Termination reason: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Active Set subset type: subvec</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>convergence tolerances: gatol=1e-08,<span class="x_Apple-converted-space"> </span><span> <span class="x_Apple-converted-space"> </span></span>steptol=0.,<span class="x_Apple-converted-space"> </span><span>
<span class="x_Apple-converted-space"> </span></span>gttol=0.</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Residual in Function/Gradient:=0.</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Objective value=0.</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of iterations=0,<span> <span class="x_Apple-converted-space"> </span></span>(max: 2000)</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Solver terminated: 0<span class="x_Apple-converted-space"> </span><span> </span></span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Rejected BFGS updates: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>CG steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Newton steps: 11</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>BFGS steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scaled gradient steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Gradient steps: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>KSP termination reasons:</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>atol: 4</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>rtol: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>ctol: 7</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>negc: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>dtol: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>iter: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>othr: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>TaoLineSearch Object: (Damage_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: more-thuente</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>maximum function evaluations=30</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of function evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of gradient evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of function/gradient evaluations=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>using variable bounds</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Termination reason: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>KSP Object: (Damage_tao_bnk_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: stcg</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>maximum iterations=10000, nonzero initial guess</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>tolerances:<span> <span class="x_Apple-converted-space"> </span></span>relative=1e-08, absolute=1e-08, divergence=1e+10</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>left preconditioning</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>using UNPRECONDITIONED norm type for convergence test</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>PC Object: (Damage_tao_bnk_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: lmvm</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Mat Object: (Damage_tao_bnk_pc_lmvm_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: lmvmbfgs</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>rows=30634, cols=30634</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scale type: DIAGONAL</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scale history: 1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scale params: alpha=1., beta=0.5, rho=1.</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Convex factors: phi=0., theta=0.125</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Max. storage: 5</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Used storage: 5</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of updates: 11</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of rejects: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of resets: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Mat Object: (Damage_tao_bnk_pc_lmvm_J0_) 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: lmvmdiagbroyden</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>rows=30634, cols=30634</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scale history: 1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Scale params: alpha=1., beta=0.5, rho=1.</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Convex factor: theta=0.125</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Max. storage: 1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Used storage: 1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of updates: 11</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of rejects: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Number of resets: 0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>linear system matrix = precond matrix:</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Mat Object: 4 MPI processes</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>type: mpiaij</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>rows=468, cols=468</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total: nonzeros=2932, allocated nonzeros=2932</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of mallocs used during MatSetValues calls=0</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>not using I-node (on process 0) routines</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total KSP iterations: 103</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Active Set subset type: subvec</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>convergence tolerances: gatol=0.0001,<span class="x_Apple-converted-space"> </span><span> <span class="x_Apple-converted-space"> </span></span>steptol=0.,<span class="x_Apple-converted-space"> </span><span> <span class="x_Apple-converted-space"> </span></span>gttol=1e-05</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Residual in Function/Gradient:=9.11153e-05</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Objective value=0.00665458</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of iterations=11,<span> <span class="x_Apple-converted-space"> </span></span>(max: 50)</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of function evaluations=17,<span> <span class="x_Apple-converted-space"> </span></span>max: -1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of gradient evaluations=13,<span> <span class="x_Apple-converted-space"> </span></span>max: -1</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>total number of Hessian evaluations=12</span></div>
<div style="margin:0px; font-style:normal; font-variant-caps:normal; font-weight:normal; font-stretch:normal; font-size:11px; line-height:normal; font-family:Menlo">
<span style="font-variant-ligatures:no-common-ligatures"><span> <span class="x_Apple-converted-space"> </span></span>Solution converged:<span> <span class="x_Apple-converted-space"> </span></span>||g(X)|| <= gatol</span></div>
<br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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; background-color:rgb(255,255,255)">Alexis</span></div>
<div id="x_x_x_m_-7317894617606137971appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt"><br>
</div>
<hr style="display:inline-block; width:1406.015625px">
<div id="x_x_x_m_-7317894617606137971divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt"><b>De :</b><span class="x_Apple-converted-space"> </span>Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>><br>
<b>Envoyé :</b><span class="x_Apple-converted-space"> </span>samedi 14 janvier 2023 01:38<br>
<b>À :</b><span class="x_Apple-converted-space"> </span>Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr">alexis.marboeuf@hotmail.fr</a>><br>
<b>Cc :</b><span class="x_Apple-converted-space"> </span><a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><span class="x_Apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
<b>Objet :</b><span class="x_Apple-converted-space"> </span>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">alexis.marboeuf@hotmail.fr</a>> wrote:<br>
</div>
<div>
<blockquote style="margin:0px 0px 0px 0.8ex; border-left-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">Hi all,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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-width:1px; border-left-style:solid; border-left-color:rgb(204,204,204); padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: ---------------------<span class="x_Apple-converted-space"> </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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: Preconditioner number of local rows<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">1122</span><span> does
not equal input vector size<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: See<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/">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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span>
GIT Date:<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: ---------------------<span class="x_Apple-converted-space"> </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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: Preconditioner number of local rows<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">1229</span><span> does
not equal input vector size<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: See<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/">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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span>
GIT Date:<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">60</span><span>) - process<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">60</span><span>) - process<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">0</span></div>
<div><span>slurmstepd:<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">error:</span><span> *** STEP<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(86,156,214)">32</span><span> seconds for job step to finish.</span></div>
<div><span>srun:<span class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">error:</span><span> bb01: tasks<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></span><span style="color:rgb(206,145,120); font-weight:bold">error:</span><span> bb01: tasks<span class="x_Apple-converted-space"> </span></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 class="x_Apple-converted-space"> </span></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; 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; background-color:rgb(255,255,255)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; 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; background-color:rgb(255,255,255)">
Regards,</div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; background-color:rgb(255,255,255)">
Alexis</div>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
--<span class="x_Apple-converted-space"> </span><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/">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>
--<span class="x_Apple-converted-space"> </span><br>
<div dir="ltr" class="x_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/">https://www.cse.buffalo.edu/~knepley/</a></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
<div>
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; overflow-wrap: break-word; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div dir="auto" style="letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none; line-break: after-white-space;">
<div>— <br>
Canada Research Chair in Mathematical and Computational Aspects of Solid Mechanics (Tier 1)</div>
<div>Professor, Department of Mathematics & Statistics<br>
Hamilton Hall room 409A, McMaster University<br>
1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada <br>
https://www.math.mcmaster.ca/bourdin | +1 (905) 525 9140 ext. 27243</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
<span id="cid:63C698A36BE9CF4D9F50B08A1CFC12CE@CANPRD01.PROD.OUTLOOK.COM"><eptorsion1.c></span></div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
<div>
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">
<div>— <br>
Canada Research Chair in Mathematical and Computational Aspects of Solid Mechanics (Tier 1)</div>
<div>Professor, Department of Mathematics & Statistics<br>
Hamilton Hall room 409A, McMaster University<br>
1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada <br>
https://www.math.mcmaster.ca/bourdin | +1 (905) 525 9140 ext. 27243</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</body>
</html>