<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);">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);"><br>
</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);">Here is the output from Petsc when I view the TAO solver:</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);"><br>
</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);">
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0">Tao Object: (Damage_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: bntr</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Tao Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: bncg</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>CG Type: ssml_bfgs</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Skipped Stepdirection Updates: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scaled gradient steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Pure gradient steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Not a descent direction: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Line search fails: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Matrix has not been preallocated yet</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>TaoLineSearch Object: (Damage_tao_bnk_cg_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: more-thuente</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>maximum function evaluations=30</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of function evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of gradient evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of function/gradient evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Termination reason: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Active Set subset type: subvec</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>convergence tolerances: gatol=1e-08, <span class="Apple-converted-space ContentPasted0">
</span>steptol=0., <span class="Apple-converted-space ContentPasted0">
</span>gttol=0.</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Residual in Function/Gradient:=0.</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Objective value=0.</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of iterations=0,<span class="Apple-converted-space ContentPasted0">
</span>(max: 2000)</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Solver terminated: 0 <span class="Apple-converted-space ContentPasted0"> </span></span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Rejected BFGS updates: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>CG steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Newton steps: 11</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>BFGS steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scaled gradient steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Gradient steps: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>KSP termination reasons:</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>atol: 4</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>rtol: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>ctol: 7</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>negc: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>dtol: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>iter: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>othr: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>TaoLineSearch Object: (Damage_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: more-thuente</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>maximum function evaluations=30</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of function evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of gradient evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of function/gradient evaluations=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>using variable bounds</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Termination reason: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>KSP Object: (Damage_tao_bnk_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: stcg</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>maximum iterations=10000, nonzero initial guess</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>tolerances:<span class="Apple-converted-space ContentPasted0"> </span>relative=1e-08, absolute=1e-08, divergence=1e+10</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>left preconditioning</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>using UNPRECONDITIONED norm type for convergence test</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>PC Object: (Damage_tao_bnk_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: lmvm</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Mat Object: (Damage_tao_bnk_pc_lmvm_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: lmvmbfgs</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>rows=30634, cols=30634</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scale type: DIAGONAL</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scale history: 1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scale params: alpha=1., beta=0.5, rho=1.</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Convex factors: phi=0., theta=0.125</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Max. storage: 5</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Used storage: 5</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of updates: 11</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of rejects: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of resets: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Mat Object: (Damage_tao_bnk_pc_lmvm_J0_) 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: lmvmdiagbroyden</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>rows=30634, cols=30634</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scale history: 1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Scale params: alpha=1., beta=0.5, rho=1.</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Convex factor: theta=0.125</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Max. storage: 1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Used storage: 1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of updates: 11</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of rejects: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Number of resets: 0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>linear system matrix = precond matrix:</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Mat Object: 4 MPI processes</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>type: mpiaij</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>rows=468, cols=468</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total: nonzeros=2932, allocated nonzeros=2932</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of mallocs used during MatSetValues calls=0</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>not using I-node (on process 0) routines</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total KSP iterations: 103</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Active Set subset type: subvec</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>convergence tolerances: gatol=0.0001, <span class="Apple-converted-space ContentPasted0">
</span>steptol=0., <span class="Apple-converted-space ContentPasted0"> </span>
gttol=1e-05</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Residual in Function/Gradient:=9.11153e-05</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Objective value=0.00665458</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of iterations=11,<span class="Apple-converted-space ContentPasted0">
</span>(max: 50)</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of function evaluations=17,<span class="Apple-converted-space ContentPasted0">
</span>max: -1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of gradient evaluations=13,<span class="Apple-converted-space ContentPasted0">
</span>max: -1</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>total number of Hessian evaluations=12</span></p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo;color:#000000"><span style="font-variant-ligatures:no-common-ligatures" class="ContentPasted0"><span class="Apple-converted-space ContentPasted0">
</span>Solution converged:<span class="Apple-converted-space ContentPasted0"> </span>
||g(X)|| <= gatol</span></p>
<br>
</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);">Thanks again for your help!</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);">Alexis</span></div>
<div id="appendonsend"></div>
<div style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</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> Matthew Knepley <knepley@gmail.com><br>
<b>Envoyé :</b> samedi 14 janvier 2023 01: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 11:22 AM Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr" data-auth="NotApplicable">alexis.marboeuf@hotmail.fr</a>> wrote:<br>
</div>
<div class="x_gmail_quote">
<blockquote class="x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_msg5485660062342652653">
<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 class="x_gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_msg5485660062342652653">
<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" class="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>
</body>
</html>