<div dir="ltr"><div dir="ltr">On Tue, Mar 23, 2021 at 4:08 PM Sepideh Kavousi <<a href="mailto:skavou1@lsu.edu">skavou1@lsu.edu</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">



<div>
<div dir="ltr">
<div>

</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you for you suggestions. In the discretization, I made sure that the main diagonal of the matrix is non-zero and I even chose Dirichlet boundary condition for pressure and velocity. And therefore, I do not get zero pivot.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
But even when setting preconditioning to be LU, the resid norm is still large:  </div></div></div></blockquote><div><br></div><div>It does not sound like this problem is well-posed. If these linear systems do not have solutions, how do you know</div><div>your problem has a solution?</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)" dir="ltr">
When I analyze the data, the velocity at regions close to the corners of solid phase becomes large. I am not sure how I could get reasonable results for one system size and totally nonsense for another.</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)" dir="ltr">
Sepideh</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
29 TS dt 3.2769e-06 time 0.00689069
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>copy!</div>
<div>Write output at step= 29!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>29 TS dt 3.2769e-06 time 0.00689069</div>
<div>copy!</div>
<div>Write output at step= 29!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<div>29 TS dt 3.2769e-06 time 0.00689069</div>
<div>copy!</div>
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>Write output at step= 29!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<div>29 TS dt 3.2769e-06 time 0.00689069</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 29!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<div>29 TS dt 3.2769e-06 time 0.00689069</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 29!</div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<div>29 TS dt 3.2769e-06 time 0.00689069</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 29!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 1.453900186218e+08 </div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>      0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03</div>
<div>      1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03</div>
<div>    1 SNES Function norm 1.298945530538e+08 </div>
<br>
</div>
<div id="gmail-m_2892199464017788175appendonsend" dir="ltr"></div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_2892199464017788175divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>><br>
<b>Sent:</b> Monday, March 22, 2021 10:48 PM<br>
<b>To:</b> Sepideh Kavousi <<a href="mailto:skavou1@lsu.edu" target="_blank">skavou1@lsu.edu</a>><br>
<b>Cc:</b> <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> Re: [petsc-users] PF+Navier stokes</font>
<div> </div>
</div>
<div style="overflow-wrap: break-word;">
<div><br>
</div>
  Each of the individual fields "converging" does not imply the overall PCFIELDSPLIT will converge.
<div><br>
</div>
<div>  This is scary,   <span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07</span><font face="Calibri, Arial, Helvetica, sans-serif" size="3"> often
 these huge differences  in true residual come from the operator having a null space that is not handled properly. Or they can come from an ILU that produces absurd pivots.</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"><br>
</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3">   Recommend using the same PCFIELDSPLIT but using a direct solver on all the splits (no ML no ILU). If that results in great overall convergence it means one of the sub-systems
 is not converging properly with ML or ILU and you need better preconditioners on some of the splits.</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"><br>
