<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div>Hi Colin,</div><div dir="ltr"><br></div><div dir="ltr"><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 19 Mar 2019 at 09:33, Cotter, Colin J <<a href="mailto:colin.cotter@imperial.ac.uk">colin.cotter@imperial.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div id="gmail-m_8482798322617626713divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Hi Dave,</p>
<div style="color:rgb(0,0,0)">
<div>
<div>
<div class="gmail-m_8482798322617626713x_gmail_quote">
<div dir="auto"><br>
</div>
<div dir="auto">>If you are doing that, then you need to tell fieldsplit to use the Amat to define the splits otherwise it will define the Schur compliment as </div>
<div dir="auto">>S = B22 - B21 inv(B11) B12</div>
<div dir="auto">>preconditiones with B22, where as what you want is </div>
<div dir="auto">>S = A22 - A21 inv(A11) A12</div>
<div dir="auto">>preconditioned with B22.</div>
<div dir="auto"><br>
</div>
<div dir="auto">>If your operators are set up this way and you didn't indicate to use Amat to define S this would definitely explain why preonly works but iterating on Schur does not.</div>
<div dir="auto"><br>
</div>
<div dir="auto">Yes, thanks - this solves it! I need pc_use_amat.</div></div></div></div></div></div></div></blockquote><div><br></div><div>Okay great. But doesn't that option eradicate your custom Schur complement object which you inserted into the Bmat in the (2,2) slot?</div><div><br></div><div>I thought you would use the option</div><div>-pc_fieldsplit_diag_use_amat<br></div><div><br></div><div>In general for fieldsplit (Schur) I found that the best way to manage user defined Schur complement preconditioners is via PCFieldSplitSetSchurPre().</div><div><br></div><div><a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCFieldSplitSetSchurPre.html#PCFieldSplitSetSchurPre">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCFieldSplitSetSchurPre.html#PCFieldSplitSetSchurPre</a><br></div><div><br></div><div>Also, for solver debugging purposes with fieldsplit and MatNest, I find it incredibly useful to attach textual names to all the matrices going to into FieldSplit. You can use PetscObjectSetName() with each of your sub-matrices in the Amat and the Bmat, and any schur complement operators. The textual names will be displayed in KSP view. In that way you have a better chance of understanding which operators are being used where. (Note that this trick is less useful with the Amat and Bmat are AIJ matrices).</div><div><br></div><div>Below is an example KSPView associated with 2x2 block system where I've attached the names Auu,Aup,Apu,App, and S* to the Amat sub-matices and the schur complement preconditioner.</div><div><br></div><div>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">PC Object:(dcy_) 1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: fieldsplit</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>FieldSplit with Schur preconditioner, factorization FULL</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Preconditioner for the Schur complement formed from Sp, an assembled approximation to S, which uses (lumped, if requested) A00's diagonal's inverse</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Split info:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Split number 0 Defined by IS</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Split number 1 Defined by IS</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP solver for A00 block</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: preonly</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>maximum iterations=10000, initial guess is zero</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerances:<span class="gmail-Apple-converted-space"> </span>relative=1e-05, absolute=1e-50, divergence=10000.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>left preconditioning</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using NONE norm type for convergence test</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>PC Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: lu</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>LU: out-of-place factorization</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerance for zero pivot 2.22045e-14</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>matrix ordering: nd</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>factor fill ratio given 0., needed 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Factored matrix follows:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object: <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=85728, cols=85728</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>package used to perform factorization: umfpack</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=0, allocated nonzeros=0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>UMFPACK run parameters:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PRL]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_STRATEGY]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_COL]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_ROW]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AMD_DENSE]: 10.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_BLOCK_SIZE]: 32.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_FIXQ]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AGGRESSIVE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PIVOT_TOLERANCE]: 0.1</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SYM_PIVOT_TOLERANCE]: 0.001</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SCALE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ALLOC_INIT]: 0.7</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DROPTOL]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_IRSTEP]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ORDERING]: AMD (not using the PETSc ordering)</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>linear system matrix = precond matrix:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>Auu<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=85728, cols=85728</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=1028736, allocated nonzeros=1028736</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using I-node routines: found 21432 nodes, limit used is 5</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP solver for S = A11 - A10 inv(A00) A01<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_p_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: fgmres</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>GMRES: happy breakdown tolerance 1e-30</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>maximum iterations=300, initial guess is zero</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerances:<span class="gmail-Apple-converted-space"> </span>relative=0.01, absolute=1e-50, divergence=10000.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>right preconditioning</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using UNPRECONDITIONED norm type for convergence test</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>PC Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_p_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: lu</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>LU: out-of-place factorization</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerance for zero pivot 2.22045e-14</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>matrix ordering: nd</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>factor fill ratio given 0., needed 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Factored matrix follows:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object: <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=32542, cols=32542</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>package used to perform factorization: umfpack</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=0, allocated nonzeros=0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>UMFPACK run parameters:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PRL]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_STRATEGY]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_COL]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_ROW]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AMD_DENSE]: 10.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_BLOCK_SIZE]: 32.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_FIXQ]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AGGRESSIVE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PIVOT_TOLERANCE]: 0.1</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SYM_PIVOT_TOLERANCE]: 0.001</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SCALE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ALLOC_INIT]: 0.7</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DROPTOL]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_IRSTEP]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ORDERING]: AMD (not using the PETSc ordering)</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>linear system matrix followed by preconditioner matrix:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_p_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: schurcomplement</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=32542, cols=32542</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Schur complement A11 - A10 inv(A00) A01</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>A11</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>App<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_p_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=32542, cols=32542</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=548482, allocated nonzeros=548482</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>A10</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>Apu <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=32542, cols=85728</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=857280, allocated nonzeros=921990</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP of A00</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>KSP Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: preonly</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>maximum iterations=10000, initial guess is zero</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerances:<span class="gmail-Apple-converted-space"> </span>relative=1e-05, absolute=1e-50, divergence=10000.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>left preconditioning</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using NONE norm type for convergence test</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>PC Object:<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: lu</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>LU: out-of-place factorization</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>tolerance for zero pivot 2.22045e-14</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>matrix ordering: nd</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>factor fill ratio given 0., needed 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Factored matrix follows:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object: <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=85728, cols=85728</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>package used to perform factorization: umfpack</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=0, allocated nonzeros=0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>UMFPACK run parameters:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PRL]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_STRATEGY]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_COL]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DENSE_ROW]: 0.2</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AMD_DENSE]: 10.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_BLOCK_SIZE]: 32.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_FIXQ]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_AGGRESSIVE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_PIVOT_TOLERANCE]: 0.1</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SYM_PIVOT_TOLERANCE]: 0.001</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_SCALE]: 1.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ALLOC_INIT]: 0.7</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_DROPTOL]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_IRSTEP]: 0.</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Control[UMFPACK_ORDERING]: AMD (not using the PETSc ordering)</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>linear system matrix = precond matrix:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>Auu<span class="gmail-Apple-converted-space"> </span>(dcy_fieldsplit_u_) <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=85728, cols=85728</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=1028736, allocated nonzeros=1028736</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using I-node routines: found 21432 nodes, limit used is 5</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>A01</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object:<span class="gmail-Apple-converted-space"> </span>Aup <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=85728, cols=32542</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=857280, allocated nonzeros=942634</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>using I-node routines: found 21432 nodes, limit used is 5</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span></span><span style="font-variant-ligatures:no-common-ligatures">Mat Object:</span><span class="gmail-Apple-converted-space" style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures">S* </span><span class="gmail-Apple-converted-space" style="font-variant-ligatures:no-common-ligatures"> </span><span style="font-variant-ligatures:no-common-ligatures">1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: seqaij</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=32542, cols=32542</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total: nonzeros=548482, allocated nonzeros=548482</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls =0</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>not using I-node routines</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>linear system matrix = precond matrix:</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Mat Object: <span class="gmail-Apple-converted-space"> </span>1 MPI processes</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type: nest</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>rows=118270, cols=118270</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>has attached null space</span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Matrix object:<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>type=nest, rows=2, cols=2<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>MatNest structure:<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>(0,0) : name="Auu", prefix="dcy_fieldsplit_u_", type=seqaij, rows=85728, cols=85728<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>(0,1) : name="Aup", type=seqaij, rows=85728, cols=32542<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>(1,0) : name="Apu", type=seqaij, rows=32542, cols=85728<span class="gmail-Apple-converted-space"> </span></span></p>
<p class="gmail-p1" style="margin:0px;font:10px Monaco;color:rgb(245,245,245);background-color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>(1,1) : name="App", prefix="dcy_fieldsplit_p_", type=seqaij, rows=32542, cols=32542<span class="gmail-Apple-converted-space"> </span></span></p></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div id="gmail-m_8482798322617626713divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols" dir="ltr"><div style="color:rgb(0,0,0)"><div><div><div class="gmail-m_8482798322617626713x_gmail_quote">
<div dir="auto"><br>
</div>
<div dir="auto">all the best</div>
<div dir="auto">--Colin</div>
<blockquote class="gmail-m_8482798322617626713x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div class="gmail-m_8482798322617626713x_gmail_quote">
<blockquote class="gmail-m_8482798322617626713x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div id="gmail-m_8482798322617626713x_m_-7508259247097921405m_5252267629230022956divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p></p>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote></div></div></div></div></div></div></div>