[petsc-users] PF+Navier stokes

Tue Mar 23 15:20:30 CDT 2021

On Tue, Mar 23, 2021 at 4:08 PM Sepideh Kavousi <skavou1 at lsu.edu> wrote:

> Thank you for you suggestions. In the discretization, I made sure that the
> main diagonal of the matrix is non-zero and I even chose Dirichlet boundary
> condition for pressure and velocity. And therefore, I do not get zero pivot.
> But even when setting preconditioning to be LU, the resid norm is still
> large:
>

It does not sound like this problem is well-posed. If these linear systems
do not have solutions, how do you know
your problem has a solution?

  Thanks,

    Matt

> When I analyze the data, the velocity at regions close to the corners of
> solid phase becomes large. I am not sure how I could get reasonable results
> for one system size and totally nonsense for another.
> Sepideh
> 29 TS dt 3.2769e-06 time 0.00689069
>       0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00
> copy!
> Write output at step= 29!
>       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 1.453900186218e+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_2_ solve converged due to CONVERGED_ITS iterations
> 1
>       2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>     1 SNES Function norm 1.298945530538e+08
>       1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
> 29 TS dt 3.2769e-06 time 0.00689069
> copy!
> Write output at step= 29!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations
> 1
>     0 SNES Function norm 1.453900186218e+08
>       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
>       2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations
> 1
>     1 SNES Function norm 1.298945530538e+08
> 29 TS dt 3.2769e-06 time 0.00689069
> copy!
>       0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00
> Write output at step= 29!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations
> 1
>     0 SNES Function norm 1.453900186218e+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_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 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00
>       1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
>       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.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
>       2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>       Linear fieldsplit_1_ solve converged due to CONVERGED_ITS iterations
> 1
>     1 SNES Function norm 1.298945530538e+08
> 29 TS dt 3.2769e-06 time 0.00689069
>       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= 29!
>       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 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00
>     0 SNES Function norm 1.453900186218e+08
>       2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations
> 1
>     1 SNES Function norm 1.298945530538e+08
> 29 TS dt 3.2769e-06 time 0.00689069
>       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= 29!
>       1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
>     0 SNES Function norm 1.453900186218e+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 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||r(i)||/||b|| 1.000000000000e+00
>       2 KSP preconditioned resid norm 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>     1 SNES Function norm 1.298945530538e+08
> 29 TS dt 3.2769e-06 time 0.00689069
> copy!
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations
> 1
>       Linear fieldsplit_0_ solve converged due to CONVERGED_ITS iterations
> 1
> Write output at step= 29!
>       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 1.453900186218e+08
>       Linear fieldsplit_2_ solve converged due to CONVERGED_ITS iterations
> 1
>       1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
>       0 KSP preconditioned resid norm 5.978616012547e+03 true resid norm
> 3.032434597885e+08 ||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_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 4.099453531745e-02 true resid norm
> 2.719463080396e+06 ||r(i)||/||b|| 8.967919975234e-03
>       1 KSP preconditioned resid norm 1.760119119843e-01 true resid norm
> 2.196480894736e+06 ||r(i)||/||b|| 7.243291895784e-03
>     1 SNES Function norm 1.298945530538e+08
>
> ------------------------------
> *From:* Barry Smith <bsmith at petsc.dev>
> *Sent:* Monday, March 22, 2021 10:48 PM
> *To:* Sepideh Kavousi <skavou1 at lsu.edu>
> *Cc:* petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
> *Subject:* Re: [petsc-users] PF+Navier stokes
>
>
>   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://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mcs.anl.gov%2Fpipermail%2Fpetsc-users%2F2015-February%2F024191.html&data=04%7C01%7Cskavou1%40lsu.edu%7Cd69125e191c8496977f008d8edae8d9d%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C1%7C637520681280339374%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=946tv7G7HkMd8Rtit8bZdBIMyftuT2nqNB4aNx2o%2FPk%3D&reserved=0>
> )
> -fieldsplit_1_fieldsplit_1_pc_type ml (based on
> https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024191.html
> <https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.mcs.anl.gov%2Fpipermail%2Fpetsc-users%2F2015-February%2F024191.html&data=04%7C01%7Cskavou1%40lsu.edu%7Cd69125e191c8496977f008d8edae8d9d%7C2d4dad3f50ae47d983a09ae2b1f466f8%7C0%7C1%7C637520681280339374%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=946tv7G7HkMd8Rtit8bZdBIMyftuT2nqNB4aNx2o%2FPk%3D&reserved=0>
> )
> -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>
> *Sent:* Monday, March 22, 2021 1:56 PM
> *To:* Sepideh Kavousi <skavou1 at lsu.edu>
> *Cc:* petsc-users at mcs.anl.gov <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> 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
>
>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210323/48afe841/attachment-0001.html>