</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3">  Barry</font></div>
<div><font face="Calibri, Arial, Helvetica, sans-serif" size="3"><br>
</font>
<div><br>
<div><br>
<blockquote type="cite">
<div>On Mar 22, 2021, at 4:05 PM, Sepideh Kavousi <<a href="mailto:skavou1@lsu.edu" target="_blank">skavou1@lsu.edu</a>> wrote:</div>
<br>
<div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I modified my BC such that on the left and right side of the interface the BC are constant value instead of Neumann(zero flux). This solves the problem but still the code has convergence problem:</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I even tried field split with the following order:<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
the block size is 9. the first block is the fields related to for PF equation, the second split block is the velocities in x and y direction and the third block is pressure.<span> </span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
-pc_type fieldsplit</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
-pc_fieldsplit_block_size 9<span> </span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
-pc_fieldsplit_0_fields 0,1  <span> </span><span style="background-color:rgb(255,255,0)">(two fields related to for the Phasefield model)<span> </span></span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
-pc_fieldsplit_1_fields 2,3 <span> </span><span style="background-color:rgb(255,255,0)">(velocity in x and y direction)</span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
-pc_fieldsplit_2_fields 4<span> </span><span style="background-color:rgb(255,255,0)">(pressure)</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255)">-fieldsplit_1_pc_fieldsplit_block_size 2<span> </span><br>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255)">-fieldsplit_1_fieldsplit_0_pc_type ml (based on<span> </span><a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mcs.anl.gov%2Fpipermail%2Fpetsc-users%2F2015-February%2F024191.html&data=04%7C01%7Cskavou1%40lsu.edu%7Cd69125e191c8496977f008d8edae8d9d%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C1%7C637520681280339374%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=946tv7G7HkMd8Rtit8bZdBIMyftuT2nqNB4aNx2o%2FPk%3D&reserved=0" id="gmail-m_2892199464017788175LPlnk191545" target="_blank">https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html</a>)<br>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255)">-fieldsplit_1_fieldsplit_1_pc_type ml<span> </span></span><span style="background-color:rgb(255,255,0)">(based on<span> </span></span><span style="background-color:rgb(255,255,255)"><a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mcs.anl.gov%2Fpipermail%2Fpetsc-users%2F2015-February%2F024191.html&data=04%7C01%7Cskavou1%40lsu.edu%7Cd69125e191c8496977f008d8edae8d9d%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C1%7C637520681280339374%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=946tv7G7HkMd8Rtit8bZdBIMyftuT2nqNB4aNx2o%2FPk%3D&reserved=0" id="gmail-m_2892199464017788175LPlnk876933" target="_blank"><span style="background-color:rgb(255,255,0)">https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html</span></a></span><span style="background-color:rgb(255,255,0)">)</span><span style="background-color:rgb(255,255,255)"><br>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255)">-fieldsplit_0_pc_type ilu<span> </span></span><span style="background-color:rgb(255,255,0)">(based on previous solutions of phase-field equations)</span><span style="background-color:rgb(255,255,255)"><br>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255)">-fieldsplit_2_pc_type ilu</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
<span style="background-color:rgb(255,255,0)"></span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I guess changing the BCs the main reason that at first few steps the code does not fail. And as time increases, true resid norm increases such that at a finite time step (~30) it reaches 1e7 and the code results non-accurate velocity calculations. Can this
 also be resulted by forward/backward discritization?</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Best,</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Sepideh<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
34 TS dt 3.12462e-07 time 0.00709097
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>copy!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 34!</div>
<div>Write output at step= 34!</div>
<div>      7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.904272255718e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>33 TS dt 1.56231e-07 time 0.00709082</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>    1 SNES Function norm 1.904272255718e+07<span> </span></div>
<div>Write output at step= 33!</div>
<div>33 TS dt 1.56231e-07 time 0.00709082</div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>copy!</div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>Write output at step= 33!</div>
<div>    0 SNES Function norm 2.464456787205e+07<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.464456787205e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      0 KSP preconditioned resid norm 2.003930340911e+01 true resid norm 6.987120567963e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      1 KSP preconditioned resid norm 1.199890501875e-02 true resid norm 1.879731143354e+07 ||r(i)||/||b|| 2.690280101896e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      2 KSP preconditioned resid norm 3.018100763012e-04 true resid norm 1.879893603977e+07 ||r(i)||/||b|| 2.690512616309e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      3 KSP preconditioned resid norm 2.835332741838e-04 true resid norm 1.879893794065e+07 ||r(i)||/||b|| 2.690512888363e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      4 KSP preconditioned resid norm 1.860011376508e-04 true resid norm 1.879893735946e+07 ||r(i)||/||b|| 2.690512805182e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.886737547995e+07<span> </span></div>
<div>31 TS dt 3.90578e-08 time 0.0070907</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 31!</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>    0 SNES Function norm 1.888557765431e+07<span> </span></div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>copy!</div>
<div>copy!</div>
<div>Write output at step= 34!</div>
<div>Write output at step= 34!</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>copy!</div>
<div>copy!</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 35!</div>
<div>Write output at step= 35!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 34!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.338586074554e-01 true resid norm 1.888557765431e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>Write output at step= 34!</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 5.451259747447e-03 true resid norm 1.887927947148e+07 ||r(i)||/||b|| 9.996665083300e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 34!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>copy!</div>
<div>      2 KSP preconditioned resid norm 9.554577345960e-04 true resid norm 1.887930135577e+07 ||r(i)||/||b|| 9.996676671129e-01</div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>copy!</div>
<div>Write output at step= 35!</div>
<div>Write output at step= 35!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>      3 KSP preconditioned resid norm 9.378991224281e-04 true resid norm 1.887930134907e+07 ||r(i)||/||b|| 9.996676667583e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 3.652611805745e-04 true resid norm 1.887930205974e+07 ||r(i)||/||b|| 9.996677043885e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 2.918222127367e-04 true resid norm 1.887930204569e+07 ||r(i)||/||b|| 9.996677036447e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 6.114488674627e-05 true resid norm 1.887930243837e+07 ||r(i)||/||b|| 9.996677244370e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 3.763532951474e-05 true resid norm 1.887930248279e+07 ||r(i)||/||b|| 9.996677267895e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 2.112644035802e-05 true resid norm 1.887930251181e+07 ||r(i)||/||b|| 9.996677283257e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 1.113068460252e-05 true resid norm 1.887930250969e+07 ||r(i)||/||b|| 9.996677282137e-01</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     10 KSP preconditioned resid norm 1.352518287887e-06 true resid norm 1.887930250333e+07 ||r(i)||/||b|| 9.996677278767e-01</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>copy!</div>
<div>copy!</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 36!</div>
<div>Write output at step= 36!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     11 KSP preconditioned resid norm 7.434707372444e-07 true resid norm 1.887930250410e+07 ||r(i)||/||b|| 9.996677279175e-01</div>
<div>    1 SNES Function norm 1.887938190335e+07<span> </span></div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>32 TS dt 7.81156e-08 time 0.00709074</div>
<div>copy!</div>
<div>Write output at step= 32!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.010558777785e+07<span> </span></div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>copy!</div>
<div>Write output at step= 35!</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>copy!</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>copy!</div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>copy!</div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>Write output at step= 35!</div>
<div>Write output at step= 36!</div>
<div>copy!</div>
<div>Write output at step= 36!</div>
<div>Write output at step= 35!</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>copy!</div>
<div>copy!</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>Write output at step= 37!</div>
<div>      0 KSP preconditioned resid norm 3.165458473526e+00 true resid norm 2.010558777785e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>Write output at step= 37!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      1 KSP preconditioned resid norm 3.655364946441e-03 true resid norm 1.904269379864e+07 ||r(i)||/||b|| 9.471343991057e-01</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 36!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>      2 KSP preconditioned resid norm 2.207564350060e-03 true resid norm 1.904265942845e+07 ||r(i)||/||b|| 9.471326896210e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 2.188447918524e-03 true resid norm 1.904266151317e+07 ||r(i)||/||b|| 9.471327933098e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 7.425314556150e-04 true resid norm 1.904265807404e+07 ||r(i)||/||b|| 9.471326222560e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 6.052794097111e-04 true resid norm 1.904265841692e+07 ||r(i)||/||b|| 9.471326393103e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 1.251197915617e-04 true resid norm 1.904266159346e+07 ||r(i)||/||b|| 9.471327973028e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>Write output at step= 36!</div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>copy!</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>    1 SNES Function norm 1.904272255718e+07<span> </span></div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>Write output at step= 36!</div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>33 TS dt 1.56231e-07 time 0.00709082</div>
<div>copy!</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>Write output at step= 33!</div>
<div>copy!</div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>copy!</div>
<div>Write output at step= 37!</div>
<div>Write output at step= 37!</div>
<div>    0 SNES Function norm 2.464456787205e+07<span> </span></div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>copy!</div>
<div>copy!</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>Write output at step= 38!</div>
<div>Write output at step= 38!</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>copy!</div>
<div>Write output at step= 37!</div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 38!</div>
<div>Write output at step= 38!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>copy!</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>Write output at step= 37!</div>
<div>      1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>copy!</div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>Write output at step= 37!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>      2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 1.938116032575e+07<span> </span></div>
<div>34 TS dt 3.12462e-07 time 0.00709097</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 34!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>    0 SNES Function norm 6.393295863037e+07<span> </span></div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>copy!</div>
<div>copy!</div>
<div>Write output at step= 39!</div>
<div>Write output at step= 39!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>copy!</div>
<div>Write output at step= 38!</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>copy!</div>
<div>Write output at step= 39!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 39!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>      1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03</div>
<div>copy!</div>
<div>Write output at step= 38!</div>
<div>    1 SNES Function norm 1.679071798524e+08<span> </span></div>
<div>38 TS dt 2.09926e-06 time 0.00709464</div>
<div>copy!</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>Write output at step= 38!</div>
<div>    0 SNES Function norm 4.969343279719e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>      1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>copy!</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 40!</div>
<div>Write output at step= 40!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 2.011100207442e+07<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>35 TS dt 6.24925e-07 time 0.00709129</div>
<div>copy!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 35!</div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>    0 SNES Function norm 2.790215258123e+08<span> </span></div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>copy!</div>
<div>Write output at step= 39!</div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>      0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>copy!</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 40!</div>
<div>Write output at step= 40!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>      1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03</div>
<div>copy!</div>
<div>Write output at step= 39!</div>
<div>    1 SNES Function norm 1.426574249724e+08<span> </span></div>
<div>39 TS dt 3.42747e-06 time 0.00709674</div>
<div>copy!</div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>Write output at step= 39!</div>
<div>    0 SNES Function norm 6.081085806316e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02</div>
<div>      1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02</div>
<div>      2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02</div>
<div>      3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>      1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>copy!</div>
<div>      4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02</div>
<div>      4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02</div>
<div>Write output at step= 40!</div>
<div>    1 SNES Function norm 3.485261902880e+07<span> </span></div>
<div>36 TS dt 1.24985e-06 time 0.00709191</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 36!</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 4.684646000717e+08<span> </span></div>
<div>      5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02</div>
<div>      5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02</div>
<div>      6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02</div>
<div>      7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02</div>
<div>      8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02</div>
<div>      9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02</div>
<div>     10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02</div>
<div>     11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02</div>
<div>      1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02</div>
<div>     12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02</div>
<div>     12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 6.245293405655e+10<span> </span></div>
<div>    1 SNES Function norm 6.245293405655e+10<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.097192603535e+10<span> </span></div>
<div>    0 SNES Function norm 2.097192603535e+10<span> </span></div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>copy!</div>
<div>      3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02</div>
<div>      3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02</div>
<div>      1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02</div>
<div>Write output at step= 40!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    1 SNES Function norm 1.906687815261e+08<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>40 TS dt 6.85494e-06 time 0.00710017</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>copy!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 40!</div>
<div>      4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02</div>
<div>      4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02</div>
<div>    0 SNES Function norm 2.386070312612e+09<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02</div>
<div>      5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02</div>
<div>      6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02</div>
<div>      7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02</div>
<div>      8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02</div>
<div>      9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02</div>
<div>     11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>     12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02</div>
<div>     12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02</div>
<div>      0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02</div>
<div>    1 SNES Function norm 6.245293405655e+10<span> </span></div>
<div>    1 SNES Function norm 6.245293405655e+10<span> </span></div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.097192603535e+10<span> </span></div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 2.097192603535e+10<span> </span></div>
<div>      1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02</div>
<div>      3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02</div>
<div>      4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02</div>
<div>    1 SNES Function norm 5.194818889106e+07<span> </span></div>
<div>37 TS dt 1.4811e-06 time 0.00709316</div>
<div>copy!</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>Write output at step= 37!</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>    0 SNES Function norm 3.786081856480e+09<span> </span></div>
<div>      5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 4.670669894784e+05 true resid norm 2.097192603535e+10 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      0 KSP preconditioned resid norm 4.670669894784e+05 true resid norm 2.097192603535e+10 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      1 KSP preconditioned resid norm 2.135766802417e+00 true resid norm 8.990660342789e+07 ||r(i)||/||b|| 4.286997926482e-03</div>
<div>      1 KSP preconditioned resid norm 2.135766802417e+00 true resid norm 8.990660342789e+07 ||r(i)||/||b|| 4.286997926482e-03</div>
<div>      Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1</div>
<div>      0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00</div>
<div>    1 SNES Function norm 7.093938024894e+08<span> </span></div>
<div>    1 SNES Function norm 7.093938024894e+08<span> </span></div>
<br>
</div>
<div id="gmail-m_2892199464017788175x_appendonsend" style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
</div>
<hr style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;display:inline-block;width:1200.5px">
<span style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;float:none;display:inline"></span>
<div id="gmail-m_2892199464017788175x_divRplyFwdMsg" dir="ltr" style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<font face="Calibri, sans-serif" style="font-size:11pt"><b>From:</b><span> </span>Barry Smith <<a href="mailto:bsmith@petsc.dev" target="_blank">bsmith@petsc.dev</a>><br>
<b>Sent:</b><span> </span>Monday, March 22, 2021 1:56 PM<br>
<b>To:</b><span> </span>Sepideh Kavousi <<a href="mailto:skavou1@lsu.edu" target="_blank">skavou1@lsu.edu</a>><br>
<b>Cc:</b><span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b><span> </span>Re: [petsc-users] PF+Navier stokes</font>
<div> </div>
</div>
<div style="font-family:Helvetica;font-size:18px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none">
<div><br>
</div>
   Singular systems come up in solving PDEs almost always due to issues related to boundary conditions. For example all Neumann (natural) boundary conditions can produce singular systems. Direct factorizations generically will eventually hit a zero pivot in
 such cases and there is no universally acceptable approaches for what to do at that point to recover. If you think your operator is singular you should start by using MatSetNullSpace(), it won't "cure" the problem but is the tool we use to manage null spaces
 in operators.
<div><br>
</div>
<div><br>
<div><br>
</div>
<div><br>
<div><br>
<blockquote type="cite">
<div>On Mar 22, 2021, at 9:04 AM, Sepideh Kavousi <<a href="mailto:skavou1@lsu.edu" target="_blank">skavou1@lsu.edu</a>> wrote:</div>
<br>
<div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Hello,</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I want to solve PF solidification+Navier stokes using Finite different method, and I have a strange problem. My code runs fine for some system sizes and fails for some of the system sizes. When I run with the following options:</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(0,255,0)">mpirun -np 2 ./one.out -ts_monitor -snes_fd_color -ts_max_snes_failures -1  -ts_type bdf -ts_bdf_adapt -pc_type bjacobi  -snes_linesearch_type l2 -snes_type ksponly -ksp_type gmres -ksp_gmres_restart 1001
 -sub_pc_type ilu -sub_ksp_type preonly -snes_monitor -ksp_monitor -snes_linesearch_monitor -ksp_monitor_true_residual -ksp_converged_reason -log_view</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,128)"><br>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,255,255);color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span><span style="background-color:rgb(255,255,128)">
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">^C    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(12,136,42)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(12,136,42)">                   PC_FAILED due to SUBPC_ERROR</span></div>
<div><span style="color:rgb(12,136,42)">    0 SNES Function norm 1.465357113711e+01</span></div>
</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="background-color:rgb(255,128,255)">Even setting pc_type to LU does not solve the problem.</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="color:rgb(237,92,87)">0 TS dt 0.0001 time 0.</span>
<div><span style="color:rgb(237,92,87)">copy!</span></div>
<div><span style="color:rgb(237,92,87)">copy!</span></div>
<div><span style="color:rgb(237,92,87)">Write output at step= 0!</span></div>
<div><span style="color:rgb(237,92,87)">Write output at step= 0!</span></div>
<div><span style="color:rgb(237,92,87)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(237,92,87)">    0 SNES Function norm 1.465357113711e+01</span></div>
<div><span style="color:rgb(237,92,87)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<div><span style="color:rgb(237,92,87)">                   PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT</span></div>
<div><span style="color:rgb(237,92,87)">    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0</span></div>
<span style="color:rgb(237,92,87)">                   PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT</span></div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
I guess the problem is that in mass conservation I used forward discretization for u (velocity in x) and for the moment in x , I used forward discretization for p (pressure) to ensure non-zero terms on the diagonal of matrix. I tried to run it with valgrind
 but it did not output anything. <span> </span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Does anyone have suggestions on how to solve this issue?<span> </span><br>
</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Best,</div>
<div style="font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">
Sepideh</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>

</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>