Thank you Hong!
I've used GMRES via
mpirun \
-n ${NP} pflotran \
-pflotranin ${INPUTFILE}.pflinput \
-flow_ksp_type gmres \
-flow_pc_type bjacobi \
-flow_sub_pc_type lu \
-flow_sub_pc_factor_nonzeros_along_diagonal \
-snes_monitor
and get:
NP 1
FLOW TS BE steps = 43 newton = 43 linear = 43 cuts
= 0
FLOW TS BE Wasted Linear Iterations = 0
FLOW TS BE SNES time = 197.0 seconds
NP 2
FLOW TS BE steps = 43 newton = 43 linear = 770 cuts
= 0
FLOW TS BE Wasted Linear Iterations = 0
FLOW TS BE SNES time = 68.7 seconds
Which looks ok to me.
Robert
On 07/07/17 15:49, h...@aspiritech.org wrote:
> What do you get with '-ksp_type gmres' or '-ksp_type bcgs' in parallel
> runs?
> Hong
>
> On Fri, Jul 7, 2017 at 6:05 AM, Robert Annewandter
> <robert.annewand...@opengosim.com
> <mailto:robert.annewand...@opengosim.com>> wrote:
>
> Yes indeed, PFLOTRAN cuts timestep after 8 failed iterations of SNES.
>
> I've rerun with -snes_monitor (attached with canonical suffix),
> their -pc_type is always PCBJACOBI + PCLU (though we'd like to try
> SUPERLU in the future, however it works only with -mat_type aij..)
>
>
> The sequential and parallel runs I did with
>
> -ksp_type preonly -pc_type lu
> -pc_factor_nonzeros_along_diagonal -snes_monitor
>
> and
>
> -ksp_type preonly -pc_type bjacobi -sub_pc_type lu
> -sub_pc_factor_nonzeros_along_diagonal -snes_monitor
>
> As expected, the sequential are bot identical and the parallel
> takes half the time compared to sequential.
>
>
>
>
> On 07/07/17 01:20, Barry Smith wrote:
>> Looks like PFLOTRAN has a maximum number of SNES iterations as 8 and
>> cuts the timestep if that fails.
>>
>> Please run with -snes_monitor I don't understand the strange densely
>> packed information that PFLOTRAN is printing.
>>
>> It looks like the linear solver is converging fine in parallel,
>> normally then there is absolutely no reason that the Newton should behave
>> different on 2 processors than 1 unless there is something wrong with the
>> Jacobian. What is the -pc_type for the two cases LU or your fancy thing?
>>
>> Please run sequential and parallel with -pc_type lu and also with
>> -snes_monitor. We need to fix all the knobs but one in order to understand
>> what is going on.
>>
>>
>> Barry
>>
>>
>>
>>> On Jul 6, 2017, at 5:11 PM, Robert Annewandter
>>> <robert.annewand...@opengosim.com>
>>> <mailto:robert.annewand...@opengosim.com> wrote:
>>>
>>> Thanks Barry!
>>>
>>> I've attached log files for np = 1 (SNES time: 218 s) and np = 2 (SNES
>>> time: 600 s). PFLOTRAN final output:
>>>
>>> NP 1
>>>
>>> FLOW TS BE steps = 43 newton = 43 linear = 43 cuts =
>>> 0
>>> FLOW TS BE Wasted Linear Iterations = 0
>>> FLOW TS BE SNES time = 218.9 seconds
>>>
>>> NP 2
>>>
>>> FLOW TS BE steps = 67 newton = 176 linear = 314 cuts =
>>> 13
>>> FLOW TS BE Wasted Linear Iterations = 208
>>> FLOW TS BE SNES time = 600.0 seconds
>>>
>>>
>>> Robert
>>>
>>> On 06/07/17 21:24, Barry Smith wrote:
>>>> So on one process the outer linear solver takes a single iteration
>>>> this is because the block Jacobi with LU and one block is a direct solver.
>>>>
>>>>
>>>>> 11 KSP preconditioned resid norm 1.131868956745e+00 true resid
>>>>> norm 1.526261825526e-05 ||r(i)||/||b|| 1.485509868409e-05
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 2.148515820410e-14 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 1.581814306485e-02 at iteration 1
>>>>> 1 KSP unpreconditioned resid norm 2.148515820410e-14 true resid
>>>>> norm 2.148698024622e-14 ||r(i)||/||b|| 1.358375642332e-12
>>>>>
>>>> On two processes the outer linear solver takes a few iterations to
>>>> solver, this is to be expected.
>>>>
>>>> But what you sent doesn't give any indication about SNES not
>>>> converging. Please turn off all inner linear solver monitoring and just
>>>> run with -ksp_monitor_true_residual -snes_monitor -snes_lineseach_monitor
>>>> -snes_converged_reason
>>>>
>>>> Barry
>>>>
>>>>
>>>>
>>>>
>>>>> On Jul 6, 2017, at 2:03 PM, Robert Annewandter
>>>>> <robert.annewand...@opengosim.com>
>>>>> <mailto:robert.annewand...@opengosim.com>
>>>>> wrote:
>>>>>
>>>>> Hi all,
>>>>>
>>>>> I like to understand why the SNES of my CPR-AMG Two-Stage
>>>>> Preconditioner (with KSPFGMRES + multipl. PCComposite (PCGalerkin with
>>>>> KSPGMRES + BoomerAMG, PCBJacobi + PCLU init) on a 24,000 x 24,000 matrix)
>>>>> struggles to converge when using two cores instead of one. Because of the
>>>>> adaptive time stepping of the Newton, this leads to severe cuts in time
>>>>> step.
>>>>>
>>>>> This is how I run it with two cores
>>>>>
>>>>> mpirun \
>>>>> -n 2 pflotran \
>>>>> -pflotranin het.pflinput \
>>>>> -ksp_monitor_true_residual \
>>>>> -flow_snes_view \
>>>>> -flow_snes_converged_reason \
>>>>> -flow_sub_1_pc_type bjacobi \
>>>>> -flow_sub_1_sub_pc_type lu \
>>>>> -flow_sub_1_sub_pc_factor_pivot_in_blocks true\
>>>>> -flow_sub_1_sub_pc_factor_nonzeros_along_diagonal \
>>>>> -options_left \
>>>>> -log_summary \
>>>>> -info
>>>>>
>>>>>
>>>>> With one core I get (after grepping the crap away from -info):
>>>>>
>>>>> Step 32 Time= 1.80000E+01
>>>>>
>>>>> [...]
>>>>>
>>>>> 0 2r: 1.58E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 7.18E-03 iu: 0.00E+00
>>>>> rsn: 0
>>>>> [0] SNESComputeJacobian(): Rebuilding preconditioner
>>>>> Residual norms for flow_ solve.
>>>>> 0 KSP unpreconditioned resid norm 1.581814306485e-02 true resid
>>>>> norm 1.581814306485e-02 ||r(i)||/||b|| 1.000000000000e+00
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 5.697603110484e+07 true resid
>>>>> norm 5.175721849125e+03 ||r(i)||/||b|| 5.037527476892e+03
>>>>> 1 KSP preconditioned resid norm 5.041509073319e+06 true resid
>>>>> norm 3.251596928176e+02 ||r(i)||/||b|| 3.164777657484e+02
>>>>> 2 KSP preconditioned resid norm 1.043761838360e+06 true resid
>>>>> norm 8.957519558348e+01 ||r(i)||/||b|| 8.718349288342e+01
>>>>> 3 KSP preconditioned resid norm 1.129189815646e+05 true resid
>>>>> norm 2.722436912053e+00 ||r(i)||/||b|| 2.649746479496e+00
>>>>> 4 KSP preconditioned resid norm 8.829637298082e+04 true resid
>>>>> norm 8.026373593492e+00 ||r(i)||/||b|| 7.812065388300e+00
>>>>> 5 KSP preconditioned resid norm 6.506021637694e+04 true resid
>>>>> norm 3.479889319880e+00 ||r(i)||/||b|| 3.386974527698e+00
>>>>> 6 KSP preconditioned resid norm 6.392263200180e+04 true resid
>>>>> norm 3.819202631980e+00 ||r(i)||/||b|| 3.717228003987e+00
>>>>> 7 KSP preconditioned resid norm 2.464946645480e+04 true resid
>>>>> norm 7.329964753388e-01 ||r(i)||/||b|| 7.134251013911e-01
>>>>> 8 KSP preconditioned resid norm 2.603879153772e+03 true resid
>>>>> norm 2.035525412004e-02 ||r(i)||/||b|| 1.981175861414e-02
>>>>> 9 KSP preconditioned resid norm 1.774410462754e+02 true resid
>>>>> norm 3.001214973121e-03 ||r(i)||/||b|| 2.921081026352e-03
>>>>> 10 KSP preconditioned resid norm 1.664227038378e+01 true resid
>>>>> norm 3.413136309181e-04 ||r(i)||/||b|| 3.322003855903e-04
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 1.131868956745e+00 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 2.067297386780e+07 at iteration 11
>>>>> 11 KSP preconditioned resid norm 1.131868956745e+00 true resid
>>>>> norm 1.526261825526e-05 ||r(i)||/||b|| 1.485509868409e-05
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 2.148515820410e-14 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 1.581814306485e-02 at iteration 1
>>>>> 1 KSP unpreconditioned resid norm 2.148515820410e-14 true resid
>>>>> norm 2.148698024622e-14 ||r(i)||/||b|| 1.358375642332e-12
>>>>> [0] SNESSolve_NEWTONLS(): iter=0, linear solve iterations=1
>>>>> [0] SNESNEWTONLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX||
>>>>> 3.590873180642e-01 near zero implies inconsistent rhs
>>>>> [0] SNESSolve_NEWTONLS(): fnorm=1.5818143064846742e-02,
>>>>> gnorm=1.0695649833687331e-02, ynorm=4.6826522561266171e+02, lssucceed=0
>>>>> [0] SNESConvergedDefault(): Converged due to small update length:
>>>>> 4.682652256127e+02 < 1.000000000000e-05 * 3.702480426117e+09
>>>>> 1 2r: 1.07E-02 2x: 3.70E+09 2u: 4.68E+02 ir: 5.05E-03 iu: 4.77E+01
>>>>> rsn: stol
>>>>> Nonlinear flow_ solve converged due to CONVERGED_SNORM_RELATIVE
>>>>> iterations 1
>>>>>
>>>>>
>>>>>
>>>>> But with two cores I get:
>>>>>
>>>>>
>>>>> Step 32 Time= 1.80000E+01
>>>>>
>>>>> [...]
>>>>>
>>>>> 0 2r: 6.16E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.63E-03 iu: 0.00E+00
>>>>> rsn: 0
>>>>> [0] SNESComputeJacobian(): Rebuilding preconditioner
>>>>>
>>>>> Residual norms for flow_ solve.
>>>>> 0 KSP unpreconditioned resid norm 6.162760088924e-03 true resid
>>>>> norm 6.162760088924e-03 ||r(i)||/||b|| 1.000000000000e+00
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 8.994949630499e+08 true resid
>>>>> norm 7.982144380936e-01 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 8.950556502615e+08 true resid
>>>>> norm 1.550138696155e+00 ||r(i)||/||b|| 1.942007839218e+00
>>>>> 2 KSP preconditioned resid norm 1.044849684205e+08 true resid
>>>>> norm 2.166193480531e+00 ||r(i)||/||b|| 2.713798920631e+00
>>>>> 3 KSP preconditioned resid norm 8.209708619718e+06 true resid
>>>>> norm 3.076045005154e-01 ||r(i)||/||b|| 3.853657436340e-01
>>>>> 4 KSP preconditioned resid norm 3.027461352422e+05 true resid
>>>>> norm 1.207731865714e-02 ||r(i)||/||b|| 1.513041869549e-02
>>>>> 5 KSP preconditioned resid norm 1.595302164817e+04 true resid
>>>>> norm 4.123713694368e-04 ||r(i)||/||b|| 5.166172769585e-04
>>>>> 6 KSP preconditioned resid norm 1.898935810797e+03 true resid
>>>>> norm 8.275885058330e-05 ||r(i)||/||b|| 1.036799719897e-04
>>>>> 7 KSP preconditioned resid norm 1.429881682558e+02 true resid
>>>>> norm 4.751240525466e-06 ||r(i)||/||b|| 5.952335987324e-06
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 8.404003313455e+00 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 8.994949630499e+08 at iteration 8
>>>>> 8 KSP preconditioned resid norm 8.404003313455e+00 true resid
>>>>> norm 3.841921844578e-07 ||r(i)||/||b|| 4.813145016211e-07
>>>>> 1 KSP unpreconditioned resid norm 6.162162548202e-03 true resid
>>>>> norm 6.162162548202e-03 ||r(i)||/||b|| 9.999030400804e-01
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 4.360556381209e+07 true resid
>>>>> norm 1.000000245433e+00 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 5.385519331932e+06 true resid
>>>>> norm 8.785183939860e-02 ||r(i)||/||b|| 8.785181783689e-02
>>>>> 2 KSP preconditioned resid norm 4.728931283459e+05 true resid
>>>>> norm 2.008708805316e-02 ||r(i)||/||b|| 2.008708312313e-02
>>>>> 3 KSP preconditioned resid norm 2.734215698319e+04 true resid
>>>>> norm 6.418720397673e-03 ||r(i)||/||b|| 6.418718822309e-03
>>>>> 4 KSP preconditioned resid norm 1.002270029334e+04 true resid
>>>>> norm 4.040289515991e-03 ||r(i)||/||b|| 4.040288524372e-03
>>>>> 5 KSP preconditioned resid norm 1.321280190971e+03 true resid
>>>>> norm 1.023292238313e-04 ||r(i)||/||b|| 1.023291987163e-04
>>>>> 6 KSP preconditioned resid norm 6.594292964815e+01 true resid
>>>>> norm 1.877106733170e-06 ||r(i)||/||b|| 1.877106272467e-06
>>>>> 7 KSP preconditioned resid norm 7.816325147216e+00 true resid
>>>>> norm 2.552611664980e-07 ||r(i)||/||b|| 2.552611038486e-07
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 6.391568446109e-01 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 4.360556381209e+07 at iteration 8
>>>>> 8 KSP preconditioned resid norm 6.391568446109e-01 true resid
>>>>> norm 1.680724939670e-08 ||r(i)||/||b|| 1.680724527166e-08
>>>>> 2 KSP unpreconditioned resid norm 4.328902922753e-07 true resid
>>>>> norm 4.328902922752e-07 ||r(i)||/||b|| 7.024292460341e-05
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 8.794597825780e+08 true resid
>>>>> norm 1.000000094566e+00 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 8.609906572102e+08 true resid
>>>>> norm 2.965044981249e+00 ||r(i)||/||b|| 2.965044700856e+00
>>>>> 2 KSP preconditioned resid norm 9.318108989314e+07 true resid
>>>>> norm 1.881262939380e+00 ||r(i)||/||b|| 1.881262761477e+00
>>>>> 3 KSP preconditioned resid norm 6.908723262483e+06 true resid
>>>>> norm 2.639592490398e-01 ||r(i)||/||b|| 2.639592240782e-01
>>>>> 4 KSP preconditioned resid norm 2.651677791227e+05 true resid
>>>>> norm 9.736480169584e-03 ||r(i)||/||b|| 9.736479248845e-03
>>>>> 5 KSP preconditioned resid norm 1.192178471172e+04 true resid
>>>>> norm 3.082839752692e-04 ||r(i)||/||b|| 3.082839461160e-04
>>>>> 6 KSP preconditioned resid norm 1.492201446262e+03 true resid
>>>>> norm 4.633866284506e-05 ||r(i)||/||b|| 4.633865846301e-05
>>>>> 7 KSP preconditioned resid norm 1.160670017241e+02 true resid
>>>>> norm 2.821157348522e-06 ||r(i)||/||b|| 2.821157081737e-06
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 6.447568262216e+00 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 8.794597825780e+08 at iteration 8
>>>>> 8 KSP preconditioned resid norm 6.447568262216e+00 true resid
>>>>> norm 1.516068561348e-07 ||r(i)||/||b|| 1.516068417980e-07
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 6.135731709822e-15 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 6.162760088924e-03 at iteration 3
>>>>> 3 KSP unpreconditioned resid norm 6.135731709822e-15 true resid
>>>>> norm 1.142020328809e-14 ||r(i)||/||b|| 1.853098793933e-12
>>>>>
>>>>> [0] SNESSolve_NEWTONLS(): iter=0, linear solve iterations=3
>>>>> [0] SNESNEWTONLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX||
>>>>> 1.998388224666e-02 near zero implies inconsistent rhs
>>>>> [0] SNESSolve_NEWTONLS(): fnorm=6.1627600889243711e-03,
>>>>> gnorm=1.0406503258190572e-02, ynorm=6.2999025681245366e+04, lssucceed=0
>>>>> 1 2r: 1.04E-02 2x: 3.70E+09 2u: 6.30E+04 ir: 6.54E-03 iu: 5.00E+04
>>>>> rsn: 0
>>>>> [0] SNESComputeJacobian(): Rebuilding preconditioner
>>>>>
>>>>> Residual norms for flow_ solve.
>>>>> 0 KSP unpreconditioned resid norm 1.040650325819e-02 true resid
>>>>> norm 1.040650325819e-02 ||r(i)||/||b|| 1.000000000000e+00
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 6.758906811264e+07 true resid
>>>>> norm 9.814998431686e-01 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 2.503922806424e+06 true resid
>>>>> norm 2.275130113021e-01 ||r(i)||/||b|| 2.318013730574e-01
>>>>> 2 KSP preconditioned resid norm 3.316753614870e+05 true resid
>>>>> norm 3.820733530238e-02 ||r(i)||/||b|| 3.892750016040e-02
>>>>> 3 KSP preconditioned resid norm 2.956751700483e+04 true resid
>>>>> norm 2.143772538677e-03 ||r(i)||/||b|| 2.184180215207e-03
>>>>> 4 KSP preconditioned resid norm 1.277067042524e+03 true resid
>>>>> norm 9.093614251311e-05 ||r(i)||/||b|| 9.265018547485e-05
>>>>> 5 KSP preconditioned resid norm 1.060996002446e+02 true resid
>>>>> norm 1.042893700050e-05 ||r(i)||/||b|| 1.062551061326e-05
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 5.058127343285e+00 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 6.758906811264e+07 at iteration 6
>>>>> 6 KSP preconditioned resid norm 5.058127343285e+00 true resid
>>>>> norm 4.054770602120e-07 ||r(i)||/||b|| 4.131198420807e-07
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 4.449606189225e-10 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 1.040650325819e-02 at iteration 1
>>>>> 1 KSP unpreconditioned resid norm 4.449606189225e-10 true resid
>>>>> norm 4.449606189353e-10 ||r(i)||/||b|| 4.275793779098e-08
>>>>>
>>>>> [0] SNESSolve_NEWTONLS(): iter=1, linear solve iterations=1
>>>>> [0] SNESNEWTONLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX||
>>>>> 4.300066663571e-02 near zero implies inconsistent rhs
>>>>> [0] SNESSolve_NEWTONLS(): fnorm=1.0406503258190572e-02,
>>>>> gnorm=7.3566280848133728e-02, ynorm=7.9500485128639993e+04, lssucceed=0
>>>>> 2 2r: 7.36E-02 2x: 3.70E+09 2u: 7.95E+04 ir: 4.62E-02 iu: 5.00E+04
>>>>> rsn: 0
>>>>> [0] SNESComputeJacobian(): Rebuilding preconditioner
>>>>>
>>>>> Residual norms for flow_ solve.
>>>>> 0 KSP unpreconditioned resid norm 7.356628084813e-02 true resid
>>>>> norm 7.356628084813e-02 ||r(i)||/||b|| 1.000000000000e+00
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 7.253424029194e+06 true resid
>>>>> norm 9.647008645250e-01 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 7.126940190688e+06 true resid
>>>>> norm 1.228009197928e+00 ||r(i)||/||b|| 1.272942984800e+00
>>>>> 2 KSP preconditioned resid norm 9.391591432635e+05 true resid
>>>>> norm 7.804929162756e-01 ||r(i)||/||b|| 8.090517433711e-01
>>>>> 3 KSP preconditioned resid norm 6.538499674761e+04 true resid
>>>>> norm 5.503467432893e-02 ||r(i)||/||b|| 5.704843475602e-02
>>>>> 4 KSP preconditioned resid norm 1.593713396575e+04 true resid
>>>>> norm 8.902701363763e-02 ||r(i)||/||b|| 9.228457951208e-02
>>>>> 5 KSP preconditioned resid norm 4.837260621464e+02 true resid
>>>>> norm 2.966772992825e-03 ||r(i)||/||b|| 3.075329464213e-03
>>>>> 6 KSP preconditioned resid norm 1.681372767335e+02 true resid
>>>>> norm 5.312467443025e-04 ||r(i)||/||b|| 5.506854651406e-04
>>>>> 7 KSP preconditioned resid norm 1.271478850717e+01 true resid
>>>>> norm 2.123810020488e-05 ||r(i)||/||b|| 2.201521838103e-05
>>>>> 8 KSP preconditioned resid norm 1.262723712696e+00 true resid
>>>>> norm 1.150572715331e-06 ||r(i)||/||b|| 1.192673042641e-06
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 9.053072585125e-02 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 7.253424029194e+06 at iteration 9
>>>>> 9 KSP preconditioned resid norm 9.053072585125e-02 true resid
>>>>> norm 9.475050575058e-08 ||r(i)||/||b|| 9.821749853747e-08
>>>>> 1 KSP unpreconditioned resid norm 8.171589173162e-03 true resid
>>>>> norm 8.171589173162e-03 ||r(i)||/||b|| 1.110779161180e-01
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 4.345765068989e+07 true resid
>>>>> norm 9.999992231691e-01 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 5.388715093466e+06 true resid
>>>>> norm 8.125387327699e-02 ||r(i)||/||b|| 8.125393639755e-02
>>>>> 2 KSP preconditioned resid norm 4.763725726436e+05 true resid
>>>>> norm 2.464285618036e-02 ||r(i)||/||b|| 2.464287532371e-02
>>>>> 3 KSP preconditioned resid norm 2.287746683380e+04 true resid
>>>>> norm 7.224823080100e-03 ||r(i)||/||b|| 7.224828692570e-03
>>>>> 4 KSP preconditioned resid norm 4.872858764091e+03 true resid
>>>>> norm 3.972261388893e-03 ||r(i)||/||b|| 3.972264474670e-03
>>>>> 5 KSP preconditioned resid norm 8.670449895323e+02 true resid
>>>>> norm 2.359005963873e-04 ||r(i)||/||b|| 2.359007796423e-04
>>>>> 6 KSP preconditioned resid norm 4.252589693890e+01 true resid
>>>>> norm 1.471904261226e-06 ||r(i)||/||b|| 1.471905404648e-06
>>>>> 7 KSP preconditioned resid norm 5.128476471782e+00 true resid
>>>>> norm 1.643725157865e-07 ||r(i)||/||b|| 1.643726434763e-07
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 4.311901915856e-01 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 4.345765068989e+07 at iteration 8
>>>>> 8 KSP preconditioned resid norm 4.311901915856e-01 true resid
>>>>> norm 1.166123921637e-08 ||r(i)||/||b|| 1.166124827519e-08
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 2.373662391739e-09 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 7.356628084813e-02 at iteration 2
>>>>> 2 KSP unpreconditioned resid norm 2.373662391739e-09 true resid
>>>>> norm 2.373662391658e-09 ||r(i)||/||b|| 3.226562990941e-08
>>>>>
>>>>> [0] SNESSolve_NEWTONLS(): iter=2, linear solve iterations=2
>>>>> [0] SNESNEWTONLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX||
>>>>> 4.343326231305e-02 near zero implies inconsistent rhs
>>>>> [0] SNESSolve_NEWTONLS(): fnorm=7.3566280848133728e-02,
>>>>> gnorm=7.2259942496422647e-02, ynorm=6.3156901950486099e+04, lssucceed=0
>>>>> 3 2r: 7.23E-02 2x: 3.70E+09 2u: 6.32E+04 ir: 4.52E-02 iu: 5.00E+04
>>>>> rsn: 0
>>>>> [0] SNESComputeJacobian(): Rebuilding preconditioner
>>>>>
>>>>> Residual norms for flow_ solve.
>>>>> 0 KSP unpreconditioned resid norm 7.225994249642e-02 true resid
>>>>> norm 7.225994249642e-02 ||r(i)||/||b|| 1.000000000000e+00
>>>>> Residual norms for flow_sub_0_galerkin_ solve.
>>>>> 0 KSP preconditioned resid norm 7.705582590638e+05 true resid
>>>>> norm 9.649751442741e-01 ||r(i)||/||b|| 1.000000000000e+00
>>>>> 1 KSP preconditioned resid norm 2.444424220392e+04 true resid
>>>>> norm 8.243110200738e-03 ||r(i)||/||b|| 8.542303135630e-03
>>>>> 2 KSP preconditioned resid norm 2.080899648412e+03 true resid
>>>>> norm 7.642343147053e-04 ||r(i)||/||b|| 7.919730567570e-04
>>>>> 3 KSP preconditioned resid norm 9.911171129874e+02 true resid
>>>>> norm 5.904182179180e-05 ||r(i)||/||b|| 6.118481096859e-05
>>>>> 4 KSP preconditioned resid norm 5.258230282482e+02 true resid
>>>>> norm 2.043366677644e-04 ||r(i)||/||b|| 2.117532964210e-04
>>>>> 5 KSP preconditioned resid norm 5.522830460456e+01 true resid
>>>>> norm 1.710780366056e-05 ||r(i)||/||b|| 1.772875059225e-05
>>>>> 6 KSP preconditioned resid norm 5.922280741715e+00 true resid
>>>>> norm 1.543198740828e-06 ||r(i)||/||b|| 1.599210870855e-06
>>>>> 7 KSP preconditioned resid norm 3.339500859115e-01 true resid
>>>>> norm 1.221335666427e-07 ||r(i)||/||b|| 1.265665414984e-07
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 3.329208597672e-02 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 7.705582590638e+05 at iteration 8
>>>>> 8 KSP preconditioned resid norm 3.329208597672e-02 true resid
>>>>> norm 9.758240835324e-09 ||r(i)||/||b|| 1.011242713683e-08
>>>>> [0] KSPConvergedDefault(): Linear solver has converged. Residual norm
>>>>> 2.697128456432e-11 is less than relative tolerance 1.000000000000e-07
>>>>> times initial right hand side norm 7.225994249642e-02 at iteration 1
>>>>> 1 KSP unpreconditioned resid norm 2.697128456432e-11 true resid
>>>>> norm 2.697128457142e-11 ||r(i)||/||b|| 3.732536124389e-10
>>>>>
>>>>> [0] SNESSolve_NEWTONLS(): iter=3, linear solve iterations=1
>>>>> [0] SNESNEWTONLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX||
>>>>> 4.329227684222e-02 near zero implies inconsistent rhs
>>>>> [0] SNESSolve_NEWTONLS(): fnorm=7.2259942496422647e-02,
>>>>> gnorm=5.4435602192925014e-01, ynorm=2.7049750229137400e+04, lssucceed=0
>>>>> [0] SNESConvergedDefault(): Converged due to small update length:
>>>>> 2.704975022914e+04 < 1.000000000000e-05 * 3.702469482296e+09
>>>>> 4 2r: 5.44E-01 2x: 3.70E+09 2u: 2.70E+04 ir: 3.84E-01 iu: 2.34E+04
>>>>> rsn: stol
>>>>> Nonlinear flow_ solve converged due to CONVERGED_SNORM_RELATIVE
>>>>> iterations 4
>>>>>
>>>>>
>>>>> As the simulation advances this behaviour leads to frequent time step
>>>>> cuts because of 8 subsequently failed Newton iterations, which brings the
>>>>> simulation practically to a halt.
>>>>>
>>>>> Is the Block Jacobi not a good choice? Better ASM with huge overlap?
>>>>> Or is there something wrong with my RHS? Maybe the SNES, SNESLS, KSP
>>>>> tolerances need better tuning?
>>>>>
>>>>> Grateful for any clarifying words!
>>>>> Robert
>>>>>
>>>>>
>>>>> My SNES_view is:
>>>>>
>>>>>
>>>>> SNES Object: (flow_) 2 MPI processes
>>>>> type: newtonls
>>>>> maximum iterations=8, maximum function evaluations=10000
>>>>> tolerances: relative=1e-05, absolute=1e-05, solution=1e-05
>>>>> total number of linear solver iterations=1
>>>>> total number of function evaluations=2
>>>>> norm schedule ALWAYS
>>>>> SNESLineSearch Object: (flow_) 2 MPI processes
>>>>> type: basic
>>>>> maxstep=1.000000e+08, minlambda=1.000000e-05
>>>>> tolerances: relative=1.000000e-05, absolute=1.000000e-05,
>>>>> lambda=1.000000e-08
>>>>> maximum iterations=40
>>>>> using user-defined precheck step
>>>>> KSP Object: (flow_) 2 MPI processes
>>>>> type: fgmres
>>>>> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt
>>>>> Orthogonalization with no iterative refinement
>>>>> GMRES: happy breakdown tolerance 1e-30
>>>>> maximum iterations=200, initial guess is zero
>>>>> tolerances: relative=1e-07, absolute=1e-50, divergence=10000.
>>>>> right preconditioning
>>>>> using UNPRECONDITIONED norm type for convergence test
>>>>> PC Object: (flow_) 2 MPI processes
>>>>> type: composite
>>>>> Composite PC type - MULTIPLICATIVE
>>>>> PCs on composite preconditioner follow
>>>>> ---------------------------------
>>>>> PC Object: (flow_sub_0_) 2 MPI processes
>>>>> type: galerkin
>>>>> Galerkin PC
>>>>> KSP on Galerkin follow
>>>>> ---------------------------------
>>>>> KSP Object: (flow_sub_0_galerkin_) 2 MPI processes
>>>>> type: gmres
>>>>> GMRES: restart=30, using Classical (unmodified)
>>>>> Gram-Schmidt Orthogonalization with no iterative refinement
>>>>> GMRES: happy breakdown tolerance 1e-30
>>>>> maximum iterations=200, initial guess is zero
>>>>> tolerances: relative=1e-07, absolute=1e-50,
>>>>> divergence=10000.
>>>>> left preconditioning
>>>>> using PRECONDITIONED norm type for convergence test
>>>>> PC Object: (flow_sub_0_galerkin_) 2 MPI processes
>>>>> type: hypre
>>>>> HYPRE BoomerAMG preconditioning
>>>>> HYPRE BoomerAMG: Cycle type V
>>>>> HYPRE BoomerAMG: Maximum number of levels 25
>>>>> HYPRE BoomerAMG: Maximum number of iterations PER hypre
>>>>> call 1
>>>>> HYPRE BoomerAMG: Convergence tolerance PER hypre call 0.
>>>>> HYPRE BoomerAMG: Threshold for strong coupling 0.25
>>>>> HYPRE BoomerAMG: Interpolation truncation factor 0.
>>>>> HYPRE BoomerAMG: Interpolation: max elements per row 0
>>>>> HYPRE BoomerAMG: Number of levels of aggressive
>>>>> coarsening 0
>>>>> HYPRE BoomerAMG: Number of paths for aggressive
>>>>> coarsening 1
>>>>> HYPRE BoomerAMG: Maximum row sums 0.9
>>>>> HYPRE BoomerAMG: Sweeps down 1
>>>>> HYPRE BoomerAMG: Sweeps up 1
>>>>> HYPRE BoomerAMG: Sweeps on coarse 1
>>>>> HYPRE BoomerAMG: Relax down symmetric-SOR/Jacobi
>>>>> HYPRE BoomerAMG: Relax up symmetric-SOR/Jacobi
>>>>> HYPRE BoomerAMG: Relax on coarse Gaussian-elimination
>>>>> HYPRE BoomerAMG: Relax weight (all) 1.
>>>>> HYPRE BoomerAMG: Outer relax weight (all) 1.
>>>>> HYPRE BoomerAMG: Using CF-relaxation
>>>>> HYPRE BoomerAMG: Not using more complex smoothers.
>>>>> HYPRE BoomerAMG: Measure type local
>>>>> HYPRE BoomerAMG: Coarsen type Falgout
>>>>> HYPRE BoomerAMG: Interpolation type classical
>>>>> linear system matrix = precond matrix:
>>>>> Mat Object: 2 MPI processes
>>>>> type: mpiaij
>>>>> rows=8000, cols=8000
>>>>> total: nonzeros=53600, allocated nonzeros=53600
>>>>> total number of mallocs used during MatSetValues calls =0
>>>>> not using I-node (on process 0) routines
>>>>> linear system matrix = precond matrix:
>>>>> Mat Object: (flow_) 2 MPI processes
>>>>> type: mpibaij
>>>>> rows=24000, cols=24000, bs=3
>>>>> total: nonzeros=482400, allocated nonzeros=482400
>>>>> total number of mallocs used during MatSetValues calls =0
>>>>> PC Object: (flow_sub_1_) 2 MPI processes
>>>>> type: bjacobi
>>>>> block Jacobi: number of blocks = 2
>>>>> Local solve is same for all blocks, in the following KSP
>>>>> and PC objects:
>>>>> KSP Object: (flow_sub_1_sub_) 1 MPI processes
>>>>> type: preonly
>>>>> maximum iterations=10000, initial guess is zero
>>>>> tolerances: relative=1e-05, absolute=1e-50,
>>>>> divergence=10000. <------ not working: -flow_sub_1_sub_ksp_rtol 1e-7
>>>>> left preconditioning
>>>>> using NONE norm type for convergence test
>>>>> PC Object: (flow_sub_1_sub_) 1 MPI processes
>>>>> type: lu
>>>>> out-of-place factorization
>>>>> tolerance for zero pivot 2.22045e-14
>>>>> matrix ordering: nd
>>>>> factor fill ratio given 5., needed 18.3108
>>>>> Factored matrix follows:
>>>>> Mat Object: 1 MPI processes
>>>>> type: seqbaij
>>>>> rows=12000, cols=12000, bs=3
>>>>> package used to perform factorization: petsc
>>>>> total: nonzeros=4350654, allocated nonzeros=4350654
>>>>> total number of mallocs used during MatSetValues
>>>>> calls =0
>>>>> block size is 3
>>>>> linear system matrix = precond matrix:
>>>>> Mat Object: (flow_) 1 MPI processes
>>>>> type: seqbaij
>>>>> rows=12000, cols=12000, bs=3
>>>>> total: nonzeros=237600, allocated nonzeros=237600
>>>>> total number of mallocs used during MatSetValues calls =0
>>>>> block size is 3
>>>>> linear system matrix = precond matrix:
>>>>> Mat Object: (flow_) 2 MPI processes
>>>>> type: mpibaij
>>>>> rows=24000, cols=24000, bs=3
>>>>> total: nonzeros=482400, allocated nonzeros=482400
>>>>> total number of mallocs used during MatSetValues calls =0
>>>>> ---------------------------------
>>>>> linear system matrix = precond matrix:
>>>>> Mat Object: (flow_) 2 MPI processes
>>>>> type: mpibaij
>>>>> rows=24000, cols=24000, bs=3
>>>>> total: nonzeros=482400, allocated nonzeros=482400
>>>>> total number of mallocs used during MatSetValues calls =0
>>>>>
>>> <het_np_1.log><het_np_2.log>
>
>
------------------------------ Provenance --------------------------------------
pflotran_compile_date_time = unknown
pflotran_compile_user = unknown
pflotran_compile_hostname = unknown
pflotran_changeset = unknown
pflotran_status = unknown
petsc_changeset = unknown
petsc_status = unknown
--------------------------------------------------------------------------------
PMTOilImsCreate()
"grid_structured_type" set to default value.
pflotran card:: DEBUG
pflotran card:: GRID
pflotran card:: TIMESTEPPER
pflotran card:: NEWTON_SOLVER
pflotran card:: LINEAR_SOLVER
pflotran card:: TIME
pflotran card:: OUTPUT
"OUTPUT,FORMAT,HDF5,# FILES" set to default value.
"OUTPUT,MASS_BALANCE,DETAILED" set to default value.
pflotran card:: MATERIAL_PROPERTY
Name :: formation
"MATERIAL_PROPERTY,rock density units" set to default value.
"MATERIAL_PROPERTY,specific heat units" set to default value.
"MATERIAL_PROPERTY,dry thermal conductivity units" set to default value.
"MATERIAL_PROPERTY,wet thermal conductivity units" set to default value.
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: CHARACTERISTIC_CURVES
Name :: ch1
pflotran card:: EOS
"EOS,OIL,VISCOSITY,CONSTANT units" set to default value.
pflotran card:: REGION
all
pflotran card:: REGION
inj_well
pflotran card:: REGION
prod_well
pflotran card:: REGION
prod_well_west
pflotran card:: REGION
south_boundary
pflotran card:: FLOW_CONDITION
initial_press
"PRESSURE UNITS" set to default value.
"OIL_SATURATION UNITS" set to default value.
"DATUM UNITS" set to default value.
"OWC UNITS" set to default value.
pflotran card:: FLOW_CONDITION
wat_injection
pflotran card:: INITIAL_CONDITION
pflotran card:: BOUNDARY_CONDITION
pflotran card:: SOURCE_SINK
pflotran card:: STRATA
pflotran card:: END_SUBSURFACE
Domain Bounds (x y z):
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
1.8250000000E+02 3.3500000000E+02 6.1000000000E+00
Testing group: Porosity
Group "Porosity" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityX
Group "PermeabilityX" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityY
Group "PermeabilityY" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityZ
Group "PermeabilityZ" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Opening hdf5 file: mat_ids_20x20x20.h5
Setting up grid cell indices
Opening group: Materials
Reading dataset: Materials/Cell Ids
0.00 Seconds to set up indices
Reading dataset: Materials/Material Ids
0.00 Seconds to read integer array.
Closing group: Materials
Closing hdf5 file: mat_ids_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityX
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityY
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityZ
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: Porosity
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Default HDF5 method is used in Initialization
Requested processors and decomposition = 1, npx,y,z= -1 -1 -1
Actual decomposition: npx,y,z= 1 1 1
Beginning setup of FLOW SNES
number of dofs = 3, number of phases = 2
Finished setting up FLOW SNES
Solver: gmres
Preconditioner: bjacobi
FLOW Time Stepper
max steps: 999999
max constant cumulative time steps: 5
max cuts: 16
ts reduction factor: 0.50000000000000000
FLOW Newton Solver
atol: 1.0000E-05
rtol: 1.0000E-05
stol: 1.0000E-05
dtol: 1.0000E+15
maxnorm: 1.0000E+20
inftolres: -9.9900E+02
inftolupd: -9.9900E+02
inftolrelupd: -9.9900E+02
inftolsclres: -9.9900E+02
max iter: 8
min iter: 1
maxf: 10000
matrix type:baij
precond. matrix type:baij
inexact newton: off
print convergence: on
print detailed convergence: off
check infinity norm: on
FLOW Linear Solver
solver: gmres
precond: bjacobi
atol: 1.0000E-50
rtol: 1.0000E-07
dtol: 1.0000E+04
max iter: 200
--> creating hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
Finished Initialization
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.168250330965e-03
0 2r: 6.17E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.207896105644e-09
1 2r: 9.21E-09 2x: 3.70E+09 2u: 1.84E+04 ir: 4.52E-09 iu: 1.69E+03 rsn: atol
Step 1 Time= 1.00000E-04 Dt= 1.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 1] linear = 1 [ 1] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 9.207896E-09 1.150987E-12 4.518315E-09
--> max chng: dpl= 1.6850E+03 dpo= 1.6850E+03 dso= 7.5593E-05 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.204655900376e-03
0 2r: 6.20E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.245052250819e-08
1 2r: 1.25E-08 2x: 3.70E+09 2u: 2.66E+04 ir: 5.54E-09 iu: 9.79E+02 rsn: atol
Step 2 Time= 3.00000E-04 Dt= 2.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 2] linear = 1 [ 2] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.245052E-08 1.556315E-12 5.543064E-09
--> max chng: dpl= 9.7897E+02 dpo= 9.7897E+02 dso= 1.5119E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205008271412e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.434590204248e-08
1 2r: 1.43E-08 2x: 3.70E+09 2u: 3.99E+04 ir: 4.04E-09 iu: 1.05E+03 rsn: atol
Step 3 Time= 7.00000E-04 Dt= 4.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 3] linear = 1 [ 3] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.434590E-08 1.793238E-12 4.040496E-09
--> max chng: dpl= 1.0535E+03 dpo= 1.0535E+03 dso= 3.0237E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205059165175e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.415138478985e-08
1 2r: 1.42E-08 2x: 3.70E+09 2u: 5.34E+04 ir: 6.20E-09 iu: 1.18E+03 rsn: atol
Step 4 Time= 1.50000E-03 Dt= 8.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 4] linear = 1 [ 4] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.415138E-08 1.768923E-12 6.197879E-09
--> max chng: dpl= 1.1776E+03 dpo= 1.1776E+03 dso= 6.0475E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205088543068e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.608901160729e-08
1 2r: 1.61E-08 2x: 3.70E+09 2u: 5.71E+04 ir: 6.30E-09 iu: 1.18E+03 rsn: atol
Step 5 Time= 3.10000E-03 Dt= 1.60000E-03 [d] snes_conv_reason: 2
newton = 1 [ 5] linear = 1 [ 5] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.608901E-08 2.011126E-12 6.296977E-09
--> max chng: dpl= 1.1798E+03 dpo= 1.1798E+03 dso= 1.2095E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205110269648e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.760879149186e-09
1 2r: 7.76E-09 2x: 3.70E+09 2u: 4.40E+04 ir: 2.22E-09 iu: 8.97E+02 rsn: atol
Step 6 Time= 6.30000E-03 Dt= 3.20000E-03 [d] snes_conv_reason: 2
newton = 1 [ 6] linear = 1 [ 6] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 7.760879E-09 9.701099E-13 2.223273E-09
--> max chng: dpl= 8.9740E+02 dpo= 8.9740E+02 dso= 2.4190E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205123874338e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.495714625814e-09
1 2r: 3.50E-09 2x: 3.70E+09 2u: 2.28E+04 ir: 1.08E-09 iu: 5.16E+02 rsn: atol
Step 7 Time= 1.27000E-02 Dt= 6.40000E-03 [d] snes_conv_reason: 2
newton = 1 [ 7] linear = 1 [ 7] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.495715E-09 4.369643E-13 1.082495E-09
--> max chng: dpl= 5.1587E+02 dpo= 5.1587E+02 dso= 4.8380E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205129605472e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.405765028609e-10
1 2r: 9.41E-10 2x: 3.70E+09 2u: 7.52E+03 ir: 2.75E-10 iu: 2.20E+02 rsn: atol
Step 8 Time= 2.55000E-02 Dt= 1.28000E-02 [d] snes_conv_reason: 2
newton = 1 [ 8] linear = 1 [ 8] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 9.405765E-10 1.175721E-13 2.748286E-10
--> max chng: dpl= 2.2049E+02 dpo= 2.2049E+02 dso= 9.6760E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131042559e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.789240999780e-10
1 2r: 1.79E-10 2x: 3.70E+09 2u: 1.45E+03 ir: 5.71E-11 iu: 4.90E+01 rsn: atol
Step 9 Time= 5.11000E-02 Dt= 2.56000E-02 [d] snes_conv_reason: 2
newton = 1 [ 9] linear = 1 [ 9] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.789241E-10 2.236551E-14 5.714697E-11
--> max chng: dpl= 4.9014E+01 dpo= 4.9014E+01 dso= 1.9352E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131242597e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.824938719594e-11
1 2r: 3.82E-11 2x: 3.70E+09 2u: 1.51E+02 ir: 3.39E-11 iu: 5.52E+00 rsn: atol
Step 10 Time= 1.02300E-01 Dt= 5.12000E-02 [d] snes_conv_reason: 2
newton = 1 [ 10] linear = 1 [ 10] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.824939E-11 4.781173E-15 3.385756E-11
--> max chng: dpl= 5.5219E+00 dpo= 5.5219E+00 dso= 3.8704E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131262406e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.855561772868e-11
1 2r: 2.86E-11 2x: 3.70E+09 2u: 8.23E+00 ir: 2.69E-11 iu: 3.12E-01 rsn: atol
Step 11 Time= 2.04700E-01 Dt= 1.02400E-01 [d] snes_conv_reason: 2
newton = 1 [ 11] linear = 1 [ 11] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 2.855562E-11 3.569452E-15 2.691689E-11
--> max chng: dpl= 3.1209E-01 dpo= 3.1209E-01 dso= 7.7408E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131263109e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.151916583270e-03
1 2r: 1.15E-03 2x: 3.70E+09 2u: 2.45E-01 ir: 1.11E-03 iu: 1.25E-01 rsn: stol
Step 12 Time= 4.09500E-01 Dt= 2.04800E-01 [d] snes_conv_reason: 4
newton = 1 [ 12] linear = 1 [ 12] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.151917E-03 1.439896E-07 1.112149E-03
--> max chng: dpl= 7.1319E-03 dpo= 7.1319E-03 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.161154480767e-03
0 2r: 6.16E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.95E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.519651739007e-03
1 2r: 3.52E-03 2x: 3.70E+09 2u: 2.00E+02 ir: 2.64E-03 iu: 4.73E+01 rsn: stol
Step 13 Time= 8.19100E-01 Dt= 4.09600E-01 [d] snes_conv_reason: 4
newton = 1 [ 13] linear = 1 [ 13] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.519652E-03 4.399565E-07 2.637128E-03
--> max chng: dpl= 4.7301E+01 dpo= 4.7301E+01 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.116654806200e-03
0 2r: 4.12E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.12E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 8.875572006760e-04
1 2r: 8.88E-04 2x: 3.70E+09 2u: 4.34E+02 ir: 6.89E-04 iu: 9.73E+01 rsn: stol
Step 14 Time= 1.00000E+00 Dt= 1.80900E-01 [d] snes_conv_reason: 4
newton = 1 [ 14] linear = 1 [ 14] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 8.875572E-04 1.109447E-07 6.891907E-04
--> max chng: dpl= 9.7334E+01 dpo= 9.7334E+01 dso= 5.0484E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.00 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 3.845220273464e-03
0 2r: 3.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.32E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.509123063551e-03
1 2r: 3.51E-03 2x: 3.70E+09 2u: 4.78E+01 ir: 2.12E-03 iu: 2.35E+01 rsn: stol
Step 15 Time= 1.81920E+00 Dt= 8.19200E-01 [d] snes_conv_reason: 4
newton = 1 [ 15] linear = 1 [ 15] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.509123E-03 4.386404E-07 2.123911E-03
--> max chng: dpl= 2.3511E+01 dpo= 2.3511E+01 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.846497651403e-03
0 2r: 4.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.73E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.523617853722e-04
1 2r: 3.52E-04 2x: 3.70E+09 2u: 6.18E+01 ir: 2.74E-04 iu: 1.82E+01 rsn: stol
Step 16 Time= 2.00000E+00 Dt= 1.80800E-01 [d] snes_conv_reason: 4
newton = 1 [ 16] linear = 1 [ 16] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.523618E-04 4.404522E-08 2.741021E-04
--> max chng: dpl= 1.8169E+01 dpo= 1.8169E+01 dso= 5.4531E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 3.858579350152e-03
0 2r: 3.86E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.80E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.972530151532e-02
1 2r: 2.97E-02 2x: 3.70E+09 2u: 4.45E+01 ir: 2.15E-02 iu: 2.59E+01 rsn: stol
Step 17 Time= 3.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 17] linear = 1 [ 17] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 2.972530E-02 3.715663E-06 2.148938E-02
--> max chng: dpl= 2.5863E+01 dpo= 2.5863E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.956666973407e-02
0 2r: 2.96E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 2.17E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.551772040002e-02
1 2r: 1.55E-02 2x: 3.70E+09 2u: 2.37E+02 ir: 1.14E-02 iu: 4.75E+01 rsn: stol
Step 18 Time= 4.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 18] linear = 1 [ 18] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.551772E-02 1.939715E-06 1.142377E-02
--> max chng: dpl= 4.7548E+01 dpo= 4.7548E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.604561374530e-02
0 2r: 1.60E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 1.18E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.723820975956e-03
1 2r: 9.72E-03 2x: 3.70E+09 2u: 2.35E+02 ir: 6.24E-03 iu: 3.69E+01 rsn: stol
Step 19 Time= 5.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 19] linear = 1 [ 19] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 9.723821E-03 1.215478E-06 6.241651E-03
--> max chng: dpl= 3.6944E+01 dpo= 3.6944E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.013385215566e-02
0 2r: 1.01E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 6.50E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.526818664055e-03
1 2r: 3.53E-03 2x: 3.70E+09 2u: 1.12E+03 ir: 1.80E-03 iu: 5.33E+02 rsn: stol
Step 20 Time= 6.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 20] linear = 1 [ 20] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.526819E-03 4.408523E-07 1.796417E-03
--> max chng: dpl= 5.3308E+02 dpo= 5.3308E+02 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.496265631168e-03
0 2r: 4.50E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.19E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.413553433289e-02
1 2r: 1.41E-02 2x: 3.70E+09 2u: 5.63E+02 ir: 9.96E-03 iu: 2.92E+02 rsn: stol
Step 21 Time= 7.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 21] linear = 1 [ 21] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.413553E-02 1.766942E-06 9.961481E-03
--> max chng: dpl= 2.9156E+02 dpo= 2.9156E+02 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.435457881275e-02
0 2r: 1.44E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 9.96E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 4.128426194700e-03
1 2r: 4.13E-03 2x: 3.70E+09 2u: 3.31E+02 ir: 1.82E-03 iu: 5.61E+01 rsn: stol
Step 22 Time= 8.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 22] linear = 1 [ 22] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 4.128426E-03 5.160533E-07 1.815179E-03
--> max chng: dpl= 5.6097E+01 dpo= 5.6097E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.541642203116e-03
0 2r: 4.54E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 1.82E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.269985747065e-03
1 2r: 1.27E-03 2x: 3.70E+09 2u: 7.56E+02 ir: 6.38E-04 iu: 8.87E+01 rsn: stol
Step 23 Time= 9.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 23] linear = 1 [ 23] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.269986E-03 1.587482E-07 6.379067E-04
--> max chng: dpl= 8.8711E+01 dpo= 8.8711E+01 dso= 1.2359E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.00 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.518756718655e-03
0 2r: 2.52E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 1.06E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.580925056447e-02
1 2r: 2.58E-02 2x: 3.70E+09 2u: 2.54E+02 ir: 1.39E-02 iu: 5.18E+01 rsn: stol
Step 24 Time= 1.00000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 24] linear = 1 [ 24] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 2.580925E-02 3.226156E-06 1.392561E-02
--> max chng: dpl= 5.1760E+01 dpo= 5.1760E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.585016910324e-02
0 2r: 2.59E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 1.39E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.576693291087e-03
1 2r: 1.58E-03 2x: 3.70E+09 2u: 3.98E+02 ir: 1.04E-03 iu: 8.85E+01 rsn: stol
Step 25 Time= 1.10000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 25] linear = 1 [ 25] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.576693E-03 1.970867E-07 1.039710E-03
--> max chng: dpl= 8.8514E+01 dpo= 8.8514E+01 dso= 1.1914E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.124880047222e-03
0 2r: 2.12E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 7.46E-04 iu: 0.00E+00 rsn: 0
1 SNES Function norm 4.708801965185e-02
1 2r: 4.71E-02 2x: 3.70E+09 2u: 4.49E+02 ir: 2.40E-02 iu: 8.07E+01 rsn: stol
Step 26 Time= 1.20000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 26] linear = 1 [ 26] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 4.708802E-02 5.886002E-06 2.399320E-02
--> max chng: dpl= 8.0735E+01 dpo= 8.0735E+01 dso= 1.0041E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.707347570687e-02
0 2r: 4.71E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 2.40E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.891387006557e-03
1 2r: 3.89E-03 2x: 3.70E+09 2u: 7.12E+02 ir: 2.25E-03 iu: 9.69E+01 rsn: stol
Step 27 Time= 1.30000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 27] linear = 1 [ 27] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 3.891387E-03 4.864234E-07 2.245873E-03
--> max chng: dpl= 9.6868E+01 dpo= 9.6868E+01 dso= 1.1124E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.406138636566e-03
0 2r: 4.41E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.25E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.354276196838e-02
1 2r: 1.35E-02 2x: 3.70E+09 2u: 2.01E+02 ir: 6.28E-03 iu: 3.51E+01 rsn: stol
Step 28 Time= 1.40000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 28] linear = 1 [ 28] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.354276E-02 1.692845E-06 6.278946E-03
--> max chng: dpl= 3.5144E+01 dpo= 3.5144E+01 dso= 1.2185E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.360174200182e-02
0 2r: 1.36E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 6.28E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.643994870466e-02
1 2r: 1.64E-02 2x: 3.70E+09 2u: 5.75E+02 ir: 8.49E-03 iu: 6.06E+01 rsn: stol
Step 29 Time= 1.50000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 29] linear = 1 [ 29] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.643995E-02 2.054994E-06 8.485369E-03
--> max chng: dpl= 6.0570E+01 dpo= 6.0570E+01 dso= 7.6407E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.649479961868e-02
0 2r: 1.65E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 8.49E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.928075369329e-02
1 2r: 1.93E-02 2x: 3.70E+09 2u: 2.64E+02 ir: 9.28E-03 iu: 5.80E+01 rsn: stol
Step 30 Time= 1.60000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 30] linear = 1 [ 30] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.928075E-02 2.410094E-06 9.276234E-03
--> max chng: dpl= 5.7996E+01 dpo= 5.7996E+01 dso= 1.0039E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.00 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.932495533810e-02
0 2r: 1.93E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 9.28E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.256062413457e-03
1 2r: 7.26E-03 2x: 3.70E+09 2u: 4.53E+02 ir: 3.54E-03 iu: 5.03E+01 rsn: stol
Step 31 Time= 1.70000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 31] linear = 1 [ 31] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 7.256062E-03 9.070078E-07 3.541330E-03
--> max chng: dpl= 5.0259E+01 dpo= 5.0259E+01 dso= 1.1286E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 7.416124927908e-03
0 2r: 7.42E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.54E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.576835310687e-02
1 2r: 1.58E-02 2x: 3.70E+09 2u: 2.19E+02 ir: 7.18E-03 iu: 3.47E+01 rsn: stol
Step 32 Time= 1.80000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 32] linear = 1 [ 32] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.576835E-02 1.971044E-06 7.175876E-03
--> max chng: dpl= 3.4726E+01 dpo= 3.4726E+01 dso= 6.5863E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.581817603399e-02
0 2r: 1.58E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 7.18E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.069563629424e-02
1 2r: 1.07E-02 2x: 3.70E+09 2u: 4.68E+02 ir: 5.05E-03 iu: 4.77E+01 rsn: stol
Step 33 Time= 1.90000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 33] linear = 1 [ 33] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.069564E-02 1.336955E-06 5.053627E-03
--> max chng: dpl= 4.7740E+01 dpo= 4.7740E+01 dso= 9.2535E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.080937873895e-02
0 2r: 1.08E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.05E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.075528185745e-02
1 2r: 1.08E-02 2x: 3.70E+09 2u: 1.60E+02 ir: 4.93E-03 iu: 3.39E+01 rsn: stol
Step 34 Time= 2.00000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 34] linear = 1 [ 34] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.075528E-02 1.344410E-06 4.933033E-03
--> max chng: dpl= 3.3911E+01 dpo= 3.3911E+01 dso= 1.0311E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.087291343571e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 4.93E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 8.884926315099e-03
1 2r: 8.88E-03 2x: 3.70E+09 2u: 5.42E+02 ir: 4.14E-03 iu: 7.88E+01 rsn: stol
Step 35 Time= 2.10000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 35] linear = 1 [ 35] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 8.884926E-03 1.110616E-06 4.136821E-03
--> max chng: dpl= 7.8832E+01 dpo= 7.8832E+01 dso= 8.5220E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 8.992846129731e-03
0 2r: 8.99E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 4.14E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.312797469131e-02
1 2r: 1.31E-02 2x: 3.70E+09 2u: 2.77E+02 ir: 7.21E-03 iu: 4.05E+01 rsn: stol
Step 36 Time= 2.20000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 36] linear = 1 [ 36] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.312797E-02 1.640997E-06 7.208984E-03
--> max chng: dpl= 4.0465E+01 dpo= 4.0465E+01 dso= 7.0028E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.318781013971e-02
0 2r: 1.32E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 7.21E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.163379421990e-03
1 2r: 7.16E-03 2x: 3.70E+09 2u: 3.78E+02 ir: 3.23E-03 iu: 4.05E+01 rsn: stol
Step 37 Time= 2.30000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 37] linear = 1 [ 37] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 7.163379E-03 8.954224E-07 3.228652E-03
--> max chng: dpl= 4.0491E+01 dpo= 4.0491E+01 dso= 8.5386E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 7.321886929978e-03
0 2r: 7.32E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.23E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.078340134516e-02
1 2r: 1.08E-02 2x: 3.70E+09 2u: 2.22E+02 ir: 5.84E-03 iu: 3.18E+01 rsn: stol
Step 38 Time= 2.40000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 38] linear = 1 [ 38] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.078340E-02 1.347925E-06 5.838300E-03
--> max chng: dpl= 3.1849E+01 dpo= 3.1849E+01 dso= 9.2871E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.087492142893e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.84E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.020760531122e-02
1 2r: 1.02E-02 2x: 3.70E+09 2u: 4.45E+02 ir: 4.65E-03 iu: 3.90E+01 rsn: stol
Step 39 Time= 2.50000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 39] linear = 1 [ 39] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.020761E-02 1.275951E-06 4.650045E-03
--> max chng: dpl= 3.9000E+01 dpo= 3.9000E+01 dso= 7.0589E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.028891472766e-02
0 2r: 1.03E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 4.65E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.777786098874e-03
1 2r: 9.78E-03 2x: 3.70E+09 2u: 2.58E+02 ir: 5.10E-03 iu: 3.62E+01 rsn: stol
Step 40 Time= 2.60000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 40] linear = 1 [ 40] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 9.777786E-03 1.222223E-06 5.097846E-03
--> max chng: dpl= 3.6234E+01 dpo= 3.6234E+01 dso= 5.6084E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 9.854784585703e-03
0 2r: 9.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.10E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.884898324061e-03
1 2r: 9.88E-03 2x: 3.70E+09 2u: 4.35E+02 ir: 4.43E-03 iu: 3.67E+01 rsn: stol
Step 41 Time= 2.70000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 41] linear = 1 [ 41] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 9.884898E-03 1.235612E-06 4.428346E-03
--> max chng: dpl= 3.6729E+01 dpo= 3.6729E+01 dso= 5.8293E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.00 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 9.969598186269e-03
0 2r: 9.97E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 4.43E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.085341335350e-02
1 2r: 1.09E-02 2x: 3.70E+09 2u: 3.00E+02 ir: 5.59E-03 iu: 3.71E+01 rsn: stol
Step 42 Time= 2.80000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 42] linear = 1 [ 42] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 1.085341E-02 1.356677E-06 5.587671E-03
--> max chng: dpl= 3.7125E+01 dpo= 3.7125E+01 dso= 6.4964E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.092752348935e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.59E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.059242556915e-04
1 2r: 7.06E-04 2x: 3.70E+09 2u: 9.45E+01 ir: 4.82E-04 iu: 1.76E+01 rsn: stol
Step 43 Time= 2.80020E+01 Dt= 2.00000E-03 [d] snes_conv_reason: 4
newton = 1 [ 43] linear = 1 [ 43] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 1 / 1
--> SNES Residual: 7.059243E-04 8.824053E-08 4.815208E-04
--> max chng: dpl= 1.7620E+01 dpo= 1.7620E+01 dso= 1.2870E-04 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
FLOW TS BE steps = 43 newton = 43 linear = 43 cuts = 0
FLOW TS BE Wasted Linear Iterations = 0
FLOW TS BE SNES time = 197.0 seconds
------------------------------ Provenance --------------------------------------
pflotran_compile_date_time = unknown
pflotran_compile_user = unknown
pflotran_compile_hostname = unknown
pflotran_changeset = unknown
pflotran_status = unknown
petsc_changeset = unknown
petsc_status = unknown
--------------------------------------------------------------------------------
PMTOilImsCreate()
PMTOilImsCreate()
"grid_structured_type" set to default value.
pflotran card:: DEBUG
pflotran card:: GRID
pflotran card:: TIMESTEPPER
pflotran card:: NEWTON_SOLVER
pflotran card:: LINEAR_SOLVER
pflotran card:: TIME
pflotran card:: OUTPUT
"OUTPUT,FORMAT,HDF5,# FILES" set to default value.
"OUTPUT,MASS_BALANCE,DETAILED" set to default value.
pflotran card:: MATERIAL_PROPERTY
Name :: formation
"MATERIAL_PROPERTY,rock density units" set to default value.
"MATERIAL_PROPERTY,specific heat units" set to default value.
"MATERIAL_PROPERTY,dry thermal conductivity units" set to default value.
"MATERIAL_PROPERTY,wet thermal conductivity units" set to default value.
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: DATASET
pflotran card:: CHARACTERISTIC_CURVES
Name :: ch1
pflotran card:: EOS
"EOS,OIL,VISCOSITY,CONSTANT units" set to default value.
pflotran card:: REGION
all
pflotran card:: REGION
inj_well
pflotran card:: REGION
prod_well
pflotran card:: REGION
prod_well_west
pflotran card:: REGION
south_boundary
pflotran card:: FLOW_CONDITION
initial_press
"PRESSURE UNITS" set to default value.
"OIL_SATURATION UNITS" set to default value.
"DATUM UNITS" set to default value.
"OWC UNITS" set to default value.
pflotran card:: FLOW_CONDITION
wat_injection
pflotran card:: INITIAL_CONDITION
pflotran card:: BOUNDARY_CONDITION
pflotran card:: SOURCE_SINK
pflotran card:: STRATA
pflotran card:: END_SUBSURFACE
Domain Bounds (x y z):
0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
1.8250000000E+02 3.3500000000E+02 6.1000000000E+00
rank= 0, nlmax= 4000, nlx,y,z= 20 20 10, nxs,e = 0 20, nys,e = 0 20, nzs,e = 0 10
rank= 1, nlmax= 4000, nlx,y,z= 20 20 10, nxs,e = 0 20, nys,e = 0 20, nzs,e = 10 20
rank= 1, ngmax= 4400, ngx,y,z= 20 20 11, ngxs,e= 0 20, ngys,e= 0 20, ngzs,e= 9 20
rank= 0, ngmax= 4400, ngx,y,z= 20 20 11, ngxs,e= 0 20, ngys,e= 0 20, ngzs,e= 0 11
Testing group: Porosity
Group "Porosity" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityX
Group "PermeabilityX" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityY
Group "PermeabilityY" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Testing group: PermeabilityZ
Group "PermeabilityZ" in HDF5 file "perm_poro_20x20x20.h5" not found in file. Therefore, assuming a cell-indexed dataset.
Opening hdf5 file: mat_ids_20x20x20.h5
Setting up grid cell indices
Opening group: Materials
Reading dataset: Materials/Cell Ids
0.00 Seconds to set up indices
Reading dataset: Materials/Material Ids
0.00 Seconds to read integer array.
Closing group: Materials
Closing hdf5 file: mat_ids_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityX
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityY
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: PermeabilityZ
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Opening hdf5 file: perm_poro_20x20x20.h5
Setting up grid cell indices
Reading dataset: Cell Ids
0.00 Seconds to set up indices
Reading dataset: Porosity
0.00 Seconds to read real array
Closing hdf5 file: perm_poro_20x20x20.h5
Default HDF5 method is used in Initialization
Requested processors and decomposition = 2, npx,y,z= -1 -1 -1
Actual decomposition: npx,y,z= 1 1 2
Beginning setup of FLOW SNES
number of dofs = 3, number of phases = 2
Finished setting up FLOW SNES
Running in CPR-AMG Two-Stage Preconditioning Mode.
Using FGMRES
Using T1: PCHYPRE AMG
Using T2: PCNONE
Solver: fgmres
Preconditioner:
FLOW Time Stepper
max steps: 999999
max constant cumulative time steps: 5
max cuts: 16
ts reduction factor: 0.50000000000000000
FLOW Newton Solver
atol: 1.0000E-05
rtol: 1.0000E-05
stol: 1.0000E-05
dtol: 1.0000E+15
maxnorm: 1.0000E+20
inftolres: -9.9900E+02
inftolupd: -9.9900E+02
inftolrelupd: -9.9900E+02
inftolsclres: -9.9900E+02
max iter: 8
min iter: 1
maxf: 10000
matrix type:baij
precond. matrix type:baij
inexact newton: off
print convergence: on
print detailed convergence: off
check infinity norm: on
FLOW Linear Solver
solver: fgmres
precond:
atol: 1.0000E-50
rtol: 1.0000E-07
dtol: 1.0000E+04
max iter: 200
--> creating hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
Finished Initialization
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.168250330965e-03
0 2r: 6.17E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.325430353299e-09
1 2r: 9.33E-09 2x: 3.70E+09 2u: 1.84E+04 ir: 4.52E-09 iu: 1.69E+03 rsn: atol
Step 1 Time= 1.00000E-04 Dt= 1.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 1] linear = 13 [ 13] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 13 / 1
--> SNES Residual: 9.325430E-09 1.165679E-12 4.523959E-09
--> max chng: dpl= 1.6850E+03 dpo= 1.6850E+03 dso= 7.5593E-05 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.204655900609e-03
0 2r: 6.20E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.243919583361e-08
1 2r: 1.24E-08 2x: 3.70E+09 2u: 2.66E+04 ir: 5.54E-09 iu: 9.79E+02 rsn: atol
Step 2 Time= 3.00000E-04 Dt= 2.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 2] linear = 14 [ 27] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 14 / 1
--> SNES Residual: 1.243920E-08 1.554899E-12 5.542073E-09
--> max chng: dpl= 9.7897E+02 dpo= 9.7897E+02 dso= 1.5119E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205008270050e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.461227665453e-08
1 2r: 1.46E-08 2x: 3.70E+09 2u: 3.99E+04 ir: 4.04E-09 iu: 1.05E+03 rsn: atol
Step 3 Time= 7.00000E-04 Dt= 4.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 3] linear = 15 [ 42] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 15 / 1
--> SNES Residual: 1.461228E-08 1.826535E-12 4.043750E-09
--> max chng: dpl= 1.0535E+03 dpo= 1.0535E+03 dso= 3.0237E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205059166499e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.642184286043e-08
1 2r: 1.64E-08 2x: 3.70E+09 2u: 5.34E+04 ir: 6.20E-09 iu: 1.18E+03 rsn: atol
Step 4 Time= 1.50000E-03 Dt= 8.00000E-04 [d] snes_conv_reason: 2
newton = 1 [ 4] linear = 16 [ 58] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 16 / 1
--> SNES Residual: 1.642184E-08 2.052730E-12 6.195889E-09
--> max chng: dpl= 1.1776E+03 dpo= 1.1776E+03 dso= 6.0475E-04 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205088543070e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.651570825449e-08
1 2r: 1.65E-08 2x: 3.70E+09 2u: 5.71E+04 ir: 6.30E-09 iu: 1.18E+03 rsn: atol
Step 5 Time= 3.10000E-03 Dt= 1.60000E-03 [d] snes_conv_reason: 2
newton = 1 [ 5] linear = 17 [ 75] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 17 / 1
--> SNES Residual: 1.651571E-08 2.064464E-12 6.295212E-09
--> max chng: dpl= 1.1798E+03 dpo= 1.1798E+03 dso= 1.2095E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205110269781e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.216408210697e-08
1 2r: 1.22E-08 2x: 3.70E+09 2u: 4.40E+04 ir: 8.94E-09 iu: 8.97E+02 rsn: atol
Step 6 Time= 6.30000E-03 Dt= 3.20000E-03 [d] snes_conv_reason: 2
newton = 1 [ 6] linear = 17 [ 92] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 17 / 1
--> SNES Residual: 1.216408E-08 1.520510E-12 8.935987E-09
--> max chng: dpl= 8.9737E+02 dpo= 8.9737E+02 dso= 2.4190E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205123875770e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 6.725955941593e-09
1 2r: 6.73E-09 2x: 3.70E+09 2u: 2.28E+04 ir: 5.89E-09 iu: 5.16E+02 rsn: atol
Step 7 Time= 1.27000E-02 Dt= 6.40000E-03 [d] snes_conv_reason: 2
newton = 1 [ 7] linear = 17 [ 109] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 17 / 1
--> SNES Residual: 6.725956E-09 8.407445E-13 5.885917E-09
--> max chng: dpl= 5.1587E+02 dpo= 5.1587E+02 dso= 4.8380E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205129604269e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.444118477586e-09
1 2r: 1.44E-09 2x: 3.70E+09 2u: 7.52E+03 ir: 1.16E-09 iu: 2.20E+02 rsn: atol
Step 8 Time= 2.55000E-02 Dt= 1.28000E-02 [d] snes_conv_reason: 2
newton = 1 [ 8] linear = 17 [ 126] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 17 / 1
--> SNES Residual: 1.444118E-09 1.805148E-13 1.162573E-09
--> max chng: dpl= 2.2050E+02 dpo= 2.2050E+02 dso= 9.6760E-03 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131043462e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.051444047816e-10
1 2r: 2.05E-10 2x: 3.70E+09 2u: 1.45E+03 ir: 5.68E-11 iu: 4.90E+01 rsn: atol
Step 9 Time= 5.11000E-02 Dt= 2.56000E-02 [d] snes_conv_reason: 2
newton = 1 [ 9] linear = 17 [ 143] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 17 / 1
--> SNES Residual: 2.051444E-10 2.564305E-14 5.678594E-11
--> max chng: dpl= 4.9019E+01 dpo= 4.9019E+01 dso= 1.9352E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131241663e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 6.903030294969e-11
1 2r: 6.90E-11 2x: 3.70E+09 2u: 1.51E+02 ir: 5.50E-11 iu: 5.52E+00 rsn: atol
Step 10 Time= 1.02300E-01 Dt= 5.12000E-02 [d] snes_conv_reason: 2
newton = 1 [ 10] linear = 16 [ 159] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 16 / 1
--> SNES Residual: 6.903030E-11 8.628788E-15 5.503347E-11
--> max chng: dpl= 5.5228E+00 dpo= 5.5228E+00 dso= 3.8704E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131263461e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 5.952716875397e-11
1 2r: 5.95E-11 2x: 3.70E+09 2u: 8.22E+00 ir: 5.75E-11 iu: 3.12E-01 rsn: atol
Step 11 Time= 2.04700E-01 Dt= 1.02400E-01 [d] snes_conv_reason: 2
newton = 1 [ 11] linear = 16 [ 175] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 16 / 1
--> SNES Residual: 5.952717E-11 7.440896E-15 5.751489E-11
--> max chng: dpl= 3.1212E-01 dpo= 3.1212E-01 dso= 7.7408E-02 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.205131263197e-03
0 2r: 6.21E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.99E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.151916567332e-03
1 2r: 1.15E-03 2x: 3.70E+09 2u: 2.44E-01 ir: 1.11E-03 iu: 1.25E-01 rsn: stol
Step 12 Time= 4.09500E-01 Dt= 2.04800E-01 [d] snes_conv_reason: 4
newton = 1 [ 12] linear = 16 [ 191] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 16 / 1
--> SNES Residual: 1.151917E-03 1.439896E-07 1.112149E-03
--> max chng: dpl= 7.1157E-03 dpo= 7.1157E-03 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 6.161154480856e-03
0 2r: 6.16E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.95E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.519653048339e-03
1 2r: 3.52E-03 2x: 3.70E+09 2u: 2.00E+02 ir: 2.64E-03 iu: 4.73E+01 rsn: stol
Step 13 Time= 8.19100E-01 Dt= 4.09600E-01 [d] snes_conv_reason: 4
newton = 1 [ 13] linear = 19 [ 210] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 3.519653E-03 4.399566E-07 2.637128E-03
--> max chng: dpl= 4.7302E+01 dpo= 4.7302E+01 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.116656042806e-03
0 2r: 4.12E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.12E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 8.875575019209e-04
1 2r: 8.88E-04 2x: 3.70E+09 2u: 4.35E+02 ir: 6.89E-04 iu: 9.73E+01 rsn: stol
Step 14 Time= 1.00000E+00 Dt= 1.80900E-01 [d] snes_conv_reason: 4
newton = 1 [ 14] linear = 19 [ 229] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 8.875575E-04 1.109447E-07 6.891909E-04
--> max chng: dpl= 9.7335E+01 dpo= 9.7335E+01 dso= 5.0484E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 3.845220304264e-03
0 2r: 3.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.32E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.509047745658e-03
1 2r: 3.51E-03 2x: 3.70E+09 2u: 4.78E+01 ir: 2.12E-03 iu: 2.35E+01 rsn: stol
Step 15 Time= 1.81920E+00 Dt= 8.19200E-01 [d] snes_conv_reason: 4
newton = 1 [ 15] linear = 19 [ 248] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 3.509048E-03 4.386310E-07 2.123847E-03
--> max chng: dpl= 2.3511E+01 dpo= 2.3511E+01 dso= 1.2500E-01 dt= 0.0000E+00
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.846442794840e-03
0 2r: 4.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.73E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.523618671110e-04
1 2r: 3.52E-04 2x: 3.70E+09 2u: 6.18E+01 ir: 2.74E-04 iu: 1.82E+01 rsn: stol
Step 16 Time= 2.00000E+00 Dt= 1.80800E-01 [d] snes_conv_reason: 4
newton = 1 [ 16] linear = 19 [ 267] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 3.523619E-04 4.404523E-08 2.741023E-04
--> max chng: dpl= 1.8170E+01 dpo= 1.8170E+01 dso= 5.4531E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 3.858579148514e-03
0 2r: 3.86E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.80E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.972528041143e-02
1 2r: 2.97E-02 2x: 3.70E+09 2u: 4.45E+01 ir: 2.15E-02 iu: 2.59E+01 rsn: stol
Step 17 Time= 3.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 17] linear = 19 [ 286] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 2.972528E-02 3.715660E-06 2.148937E-02
--> max chng: dpl= 2.5863E+01 dpo= 2.5863E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.956664877583e-02
0 2r: 2.96E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 2.17E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.551770937239e-02
1 2r: 1.55E-02 2x: 3.70E+09 2u: 2.37E+02 ir: 1.14E-02 iu: 4.75E+01 rsn: stol
Step 18 Time= 4.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 18] linear = 19 [ 305] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.551771E-02 1.939714E-06 1.142376E-02
--> max chng: dpl= 4.7548E+01 dpo= 4.7548E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.604560280156e-02
0 2r: 1.60E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 1.18E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.723562703706e-03
1 2r: 9.72E-03 2x: 3.70E+09 2u: 2.35E+02 ir: 6.24E-03 iu: 3.69E+01 rsn: stol
Step 19 Time= 5.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 19] linear = 19 [ 324] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 9.723563E-03 1.215445E-06 6.241646E-03
--> max chng: dpl= 3.6944E+01 dpo= 3.6944E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.013360787586e-02
0 2r: 1.01E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 6.50E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.526798894974e-03
1 2r: 3.53E-03 2x: 3.70E+09 2u: 1.12E+03 ir: 1.80E-03 iu: 5.33E+02 rsn: stol
Step 20 Time= 6.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 20] linear = 19 [ 343] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 3.526799E-03 4.408499E-07 1.796413E-03
--> max chng: dpl= 5.3308E+02 dpo= 5.3308E+02 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.496250350068e-03
0 2r: 4.50E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.19E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.413595302532e-02
1 2r: 1.41E-02 2x: 3.70E+09 2u: 5.63E+02 ir: 9.96E-03 iu: 2.92E+02 rsn: stol
Step 21 Time= 7.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 21] linear = 18 [ 361] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.413595E-02 1.766994E-06 9.961779E-03
--> max chng: dpl= 2.9156E+02 dpo= 2.9156E+02 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.435499001253e-02
0 2r: 1.44E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 9.96E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 4.128435499262e-03
1 2r: 4.13E-03 2x: 3.70E+09 2u: 3.31E+02 ir: 1.81E-03 iu: 5.61E+01 rsn: stol
Step 22 Time= 8.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 22] linear = 19 [ 380] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 4.128435E-03 5.160544E-07 1.814999E-03
--> max chng: dpl= 5.6095E+01 dpo= 5.6095E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.541651828913e-03
0 2r: 4.54E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 1.82E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.270751427973e-03
1 2r: 1.27E-03 2x: 3.70E+09 2u: 7.56E+02 ir: 6.38E-04 iu: 8.87E+01 rsn: stol
Step 23 Time= 9.00000E+00 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 23] linear = 19 [ 399] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.270751E-03 1.588439E-07 6.379279E-04
--> max chng: dpl= 8.8712E+01 dpo= 8.8712E+01 dso= 1.2359E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.519109846080e-03
0 2r: 2.52E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 1.06E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 2.580892914792e-02
1 2r: 2.58E-02 2x: 3.70E+09 2u: 2.54E+02 ir: 1.39E-02 iu: 5.18E+01 rsn: stol
Step 24 Time= 1.00000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 24] linear = 19 [ 418] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 2.580893E-02 3.226116E-06 1.392543E-02
--> max chng: dpl= 5.1761E+01 dpo= 5.1761E+01 dso= 1.2500E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.584984945943e-02
0 2r: 2.58E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 1.39E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.576652114662e-03
1 2r: 1.58E-03 2x: 3.70E+09 2u: 3.98E+02 ir: 1.04E-03 iu: 8.85E+01 rsn: stol
Step 25 Time= 1.10000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 25] linear = 19 [ 437] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.576652E-03 1.970815E-07 1.039702E-03
--> max chng: dpl= 8.8513E+01 dpo= 8.8513E+01 dso= 1.1914E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 2.124856560488e-03
0 2r: 2.12E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 7.46E-04 iu: 0.00E+00 rsn: 0
1 SNES Function norm 4.708769524644e-02
1 2r: 4.71E-02 2x: 3.70E+09 2u: 4.49E+02 ir: 2.40E-02 iu: 8.07E+01 rsn: stol
Step 26 Time= 1.20000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 26] linear = 19 [ 456] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 4.708770E-02 5.885962E-06 2.399306E-02
--> max chng: dpl= 8.0734E+01 dpo= 8.0734E+01 dso= 1.0041E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.707315176730e-02
0 2r: 4.71E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 2.40E-02 iu: 0.00E+00 rsn: 0
1 SNES Function norm 3.891651360888e-03
1 2r: 3.89E-03 2x: 3.70E+09 2u: 7.12E+02 ir: 2.25E-03 iu: 9.69E+01 rsn: stol
Step 27 Time= 1.30000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 27] linear = 19 [ 475] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 3.891651E-03 4.864564E-07 2.246017E-03
--> max chng: dpl= 9.6869E+01 dpo= 9.6869E+01 dso= 1.1124E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 4.406368552447e-03
0 2r: 4.41E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 2.25E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.354255478633e-02
1 2r: 1.35E-02 2x: 3.70E+09 2u: 2.01E+02 ir: 6.28E-03 iu: 3.51E+01 rsn: stol
Step 28 Time= 1.40000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 28] linear = 18 [ 493] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.354255E-02 1.692819E-06 6.278871E-03
--> max chng: dpl= 3.5144E+01 dpo= 3.5144E+01 dso= 1.2185E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.360153647921e-02
0 2r: 1.36E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 6.28E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.643998042052e-02
1 2r: 1.64E-02 2x: 3.70E+09 2u: 5.75E+02 ir: 8.49E-03 iu: 6.06E+01 rsn: stol
Step 29 Time= 1.50000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 29] linear = 19 [ 512] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.643998E-02 2.054998E-06 8.485402E-03
--> max chng: dpl= 6.0570E+01 dpo= 6.0570E+01 dso= 7.6407E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.649483162958e-02
0 2r: 1.65E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 8.49E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.928108036313e-02
1 2r: 1.93E-02 2x: 3.70E+09 2u: 2.64E+02 ir: 9.28E-03 iu: 5.80E+01 rsn: stol
Step 30 Time= 1.60000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 30] linear = 18 [ 530] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.928108E-02 2.410135E-06 9.276394E-03
--> max chng: dpl= 5.7996E+01 dpo= 5.7996E+01 dso= 1.0039E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.932528036104e-02
0 2r: 1.93E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 9.28E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.255733981009e-03
1 2r: 7.26E-03 2x: 3.70E+09 2u: 4.53E+02 ir: 3.54E-03 iu: 5.03E+01 rsn: stol
Step 31 Time= 1.70000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 31] linear = 19 [ 549] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 7.255734E-03 9.069667E-07 3.541171E-03
--> max chng: dpl= 5.0259E+01 dpo= 5.0259E+01 dso= 1.1286E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 7.415803723865e-03
0 2r: 7.42E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.54E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.576833677600e-02
1 2r: 1.58E-02 2x: 3.70E+09 2u: 2.19E+02 ir: 7.18E-03 iu: 3.47E+01 rsn: stol
Step 32 Time= 1.80000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 32] linear = 18 [ 567] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.576834E-02 1.971042E-06 7.175858E-03
--> max chng: dpl= 3.4725E+01 dpo= 3.4725E+01 dso= 6.5864E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.581816001534e-02
0 2r: 1.58E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 7.18E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.069571155903e-02
1 2r: 1.07E-02 2x: 3.70E+09 2u: 4.68E+02 ir: 5.05E-03 iu: 4.77E+01 rsn: stol
Step 33 Time= 1.90000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 33] linear = 19 [ 586] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.069571E-02 1.336964E-06 5.053663E-03
--> max chng: dpl= 4.7740E+01 dpo= 4.7740E+01 dso= 9.2535E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.080945314680e-02
0 2r: 1.08E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.05E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.075524867245e-02
1 2r: 1.08E-02 2x: 3.70E+09 2u: 1.60E+02 ir: 4.93E-03 iu: 3.39E+01 rsn: stol
Step 34 Time= 2.00000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 34] linear = 18 [ 604] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.075525E-02 1.344406E-06 4.933020E-03
--> max chng: dpl= 3.3912E+01 dpo= 3.3912E+01 dso= 1.0311E-01 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.087288017689e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 4.93E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 8.885129226734e-03
1 2r: 8.89E-03 2x: 3.70E+09 2u: 5.42E+02 ir: 4.14E-03 iu: 7.88E+01 rsn: stol
Step 35 Time= 2.10000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 35] linear = 19 [ 623] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 8.885129E-03 1.110641E-06 4.136924E-03
--> max chng: dpl= 7.8833E+01 dpo= 7.8833E+01 dso= 8.5220E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 8.993045622162e-03
0 2r: 8.99E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 4.14E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.312805097768e-02
1 2r: 1.31E-02 2x: 3.70E+09 2u: 2.77E+02 ir: 7.21E-03 iu: 4.05E+01 rsn: stol
Step 36 Time= 2.20000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 36] linear = 18 [ 641] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.312805E-02 1.641006E-06 7.209029E-03
--> max chng: dpl= 4.0465E+01 dpo= 4.0465E+01 dso= 7.0028E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.318788580896e-02
0 2r: 1.32E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 7.21E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.163073827835e-03
1 2r: 7.16E-03 2x: 3.70E+09 2u: 3.78E+02 ir: 3.23E-03 iu: 4.05E+01 rsn: stol
Step 37 Time= 2.30000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 37] linear = 19 [ 660] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 7.163074E-03 8.953842E-07 3.228516E-03
--> max chng: dpl= 4.0490E+01 dpo= 4.0490E+01 dso= 8.5386E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 7.321589835161e-03
0 2r: 7.32E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 3.23E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.078365787178e-02
1 2r: 1.08E-02 2x: 3.70E+09 2u: 2.22E+02 ir: 5.84E-03 iu: 3.18E+01 rsn: stol
Step 38 Time= 2.40000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 38] linear = 18 [ 678] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.078366E-02 1.347957E-06 5.838427E-03
--> max chng: dpl= 3.1849E+01 dpo= 3.1849E+01 dso= 9.2871E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.087517466033e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.84E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.020732554613e-02
1 2r: 1.02E-02 2x: 3.70E+09 2u: 4.45E+02 ir: 4.65E-03 iu: 3.90E+01 rsn: stol
Step 39 Time= 2.50000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 39] linear = 19 [ 697] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 1.020733E-02 1.275916E-06 4.649910E-03
--> max chng: dpl= 3.9000E+01 dpo= 3.9000E+01 dso= 7.0589E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.028863831306e-02
0 2r: 1.03E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 4.65E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.777863141876e-03
1 2r: 9.78E-03 2x: 3.70E+09 2u: 2.58E+02 ir: 5.10E-03 iu: 3.62E+01 rsn: stol
Step 40 Time= 2.60000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 40] linear = 18 [ 715] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 9.777863E-03 1.222233E-06 5.097880E-03
--> max chng: dpl= 3.6234E+01 dpo= 3.6234E+01 dso= 5.6084E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 9.854860649931e-03
0 2r: 9.85E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 5.10E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 9.884825033543e-03
1 2r: 9.88E-03 2x: 3.70E+09 2u: 4.35E+02 ir: 4.43E-03 iu: 3.67E+01 rsn: stol
Step 41 Time= 2.70000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 41] linear = 19 [ 734] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 19 / 1
--> SNES Residual: 9.884825E-03 1.235603E-06 4.428311E-03
--> max chng: dpl= 3.6729E+01 dpo= 3.6729E+01 dso= 5.8293E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 9.969525868361e-03
0 2r: 9.97E-03 2x: 0.00E+00 2u: 0.00E+00 ir: 4.43E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 1.085371422476e-02
1 2r: 1.09E-02 2x: 3.70E+09 2u: 3.00E+02 ir: 5.59E-03 iu: 3.71E+01 rsn: stol
Step 42 Time= 2.80000E+01 Dt= 1.00000E+00 [d] snes_conv_reason: 4
newton = 1 [ 42] linear = 18 [ 752] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 1.085371E-02 1.356714E-06 5.587821E-03
--> max chng: dpl= 3.7125E+01 dpo= 3.7125E+01 dso= 6.4964E-02 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
== TOIL_IMS FLOW ================================================================
0 SNES Function norm 1.092782102539e-02
0 2r: 1.09E-02 2x: 0.00E+00 2u: 0.00E+00 ir: 5.59E-03 iu: 0.00E+00 rsn: 0
1 SNES Function norm 7.059345949085e-04
1 2r: 7.06E-04 2x: 3.70E+09 2u: 9.45E+01 ir: 4.82E-04 iu: 1.76E+01 rsn: stol
Step 43 Time= 2.80020E+01 Dt= 2.00000E-03 [d] snes_conv_reason: 4
newton = 1 [ 43] linear = 18 [ 770] cuts = 0 [ 0]
--> SNES Linear/Non-Linear Iterations = 18 / 1
--> SNES Residual: 7.059346E-04 8.824182E-08 4.815269E-04
--> max chng: dpl= 1.7620E+01 dpo= 1.7620E+01 dso= 1.2871E-04 dt= 0.0000E+00
--> appending to hdf5 output file: het.h5
0.01 Seconds to write HDF5 file.
FLOW TS BE steps = 43 newton = 43 linear = 770 cuts = 0
FLOW TS BE Wasted Linear Iterations = 0
FLOW TS BE SNES time = 68.7 seconds
