[petsc-users] PF+Navier stokes

Sepideh Kavousi skavou1 at lsu.edu
Tue Mar 23 15:07:56 CDT 2021


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:
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<mailto: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<mailto:bsmith at petsc.dev>>
Sent: Monday, March 22, 2021 1:56 PM
To: Sepideh Kavousi <skavou1 at lsu.edu<mailto:skavou1 at lsu.edu>>
Cc: petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov> <petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov>>
Subject: Re: [petsc-users] PF+Navier stokes


   Singular systems come up in solving PDEs almost always due to issues related to boundary conditions. For example all Neumann (natural) boundary conditions can produce singular systems. Direct factorizations generically will eventually hit a zero pivot in such cases and there is no universally acceptable approaches for what to do at that point to recover. If you think your operator is singular you should start by using MatSetNullSpace(), it won't "cure" the problem but is the tool we use to manage null spaces in operators.





On Mar 22, 2021, at 9:04 AM, Sepideh Kavousi <skavou1 at lsu.edu<mailto:skavou1 at lsu.edu>> wrote:

Hello,
I want to solve PF solidification+Navier stokes using Finite different method, and I have a strange problem. My code runs fine for some system sizes and fails for some of the system sizes. When I run with the following options:
mpirun -np 2 ./one.out -ts_monitor -snes_fd_color -ts_max_snes_failures -1  -ts_type bdf -ts_bdf_adapt -pc_type bjacobi  -snes_linesearch_type l2 -snes_type ksponly -ksp_type gmres -ksp_gmres_restart 1001 -sub_pc_type ilu -sub_ksp_type preonly -snes_monitor -ksp_monitor -snes_linesearch_monitor -ksp_monitor_true_residual -ksp_converged_reason -log_view

    0 SNES Function norm 1.465357113711e+01
    0 SNES Function norm 1.465357113711e+01
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    0 SNES Function norm 1.465357113711e+01
    0 SNES Function norm 1.465357113711e+01
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    0 SNES Function norm 1.465357113711e+01
    0 SNES Function norm 1.465357113711e+01
^C    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    0 SNES Function norm 1.465357113711e+01
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to SUBPC_ERROR
    0 SNES Function norm 1.465357113711e+01

Even setting pc_type to LU does not solve the problem.
0 TS dt 0.0001 time 0.
copy!
copy!
Write output at step= 0!
Write output at step= 0!
    0 SNES Function norm 1.465357113711e+01
    0 SNES Function norm 1.465357113711e+01
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC_FAILED due to FACTOR_NUMERIC_ZEROPIVOT

I guess the problem is that in mass conservation I used forward discretization for u (velocity in x) and for the moment in x , I used forward discretization for p (pressure) to ensure non-zero terms on the diagonal of matrix. I tried to run it with valgrind but it did not output anything.

Does anyone have suggestions on how to solve this issue?
Best,
Sepideh

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


More information about the petsc-users mailing list