[petsc-users] PF+Navier stokes

Mon Mar 22 22:48:39 CDT 2021

  Each of the individual fields "converging" does not imply the overall PCFIELDSPLIT will converge.

  This is scary,   7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 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.

   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.

  Barry

> On Mar 22, 2021, at 4:05 PM, Sepideh Kavousi <skavou1 at lsu.edu> wrote:
> 
> 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:
> I even tried field split with the following order:
> 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. 
> 
> -pc_type fieldsplit
> -pc_fieldsplit_block_size 9 
> -pc_fieldsplit_0_fields 0,1   (two fields related to for the Phasefield model) 
> -pc_fieldsplit_1_fields 2,3  (velocity in x and y direction)
> -pc_fieldsplit_2_fields 4 (pressure)
> -fieldsplit_1_pc_fieldsplit_block_size 2 
> -fieldsplit_1_fieldsplit_0_pc_type ml (based on https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html <https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html>)
> -fieldsplit_1_fieldsplit_1_pc_type ml (based on https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html <https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html>)
> -fieldsplit_0_pc_type ilu (based on previous solutions of phase-field equations)
> -fieldsplit_2_pc_type ilu
> 
> 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?
> Best,
> Sepideh
> 
> 
> 34 TS dt 3.12462e-07 time 0.00709097
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> copy!
> copy!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 34!
> Write output at step= 34!
>       7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01
>     0 SNES Function norm 6.393295863037e+07 
>     0 SNES Function norm 6.393295863037e+07 
>       7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.904272255718e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> 33 TS dt 1.56231e-07 time 0.00709082
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>     1 SNES Function norm 1.904272255718e+07 
> Write output at step= 33!
> 33 TS dt 1.56231e-07 time 0.00709082
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
> copy!
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
> Write output at step= 33!
>     0 SNES Function norm 2.464456787205e+07 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.464456787205e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       0 KSP preconditioned resid norm 2.003930340911e+01 true resid norm 6.987120567963e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       1 KSP preconditioned resid norm 1.199890501875e-02 true resid norm 1.879731143354e+07 ||r(i)||/||b|| 2.690280101896e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       2 KSP preconditioned resid norm 3.018100763012e-04 true resid norm 1.879893603977e+07 ||r(i)||/||b|| 2.690512616309e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       3 KSP preconditioned resid norm 2.835332741838e-04 true resid norm 1.879893794065e+07 ||r(i)||/||b|| 2.690512888363e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       4 KSP preconditioned resid norm 1.860011376508e-04 true resid norm 1.879893735946e+07 ||r(i)||/||b|| 2.690512805182e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.886737547995e+07 
> 31 TS dt 3.90578e-08 time 0.0070907
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 31!
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>     0 SNES Function norm 1.888557765431e+07 
>     1 SNES Function norm 1.938116032575e+07 
>     1 SNES Function norm 1.938116032575e+07 
> 34 TS dt 3.12462e-07 time 0.00709097
> 34 TS dt 3.12462e-07 time 0.00709097
> copy!
> copy!
> Write output at step= 34!
> Write output at step= 34!
>     0 SNES Function norm 6.393295863037e+07 
>     0 SNES Function norm 6.393295863037e+07 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 2.011100207442e+07 
>     1 SNES Function norm 2.011100207442e+07 
> 35 TS dt 6.24925e-07 time 0.00709129
> 35 TS dt 6.24925e-07 time 0.00709129
> copy!
> copy!
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 35!
> Write output at step= 35!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.790215258123e+08 
>     0 SNES Function norm 2.790215258123e+08 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 1.938116032575e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> 34 TS dt 3.12462e-07 time 0.00709097
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 34!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
>     0 SNES Function norm 6.393295863037e+07 
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.338586074554e-01 true resid norm 1.888557765431e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.938116032575e+07 
> 34 TS dt 3.12462e-07 time 0.00709097
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>     1 SNES Function norm 1.938116032575e+07 
> 34 TS dt 3.12462e-07 time 0.00709097
> Write output at step= 34!
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 5.451259747447e-03 true resid norm 1.887927947148e+07 ||r(i)||/||b|| 9.996665083300e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 34!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 6.393295863037e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 6.393295863037e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>     1 SNES Function norm 2.011100207442e+07 
> 35 TS dt 6.24925e-07 time 0.00709129
>     1 SNES Function norm 2.011100207442e+07 
> copy!
>       2 KSP preconditioned resid norm 9.554577345960e-04 true resid norm 1.887930135577e+07 ||r(i)||/||b|| 9.996676671129e-01
> 35 TS dt 6.24925e-07 time 0.00709129
> copy!
> Write output at step= 35!
> Write output at step= 35!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.790215258123e+08 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.790215258123e+08 
>       3 KSP preconditioned resid norm 9.378991224281e-04 true resid norm 1.887930134907e+07 ||r(i)||/||b|| 9.996676667583e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 3.652611805745e-04 true resid norm 1.887930205974e+07 ||r(i)||/||b|| 9.996677043885e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 2.918222127367e-04 true resid norm 1.887930204569e+07 ||r(i)||/||b|| 9.996677036447e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 6.114488674627e-05 true resid norm 1.887930243837e+07 ||r(i)||/||b|| 9.996677244370e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 3.763532951474e-05 true resid norm 1.887930248279e+07 ||r(i)||/||b|| 9.996677267895e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 2.112644035802e-05 true resid norm 1.887930251181e+07 ||r(i)||/||b|| 9.996677283257e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 1.113068460252e-05 true resid norm 1.887930250969e+07 ||r(i)||/||b|| 9.996677282137e-01
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>      10 KSP preconditioned resid norm 1.352518287887e-06 true resid norm 1.887930250333e+07 ||r(i)||/||b|| 9.996677278767e-01
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 3.485261902880e+07 
>     1 SNES Function norm 3.485261902880e+07 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> 36 TS dt 1.24985e-06 time 0.00709191
> 36 TS dt 1.24985e-06 time 0.00709191
> copy!
> copy!
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 36!
> Write output at step= 36!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 4.684646000717e+08 
>     0 SNES Function norm 4.684646000717e+08 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>      11 KSP preconditioned resid norm 7.434707372444e-07 true resid norm 1.887930250410e+07 ||r(i)||/||b|| 9.996677279175e-01
>     1 SNES Function norm 1.887938190335e+07 
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
> 32 TS dt 7.81156e-08 time 0.00709074
> copy!
> Write output at step= 32!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.010558777785e+07 
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>     1 SNES Function norm 2.011100207442e+07 
> 35 TS dt 6.24925e-07 time 0.00709129
> copy!
> Write output at step= 35!
>     0 SNES Function norm 2.790215258123e+08 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>     1 SNES Function norm 2.011100207442e+07 
>     1 SNES Function norm 3.485261902880e+07 
> 36 TS dt 1.24985e-06 time 0.00709191
> 35 TS dt 6.24925e-07 time 0.00709129
> copy!
>     1 SNES Function norm 3.485261902880e+07 
> copy!
> 36 TS dt 1.24985e-06 time 0.00709191
>     1 SNES Function norm 2.011100207442e+07 
> copy!
> 35 TS dt 6.24925e-07 time 0.00709129
> Write output at step= 35!
> Write output at step= 36!
> copy!
> Write output at step= 36!
> Write output at step= 35!
>     0 SNES Function norm 2.790215258123e+08 
>     0 SNES Function norm 4.684646000717e+08 
>     0 SNES Function norm 4.684646000717e+08 
>     0 SNES Function norm 2.790215258123e+08 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 5.194818889106e+07 
>     1 SNES Function norm 5.194818889106e+07 
> 37 TS dt 1.4811e-06 time 0.00709316
> 37 TS dt 1.4811e-06 time 0.00709316
> copy!
> copy!
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
> Write output at step= 37!
>       0 KSP preconditioned resid norm 3.165458473526e+00 true resid norm 2.010558777785e+07 ||r(i)||/||b|| 1.000000000000e+00
> Write output at step= 37!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 3.786081856480e+09 
>     0 SNES Function norm 3.786081856480e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       1 KSP preconditioned resid norm 3.655364946441e-03 true resid norm 1.904269379864e+07 ||r(i)||/||b|| 9.471343991057e-01
>     1 SNES Function norm 3.485261902880e+07 
> 36 TS dt 1.24985e-06 time 0.00709191
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 36!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 4.684646000717e+08 
>       2 KSP preconditioned resid norm 2.207564350060e-03 true resid norm 1.904265942845e+07 ||r(i)||/||b|| 9.471326896210e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 2.188447918524e-03 true resid norm 1.904266151317e+07 ||r(i)||/||b|| 9.471327933098e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 7.425314556150e-04 true resid norm 1.904265807404e+07 ||r(i)||/||b|| 9.471326222560e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 6.052794097111e-04 true resid norm 1.904265841692e+07 ||r(i)||/||b|| 9.471326393103e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 1.251197915617e-04 true resid norm 1.904266159346e+07 ||r(i)||/||b|| 9.471327973028e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 7.989463672336e-05 true resid norm 1.904266194607e+07 ||r(i)||/||b|| 9.471328148409e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 4.768163844493e-05 true resid norm 1.904266226755e+07 ||r(i)||/||b|| 9.471328308303e-01
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
>     1 SNES Function norm 3.485261902880e+07 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> 36 TS dt 1.24985e-06 time 0.00709191
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 2.073486308871e-05 true resid norm 1.904266231364e+07 ||r(i)||/||b|| 9.471328331228e-01
>     1 SNES Function norm 3.485261902880e+07 
> Write output at step= 36!
> 36 TS dt 1.24985e-06 time 0.00709191
> copy!
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>     1 SNES Function norm 1.904272255718e+07 
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
> Write output at step= 36!
>     0 SNES Function norm 4.684646000717e+08 
> 33 TS dt 1.56231e-07 time 0.00709082
> copy!
>     1 SNES Function norm 5.194818889106e+07 
> 37 TS dt 1.4811e-06 time 0.00709316
>     1 SNES Function norm 5.194818889106e+07 
> Write output at step= 33!
> copy!
>     0 SNES Function norm 4.684646000717e+08 
> 37 TS dt 1.4811e-06 time 0.00709316
> copy!
> Write output at step= 37!
> Write output at step= 37!
>     0 SNES Function norm 2.464456787205e+07 
>     0 SNES Function norm 3.786081856480e+09 
>     0 SNES Function norm 3.786081856480e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.679071798524e+08 
>     1 SNES Function norm 1.679071798524e+08 
> 38 TS dt 2.09926e-06 time 0.00709464
> 38 TS dt 2.09926e-06 time 0.00709464
> copy!
> copy!
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
> Write output at step= 38!
> Write output at step= 38!
>     0 SNES Function norm 4.969343279719e+09 
>     0 SNES Function norm 4.969343279719e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>     1 SNES Function norm 5.194818889106e+07 
> 37 TS dt 1.4811e-06 time 0.00709316
> copy!
> Write output at step= 37!
>     0 SNES Function norm 3.786081856480e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.502321597276e+01 true resid norm 2.464456787205e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.679071798524e+08 
> 38 TS dt 2.09926e-06 time 0.00709464
>     1 SNES Function norm 1.679071798524e+08 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> copy!
> 38 TS dt 2.09926e-06 time 0.00709464
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 38!
> Write output at step= 38!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 5.194818889106e+07 
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
> 37 TS dt 1.4811e-06 time 0.00709316
>     0 SNES Function norm 4.969343279719e+09 
> copy!
>     0 SNES Function norm 4.969343279719e+09 
> Write output at step= 37!
>       1 KSP preconditioned resid norm 8.755672227974e-03 true resid norm 1.938151775153e+07 ||r(i)||/||b|| 7.864417770341e-01
>     1 SNES Function norm 5.194818889106e+07 
> 37 TS dt 1.4811e-06 time 0.00709316
> copy!
>     0 SNES Function norm 3.786081856480e+09 
> Write output at step= 37!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 3.786081856480e+09 
>       2 KSP preconditioned resid norm 1.508190120513e-03 true resid norm 1.938016275793e+07 ||r(i)||/||b|| 7.863867956035e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 1.422130457598e-03 true resid norm 1.938011909055e+07 ||r(i)||/||b|| 7.863850237166e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       4 KSP preconditioned resid norm 4.957954730602e-04 true resid norm 1.938028047800e+07 ||r(i)||/||b|| 7.863915723180e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 3.552719096462e-04 true resid norm 1.938029700510e+07 ||r(i)||/||b|| 7.863922429363e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 1.508455316659e-04 true resid norm 1.938039640595e+07 ||r(i)||/||b|| 7.863962763140e-01
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 7.631525161839e-05 true resid norm 1.938041339557e+07 ||r(i)||/||b|| 7.863969657001e-01
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 1.938116032575e+07 
> 34 TS dt 3.12462e-07 time 0.00709097
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 34!
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>     0 SNES Function norm 6.393295863037e+07 
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>     1 SNES Function norm 1.426574249724e+08 
>     1 SNES Function norm 1.426574249724e+08 
> 39 TS dt 3.42747e-06 time 0.00709674
> 39 TS dt 3.42747e-06 time 0.00709674
> copy!
> copy!
> Write output at step= 39!
> Write output at step= 39!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 6.081085806316e+09 
>     0 SNES Function norm 6.081085806316e+09 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>     1 SNES Function norm 1.679071798524e+08 
> 38 TS dt 2.09926e-06 time 0.00709464
> copy!
> Write output at step= 38!
>     0 SNES Function norm 4.969343279719e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.626907080947e+04 true resid norm 3.786081856480e+09 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 1.426574249724e+08 
> 39 TS dt 3.42747e-06 time 0.00709674
>     1 SNES Function norm 1.426574249724e+08 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> 39 TS dt 3.42747e-06 time 0.00709674
> copy!
> Write output at step= 39!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 39!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 6.081085806316e+09 
>     0 SNES Function norm 6.081085806316e+09 
>     1 SNES Function norm 1.679071798524e+08 
> 38 TS dt 2.09926e-06 time 0.00709464
>       1 KSP preconditioned resid norm 1.411143979706e-01 true resid norm 3.252255320579e+07 ||r(i)||/||b|| 8.590029069268e-03
> copy!
> Write output at step= 38!
>     1 SNES Function norm 1.679071798524e+08 
> 38 TS dt 2.09926e-06 time 0.00709464
> copy!
>     0 SNES Function norm 4.969343279719e+09 
> Write output at step= 38!
>     0 SNES Function norm 4.969343279719e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.212273553655e+02 true resid norm 6.393295863037e+07 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>     1 SNES Function norm 1.906687815261e+08 
>     1 SNES Function norm 1.906687815261e+08 
>       1 KSP preconditioned resid norm 2.328823680294e-02 true resid norm 2.016516515602e+07 ||r(i)||/||b|| 3.154111054458e-01
> 40 TS dt 6.85494e-06 time 0.00710017
> 40 TS dt 6.85494e-06 time 0.00710017
> copy!
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 40!
> Write output at step= 40!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.386070312612e+09 
>     0 SNES Function norm 2.386070312612e+09 
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 2.757087112828e-04 true resid norm 2.010380358203e+07 ||r(i)||/||b|| 3.144513254621e-01
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 2.011100207442e+07 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> 35 TS dt 6.24925e-07 time 0.00709129
> copy!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 35!
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>     0 SNES Function norm 2.790215258123e+08 
>     1 SNES Function norm 1.426574249724e+08 
> 39 TS dt 3.42747e-06 time 0.00709674
> copy!
> Write output at step= 39!
>     0 SNES Function norm 6.081085806316e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>     1 SNES Function norm 1.906687815261e+08 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> 40 TS dt 6.85494e-06 time 0.00710017
>     1 SNES Function norm 1.906687815261e+08 
>       0 KSP preconditioned resid norm 8.865880769163e+04 true resid norm 4.969343279719e+09 ||r(i)||/||b|| 1.000000000000e+00
> copy!
> 40 TS dt 6.85494e-06 time 0.00710017
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 40!
> Write output at step= 40!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.386070312612e+09 
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.386070312612e+09 
>     1 SNES Function norm 1.426574249724e+08 
> 39 TS dt 3.42747e-06 time 0.00709674
>       1 KSP preconditioned resid norm 4.636177252582e-01 true resid norm 4.472959395106e+07 ||r(i)||/||b|| 9.001107678275e-03
> copy!
> Write output at step= 39!
>     1 SNES Function norm 1.426574249724e+08 
> 39 TS dt 3.42747e-06 time 0.00709674
> copy!
>     0 SNES Function norm 6.081085806316e+09 
> Write output at step= 39!
>     0 SNES Function norm 6.081085806316e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02
>       1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02
>       2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.143120146053e+03 true resid norm 2.790215258123e+08 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02
>       3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.906687815261e+08 
>       1 KSP preconditioned resid norm 1.086405760770e-02 true resid norm 2.165319343001e+07 ||r(i)||/||b|| 7.760402487575e-02
> 40 TS dt 6.85494e-06 time 0.00710017
> copy!
>       4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02
>       4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02
> Write output at step= 40!
>     1 SNES Function norm 3.485261902880e+07 
> 36 TS dt 1.24985e-06 time 0.00709191
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.386070312612e+09 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 36!
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 4.684646000717e+08 
>       5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02
>       5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02
>       6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02
>       7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02
>       8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02
>       9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>      10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02
>      10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>      11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02
>      11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02
>       1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02
>      12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02
>      12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.932453778489e+05 true resid norm 6.081085806316e+09 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 6.245293405655e+10 
>     1 SNES Function norm 6.245293405655e+10 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.097192603535e+10 
>     0 SNES Function norm 2.097192603535e+10 
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.906687815261e+08 
> 40 TS dt 6.85494e-06 time 0.00710017
> copy!
>       3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02
>       3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02
>       1 KSP preconditioned resid norm 1.434243067909e+00 true resid norm 7.171981611988e+07 ||r(i)||/||b|| 1.179391615316e-02
> Write output at step= 40!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     1 SNES Function norm 1.906687815261e+08 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.386070312612e+09 
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
> 40 TS dt 6.85494e-06 time 0.00710017
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> copy!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 40!
>       4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02
>       4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02
>     0 SNES Function norm 2.386070312612e+09 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02
>       5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02
>       6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02
>       7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02
>       8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02
>       9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>      10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>      10 KSP preconditioned resid norm 1.779870885679e+00 true resid norm 4.269989357798e+07 ||r(i)||/||b|| 1.789548839038e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>      11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02
>      11 KSP preconditioned resid norm 1.701919459227e+00 true resid norm 4.102844905926e+07 ||r(i)||/||b|| 1.719498744124e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 5.030869192505e+00 true resid norm 1.519850942875e+08 ||r(i)||/||b|| 6.369682128989e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>      12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02
>      12 KSP preconditioned resid norm 1.633671715659e+00 true resid norm 3.959508749664e+07 ||r(i)||/||b|| 1.659426685264e-02
>       0 KSP preconditioned resid norm 4.028400414901e+03 true resid norm 4.684646000717e+08 ||r(i)||/||b|| 1.000000000000e+00
>       2 KSP preconditioned resid norm 3.859648344440e+00 true resid norm 9.578310406856e+07 ||r(i)||/||b|| 4.014261589958e-02
>     1 SNES Function norm 6.245293405655e+10 
>     1 SNES Function norm 6.245293405655e+10 
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.097192603535e+10 
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 2.097192603535e+10 
>       1 KSP preconditioned resid norm 7.513114811690e-02 true resid norm 2.515581080270e+07 ||r(i)||/||b|| 5.369842416878e-02
>       3 KSP preconditioned resid norm 3.133252984070e+00 true resid norm 7.649343929054e+07 ||r(i)||/||b|| 3.205833410953e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       2 KSP preconditioned resid norm 3.361950980934e-02 true resid norm 2.589885319795e+07 ||r(i)||/||b|| 5.528454699457e-02
>       4 KSP preconditioned resid norm 2.753196206679e+00 true resid norm 6.468496966917e+07 ||r(i)||/||b|| 2.710941472566e-02
>     1 SNES Function norm 5.194818889106e+07 
> 37 TS dt 1.4811e-06 time 0.00709316
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
> Write output at step= 37!
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>     0 SNES Function norm 3.786081856480e+09 
>       5 KSP preconditioned resid norm 2.500754637621e+00 true resid norm 5.655851544362e+07 ||r(i)||/||b|| 2.370362480296e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       6 KSP preconditioned resid norm 2.263885927226e+00 true resid norm 5.368792420993e+07 ||r(i)||/||b|| 2.250056250487e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       7 KSP preconditioned resid norm 2.104083158879e+00 true resid norm 5.002745241660e+07 ||r(i)||/||b|| 2.096646194883e-02
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       8 KSP preconditioned resid norm 1.976241956139e+00 true resid norm 4.707835526096e+07 ||r(i)||/||b|| 1.973049788689e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 4.670669894784e+05 true resid norm 2.097192603535e+10 ||r(i)||/||b|| 1.000000000000e+00
>       0 KSP preconditioned resid norm 4.670669894784e+05 true resid norm 2.097192603535e+10 ||r(i)||/||b|| 1.000000000000e+00
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       9 KSP preconditioned resid norm 1.870078692623e+00 true resid norm 4.468270947214e+07 ||r(i)||/||b|| 1.872648481311e-02
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       1 KSP preconditioned resid norm 2.135766802417e+00 true resid norm 8.990660342789e+07 ||r(i)||/||b|| 4.286997926482e-03
>       1 KSP preconditioned resid norm 2.135766802417e+00 true resid norm 8.990660342789e+07 ||r(i)||/||b|| 4.286997926482e-03
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations 1
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations 1
>       0 KSP preconditioned resid norm 1.698771295747e+05 true resid norm 2.386070312612e+09 ||r(i)||/||b|| 1.000000000000e+00
>     1 SNES Function norm 7.093938024894e+08 
>     1 SNES Function norm 7.093938024894e+08 
> 
> From: Barry Smith <bsmith at petsc.dev <mailto:bsmith at petsc.dev>>
> Sent: Monday, March 22, 2021 1:56 PM
> To: Sepideh Kavousi <skavou1 at lsu.edu <mailto:skavou1 at lsu.edu>>
> Cc: petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov> <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov>>
> Subject: Re: [petsc-users] PF+Navier stokes
>  
> 
>    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.
> 
> 
> 
> 
> 
>> On Mar 22, 2021, at 9:04 AM, Sepideh Kavousi <skavou1 at lsu.edu <mailto:skavou1 at lsu.edu>> wrote:
>> 
>> Hello,
>> 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:
>> 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
>> 
>>     0 SNES Function norm 1.465357113711e+01
>>     0 SNES Function norm 1.465357113711e+01
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     0 SNES Function norm 1.465357113711e+01
>>     0 SNES Function norm 1.465357113711e+01
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     0 SNES Function norm 1.465357113711e+01
>>     0 SNES Function norm 1.465357113711e+01
>> ^C    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     0 SNES Function norm 1.465357113711e+01
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to SUBPC_ERROR
>>     0 SNES Function norm 1.465357113711e+01
>> 
>> Even setting pc_type to LU does not solve the problem.
>> 0 TS dt 0.0001 time 0.
>> copy!
>> copy!
>> Write output at step= 0!
>> Write output at step= 0!
>>     0 SNES Function norm 1.465357113711e+01
>>     0 SNES Function norm 1.465357113711e+01
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT
>>     Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
>>                    PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT
>> 
>> 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.  
>> 
>> Does anyone have suggestions on how to solve this issue? 
>> Best,
>> Sepideh

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210322/6af9d199/attachment-0001.html>