# [petsc-users] A number of questions about DMDA with SNES and Quasi-Newton methods

Barry Smith bsmith at mcs.anl.gov
Sun Aug 27 20:33:49 CDT 2017

```> On Aug 27, 2017, at 6:51 PM, zakaryah . <zakaryah at gmail.com> wrote:
>
> OK, I'm with you.  Here's my call to DMDACreate:
>
> ierr = DMDACreate3d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE, DMDA_STENCIL_BOX, mp->W, mp->H, mp->D,PETSC_DECIDE, PETSC_DECIDE, PETSC_DECIDE, 3, 1, NULL, NULL, NULL, &(mp->PETSc_da));CHKERRQ(ierr);
>
> The grid dimensions are mp->W x mp->H x mp->D, and there are three degrees of freedom at each grid location (x, y, z terms of field).  In my Jacobian routine, at the row representing grid coordinate ix, iy, iz, and each of the 3 dof, I set values in columns corresponding to ix+/-1, iy+/-1, iz+/-1, but not all 27 such coordinates (the Jacobian is zero at the 8 corners of the cube).  The code for the Jacobian is explicit at the boundaries.  In other words, the ix-1 terms are only set when ix>0.
>
> I'm not sure how the code can be detecting a new non-zero element in the Jacobian that early in the code.  As I said before, I do not explicitly allocate or set the Jacobian anywhere, following example 48 which also uses a DMDA.

The DMDA preallocates the matrix based on the stencil information you provide. It seems your matrix entries are within the preallocated space; maybe the boundary condition has a wider stencil?

Anyways you need to use the debugger to find out why the code is allocating into a space it shouldn't (based on the DMDA).

Barry

>
>
>
>
> On Sun, Aug 27, 2017 at 7:29 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
> > On Aug 27, 2017, at 6:19 PM, zakaryah . <zakaryah at gmail.com> wrote:
> >
> > Heh, I guess it's at least a sign that the two fd methods agreed.  This takes me back to my earlier question - when I call MatSetOption(jac,MAT_NEW_NONZERO_LOCATION_ERR,PETSC_FALSE), what is jac?  Do I need to preallocate it?  Do I need to get it from the SNES?  Do I need to set that matrix as the SNES Jacobian?  Do I allocate the actual number of non-zero rows, or the entirety of the stencil?  I notice that -snes_test_display reports values for the entire stencil, even though of the 27 points in the stencil, only 19 can be non-zero.
>
>   Back up a minute I didn't read your previous email well enough, sorry. See below
>
> >
> > On Sun, Aug 27, 2017 at 5:51 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> >
> >    Sorry this a flaw with the current release. Call
> >
> > MatSetOption(jac,MAT_NEW_NONZERO_LOCATION_ERR,PETSC_FALSE)
> >
> > before the call to the solver.
> >
> >
> > > On Aug 27, 2017, at 12:54 PM, zakaryah . <zakaryah at gmail.com> wrote:
> > >
> > > Is it suspicious that the KSP converges so quickly?  I have no experience with how the solvers behave on such small grids.
> > >
> > > Also, I ran the code with -snes_type test on a very small grid (1530 elements).  The simpler version of the PDE ran fine, and showed good agreement of the matrices.  However, the more complicated version crashes with the following errors:
> > >
> > > Testing hand-coded Jacobian, if the ratio is
> > >
> > > O(1.e-8), the hand-coded Jacobian is probably correct.
> > >
> > > Run with -snes_test_display to show difference
> > >
> > > of hand-coded and finite difference Jacobian.
> > >
> > > [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> > >
> > > [0]PETSC ERROR: Argument out of range
> > >
> > > [0]PETSC ERROR: Inserting a new nonzero at (7,0) in the matrix
> > >
> > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> > >
> > > [0]PETSC ERROR: Petsc Release Version 3.7.3, Jul, 24, 2016
> > >
> > > [0]PETSC ERROR: #1 MatSetValues_SeqAIJ() line 484 in <PETSC path>/PETSc/build/petsc-3.7.3/src/mat/impls/aij/seq/aij.c
> > >
> > > [0]PETSC ERROR: #2 MatSetValues() line 1190 in <PETSC path>/PETSc/build/petsc-3.7.3/src/mat/interface/matrix.c
> > >
> > > [0]PETSC ERROR: #3 MatSetValuesLocal() line 2053 in <PETSC path>/PETSc/build/petsc-3.7.3/src/mat/interface/matrix.c
> > >
> > > [0]PETSC ERROR: #4 MatSetValuesStencil() line 1447 in <PETSC path>/PETSc/build/petsc-3.7.3/src/mat/interface/matrix.c
> > >
> > > Does this indicate a bug in the Jacobian calculation?  I don't think I set values outside of the stencil...
>
>      Yes something is definitely wrong, you are inserting a value in a location the matrix was not expecting. What are your arguments to DMDACreate()? For stencil width and stencil type? You can run in the debugger to see why the stencil value doesn't "fit".
>
>   Barry
>
>
>
> > >
> > >
> > > On Sun, Aug 27, 2017 at 12:14 AM, zakaryah . <zakaryah at gmail.com> wrote:
> > >
> > > mpiexec -n 1 <exec name> <exec params> -snes_fd_color -snes_type newtonls -snes_monitor -snes_linesearch_monitor -ksp_monitor -pc_type lu
> > >
> > > On Sun, Aug 27, 2017 at 12:11 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > >
> > > > On Aug 26, 2017, at 10:55 PM, zakaryah . <zakaryah at gmail.com> wrote:
> > > >
> > > > Yes, except I changed "-snes_type ls" to "-snes_type newtonls" because PETSc complained that ls wasn't a known method.
> > >
> > >   Sorry, my memory is not as good as it used to be; but what other options did you use? Send the exact command line.
> > >
> > >
> > >    Barry
> > > '
> > > >
> > > > On Sat, Aug 26, 2017 at 11:52 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > > >
> > > >   My exact options did you run with?
> > > >
> > > > > On Aug 26, 2017, at 10:48 PM, zakaryah . <zakaryah at gmail.com> wrote:
> > > > >
> > > > > Here's the output, from the data set which does not converge with l-bfgs:
> > > > >
> > > > >   0 SNES Function norm 1.370293318432e+04
> > > > >     0 KSP Residual norm 7.457506389218e+02
> > > > >     1 KSP Residual norm 2.244764079994e-10
> > > > >       Line search: gnorm after quadratic fit 6.541298279196e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.963717927372e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.788909416707e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.340565762719e+04 lambda=1.2500000000000001e-02
> > > > >   1 SNES Function norm 1.340565762719e+04
> > > > >     0 KSP Residual norm 4.096066553372e+02
> > > > >     1 KSP Residual norm 3.313671167444e-11
> > > > >       Line search: gnorm after quadratic fit 1.113749573695e+06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.406390140546e+05 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.126781917443e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.272988597973e+04 lambda=1.2500000000000001e-02
> > > > >   2 SNES Function norm 1.272988597973e+04
> > > > >     0 KSP Residual norm 8.291344597817e+01
> > > > >     1 KSP Residual norm 7.182258362182e-12
> > > > >       Line search: gnorm after quadratic fit 1.121913743889e+04
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   3 SNES Function norm 1.121913743889e+04
> > > > >     0 KSP Residual norm 6.830535877014e+01
> > > > >     1 KSP Residual norm 3.629826411580e-12
> > > > >       Line search: gnorm after quadratic fit 9.966344973786e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   4 SNES Function norm 9.966344973786e+03
> > > > >     0 KSP Residual norm 5.137691531345e+01
> > > > >     1 KSP Residual norm 1.954502098181e-12
> > > > >       Line search: gnorm after quadratic fit 8.905508641232e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   5 SNES Function norm 8.905508641232e+03
> > > > >     0 KSP Residual norm 4.295019262963e+01
> > > > >     1 KSP Residual norm 1.581527994925e-12
> > > > >       Line search: gnorm after quadratic fit 7.599173694189e+03
> > > > >       Line search: Quadratically determined step, lambda=1.4157226463489950e-01
> > > > >   6 SNES Function norm 7.599173694189e+03
> > > > >     0 KSP Residual norm 3.520472291520e+01
> > > > >     1 KSP Residual norm 1.169994130568e-12
> > > > >       Line search: gnorm after quadratic fit 6.797214843928e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   7 SNES Function norm 6.797214843928e+03
> > > > >     0 KSP Residual norm 2.683523206056e+01
> > > > >     1 KSP Residual norm 9.039858014119e-13
> > > > >       Line search: gnorm after quadratic fit 6.098038917364e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   8 SNES Function norm 6.098038917364e+03
> > > > >     0 KSP Residual norm 2.300710560681e+01
> > > > >     1 KSP Residual norm 1.010402303464e-12
> > > > >       Line search: gnorm after quadratic fit 5.486151385552e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   9 SNES Function norm 5.486151385552e+03
> > > > >     0 KSP Residual norm 2.824628827619e+01
> > > > >     1 KSP Residual norm 1.009589866569e-12
> > > > >       Line search: gnorm after quadratic fit 4.964991021247e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  10 SNES Function norm 4.964991021247e+03
> > > > >     0 KSP Residual norm 6.285926028595e+01
> > > > >     1 KSP Residual norm 2.833546214180e-12
> > > > >       Line search: gnorm after quadratic fit 2.100041391472e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.804051080767e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.893296169579e+03 lambda=2.5000000000000001e-02
> > > > >  11 SNES Function norm 4.893296169579e+03
> > > > >     0 KSP Residual norm 1.688978282804e+01
> > > > >     1 KSP Residual norm 5.927677470786e-13
> > > > >       Line search: gnorm after quadratic fit 4.400894622431e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  12 SNES Function norm 4.400894622431e+03
> > > > >     0 KSP Residual norm 1.481197699600e+01
> > > > >     1 KSP Residual norm 4.522572128105e-13
> > > > >       Line search: Using full step: fnorm 4.400894622431e+03 gnorm 2.094971849007e+03
> > > > >  13 SNES Function norm 2.094971849007e+03
> > > > >     0 KSP Residual norm 3.514164563318e+00
> > > > >     1 KSP Residual norm 3.022385947499e-13
> > > > >       Line search: gnorm after quadratic fit 1.687641025382e+03
> > > > >       Line search: Quadratically determined step, lambda=2.0307116693368590e-01
> > > > >  14 SNES Function norm 1.687641025382e+03
> > > > >     0 KSP Residual norm 2.138591884686e+00
> > > > >     1 KSP Residual norm 4.023500814078e-14
> > > > >       Line search: Using full step: fnorm 1.687641025382e+03 gnorm 1.131934218470e+03
> > > > >  15 SNES Function norm 1.131934218470e+03
> > > > >     0 KSP Residual norm 3.245035110327e+00
> > > > >     1 KSP Residual norm 1.293626215269e-13
> > > > >       Line search: Using full step: fnorm 1.131934218470e+03 gnorm 7.340573607307e+02
> > > > >  16 SNES Function norm 7.340573607307e+02
> > > > >     0 KSP Residual norm 1.957698957904e+00
> > > > >     1 KSP Residual norm 3.444648967078e-13
> > > > >       Line search: Using full step: fnorm 7.340573607307e+02 gnorm 5.024723505168e+02
> > > > >  17 SNES Function norm 5.024723505168e+02
> > > > >     0 KSP Residual norm 2.510538703839e+00
> > > > >     1 KSP Residual norm 8.570440675253e-14
> > > > >       Line search: gnorm after quadratic fit 4.548776456608e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  18 SNES Function norm 4.548776456608e+02
> > > > >     0 KSP Residual norm 6.741705718178e-01
> > > > >     1 KSP Residual norm 1.650556120809e-14
> > > > >       Line search: Using full step: fnorm 4.548776456608e+02 gnorm 1.351189086642e+02
> > > > >  19 SNES Function norm 1.351189086642e+02
> > > > >     0 KSP Residual norm 7.653741241950e+00
> > > > >     1 KSP Residual norm 8.326848236950e-14
> > > > >       Line search: gnorm after quadratic fit 4.215455105006e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.889750369356e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.358691836999e+02 lambda=1.0303015930243279e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.348911963390e+02 lambda=3.5279526770504110e-03
> > > > >  20 SNES Function norm 1.348911963390e+02
> > > > >     0 KSP Residual norm 4.209281176918e+00
> > > > >     1 KSP Residual norm 1.233232881375e-13
> > > > >       Line search: gnorm after quadratic fit 1.860023264470e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.413911132196e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.339597869053e+02 lambda=1.5787957598629023e-02
> > > > >  21 SNES Function norm 1.339597869053e+02
> > > > >     0 KSP Residual norm 1.718484765841e+00
> > > > >     1 KSP Residual norm 6.666016296543e-14
> > > > >       Line search: gnorm after quadratic fit 1.326878521472e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  22 SNES Function norm 1.326878521472e+02
> > > > >     0 KSP Residual norm 1.515373499919e+00
> > > > >     1 KSP Residual norm 2.259154130795e-12
> > > > >       Line search: gnorm after quadratic fit 1.226907316391e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  23 SNES Function norm 1.226907316391e+02
> > > > >     0 KSP Residual norm 1.080252936903e+00
> > > > >     1 KSP Residual norm 1.847020526683e-14
> > > > >       Line search: gnorm after quadratic fit 1.163632111851e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  24 SNES Function norm 1.163632111851e+02
> > > > >     0 KSP Residual norm 2.491746958382e+00
> > > > >     1 KSP Residual norm 6.389134193637e-14
> > > > >       Line search: gnorm after quadratic fit 1.345487153563e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.169738363622e+02 lambda=4.4108900954515480e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.152177638108e+02 lambda=1.9997442889243631e-02
> > > > >  25 SNES Function norm 1.152177638108e+02
> > > > >     0 KSP Residual norm 5.364385501004e+00
> > > > >     1 KSP Residual norm 8.457149562463e-14
> > > > >       Line search: gnorm after quadratic fit 2.722095758964e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.480946353374e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.150680255994e+02 lambda=6.3560128131639167e-03
> > > > >  26 SNES Function norm 1.150680255994e+02
> > > > >     0 KSP Residual norm 4.302025268263e+00
> > > > >     1 KSP Residual norm 1.017526937941e-13
> > > > >       Line search: gnorm after quadratic fit 1.956300983573e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.318600522939e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.147123505014e+02 lambda=7.6060788653548803e-03
> > > > >  27 SNES Function norm 1.147123505014e+02
> > > > >     0 KSP Residual norm 6.195463111324e+00
> > > > >     1 KSP Residual norm 1.235283250361e-13
> > > > >       Line search: gnorm after quadratic fit 3.324821547049e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.612593208504e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.147389525772e+02 lambda=6.1750939181374294e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.145411432123e+02 lambda=3.0078592586792975e-03
> > > > >  28 SNES Function norm 1.145411432123e+02
> > > > >     0 KSP Residual norm 1.454704250187e+01
> > > > >     1 KSP Residual norm 2.368485174123e-13
> > > > >       Line search: gnorm after quadratic fit 1.216293873116e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.796524621727e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.359314163651e+02 lambda=1.4867675803955788e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.146009802557e+02 lambda=1.4867675803955788e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.145096815033e+02 lambda=5.5250912472932839e-04
> > > > >  29 SNES Function norm 1.145096815033e+02
> > > > >     0 KSP Residual norm 3.430451345093e+01
> > > > >     1 KSP Residual norm 1.069396165938e-12
> > > > >       Line search: gnorm after quadratic fit 1.161329407098e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.260690718050e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.696806489042e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.960149915648e+02 lambda=1.1276129246469476e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.158113641658e+02 lambda=1.5873910451430554e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.145062232916e+02 lambda=1.5873910451430555e-04
> > > > >  30 SNES Function norm 1.145062232916e+02
> > > > >     0 KSP Residual norm 2.662578971526e+01
> > > > >     1 KSP Residual norm 5.464011789728e-13
> > > > >       Line search: gnorm after quadratic fit 4.375119254490e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.038018103928e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.169328002874e+02 lambda=2.3789161387159519e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.262835432714e+02 lambda=5.9737415571960795e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.145626045197e+02 lambda=5.9737415571960802e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144968604079e+02 lambda=1.6421773527706333e-04
> > > > >  31 SNES Function norm 1.144968604079e+02
> > > > >     0 KSP Residual norm 6.322033697338e+01
> > > > >     1 KSP Residual norm 1.507157991448e-12
> > > > >       Line search: gnorm after quadratic fit 5.809243699277e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.460487584142e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.893183483197e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.960337007072e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.771527792790e+02 lambda=5.3912931685137378e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.150137639707e+02 lambda=5.3912931685137376e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144964484571e+02 lambda=5.3912931685137380e-05
> > > > >  32 SNES Function norm 1.144964484571e+02
> > > > >     0 KSP Residual norm 3.859510930996e+01
> > > > >     1 KSP Residual norm 8.069118044382e-13
> > > > >       Line search: gnorm after quadratic fit 1.106329930525e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.155884463041e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.933963819094e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.849629797377e+02 lambda=9.7744286136270241e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.150857687050e+02 lambda=9.7744286136270254e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144922906340e+02 lambda=9.7744286136270254e-05
> > > > >  33 SNES Function norm 1.144922906340e+02
> > > > >     0 KSP Residual norm 4.952038022394e+01
> > > > >     1 KSP Residual norm 7.460263245424e-13
> > > > >       Line search: gnorm after quadratic fit 3.005091127220e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.229308416025e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.138600794682e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.390856262238e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.394997214983e+02 lambda=4.4582997750629580e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.146836358799e+02 lambda=4.4582997750629581e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144895966587e+02 lambda=4.7644082443915877e-05
> > > > >  34 SNES Function norm 1.144895966587e+02
> > > > >     0 KSP Residual norm 1.146914786457e+02
> > > > >     1 KSP Residual norm 4.791627042400e-12
> > > > >       Line search: gnorm after quadratic fit 2.601294145944e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.324148513778e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.247478406023e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.200340900661e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.696223112300e+02 lambda=6.1514413845115794e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.342365431983e+02 lambda=1.7502929694284564e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.146686063035e+02 lambda=1.7502929694284566e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144895868489e+02 lambda=1.7502929694284566e-05
> > > > >  35 SNES Function norm 1.144895868489e+02
> > > > >     0 KSP Residual norm 6.311017209658e+01
> > > > >     1 KSP Residual norm 8.384445939117e-13
> > > > >       Line search: gnorm after quadratic fit 5.781533400572e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.419557746031e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.886081770030e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.945810143740e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.767820854837e+02 lambda=5.3859001959289735e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.150034040559e+02 lambda=5.3859001959289732e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144891501665e+02 lambda=5.3859001959289732e-05
> > > > >  36 SNES Function norm 1.144891501665e+02
> > > > >     0 KSP Residual norm 3.880748384332e+01
> > > > >     1 KSP Residual norm 6.400339290453e-13
> > > > >       Line search: gnorm after quadratic fit 1.122504184426e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.180728237366e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.987870621566e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.863711604331e+02 lambda=9.8150771384688824e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.150916783781e+02 lambda=9.8150771384688824e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144850837411e+02 lambda=9.8150771384688832e-05
> > > > >  37 SNES Function norm 1.144850837411e+02
> > > > >     0 KSP Residual norm 4.811537498557e+01
> > > > >     1 KSP Residual norm 9.738167670242e-13
> > > > >       Line search: gnorm after quadratic fit 2.784098355218e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.885118868254e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.074843354348e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.255202820836e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.365202996130e+02 lambda=4.3204204582016374e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.146504919412e+02 lambda=4.3204204582016374e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144822306241e+02 lambda=5.0376546850227848e-05
> > > > >  38 SNES Function norm 1.144822306241e+02
> > > > >     0 KSP Residual norm 1.120914002482e+02
> > > > >     1 KSP Residual norm 5.452125292284e-12
> > > > >       Line search: gnorm after quadratic fit 2.427186680567e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.112196656857e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.967397366072e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.149551659770e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.529630886791e+02 lambda=6.0885216927030559e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.315267519231e+02 lambda=1.6656233536183130e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.146353659250e+02 lambda=1.6656233536183130e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.144820487391e+02 lambda=1.6656233536183130e-05
> > > > >  39 SNES Function norm 1.144820487391e+02
> > > > >     0 KSP Residual norm 7.182416170980e+01
> > > > >     1 KSP Residual norm 1.292147133333e-12
> > > > >       Line search: gnorm after quadratic fit 8.255331293322e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.302939895170e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.501228089911e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.185010378583e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.088179821424e+02 lambda=5.7489510178867142e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.168272901121e+02 lambda=9.7452649444612811e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144951972300e+02 lambda=9.7452649444612822e-05
> > > > >       Line search: Cubically determined step, current gnorm 1.144807676687e+02 lambda=2.2399506502463798e-05
> > > > >  40 SNES Function norm 1.144807676687e+02
> > > > >     0 KSP Residual norm 1.728679810285e+02
> > > > >     1 KSP Residual norm 3.878068639716e-12
> > > > >       Line search: gnorm after quadratic fit 9.011020578535e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.111818830011e+05 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.504789858103e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.754312133162e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.212492111975e+02 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.199969180537e+02 lambda=2.6424184504193135e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.154794639440e+02 lambda=2.6424184504193136e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144880673342e+02 lambda=2.6424184504193136e-05
> > > > >       Line search: Cubically determined step, current gnorm 1.144805461570e+02 lambda=3.8712107591125690e-06
> > > > >  41 SNES Function norm 1.144805461570e+02
> > > > >     0 KSP Residual norm 4.162780846202e+02
> > > > >     1 KSP Residual norm 1.732341544934e-11
> > > > >       Line search: gnorm after quadratic fit 1.355673864369e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.747523002884e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.342205048417e+05 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.425635699680e+04 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.882150659178e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.259664786336e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.653670676209e+02 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.453655830144e+02 lambda=5.8237455565260388e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.147659525018e+02 lambda=5.8237455565260393e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144827915995e+02 lambda=5.8237455565260400e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.144805080017e+02 lambda=6.6689295146509104e-07
> > > > >  42 SNES Function norm 1.144805080017e+02
> > > > >     0 KSP Residual norm 1.004135512581e+03
> > > > >     1 KSP Residual norm 3.867626041000e-09
> > > > >       Line search: gnorm after quadratic fit 1.832578994882e+08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.269698278878e+07 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.792476589915e+06 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.420660371083e+05 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.321686926168e+04 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.548358078234e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.429504802276e+03 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.317723385004e+02 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.460079723095e+02 lambda=2.5078917343692407e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.147908977725e+02 lambda=2.5078917343692408e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144833604238e+02 lambda=2.5078917343692412e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144805106824e+02 lambda=2.5078917343692411e-07
> > > > >       Line search: Cubically determined step, current gnorm 1.144805014337e+02 lambda=1.1468707119027746e-07
> > > > >  43 SNES Function norm 1.144805014337e+02
> > > > >     0 KSP Residual norm 2.418614573592e+03
> > > > >     1 KSP Residual norm 6.085078742847e-10
> > > > >       Line search: gnorm after quadratic fit 2.598476259775e+09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.262759415873e+08 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.116845970403e+07 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.250465130250e+06 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.863493884574e+05 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.501468163771e+04 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.482658980483e+04 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.802395557883e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.788149582374e+02 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.248828337038e+02 lambda=1.8333058767960295e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.186285836222e+02 lambda=3.7550993478654105e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.145209710139e+02 lambda=3.7550993478654107e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144808670668e+02 lambda=3.7550993478654111e-07
> > > > >       Line search: Cubically determined step, current gnorm 1.144805012252e+02 lambda=3.7550993478654114e-08
> > > > >  44 SNES Function norm 1.144805012252e+02
> > > > >     0 KSP Residual norm 1.432476672735e+03
> > > > >     1 KSP Residual norm 7.850524783892e-11
> > > > >       Line search: gnorm after quadratic fit 5.336191593632e+08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.625576866625e+07 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.181408387845e+06 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.003310644422e+06 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.236384273292e+05 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.658430638069e+04 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.977463068161e+03 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.674238499253e+02 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.333616820791e+02 lambda=3.3764776729281175e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.161349153752e+02 lambda=4.0497161764116874e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144966925969e+02 lambda=4.0497161764116874e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144806214839e+02 lambda=4.0497161764116878e-07
> > > > >       Line search: Cubically determined step, current gnorm 1.144804979966e+02 lambda=5.6339112427782378e-08
> > > > >  45 SNES Function norm 1.144804979966e+02
> > > > >     0 KSP Residual norm 3.453401579761e+03
> > > > >     1 KSP Residual norm 4.197968028126e-09
> > > > >       Line search: gnorm after quadratic fit 7.554006183767e+09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.471974940372e+08 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.191613865049e+08 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.509784334249e+07 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.943752478560e+06 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.597601716755e+05 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.775572331075e+04 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.418045407608e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.357066758731e+03 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.859267839080e+02 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.501494912160e+02 lambda=7.5160826909319168e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.148144926477e+02 lambda=7.5160826909319171e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144837498478e+02 lambda=7.5160826909319179e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144805227746e+02 lambda=7.5160826909319182e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.144804974438e+02 lambda=9.6864021663644795e-09
> > > > >  46 SNES Function norm 1.144804974438e+02
> > > > >     0 KSP Residual norm 8.345139428850e+03
> > > > >     1 KSP Residual norm 1.027819754739e-09
> > > > >       Line search: gnorm after quadratic fit 1.061319903906e+11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.325012985379e+10 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.652236226640e+09 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.055533971630e+08 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.546607716763e+07 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.134442858835e+06 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.839168570687e+05 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.830568178626e+04 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.201776786081e+03 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.540520468013e+03 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.573504890506e+02 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.516203908013e+02 lambda=3.2703670177372021e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.148480075676e+02 lambda=3.2703670177372022e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144841475328e+02 lambda=3.2703670177372023e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144805305720e+02 lambda=3.2703670177372021e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.144804974367e+02 lambda=3.2703670177372025e-09
> > > > >  47 SNES Function norm 1.144804974367e+02
> > > > >     0 KSP Residual norm 4.668983500382e+03
> > > > >     1 KSP Residual norm 1.398997665317e-09
> > > > >       Line search: gnorm after quadratic fit 1.865309565532e+10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.336975299727e+09 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.934905088739e+08 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.704520478641e+07 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.728477809448e+06 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.193001670096e+05 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.610673238239e+04 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.354711683605e+04 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.589557246433e+03 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.367851970709e+02 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.136930269795e+02 lambda=9.0316845802227864e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.172186635069e+02 lambda=1.5831613287181393e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.145074017254e+02 lambda=1.5831613287181394e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144807499816e+02 lambda=1.5831613287181396e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144804983345e+02 lambda=1.5831613287181397e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.144804971346e+02 lambda=5.2932921109871310e-09
> > > > >  48 SNES Function norm 1.144804971346e+02
> > > > >     0 KSP Residual norm 1.132678078317e+04
> > > > >     1 KSP Residual norm 2.477518733314e-10
> > > > >       Line search: gnorm after quadratic fit 2.654659022365e+11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.315296223725e+10 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.136635677074e+09 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.152507060235e+08 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.397060094478e+07 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.898362012287e+06 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.685179520906e+05 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.194040779842e+05 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.606415730227e+04 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.902698091972e+03 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.521549384652e+02 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.288999778994e+02 lambda=4.1899746703195281e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.157880829786e+02 lambda=4.5470218451893824e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144935741206e+02 lambda=4.5470218451893828e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144806232638e+02 lambda=4.5470218451893831e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144804979250e+02 lambda=4.5470218451893835e-09
> > > > >       Line search: Cubically determined step, current gnorm 1.144804970825e+02 lambda=9.0293679903918216e-10
> > > > >  49 SNES Function norm 1.144804970825e+02
> > > > >     0 KSP Residual norm 2.712544829144e+04
> > > > >     1 KSP Residual norm 4.949246309677e-09
> > > > >       Line search: gnorm after quadratic fit 3.650846005745e+12
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.565321055911e+11 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.711080201118e+10 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.150022242308e+09 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.965954117985e+08 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.128096393645e+08 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.429702091584e+07 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.841819946536e+06 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.465032956946e+05 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.594210253794e+04 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.141110278944e+03 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.306952279045e+03 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.754164864338e+02 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.476861367158e+02 lambda=9.2451761406722810e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.147929263780e+02 lambda=9.2451761406722810e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144836022952e+02 lambda=9.2451761406722821e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144805271860e+02 lambda=9.2451761406722831e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.144804972895e+02 lambda=9.2451761406722839e-10
> > > > >       Line search: Cubically determined step, current gnorm 1.144804970737e+02 lambda=1.5633616333063305e-10
> > > > >  50 SNES Function norm 1.144804970737e+02
> > > > >   0 SNES Function norm 2.308693894796e+03
> > > > >     0 KSP Residual norm 8.981999720532e+00
> > > > >     1 KSP Residual norm 2.363170936183e-13
> > > > >       Line search: Using full step: fnorm 2.308693894796e+03 gnorm 7.571661187318e+02
> > > > >   1 SNES Function norm 7.571661187318e+02
> > > > >     0 KSP Residual norm 2.149614048903e+00
> > > > >     1 KSP Residual norm 9.511247057888e-13
> > > > >       Line search: Using full step: fnorm 7.571661187318e+02 gnorm 4.450081004997e+02
> > > > >   2 SNES Function norm 4.450081004997e+02
> > > > >     0 KSP Residual norm 1.706469075123e+01
> > > > >     1 KSP Residual norm 5.803815175472e-13
> > > > >       Line search: gnorm after quadratic fit 1.510518198899e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.850796532980e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.515983139755e+02 lambda=2.0602556887689423e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.434800867235e+02 lambda=8.2396279492937211e-03
> > > > >   3 SNES Function norm 4.434800867235e+02
> > > > >     0 KSP Residual norm 3.208587171626e+00
> > > > >     1 KSP Residual norm 1.099072610659e-13
> > > > >       Line search: gnorm after quadratic fit 4.080637194736e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   4 SNES Function norm 4.080637194736e+02
> > > > >     0 KSP Residual norm 8.136557176540e+00
> > > > >     1 KSP Residual norm 8.800137844173e-13
> > > > >       Line search: gnorm after quadratic fit 5.261656549654e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.131114070934e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.031614746044e+02 lambda=2.4674842039461482e-02
> > > > >   5 SNES Function norm 4.031614746044e+02
> > > > >     0 KSP Residual norm 7.609918907567e+00
> > > > >     1 KSP Residual norm 1.613159932655e-13
> > > > >       Line search: gnorm after quadratic fit 5.353908064984e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.110480554132e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 3.991168240716e+02 lambda=2.3079500036735322e-02
> > > > >   6 SNES Function norm 3.991168240716e+02
> > > > >     0 KSP Residual norm 3.800845472041e+00
> > > > >     1 KSP Residual norm 4.841682188278e-13
> > > > >       Line search: gnorm after quadratic fit 3.787886582952e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   7 SNES Function norm 3.787886582952e+02
> > > > >     0 KSP Residual norm 1.892621919219e+00
> > > > >     1 KSP Residual norm 3.576655371800e-12
> > > > >       Line search: gnorm after quadratic fit 3.440181918909e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   8 SNES Function norm 3.440181918909e+02
> > > > >     0 KSP Residual norm 3.559759789147e+00
> > > > >     1 KSP Residual norm 8.347521826622e-13
> > > > >       Line search: gnorm after quadratic fit 3.287993515195e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   9 SNES Function norm 3.287993515195e+02
> > > > >     0 KSP Residual norm 1.825073540507e+00
> > > > >     1 KSP Residual norm 8.352745008123e-14
> > > > >       Line search: Using full step: fnorm 3.287993515195e+02 gnorm 2.710650507416e+02
> > > > >  10 SNES Function norm 2.710650507416e+02
> > > > >     0 KSP Residual norm 7.568432945608e-01
> > > > >     1 KSP Residual norm 2.780715252274e-14
> > > > >       Line search: Using full step: fnorm 2.710650507416e+02 gnorm 1.513359047342e+02
> > > > >  11 SNES Function norm 1.513359047342e+02
> > > > >     0 KSP Residual norm 9.387460484575e+00
> > > > >     1 KSP Residual norm 7.091233040676e-13
> > > > >       Line search: gnorm after quadratic fit 3.998857979851e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.060972049714e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.512140169420e+02 lambda=5.4338709720680610e-03
> > > > >  12 SNES Function norm 1.512140169420e+02
> > > > >     0 KSP Residual norm 4.399424360499e+00
> > > > >     1 KSP Residual norm 1.614362118182e-13
> > > > >       Line search: gnorm after quadratic fit 1.913552832643e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.559719302467e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.499985374733e+02 lambda=1.7102027220313589e-02
> > > > >  13 SNES Function norm 1.499985374733e+02
> > > > >     0 KSP Residual norm 3.740397849742e+00
> > > > >     1 KSP Residual norm 8.986775577006e-13
> > > > >       Line search: gnorm after quadratic fit 1.764118876632e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.522381536119e+02 lambda=4.8196830215713270e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.486175880390e+02 lambda=1.8881300948189364e-02
> > > > >  14 SNES Function norm 1.486175880390e+02
> > > > >     0 KSP Residual norm 6.067973818528e+00
> > > > >     1 KSP Residual norm 4.527845229549e-13
> > > > >       Line search: gnorm after quadratic fit 2.514021299115e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.679972156573e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.481015724304e+02 lambda=9.0116621678703428e-03
> > > > >  15 SNES Function norm 1.481015724304e+02
> > > > >     0 KSP Residual norm 6.108754994263e+00
> > > > >     1 KSP Residual norm 2.557635976440e-13
> > > > >       Line search: gnorm after quadratic fit 2.458081150104e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.681837029397e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.476148340442e+02 lambda=7.9580911933571953e-03
> > > > >  16 SNES Function norm 1.476148340442e+02
> > > > >     0 KSP Residual norm 7.914946163932e+00
> > > > >     1 KSP Residual norm 1.512066885699e-13
> > > > >       Line search: gnorm after quadratic fit 3.458938604598e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.883018868359e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.474286108114e+02 lambda=6.7262521208543979e-03
> > > > >  17 SNES Function norm 1.474286108114e+02
> > > > >     0 KSP Residual norm 5.285449532689e+00
> > > > >     1 KSP Residual norm 6.693733977456e-12
> > > > >       Line search: gnorm after quadratic fit 2.170263102661e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.607560548008e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.467785507822e+02 lambda=9.9244383205188240e-03
> > > > >  18 SNES Function norm 1.467785507822e+02
> > > > >     0 KSP Residual norm 8.124903695635e+00
> > > > >     1 KSP Residual norm 2.322439601127e-11
> > > > >       Line search: gnorm after quadratic fit 3.589716592364e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.907315184877e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.466341888638e+02 lambda=6.5538271222582893e-03
> > > > >  19 SNES Function norm 1.466341888638e+02
> > > > >     0 KSP Residual norm 5.253550718343e+00
> > > > >     1 KSP Residual norm 1.575657996883e-12
> > > > >       Line search: gnorm after quadratic fit 2.155795776725e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.598564001687e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.459869404602e+02 lambda=9.9195822632931301e-03
> > > > >  20 SNES Function norm 1.459869404602e+02
> > > > >     0 KSP Residual norm 8.405987790193e+00
> > > > >     1 KSP Residual norm 2.046159481178e-13
> > > > >       Line search: gnorm after quadratic fit 3.767908298426e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.941885062002e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.458995232755e+02 lambda=6.4359684554015136e-03
> > > > >  21 SNES Function norm 1.458995232755e+02
> > > > >     0 KSP Residual norm 5.036348246593e+00
> > > > >     1 KSP Residual norm 3.645081669075e-13
> > > > >       Line search: gnorm after quadratic fit 2.083222977519e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.575737508694e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.452039802458e+02 lambda=1.0554092389542354e-02
> > > > >  22 SNES Function norm 1.452039802458e+02
> > > > >     0 KSP Residual norm 8.562292355998e+00
> > > > >     1 KSP Residual norm 3.256332645483e-13
> > > > >       Line search: gnorm after quadratic fit 3.869470823355e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.959483128249e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.451523830598e+02 lambda=6.3780526507799607e-03
> > > > >  23 SNES Function norm 1.451523830598e+02
> > > > >     0 KSP Residual norm 4.919667503862e+00
> > > > >     1 KSP Residual norm 4.823931177115e-13
> > > > >       Line search: gnorm after quadratic fit 2.042891039525e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.560623102934e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.444331916523e+02 lambda=1.0890355421223689e-02
> > > > >  24 SNES Function norm 1.444331916523e+02
> > > > >     0 KSP Residual norm 8.727282395369e+00
> > > > >     1 KSP Residual norm 7.090683582186e-13
> > > > >       Line search: gnorm after quadratic fit 3.977929422821e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.978556223200e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.444215965394e+02 lambda=6.3477766966048947e-03
> > > > >  25 SNES Function norm 1.444215965394e+02
> > > > >     0 KSP Residual norm 4.759732971179e+00
> > > > >     1 KSP Residual norm 7.539889498211e-13
> > > > >       Line search: gnorm after quadratic fit 1.990589649135e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.542648241555e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.436629416941e+02 lambda=1.1434720389464151e-02
> > > > >  26 SNES Function norm 1.436629416941e+02
> > > > >     0 KSP Residual norm 8.844861828477e+00
> > > > >     1 KSP Residual norm 4.372001279689e-13
> > > > >       Line search: gnorm after quadratic fit 4.055598972133e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.990820671396e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.436827663113e+02 lambda=6.3302081701296859e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.434397789472e+02 lambda=3.1285674619640730e-03
> > > > >  27 SNES Function norm 1.434397789472e+02
> > > > >     0 KSP Residual norm 2.168630760128e+01
> > > > >     1 KSP Residual norm 3.975838928402e-13
> > > > >       Line search: gnorm after quadratic fit 1.655681371309e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.972331169594e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.804118272840e+02 lambda=1.6710557456633125e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.435928635079e+02 lambda=1.6710557456633126e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.434032742903e+02 lambda=5.1357350159978810e-04
> > > > >  28 SNES Function norm 1.434032742903e+02
> > > > >     0 KSP Residual norm 5.259508821923e+01
> > > > >     1 KSP Residual norm 1.358381204928e-11
> > > > >       Line search: gnorm after quadratic fit 1.908330224731e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.431279702545e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.060945045253e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.791318323447e+02 lambda=1.2124172979783788e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.491140019897e+02 lambda=2.6952165802905347e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434246250245e+02 lambda=2.6952165802905350e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.433970449422e+02 lambda=8.6980371578986910e-05
> > > > >  29 SNES Function norm 1.433970449422e+02
> > > > >     0 KSP Residual norm 1.326287161615e+02
> > > > >     1 KSP Residual norm 5.692176796134e-12
> > > > >       Line search: gnorm after quadratic fit 2.077891749832e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.654187989332e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.213029457851e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.001385334742e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.323976827250e+02 lambda=5.9607661619885998e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.489839529892e+02 lambda=1.0476257410701045e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434396736634e+02 lambda=1.0476257410701046e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.433960671821e+02 lambda=1.3664867646895530e-05
> > > > >  30 SNES Function norm 1.433960671821e+02
> > > > >     0 KSP Residual norm 3.320710360189e+02
> > > > >     1 KSP Residual norm 1.365660123102e-11
> > > > >       Line search: gnorm after quadratic fit 3.597335441013e+06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.714766420450e+05 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.566809766217e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.038660363150e+04 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.026997416957e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.382653551072e+02 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.071747386385e+02 lambda=1.3384948046761360e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.439693866777e+02 lambda=1.3384948046761360e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434000505249e+02 lambda=1.3384948046761361e-05
> > > > >       Line search: Cubically determined step, current gnorm 1.433959111179e+02 lambda=2.1771867000079373e-06
> > > > >  31 SNES Function norm 1.433959111179e+02
> > > > >     0 KSP Residual norm 8.385024273664e+02
> > > > >     1 KSP Residual norm 2.688330817732e-11
> > > > >       Line search: gnorm after quadratic fit 5.487352481970e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.776642694838e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.310551423141e+05 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.023322644608e+05 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.368662233379e+04 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.472194884015e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.718801059897e+02 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.274234972424e+02 lambda=6.3064412238487922e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.442194384053e+02 lambda=6.3064412238487925e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434033578642e+02 lambda=6.3064412238487927e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.433959042208e+02 lambda=6.3064412238487929e-07
> > > > >  32 SNES Function norm 1.433959042208e+02
> > > > >     0 KSP Residual norm 5.310191626867e+02
> > > > >     1 KSP Residual norm 2.270033897601e-10
> > > > >       Line search: gnorm after quadratic fit 1.447633783740e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.859761774309e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.472992503503e+05 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.557594026810e+04 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.969012939380e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.269212161982e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.859005100842e+02 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.695095262046e+02 lambda=5.4509037994531233e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.436389958907e+02 lambda=5.4509037994531237e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433976258223e+02 lambda=5.4509037994531242e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.433958431980e+02 lambda=8.5124820362933632e-07
> > > > >  33 SNES Function norm 1.433958431980e+02
> > > > >     0 KSP Residual norm 1.341142541825e+03
> > > > >     1 KSP Residual norm 4.194815344365e-11
> > > > >       Line search: gnorm after quadratic fit 2.256410317099e+08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.797696259877e+07 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.447237377751e+06 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.221898175722e+05 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.260244723327e+04 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.539148524881e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.558034531173e+03 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.788188892302e+02 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.753621043454e+02 lambda=2.4427125599600652e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.437122397600e+02 lambda=2.4427125599600654e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433986985128e+02 lambda=2.4427125599600655e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.433958402293e+02 lambda=2.4427125599600656e-07
> > > > >  34 SNES Function norm 1.433958402293e+02
> > > > >     0 KSP Residual norm 8.620962950418e+02
> > > > >     1 KSP Residual norm 4.517777375659e-11
> > > > >       Line search: gnorm after quadratic fit 6.134442313460e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.790495969255e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.008365220843e+06 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.364572513478e+05 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.037891335918e+04 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.640571521199e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.472549649319e+02 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.906041216274e+02 lambda=7.6348458761785112e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.504925859329e+02 lambda=1.7740973415567094e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434632637071e+02 lambda=1.7740973415567094e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433962847027e+02 lambda=1.7740973415567095e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.433958170795e+02 lambda=3.2293255704969003e-07
> > > > >  35 SNES Function norm 1.433958170795e+02
> > > > >     0 KSP Residual norm 2.177497039945e+03
> > > > >     1 KSP Residual norm 8.546235181505e-11
> > > > >       Line search: gnorm after quadratic fit 9.689178330938e+08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.204850731541e+08 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.491460480376e+07 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.833699292784e+06 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.246639021599e+05 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.860166571872e+04 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.484329051490e+03 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.050411687443e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.488366972009e+02 lambda=3.7767265396089055e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.508326035050e+02 lambda=7.2725649346261757e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434696063559e+02 lambda=7.2725649346261764e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433964618996e+02 lambda=7.2725649346261768e-07
> > > > >       Line search: Cubically determined step, current gnorm 1.433958141405e+02 lambda=7.2725649346261771e-08
> > > > >  36 SNES Function norm 1.433958141405e+02
> > > > >     0 KSP Residual norm 2.164813477994e+03
> > > > >     1 KSP Residual norm 1.148881458292e-10
> > > > >       Line search: gnorm after quadratic fit 9.626940870744e+08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.210459267934e+08 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.531855101756e+07 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.966826097072e+06 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.611808255272e+05 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.746062783262e+04 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.252259871244e+03 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.319447372021e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.963055448218e+02 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.719034797105e+02 lambda=1.3935000445038475e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.436664111092e+02 lambda=1.3935000445038476e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433983336239e+02 lambda=1.3935000445038476e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958213497e+02 lambda=1.3935000445038476e-07
> > > > >       Line search: Cubically determined step, current gnorm 1.433958104705e+02 lambda=5.1202466521517630e-08
> > > > >  37 SNES Function norm 1.433958104705e+02
> > > > >     0 KSP Residual norm 5.470482746849e+03
> > > > >     1 KSP Residual norm 2.077456833170e-10
> > > > >       Line search: gnorm after quadratic fit 1.541335606193e+10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.922526157954e+09 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.393079394383e+08 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.967573549050e+07 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.657319925168e+06 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.479516204895e+05 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.573399122917e+04 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.931177424479e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.619710606187e+03 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.924551713717e+02 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.786095404459e+02 lambda=6.2808469554243522e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.437467476469e+02 lambda=6.2808469554243527e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433992463439e+02 lambda=6.2808469554243527e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958367271e+02 lambda=6.2808469554243525e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.433958098948e+02 lambda=8.0208105170108588e-09
> > > > >  38 SNES Function norm 1.433958098948e+02
> > > > >     0 KSP Residual norm 1.380568582849e+04
> > > > >     1 KSP Residual norm 3.830866223513e-10
> > > > >       Line search: gnorm after quadratic fit 2.484973518314e+11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.108953896984e+10 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.893104137528e+09 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.884003910853e+08 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.150765563542e+07 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.811161146103e+06 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.011017986781e+06 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.368084343451e+05 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.042876665700e+04 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.648735196752e+03 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.489051252413e+02 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.910677756717e+02 lambda=4.7728537787817724e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.505452645186e+02 lambda=1.1099980464343249e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434658984864e+02 lambda=1.1099980464343249e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433964956476e+02 lambda=1.1099980464343249e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958153212e+02 lambda=1.1099980464343250e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.433958098048e+02 lambda=1.2587012037197021e-09
> > > > >  39 SNES Function norm 1.433958098048e+02
> > > > >     0 KSP Residual norm 3.492284741552e+04
> > > > >     1 KSP Residual norm 2.459788921244e-09
> > > > >       Line search: gnorm after quadratic fit 4.017424324975e+12
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.020053801097e+11 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.270768402853e+10 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.827801904036e+09 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.758550488105e+08 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.213492627852e+08 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.502191365805e+07 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.846947104704e+06 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.262850216093e+05 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.879926646684e+04 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.510203332079e+03 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.055028679619e+03 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.503903269554e+02 lambda=2.3633949212111800e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.510245991472e+02 lambda=4.5897967908360529e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434724062782e+02 lambda=4.5897967908360533e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433965706606e+02 lambda=4.5897967908360533e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958168196e+02 lambda=4.5897967908360538e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958098155e+02 lambda=4.5897967908360540e-10
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097906e+02 lambda=1.9745369723997599e-10
> > > > >  40 SNES Function norm 1.433958097906e+02
> > > > >     0 KSP Residual norm 8.722495521587e+04
> > > > >     1 KSP Residual norm 3.742695919275e-09
> > > > >       Line search: gnorm after quadratic fit 6.262538457180e+13
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.829256308992e+12 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.789282877890e+11 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.224340679566e+11 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.532137613500e+10 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.919506047737e+09 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.410488718338e+08 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.042211373104e+07 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.881986305609e+06 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.080857001861e+05 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.054449734307e+04 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.109051581397e+04 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.144985497099e+03 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.617627947761e+02 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.132828960986e+02 lambda=5.3137869181552773e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.440403409760e+02 lambda=5.3137869181552777e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434022226216e+02 lambda=5.3137869181552781e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958732177e+02 lambda=5.3137869181552781e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958103569e+02 lambda=5.3137869181552779e-10
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097894e+02 lambda=5.3137869181552781e-11
> > > > >  41 SNES Function norm 1.433958097894e+02
> > > > >     0 KSP Residual norm 6.453169081212e+04
> > > > >     1 KSP Residual norm 2.762540688686e-09
> > > > >       Line search: gnorm after quadratic fit 2.535166112592e+13
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.168366990040e+12 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.958985371462e+11 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.945064622488e+10 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.172244865713e+09 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.693002384290e+08 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.562568994785e+07 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.182957296890e+07 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.453260254263e+06 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.781984147743e+05 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.294934468515e+04 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.740461720782e+03 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.162621855154e+02 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.042417294890e+02 lambda=1.1290474210136388e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.458920680477e+02 lambda=1.4201366604660443e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434208620254e+02 lambda=1.4201366604660444e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433960586269e+02 lambda=1.4201366604660445e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958120934e+02 lambda=1.4201366604660445e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958097940e+02 lambda=1.4201366604660446e-10
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097852e+02 lambda=5.7981795207229486e-11
> > > > >  42 SNES Function norm 1.433958097852e+02
> > > > >     0 KSP Residual norm 1.596625793747e+05
> > > > >     1 KSP Residual norm 6.465899717071e-09
> > > > >       Line search: gnorm after quadratic fit 3.840699863976e+14
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.801237514773e+13 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.002454410823e+12 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.505340831651e+11 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.387378188418e+10 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.174857842415e+10 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.472211181784e+09 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.849608933788e+08 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.336592011613e+07 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.988079805895e+06 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.930858080425e+05 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.521190722497e+04 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.872153015323e+03 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.772337662956e+03 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.879470852054e+02 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.942774855488e+02 lambda=2.4957912736226801e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.438719078958e+02 lambda=2.4957912736226800e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434005516391e+02 lambda=2.4957912736226800e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958568732e+02 lambda=2.4957912736226803e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958102244e+02 lambda=2.4957912736226804e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958097865e+02 lambda=2.4957912736226805e-11
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097846e+02 lambda=9.2900433293108317e-12
> > > > >  43 SNES Function norm 1.433958097846e+02
> > > > >     0 KSP Residual norm 4.230869747157e+05
> > > > >     1 KSP Residual norm 2.238707423027e-08
> > > > >       Line search: gnorm after quadratic fit 7.145700515048e+15
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.931871281700e+14 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.116420341041e+14 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.395366610168e+13 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.743811756566e+12 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.178776103292e+11 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.721012032848e+10 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.395186924428e+09 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.229125978585e+08 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.250971135953e+07 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.483856203094e+06 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.950671355228e+05 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.795408218555e+04 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.315025696280e+04 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.395961070555e+03 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.565070307576e+02 lambda=3.0517578125000002e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.228949713871e+02 lambda=1.2154989071946628e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.441831998014e+02 lambda=1.2154989071946629e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434037055996e+02 lambda=1.2154989071946630e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958886074e+02 lambda=1.2154989071946630e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958105564e+02 lambda=1.2154989071946630e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958097907e+02 lambda=1.2154989071946631e-11
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097845e+02 lambda=1.3559489351735629e-12
> > > > >  44 SNES Function norm 1.433958097845e+02
> > > > >     0 KSP Residual norm 1.017589039922e+06
> > > > >     1 KSP Residual norm 5.483283808994e-08
> > > > >       Line search: gnorm after quadratic fit 9.942386816509e+16
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.242813073279e+16 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.553553149772e+15 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.942033483320e+14 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.427772097758e+13 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.035291372732e+12 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.795558047824e+11 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.748073147307e+10 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.944235240368e+09 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.453550734860e+08 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.377045707707e+07 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.188119937132e+07 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.529730353136e+06 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.044646611329e+05 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.974477616118e+04 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.087923877078e+03 lambda=3.0517578125000002e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.112257173311e+03 lambda=1.5258789062500001e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.536190077033e+02 lambda=7.6293945312500004e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.622567424729e+02 lambda=2.4244966960965791e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.435780438142e+02 lambda=2.4244966960965793e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433976282603e+02 lambda=2.4244966960965796e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958279382e+02 lambda=2.4244966960965797e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958099632e+02 lambda=2.4244966960965799e-11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958097860e+02 lambda=2.4244966960965801e-12
> > > > >       Line search: Cubically determined step, current gnorm 1.433958097845e+02 lambda=2.4244966960965801e-13
> > > > >  45 SNES Function norm 1.433958097845e+02
> > > > >     0 KSP Residual norm 2.206687220910e+06
> > > > >     1 KSP Residual norm 4.897651438902e-08
> > > > >       Line search: gnorm after quadratic fit 1.013883843512e+18
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.267347883019e+17 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.584167551460e+16 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.980166189110e+15 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.475099638694e+14 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.093604443378e+13 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.866330988164e+12 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.831230804145e+11 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.034848618023e+10 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.533173457427e+09 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.390937545910e+08 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.167706366282e+08 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.445357932595e+07 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.776833600920e+06 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.177171496134e+05 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.775750596740e+04 lambda=3.0517578125000002e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.374283858049e+03 lambda=1.5258789062500001e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.030949543491e+03 lambda=7.6293945312500004e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.423064811256e+02 lambda=3.6673216108643130e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.499924205417e+02 lambda=6.7541706529544844e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.434620968378e+02 lambda=6.7541706529544851e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433964732224e+02 lambda=6.7541706529544853e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958164087e+02 lambda=6.7541706529544856e-11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958098496e+02 lambda=6.7541706529544860e-12
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.433958097850e+02 lambda=6.7541706529544869e-13
> > > > >       Line search: unable to find good step length! After 24 tries
> > > > >       Line search: fnorm=1.4339580978447020e+02, gnorm=1.4339580978501337e+02, ynorm=2.2066872446260620e+06, minlambda=9.9999999999999998e-13, lambda=6.7541706529544869e-13, initial slope=-2.0562359199040133e+04
> > > > >   0 SNES Function norm 7.494832241120e+03
> > > > >     0 KSP Residual norm 2.096735360816e+01
> > > > >     1 KSP Residual norm 1.310027756820e-12
> > > > >       Line search: gnorm after quadratic fit 6.728375857515e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   1 SNES Function norm 6.728375857515e+03
> > > > >     0 KSP Residual norm 2.051806116405e+01
> > > > >     1 KSP Residual norm 4.624620138831e-13
> > > > >       Line search: gnorm after quadratic fit 6.028616261650e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   2 SNES Function norm 6.028616261650e+03
> > > > >     0 KSP Residual norm 2.416503160016e+01
> > > > >     1 KSP Residual norm 8.456041680630e-13
> > > > >       Line search: gnorm after quadratic fit 5.407473517447e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   3 SNES Function norm 5.407473517447e+03
> > > > >     0 KSP Residual norm 1.429873750399e+01
> > > > >     1 KSP Residual norm 2.956316145819e-13
> > > > >       Line search: Using full step: fnorm 5.407473517447e+03 gnorm 1.894530516076e+03
> > > > >   4 SNES Function norm 1.894530516076e+03
> > > > >     0 KSP Residual norm 2.898580554597e+00
> > > > >     1 KSP Residual norm 6.622560162623e-14
> > > > >       Line search: Using full step: fnorm 1.894530516076e+03 gnorm 1.794029421846e+03
> > > > >   5 SNES Function norm 1.794029421846e+03
> > > > >     0 KSP Residual norm 1.749678611959e+01
> > > > >     1 KSP Residual norm 8.138905825078e-12
> > > > >       Line search: gnorm after quadratic fit 1.767361408940e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   6 SNES Function norm 1.767361408940e+03
> > > > >     0 KSP Residual norm 1.094404109423e+02
> > > > >     1 KSP Residual norm 7.696691527700e-12
> > > > >       Line search: gnorm after quadratic fit 3.545750923895e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.153479540314e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.205683629874e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.802976525485e+03 lambda=1.2441094586006883e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.765063299263e+03 lambda=4.9240328094708914e-03
> > > > >   7 SNES Function norm 1.765063299263e+03
> > > > >     0 KSP Residual norm 1.778095975189e+01
> > > > >     1 KSP Residual norm 2.814661813934e-12
> > > > >       Line search: gnorm after quadratic fit 2.620761680117e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.793045753271e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.740299227935e+03 lambda=2.5000000000000001e-02
> > > > >   8 SNES Function norm 1.740299227935e+03
> > > > >     0 KSP Residual norm 3.284236894535e+02
> > > > >     1 KSP Residual norm 5.172103783952e-10
> > > > >       Line search: gnorm after quadratic fit 2.857820523902e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.063993491139e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.588139591650e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.491163684413e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.802614760566e+03 lambda=5.7977212395253904e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.741444490992e+03 lambda=1.8448315876950813e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.739663656043e+03 lambda=7.7468345598227730e-04
> > > > >   9 SNES Function norm 1.739663656043e+03
> > > > >     0 KSP Residual norm 4.581839103558e+02
> > > > >     1 KSP Residual norm 2.627035885943e-11
> > > > >       Line search: gnorm after quadratic fit 8.199478499066e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.127270076812e+05 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.816774161281e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.053723877333e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.971814513392e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.755066208174e+03 lambda=2.7232010924558834e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.739559834479e+03 lambda=8.1458417855297680e-04
> > > > >  10 SNES Function norm 1.739559834479e+03
> > > > >     0 KSP Residual norm 1.463056810139e+02
> > > > >     1 KSP Residual norm 5.383584598825e-12
> > > > >       Line search: gnorm after quadratic fit 2.885699620595e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.847717401708e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.244464170034e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.769157604338e+03 lambda=1.0857209099577540e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.737356225452e+03 lambda=4.0312937622950231e-03
> > > > >  11 SNES Function norm 1.737356225452e+03
> > > > >     0 KSP Residual norm 1.564105677925e+02
> > > > >     1 KSP Residual norm 6.671164745832e-11
> > > > >       Line search: gnorm after quadratic fit 3.815800291873e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.197027796134e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.373151605545e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.783976520867e+03 lambda=1.2093374970461970e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.735737929915e+03 lambda=4.9529698477569920e-03
> > > > >  12 SNES Function norm 1.735737929915e+03
> > > > >     0 KSP Residual norm 5.030757139645e+01
> > > > >     1 KSP Residual norm 1.157542730692e-12
> > > > >       Line search: gnorm after quadratic fit 3.004544441458e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.850403239750e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.722893188305e+03 lambda=2.0881992439921702e-02
> > > > >  13 SNES Function norm 1.722893188305e+03
> > > > >     0 KSP Residual norm 7.123428853845e+01
> > > > >     1 KSP Residual norm 8.671726774894e-12
> > > > >       Line search: gnorm after quadratic fit 4.361041499279e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.032697222107e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.726212330814e+03 lambda=1.9909470348267504e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.714127356920e+03 lambda=9.9547351741337518e-03
> > > > >  14 SNES Function norm 1.714127356920e+03
> > > > >     0 KSP Residual norm 8.304557447257e+00
> > > > >     1 KSP Residual norm 6.268295862688e-13
> > > > >       Line search: gnorm after quadratic fit 1.558163002487e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  15 SNES Function norm 1.558163002487e+03
> > > > >     0 KSP Residual norm 2.820803287164e+00
> > > > >     1 KSP Residual norm 5.477413853752e-13
> > > > >       Line search: gnorm after quadratic fit 1.065673478530e+03
> > > > >       Line search: Quadratically determined step, lambda=3.8943438938591052e-01
> > > > >  16 SNES Function norm 1.065673478530e+03
> > > > >     0 KSP Residual norm 2.296739120664e+00
> > > > >     1 KSP Residual norm 6.273731899885e-14
> > > > >       Line search: gnorm after quadratic fit 8.964342150621e+02
> > > > >       Line search: Quadratically determined step, lambda=1.8740736900935545e-01
> > > > >  17 SNES Function norm 8.964342150621e+02
> > > > >     0 KSP Residual norm 5.567568505830e+00
> > > > >     1 KSP Residual norm 8.649083600764e-12
> > > > >       Line search: gnorm after quadratic fit 8.322514858992e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  18 SNES Function norm 8.322514858992e+02
> > > > >     0 KSP Residual norm 1.164393577210e+00
> > > > >     1 KSP Residual norm 2.019248309015e-14
> > > > >       Line search: Using full step: fnorm 8.322514858992e+02 gnorm 3.179750274746e+02
> > > > >  19 SNES Function norm 3.179750274746e+02
> > > > >     0 KSP Residual norm 5.339294310641e+00
> > > > >     1 KSP Residual norm 2.070321587238e-13
> > > > >       Line search: gnorm after quadratic fit 3.433012316833e+02
> > > > >       Line search: Cubically determined step, current gnorm 3.140305403821e+02 lambda=5.0000000000000003e-02
> > > > >  20 SNES Function norm 3.140305403821e+02
> > > > >     0 KSP Residual norm 1.177016565578e+01
> > > > >     1 KSP Residual norm 2.413027540446e-13
> > > > >       Line search: gnorm after quadratic fit 7.377851778409e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.896721016582e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 3.136619309181e+02 lambda=9.7219892278716195e-03
> > > > >  21 SNES Function norm 3.136619309181e+02
> > > > >     0 KSP Residual norm 3.973565083977e+00
> > > > >     1 KSP Residual norm 9.726786674410e-14
> > > > >       Line search: gnorm after quadratic fit 3.098277326090e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  22 SNES Function norm 3.098277326090e+02
> > > > >     0 KSP Residual norm 9.389290683519e+00
> > > > >     1 KSP Residual norm 1.651497163724e-13
> > > > >       Line search: gnorm after quadratic fit 6.887970540455e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.588146381477e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.099274823603e+02 lambda=1.6890426151283184e-02
> > > > >       Line search: Cubically determined step, current gnorm 3.084890760827e+02 lambda=8.4452130756415920e-03
> > > > >  23 SNES Function norm 3.084890760827e+02
> > > > >     0 KSP Residual norm 2.381130479641e+00
> > > > >     1 KSP Residual norm 4.694489283156e-14
> > > > >       Line search: gnorm after quadratic fit 2.872151876535e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  24 SNES Function norm 2.872151876535e+02
> > > > >     0 KSP Residual norm 2.405498502269e+00
> > > > >     1 KSP Residual norm 3.316846070538e-13
> > > > >       Line search: gnorm after quadratic fit 2.686789761232e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  25 SNES Function norm 2.686789761232e+02
> > > > >     0 KSP Residual norm 1.728772970554e+00
> > > > >     1 KSP Residual norm 4.132974383339e-14
> > > > >       Line search: gnorm after quadratic fit 2.365009918229e+02
> > > > >       Line search: Quadratically determined step, lambda=1.7273357529379074e-01
> > > > >  26 SNES Function norm 2.365009918229e+02
> > > > >     0 KSP Residual norm 4.413764759374e+00
> > > > >     1 KSP Residual norm 5.210690816917e-14
> > > > >       Line search: gnorm after quadratic fit 2.756730968257e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.372321757171e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.335020811250e+02 lambda=2.5000000000000001e-02
> > > > >  27 SNES Function norm 2.335020811250e+02
> > > > >     0 KSP Residual norm 1.606246507553e+00
> > > > >     1 KSP Residual norm 3.564847846678e-14
> > > > >       Line search: gnorm after quadratic fit 2.078263630653e+02
> > > > >       Line search: Quadratically determined step, lambda=1.5913860995949433e-01
> > > > >  28 SNES Function norm 2.078263630653e+02
> > > > >     0 KSP Residual norm 3.700632954873e+00
> > > > >     1 KSP Residual norm 5.113863400264e-13
> > > > >       Line search: gnorm after quadratic fit 2.142503514807e+02
> > > > >       Line search: Cubically determined step, current gnorm 2.021095009118e+02 lambda=5.0000000000000003e-02
> > > > >  29 SNES Function norm 2.021095009118e+02
> > > > >     0 KSP Residual norm 3.056449560135e+00
> > > > >     1 KSP Residual norm 3.207987681334e-14
> > > > >       Line search: gnorm after quadratic fit 2.064252530802e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.978069081639e+02 lambda=5.0000000000000003e-02
> > > > >  30 SNES Function norm 1.978069081639e+02
> > > > >     0 KSP Residual norm 2.024695620703e+00
> > > > >     1 KSP Residual norm 5.460360737995e-14
> > > > >       Line search: gnorm after quadratic fit 1.839556530796e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  31 SNES Function norm 1.839556530796e+02
> > > > >     0 KSP Residual norm 1.610676619931e+01
> > > > >     1 KSP Residual norm 6.703552136307e-13
> > > > >       Line search: gnorm after quadratic fit 7.186735314913e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.905672430097e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.856766286600e+02 lambda=1.0830798014298563e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.836818937131e+02 lambda=3.3207362501459785e-03
> > > > >  32 SNES Function norm 1.836818937131e+02
> > > > >     0 KSP Residual norm 7.471722173131e+01
> > > > >     1 KSP Residual norm 2.671534648829e-12
> > > > >       Line search: gnorm after quadratic fit 9.613379909411e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.509491298762e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.808861367705e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.767283758299e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.665421328348e+02 lambda=6.0369453941238101e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.869726717041e+02 lambda=1.4394327006264963e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.836897704575e+02 lambda=1.4394327006264963e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.836767951408e+02 lambda=5.5567420418543164e-05
> > > > >  33 SNES Function norm 1.836767951408e+02
> > > > >     0 KSP Residual norm 2.702367712138e+01
> > > > >     1 KSP Residual norm 2.731656687850e-12
> > > > >       Line search: gnorm after quadratic fit 3.331563146098e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.723694790688e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.941243220944e+02 lambda=2.1443568774664485e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.847231733910e+02 lambda=2.7571330376860012e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.836356729409e+02 lambda=4.7841280418270480e-04
> > > > >  34 SNES Function norm 1.836356729409e+02
> > > > >     0 KSP Residual norm 1.228333177997e+01
> > > > >     1 KSP Residual norm 2.681127387880e-13
> > > > >       Line search: gnorm after quadratic fit 6.936329223475e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.749327218775e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.871557327742e+02 lambda=1.4160505301999991e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.833523080682e+02 lambda=3.5196326548579192e-03
> > > > >  35 SNES Function norm 1.833523080682e+02
> > > > >     0 KSP Residual norm 3.531886257396e+02
> > > > >     1 KSP Residual norm 2.364634092895e-11
> > > > >       Line search: gnorm after quadratic fit 3.994783047929e+06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.753715960726e+05 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.516669898185e+04 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.918985942348e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.363599250418e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.344906818752e+02 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.991088183980e+02 lambda=1.0108094710462592e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.834621681844e+02 lambda=1.0108094710462593e-04
> > > > >       Line search: Cubically determined step, current gnorm 1.833517377324e+02 lambda=1.0108094710462594e-05
> > > > >  36 SNES Function norm 1.833517377324e+02
> > > > >     0 KSP Residual norm 9.005833507118e+01
> > > > >     1 KSP Residual norm 6.116387880629e-12
> > > > >       Line search: gnorm after quadratic fit 8.899847103226e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.388111536526e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.532652074749e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.864912734009e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.421068789960e+02 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.873498120960e+02 lambda=2.1924025345283278e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.833506987706e+02 lambda=2.1924025345283280e-04
> > > > >  37 SNES Function norm 1.833506987706e+02
> > > > >     0 KSP Residual norm 1.548426525207e+01
> > > > >     1 KSP Residual norm 1.038038773942e-12
> > > > >       Line search: gnorm after quadratic fit 6.702848665441e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.789880616078e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.846521504351e+02 lambda=1.0567302644323170e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.830548481815e+02 lambda=3.5274490352771972e-03
> > > > >  38 SNES Function norm 1.830548481815e+02
> > > > >     0 KSP Residual norm 1.523095478807e+01
> > > > >     1 KSP Residual norm 1.963823596119e-12
> > > > >       Line search: gnorm after quadratic fit 9.255727478491e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.496757253461e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.863126241242e+02 lambda=9.3143291632966311e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.829091761403e+02 lambda=1.8107234819462800e-03
> > > > >  39 SNES Function norm 1.829091761403e+02
> > > > >     0 KSP Residual norm 1.311320726934e+01
> > > > >     1 KSP Residual norm 9.084170520902e-13
> > > > >       Line search: gnorm after quadratic fit 7.488610069188e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.691148243365e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.885558228772e+02 lambda=1.9439894235886539e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.825540619134e+02 lambda=5.4967824488704169e-03
> > > > >  40 SNES Function norm 1.825540619134e+02
> > > > >     0 KSP Residual norm 1.218635557505e+01
> > > > >     1 KSP Residual norm 7.760485728617e-13
> > > > >       Line search: gnorm after quadratic fit 6.636453826807e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.880828006022e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.830493790584e+02 lambda=6.6234802648049863e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.823397545268e+02 lambda=2.4308217464828865e-03
> > > > >  41 SNES Function norm 1.823397545268e+02
> > > > >     0 KSP Residual norm 9.456543593865e+00
> > > > >     1 KSP Residual norm 7.046593270309e-13
> > > > >       Line search: gnorm after quadratic fit 3.369769163110e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.105230957789e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.817897533119e+02 lambda=9.1612481372917148e-03
> > > > >  42 SNES Function norm 1.817897533119e+02
> > > > >     0 KSP Residual norm 7.290035145805e+00
> > > > >     1 KSP Residual norm 3.198962158141e-12
> > > > >       Line search: gnorm after quadratic fit 2.858311567415e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.965004842981e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.808845577764e+02 lambda=1.3826421990147811e-02
> > > > >  43 SNES Function norm 1.808845577764e+02
> > > > >     0 KSP Residual norm 7.256785036157e+00
> > > > >     1 KSP Residual norm 9.243179217625e-14
> > > > >       Line search: gnorm after quadratic fit 2.534388176390e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.916726818893e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.797953125543e+02 lambda=1.3677747819164022e-02
> > > > >  44 SNES Function norm 1.797953125543e+02
> > > > >     0 KSP Residual norm 1.716201096352e+01
> > > > >     1 KSP Residual norm 2.701418800808e-13
> > > > >       Line search: gnorm after quadratic fit 1.331064301474e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.461842246267e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.920668830294e+02 lambda=1.2856504859057132e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.797121460957e+02 lambda=1.4396611381500967e-03
> > > > >  45 SNES Function norm 1.797121460957e+02
> > > > >     0 KSP Residual norm 6.613432967985e+00
> > > > >     1 KSP Residual norm 1.357752977076e-13
> > > > >       Line search: gnorm after quadratic fit 2.721288927858e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.939638282615e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.787985482018e+02 lambda=1.2647907914070569e-02
> > > > >  46 SNES Function norm 1.787985482018e+02
> > > > >     0 KSP Residual norm 1.225736349281e+01
> > > > >     1 KSP Residual norm 3.819378790812e-13
> > > > >       Line search: gnorm after quadratic fit 4.548006065036e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.304803559060e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.787344729174e+02 lambda=8.9002284237256531e-03
> > > > >  47 SNES Function norm 1.787344729174e+02
> > > > >     0 KSP Residual norm 6.302465402340e+00
> > > > >     1 KSP Residual norm 6.980931947122e-14
> > > > >       Line search: gnorm after quadratic fit 2.879477700004e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.980806946742e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.780128006935e+02 lambda=1.0223293444602008e-02
> > > > >  48 SNES Function norm 1.780128006935e+02
> > > > >     0 KSP Residual norm 4.323829277968e+00
> > > > >     1 KSP Residual norm 5.223881369308e-13
> > > > >       Line search: gnorm after quadratic fit 1.957928151218e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.768899805640e+02 lambda=5.0000000000000003e-02
> > > > >  49 SNES Function norm 1.768899805640e+02
> > > > >     0 KSP Residual norm 2.249256107033e+00
> > > > >     1 KSP Residual norm 3.624400957270e-14
> > > > >       Line search: gnorm after quadratic fit 1.704782114773e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  50 SNES Function norm 1.704782114773e+02
> > > > >   0 SNES Function norm 3.513783397332e+03
> > > > >     0 KSP Residual norm 1.650214950648e+01
> > > > >     1 KSP Residual norm 3.574269752690e-13
> > > > >       Line search: gnorm after quadratic fit 3.155830404050e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   1 SNES Function norm 3.155830404050e+03
> > > > >     0 KSP Residual norm 1.443734018042e+01
> > > > >     1 KSP Residual norm 3.685565191058e-13
> > > > >       Line search: Using full step: fnorm 3.155830404050e+03 gnorm 8.581212204155e+02
> > > > >   2 SNES Function norm 8.581212204155e+02
> > > > >     0 KSP Residual norm 2.492601775565e+00
> > > > >     1 KSP Residual norm 9.698440210389e-14
> > > > >       Line search: Using full step: fnorm 8.581212204155e+02 gnorm 7.018609898542e+02
> > > > >   3 SNES Function norm 7.018609898542e+02
> > > > >     0 KSP Residual norm 1.740320986421e+00
> > > > >     1 KSP Residual norm 2.868780682435e-14
> > > > >       Line search: Using full step: fnorm 7.018609898542e+02 gnorm 2.135824235887e+02
> > > > >   4 SNES Function norm 2.135824235887e+02
> > > > >     0 KSP Residual norm 1.647413497077e+00
> > > > >     1 KSP Residual norm 2.731226225666e-14
> > > > >       Line search: gnorm after quadratic fit 1.936469032649e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   5 SNES Function norm 1.936469032649e+02
> > > > >     0 KSP Residual norm 1.406615972124e+00
> > > > >     1 KSP Residual norm 3.167179626699e-14
> > > > >       Line search: gnorm after quadratic fit 1.755362571467e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   6 SNES Function norm 1.755362571467e+02
> > > > >     0 KSP Residual norm 1.480681706594e+00
> > > > >     1 KSP Residual norm 2.935210968935e-14
> > > > >       Line search: gnorm after quadratic fit 1.597659506616e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   7 SNES Function norm 1.597659506616e+02
> > > > >     0 KSP Residual norm 4.013154698097e+00
> > > > >     1 KSP Residual norm 1.021628704832e-13
> > > > >       Line search: gnorm after quadratic fit 1.728408060966e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.570815760721e+02 lambda=5.0000000000000003e-02
> > > > >   8 SNES Function norm 1.570815760721e+02
> > > > >     0 KSP Residual norm 1.510335662636e+00
> > > > >     1 KSP Residual norm 2.728243244724e-14
> > > > >       Line search: gnorm after quadratic fit 1.450162880692e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   9 SNES Function norm 1.450162880692e+02
> > > > >     0 KSP Residual norm 1.923349271077e+01
> > > > >     1 KSP Residual norm 1.640711956406e-11
> > > > >       Line search: gnorm after quadratic fit 1.016537223115e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.584263092048e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.490729346583e+02 lambda=9.3725990358703350e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.449349273006e+02 lambda=1.5464889706033082e-03
> > > > >  10 SNES Function norm 1.449349273006e+02
> > > > >     0 KSP Residual norm 1.154999084194e+01
> > > > >     1 KSP Residual norm 1.292167633338e-11
> > > > >       Line search: gnorm after quadratic fit 6.060386918705e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.236630526035e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.491201927819e+02 lambda=1.6816056300124185e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.447023047013e+02 lambda=4.1484374564153808e-03
> > > > >  11 SNES Function norm 1.447023047013e+02
> > > > >     0 KSP Residual norm 7.278906180502e+00
> > > > >     1 KSP Residual norm 2.815632651924e-12
> > > > >       Line search: gnorm after quadratic fit 2.512278711773e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.624350957814e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.441656049510e+02 lambda=1.0879389002513012e-02
> > > > >  12 SNES Function norm 1.441656049510e+02
> > > > >     0 KSP Residual norm 4.449836480267e+00
> > > > >     1 KSP Residual norm 3.152365624813e-13
> > > > >       Line search: gnorm after quadratic fit 1.749680739878e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.458763435996e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.425715025152e+02 lambda=2.2766809271842225e-02
> > > > >  13 SNES Function norm 1.425715025152e+02
> > > > >     0 KSP Residual norm 4.218495037726e+00
> > > > >     1 KSP Residual norm 1.398472536142e-12
> > > > >       Line search: gnorm after quadratic fit 1.645159536467e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.421991559212e+02 lambda=4.1883769431247941e-02
> > > > >  14 SNES Function norm 1.421991559212e+02
> > > > >     0 KSP Residual norm 1.968853538608e+00
> > > > >     1 KSP Residual norm 7.594023450519e-12
> > > > >       Line search: gnorm after quadratic fit 1.350960142124e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  15 SNES Function norm 1.350960142124e+02
> > > > >     0 KSP Residual norm 3.444721686085e+00
> > > > >     1 KSP Residual norm 5.312609674019e-14
> > > > >       Line search: gnorm after quadratic fit 1.472880565107e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.336714620145e+02 lambda=4.1341550820829437e-02
> > > > >  16 SNES Function norm 1.336714620145e+02
> > > > >     0 KSP Residual norm 3.276288899397e+00
> > > > >     1 KSP Residual norm 8.394674946715e-14
> > > > >       Line search: gnorm after quadratic fit 1.459274497150e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.327618539486e+02 lambda=5.0000000000000003e-02
> > > > >  17 SNES Function norm 1.327618539486e+02
> > > > >     0 KSP Residual norm 2.630052244303e+00
> > > > >     1 KSP Residual norm 4.351283410507e-14
> > > > >       Line search: gnorm after quadratic fit 1.350378060116e+02
> > > > >       Line search: Cubically determined step, current gnorm 1.299403903934e+02 lambda=4.9913529436269283e-02
> > > > >  18 SNES Function norm 1.299403903934e+02
> > > > >     0 KSP Residual norm 6.124953430138e+00
> > > > >     1 KSP Residual norm 2.295381352938e-13
> > > > >       Line search: gnorm after quadratic fit 2.275621676963e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.469611730657e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.294890694898e+02 lambda=1.0071939100661648e-02
> > > > >  19 SNES Function norm 1.294890694898e+02
> > > > >     0 KSP Residual norm 1.036733907011e+01
> > > > >     1 KSP Residual norm 1.800941381283e-13
> > > > >       Line search: gnorm after quadratic fit 3.844879463595e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.875071148504e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 1.294494430642e+02 lambda=5.0000000000000010e-03
> > > > >  20 SNES Function norm 1.294494430642e+02
> > > > >     0 KSP Residual norm 8.224001595443e+00
> > > > >     1 KSP Residual norm 3.337810156182e-13
> > > > >       Line search: gnorm after quadratic fit 3.332325263497e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.682441841234e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.294720538860e+02 lambda=8.5401597581228530e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.291762382985e+02 lambda=4.2555418164960989e-03
> > > > >  21 SNES Function norm 1.291762382985e+02
> > > > >     0 KSP Residual norm 3.920749219860e+01
> > > > >     1 KSP Residual norm 1.610902886435e-12
> > > > >       Line search: gnorm after quadratic fit 3.472422440640e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.023065712859e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.228184088700e+02 lambda=1.6004598823093592e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.298607525058e+02 lambda=1.6004598823093593e-03
> > > > >       Line search: Cubically determined step, current gnorm 1.291643460998e+02 lambda=1.9319076533745672e-04
> > > > >  22 SNES Function norm 1.291643460998e+02
> > > > >     0 KSP Residual norm 1.520092314848e+02
> > > > >     1 KSP Residual norm 7.089159192434e-11
> > > > >       Line search: gnorm after quadratic fit 2.752539686422e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.237188995536e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.484370498858e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.571039994484e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.280898858362e+02 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.708927009106e+02 lambda=2.6114913411793973e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.294919567784e+02 lambda=2.6114913411793975e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291645609810e+02 lambda=2.6114913411793977e-05
> > > > >       Line search: Cubically determined step, current gnorm 1.291635530596e+02 lambda=1.2278773895355162e-05
> > > > >  23 SNES Function norm 1.291635530596e+02
> > > > >     0 KSP Residual norm 7.569206058965e+02
> > > > >     1 KSP Residual norm 3.454693729751e-11
> > > > >       Line search: gnorm after quadratic fit 2.570888738395e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.064965025187e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.498514098182e+05 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.807487005202e+04 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.233167830543e+03 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.396736912757e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.547572543330e+02 lambda=1.2623406091699435e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.335889024473e+02 lambda=1.8542137907683829e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.292059042479e+02 lambda=1.8542137907683831e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291637618568e+02 lambda=1.8542137907683832e-06
> > > > >       Line search: Cubically determined step, current gnorm 1.291635210734e+02 lambda=4.9524172913414387e-07
> > > > >  24 SNES Function norm 1.291635210734e+02
> > > > >     0 KSP Residual norm 3.761913422861e+03
> > > > >     1 KSP Residual norm 6.181718776929e-10
> > > > >       Line search: gnorm after quadratic fit 3.342370445037e+09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.218326646111e+08 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.375128803204e+07 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.980626157021e+06 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.407418671219e+05 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.357321025399e+05 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.188948059047e+04 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.105117929713e+03 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.331054736818e+02 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.942479792843e+02 lambda=1.9412500272460620e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.436944624694e+02 lambda=6.4483223307578726e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.292972287211e+02 lambda=6.4483223307578726e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291647778160e+02 lambda=6.4483223307578730e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635261433e+02 lambda=6.4483223307578730e-08
> > > > >       Line search: Cubically determined step, current gnorm 1.291635197798e+02 lambda=2.0042115918485846e-08
> > > > >  25 SNES Function norm 1.291635197798e+02
> > > > >     0 KSP Residual norm 1.874745766083e+04
> > > > >     1 KSP Residual norm 1.476978150334e-09
> > > > >       Line search: gnorm after quadratic fit 4.089262111368e+11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.101717203770e+10 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.352565940292e+09 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.879613196620e+08 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.698613558125e+07 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.175566265039e+07 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.382919964952e+06 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.545431975334e+05 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.713561369103e+04 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.031789308419e+03 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.205847020738e+02 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.957203153158e+02 lambda=2.8132879163728503e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.297923302791e+02 lambda=2.8132879163728504e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291698100579e+02 lambda=2.8132879163728504e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635794525e+02 lambda=2.8132879163728506e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635200489e+02 lambda=2.8132879163728508e-09
> > > > >       Line search: Cubically determined step, current gnorm 1.291635197275e+02 lambda=8.0832061492461345e-10
> > > > >  26 SNES Function norm 1.291635197275e+02
> > > > >     0 KSP Residual norm 9.248381077528e+04
> > > > >     1 KSP Residual norm 7.262307068447e-09
> > > > >       Line search: gnorm after quadratic fit 4.920647210624e+13
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.153216818065e+12 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.697543960375e+11 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.637004379805e+10 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.208402653149e+10 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.519988135798e+09 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.923903214787e+08 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.465663001976e+07 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.238603967195e+06 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.459179143177e+05 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.676301042792e+04 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.135808875752e+04 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.275714918657e+03 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.719037747616e+02 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.039826156716e+02 lambda=5.5609199181402721e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.308092967809e+02 lambda=9.1057337296683134e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291796742824e+02 lambda=9.1057337296683134e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291636800174e+02 lambda=9.1057337296683134e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635212255e+02 lambda=9.1057337296683134e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635197320e+02 lambda=9.1057337296683139e-11
> > > > >       Line search: Cubically determined step, current gnorm 1.291635197254e+02 lambda=3.2898162617097329e-11
> > > > >  27 SNES Function norm 1.291635197254e+02
> > > > >     0 KSP Residual norm 4.818745996826e+05
> > > > >     1 KSP Residual norm 3.300979345194e-08
> > > > >       Line search: gnorm after quadratic fit 6.957046028867e+15
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.695654311477e+14 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.086793500688e+14 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.358083744813e+13 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.696584801696e+12 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.118183549773e+11 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.641372080976e+10 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.285878535769e+09 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.068045470276e+08 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.988290932539e+07 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.001404304764e+06 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.962234585736e+05 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.648190078115e+04 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.098288624714e+03 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.025132922751e+03 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.536770142392e+02 lambda=3.0517578125000002e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.529553964021e+02 lambda=6.6615844706846272e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.293970371303e+02 lambda=6.6615844706846277e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291658631829e+02 lambda=6.6615844706846284e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635430890e+02 lambda=6.6615844706846284e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635199508e+02 lambda=6.6615844706846284e-11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635197268e+02 lambda=6.6615844706846284e-12
> > > > >       Line search: Cubically determined step, current gnorm 1.291635197253e+02 lambda=1.2496728806533799e-12
> > > > >  28 SNES Function norm 1.291635197253e+02
> > > > >     0 KSP Residual norm 2.124622394867e+06
> > > > >     1 KSP Residual norm 2.766225333896e-07
> > > > >       Line search: gnorm after quadratic fit 5.963587884154e+17
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.454611858824e+16 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.318582340500e+15 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.164902175749e+15 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.456326197371e+14 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.820904039469e+13 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.277371273495e+12 lambda=1.5625000000000001e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.849819609184e+11 lambda=7.8125000000000004e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.570050545145e+10 lambda=3.9062500000000002e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.482064028328e+09 lambda=1.9531250000000001e-04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.651631635063e+08 lambda=9.7656250000000005e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.188623828904e+07 lambda=4.8828125000000003e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.302868645305e+06 lambda=2.4414062500000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.245214424313e+06 lambda=1.2207031250000001e-05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.775004618570e+05 lambda=6.1035156250000003e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.810228834086e+04 lambda=3.0517578125000002e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.148910945566e+03 lambda=1.5258789062500001e-06
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.133679879416e+03 lambda=7.6293945312500004e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.382468615011e+02 lambda=3.8146972656250002e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.519773483633e+02 lambda=1.4103864371136453e-07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.293690599792e+02 lambda=1.4103864371136454e-08
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291655645282e+02 lambda=1.4103864371136455e-09
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635401518e+02 lambda=1.4103864371136456e-10
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635199283e+02 lambda=1.4103864371136456e-11
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635197272e+02 lambda=1.4103864371136456e-12
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.291635197253e+02 lambda=1.4103864371136457e-13
> > > > >       Line search: unable to find good step length! After 25 tries
> > > > >       Line search: fnorm=1.2916351972528861e+02, gnorm=1.2916351972529532e+02, ynorm=2.1246218291095798e+06, minlambda=9.9999999999999998e-13, lambda=1.4103864371136457e-13, initial slope=-1.6683210999382733e+04
> > > > >   0 SNES Function norm 7.158176401507e+03
> > > > >     0 KSP Residual norm 7.545626598553e+02
> > > > >     1 KSP Residual norm 2.192822940623e-11
> > > > >       Line search: gnorm after quadratic fit 6.003975664723e+07
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.501930637484e+06 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 9.380142516806e+05 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.179837352860e+05 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.656902039419e+04 lambda=6.2500000000000003e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.340091686120e+03 lambda=3.1250000000000002e-03
> > > > >       Line search: Cubically determined step, current gnorm 7.127478647243e+03 lambda=1.5625000000000001e-03
> > > > >   1 SNES Function norm 7.127478647243e+03
> > > > >     0 KSP Residual norm 3.139792512249e+01
> > > > >     1 KSP Residual norm 5.765552480238e-13
> > > > >       Line search: gnorm after quadratic fit 7.131728282938e+03
> > > > >       Line search: Cubically determined step, current gnorm 6.730387269330e+03 lambda=5.0000000000000003e-02
> > > > >   2 SNES Function norm 6.730387269330e+03
> > > > >     0 KSP Residual norm 2.247436817553e+01
> > > > >     1 KSP Residual norm 4.129717651221e-13
> > > > >       Line search: gnorm after quadratic fit 6.018304311910e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   3 SNES Function norm 6.018304311910e+03
> > > > >     0 KSP Residual norm 1.605973421081e+01
> > > > >     1 KSP Residual norm 3.614891198134e-13
> > > > >       Line search: gnorm after quadratic fit 5.394599276028e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   4 SNES Function norm 5.394599276028e+03
> > > > >     0 KSP Residual norm 1.491002227850e+01
> > > > >     1 KSP Residual norm 5.905756964447e-13
> > > > >       Line search: gnorm after quadratic fit 4.845108544722e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   5 SNES Function norm 4.845108544722e+03
> > > > >     0 KSP Residual norm 1.562997766581e+02
> > > > >     1 KSP Residual norm 5.722461855805e-12
> > > > >       Line search: gnorm after quadratic fit 1.663505509947e+05
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.368547472010e+04 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.067825049920e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.852173464442e+03 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.819549937425e+03 lambda=6.2500000000000003e-03
> > > > >   6 SNES Function norm 4.819549937425e+03
> > > > >     0 KSP Residual norm 1.652787385042e+01
> > > > >     1 KSP Residual norm 7.774483442550e-13
> > > > >       Line search: gnorm after quadratic fit 2.868840270198e+03
> > > > >       Line search: Quadratically determined step, lambda=3.9586658383856743e-01
> > > > >   7 SNES Function norm 2.868840270198e+03
> > > > >     0 KSP Residual norm 1.220061851091e+01
> > > > >     1 KSP Residual norm 4.785407437718e-13
> > > > >       Line search: gnorm after quadratic fit 2.616720732407e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >   8 SNES Function norm 2.616720732407e+03
> > > > >     0 KSP Residual norm 2.884722153958e+01
> > > > >     1 KSP Residual norm 5.474553921306e-13
> > > > >       Line search: gnorm after quadratic fit 3.025386805317e+03
> > > > >       Line search: Cubically determined step, current gnorm 2.572110173067e+03 lambda=5.0000000000000003e-02
> > > > >   9 SNES Function norm 2.572110173067e+03
> > > > >     0 KSP Residual norm 5.239447686736e+01
> > > > >     1 KSP Residual norm 1.006906008399e-12
> > > > >       Line search: gnorm after quadratic fit 5.237562098046e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.722529781980e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.553072891461e+03 lambda=2.5000000000000001e-02
> > > > >  10 SNES Function norm 2.553072891461e+03
> > > > >     0 KSP Residual norm 6.511316788880e+00
> > > > >     1 KSP Residual norm 1.340296659008e-13
> > > > >       Line search: gnorm after quadratic fit 2.299704119080e+03
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  11 SNES Function norm 2.299704119080e+03
> > > > >     0 KSP Residual norm 5.268131426587e+00
> > > > >     1 KSP Residual norm 1.017563310127e-13
> > > > >       Line search: Using full step: fnorm 2.299704119080e+03 gnorm 7.642690039388e+02
> > > > >  12 SNES Function norm 7.642690039388e+02
> > > > >     0 KSP Residual norm 1.467735721574e+01
> > > > >     1 KSP Residual norm 1.090242963543e-12
> > > > >       Line search: gnorm after quadratic fit 1.042766084830e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.924803860137e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 7.581126971182e+02 lambda=1.8508848315139208e-02
> > > > >  13 SNES Function norm 7.581126971182e+02
> > > > >     0 KSP Residual norm 2.222609581896e+00
> > > > >     1 KSP Residual norm 5.661437314341e-14
> > > > >       Line search: gnorm after quadratic fit 6.235337602260e+02
> > > > >       Line search: Quadratically determined step, lambda=2.1172883827013136e-01
> > > > >  14 SNES Function norm 6.235337602260e+02
> > > > >     0 KSP Residual norm 3.080209609798e+00
> > > > >     1 KSP Residual norm 6.073476899313e-14
> > > > >       Line search: gnorm after quadratic fit 5.704745248520e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0142301129466913e-01
> > > > >  15 SNES Function norm 5.704745248520e+02
> > > > >     0 KSP Residual norm 5.628316803979e+00
> > > > >     1 KSP Residual norm 9.930477041779e-14
> > > > >       Line search: gnorm after quadratic fit 5.401354794776e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  16 SNES Function norm 5.401354794776e+02
> > > > >     0 KSP Residual norm 2.052541406719e+01
> > > > >     1 KSP Residual norm 4.328836834239e-13
> > > > >       Line search: gnorm after quadratic fit 1.709850100987e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.717966555703e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.414032576976e+02 lambda=8.7805694170804971e-03
> > > > >       Line search: Cubically determined step, current gnorm 5.391967144330e+02 lambda=3.6341472475199424e-03
> > > > >  17 SNES Function norm 5.391967144330e+02
> > > > >     0 KSP Residual norm 2.246993769488e+00
> > > > >     1 KSP Residual norm 5.078373642913e-14
> > > > >       Line search: gnorm after quadratic fit 4.939618639951e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  18 SNES Function norm 4.939618639951e+02
> > > > >     0 KSP Residual norm 2.155867648564e+00
> > > > >     1 KSP Residual norm 4.438622106380e-14
> > > > >       Line search: gnorm after quadratic fit 4.334377322188e+02
> > > > >       Line search: Quadratically determined step, lambda=1.5092486954829210e-01
> > > > >  19 SNES Function norm 4.334377322188e+02
> > > > >     0 KSP Residual norm 8.318284791610e+00
> > > > >     1 KSP Residual norm 6.910116607023e-13
> > > > >       Line search: gnorm after quadratic fit 6.148799641401e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.502386288041e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.297879862602e+02 lambda=1.9565738732695813e-02
> > > > >  20 SNES Function norm 4.297879862602e+02
> > > > >     0 KSP Residual norm 7.593855254976e+01
> > > > >     1 KSP Residual norm 1.300639045149e-12
> > > > >       Line search: gnorm after quadratic fit 7.318016560287e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 7.832525429452e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.543595712952e+03 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.752901058720e+02 lambda=1.2500000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.309726315570e+02 lambda=1.2500000000000002e-03
> > > > >       Line search: Cubically determined step, current gnorm 4.297464967378e+02 lambda=2.0986409143217158e-04
> > > > >  21 SNES Function norm 4.297464967378e+02
> > > > >     0 KSP Residual norm 8.492609701850e+00
> > > > >     1 KSP Residual norm 2.830329105288e-13
> > > > >       Line search: gnorm after quadratic fit 6.575200413213e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.532760802544e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 4.268212296998e+02 lambda=1.8460064973695799e-02
> > > > >  22 SNES Function norm 4.268212296998e+02
> > > > >     0 KSP Residual norm 2.527155905469e+00
> > > > >     1 KSP Residual norm 2.235724978394e-13
> > > > >       Line search: gnorm after quadratic fit 3.958132075117e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  23 SNES Function norm 3.958132075117e+02
> > > > >     0 KSP Residual norm 3.425644398425e+00
> > > > >     1 KSP Residual norm 7.166017790799e-14
> > > > >       Line search: gnorm after quadratic fit 3.803199338641e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  24 SNES Function norm 3.803199338641e+02
> > > > >     0 KSP Residual norm 3.581435186688e+00
> > > > >     1 KSP Residual norm 1.730091027109e-13
> > > > >       Line search: gnorm after quadratic fit 3.678433353709e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  25 SNES Function norm 3.678433353709e+02
> > > > >     0 KSP Residual norm 2.817439291614e+00
> > > > >     1 KSP Residual norm 8.592333136485e-13
> > > > >       Line search: gnorm after quadratic fit 3.472671313564e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  26 SNES Function norm 3.472671313564e+02
> > > > >     0 KSP Residual norm 1.830066839089e+00
> > > > >     1 KSP Residual norm 3.248934231575e-14
> > > > >       Line search: gnorm after quadratic fit 3.195079998571e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  27 SNES Function norm 3.195079998571e+02
> > > > >     0 KSP Residual norm 3.775823654589e+00
> > > > >     1 KSP Residual norm 1.316233059492e-13
> > > > >       Line search: gnorm after quadratic fit 3.252874639934e+02
> > > > >       Line search: Cubically determined step, current gnorm 3.125168318399e+02 lambda=5.0000000000000003e-02
> > > > >  28 SNES Function norm 3.125168318399e+02
> > > > >     0 KSP Residual norm 3.892613775622e+00
> > > > >     1 KSP Residual norm 7.093200713216e-11
> > > > >       Line search: gnorm after quadratic fit 3.137590389484e+02
> > > > >       Line search: Cubically determined step, current gnorm 3.045559021319e+02 lambda=5.0000000000000003e-02
> > > > >  29 SNES Function norm 3.045559021319e+02
> > > > >     0 KSP Residual norm 2.364531179390e+00
> > > > >     1 KSP Residual norm 1.615423543115e-12
> > > > >       Line search: gnorm after quadratic fit 2.864449141431e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  30 SNES Function norm 2.864449141431e+02
> > > > >     0 KSP Residual norm 5.187063704081e+00
> > > > >     1 KSP Residual norm 4.254799504045e-13
> > > > >       Line search: gnorm after quadratic fit 3.171238299438e+02
> > > > >       Line search: Cubically determined step, current gnorm 2.859321807369e+02 lambda=4.9550691080557742e-02
> > > > >  31 SNES Function norm 2.859321807369e+02
> > > > >     0 KSP Residual norm 2.400449127985e+01
> > > > >     1 KSP Residual norm 5.221032480453e-13
> > > > >       Line search: gnorm after quadratic fit 1.812538817802e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.647888939229e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.905120335847e+02 lambda=7.3371687404129469e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.857698450683e+02 lambda=1.3231746793531958e-03
> > > > >  32 SNES Function norm 2.857698450683e+02
> > > > >     0 KSP Residual norm 7.438521293559e+00
> > > > >     1 KSP Residual norm 1.293652874831e-12
> > > > >       Line search: gnorm after quadratic fit 3.600694148787e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.932936811991e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.832774969882e+02 lambda=1.8636088141277370e-02
> > > > >  33 SNES Function norm 2.832774969882e+02
> > > > >     0 KSP Residual norm 2.676138557891e+00
> > > > >     1 KSP Residual norm 5.204105674042e-12
> > > > >       Line search: gnorm after quadratic fit 2.698470565576e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  34 SNES Function norm 2.698470565576e+02
> > > > >     0 KSP Residual norm 1.009562156863e+01
> > > > >     1 KSP Residual norm 3.544140587695e-13
> > > > >       Line search: gnorm after quadratic fit 5.672695948559e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.197216154647e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.694068135292e+02 lambda=1.0762403784106927e-02
> > > > >  35 SNES Function norm 2.694068135292e+02
> > > > >     0 KSP Residual norm 3.927549314525e+00
> > > > >     1 KSP Residual norm 2.619134786598e-13
> > > > >       Line search: gnorm after quadratic fit 2.646675787349e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  36 SNES Function norm 2.646675787349e+02
> > > > >     0 KSP Residual norm 3.630219417922e+01
> > > > >     1 KSP Residual norm 1.546302717349e-12
> > > > >       Line search: gnorm after quadratic fit 1.827432666108e+04
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.342968941151e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.008633273369e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.490753494196e+02 lambda=1.2345129233072847e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.682474011960e+02 lambda=3.3291959087019770e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.646227701989e+02 lambda=4.0529212866107715e-04
> > > > >  37 SNES Function norm 2.646227701989e+02
> > > > >     0 KSP Residual norm 8.122774697994e+00
> > > > >     1 KSP Residual norm 1.316362046092e-13
> > > > >       Line search: gnorm after quadratic fit 4.731092571363e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.030088374328e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.637462652877e+02 lambda=9.1057248517509397e-03
> > > > >  38 SNES Function norm 2.637462652877e+02
> > > > >     0 KSP Residual norm 2.601520363063e+00
> > > > >     1 KSP Residual norm 1.060764270007e-12
> > > > >       Line search: gnorm after quadratic fit 2.536856446094e+02
> > > > >       Line search: Quadratically determined step, lambda=1.0000000000000001e-01
> > > > >  39 SNES Function norm 2.536856446094e+02
> > > > >     0 KSP Residual norm 5.447955134327e+01
> > > > >     1 KSP Residual norm 2.556989975730e-12
> > > > >       Line search: gnorm after quadratic fit 9.199250935618e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.998238940105e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 6.227083769291e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.852215674478e+02 lambda=8.5024965111563117e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.537821339347e+02 lambda=8.5024965111563122e-04
> > > > >       Line search: Cubically determined step, current gnorm 2.536484146652e+02 lambda=2.9594242443314720e-04
> > > > >  40 SNES Function norm 2.536484146652e+02
> > > > >     0 KSP Residual norm 3.357359266965e+01
> > > > >     1 KSP Residual norm 8.485166742752e-11
> > > > >       Line search: gnorm after quadratic fit 3.327150899857e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 8.660943990394e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.428891921377e+02 lambda=2.2262238648584176e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.550059920785e+02 lambda=3.9120439672600755e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.535425543202e+02 lambda=8.8078244093807846e-04
> > > > >  41 SNES Function norm 2.535425543202e+02
> > > > >     0 KSP Residual norm 1.982789391732e+01
> > > > >     1 KSP Residual norm 1.156473750758e-11
> > > > >       Line search: gnorm after quadratic fit 9.648483369708e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.023504841042e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.554919970151e+02 lambda=8.7092873850291869e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.532490636747e+02 lambda=2.4564558988847251e-03
> > > > >  42 SNES Function norm 2.532490636747e+02
> > > > >     0 KSP Residual norm 1.762994613361e+01
> > > > >     1 KSP Residual norm 4.550684860593e-12
> > > > >       Line search: gnorm after quadratic fit 6.772107527838e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.370541870778e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.532869639177e+02 lambda=7.7662700854918328e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.527651129995e+02 lambda=3.8686733161537824e-03
> > > > >  43 SNES Function norm 2.527651129995e+02
> > > > >     0 KSP Residual norm 1.559278923951e+01
> > > > >     1 KSP Residual norm 1.060470887103e-11
> > > > >       Line search: gnorm after quadratic fit 7.344361373020e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.498001534426e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.530518382567e+02 lambda=7.6730113050967469e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.523423548732e+02 lambda=3.4156098513217236e-03
> > > > >  44 SNES Function norm 2.523423548732e+02
> > > > >     0 KSP Residual norm 1.190500639095e+01
> > > > >     1 KSP Residual norm 6.311739265577e-12
> > > > >       Line search: gnorm after quadratic fit 4.875196633557e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.933215922947e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.516782376717e+02 lambda=9.8878729665301743e-03
> > > > >  45 SNES Function norm 2.516782376717e+02
> > > > >     0 KSP Residual norm 1.003632001309e+01
> > > > >     1 KSP Residual norm 7.039473047666e-13
> > > > >       Line search: gnorm after quadratic fit 3.417133560813e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.636443273929e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.499729768042e+02 lambda=1.5332495922160540e-02
> > > > >  46 SNES Function norm 2.499729768042e+02
> > > > >     0 KSP Residual norm 5.252587112411e+01
> > > > >     1 KSP Residual norm 2.072629114336e-12
> > > > >       Line search: gnorm after quadratic fit 7.891803681498e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.819076698717e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 5.911876424956e+02 lambda=2.4538865368827514e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.688195882303e+02 lambda=6.5764457912135515e-03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.500076295129e+02 lambda=6.5764457912135523e-04
> > > > >       Line search: Cubically determined step, current gnorm 2.499390700783e+02 lambda=2.7231956365922875e-04
> > > > >  47 SNES Function norm 2.499390700783e+02
> > > > >     0 KSP Residual norm 4.007648116930e+01
> > > > >     1 KSP Residual norm 5.646654234336e-11
> > > > >       Line search: gnorm after quadratic fit 6.276152857499e+03
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 1.418274035925e+03 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 4.785345842563e+02 lambda=2.5000000000000001e-02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 2.665412152699e+02 lambda=8.0607289886479045e-03
> > > > >       Line search: Cubically determined step, current gnorm 2.499109739904e+02 lambda=8.0607289886479045e-04
> > > > >  48 SNES Function norm 2.499109739904e+02
> > > > >     0 KSP Residual norm 1.339756751889e+01
> > > > >     1 KSP Residual norm 2.559175980945e-13
> > > > >       Line search: gnorm after quadratic fit 6.052259901869e+02
> > > > >       Line search: Cubic step no good, shrinking lambda, current gnorm 3.217431518450e+02 lambda=5.0000000000000003e-02
> > > > >       Line search: Cubically determined step, current gnorm 2.496541765916e+02 lambda=7.0239599632283649e-03
> > > > >  49 SNES Function norm 2.496541765916e+02
> > > > >     0 KSP Residual norm 5.340771873687e+00
> > > > >     1 KSP Residual norm 1.207454778077e-13
> > > > >       Line search: gnorm after quadratic fit 2.760767095070e+02
> > > > >       Line search: Cubically determined step, current gnorm 2.491630276458e+02 lambda=5.0000000000000003e-02
> > > > >  50 SNES Function norm 2.491630276458e+02
> > > > >
> > > > >
> > > > > On Sat, Aug 26, 2017 at 11:45 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > > > >
> > > > > > On Aug 26, 2017, at 10:43 PM, zakaryah . <zakaryah at gmail.com> wrote:
> > > > > >
> > > > > > Hi Barry - many thanks for taking the time to understand my many problems and providing so much help.
> > > > > >
> > > > > > The reason I was concerned that I could not alter the linesearch was when I tried to use bt instead of the L-BFGS default, cp, the code crashed with an error like "Could not get Jacobian".  Maybe this is an incompatibility like you say, since L-BFGS only uses the initial Jacobian and I never tried setting the scale type.
> > > > > >
> > > > > > I took your advice and tried to shrink the problem.  First I tried shrinking by a factor 1000 but this converged very quickly with all test data I could provide, including the data which was problematic with the large grid.  So I settled for a reduction in size by a factor 125.  The grid size is 13,230.  This is a decent test case because the solver fails to converge with the options I was using before, and it is small enough that I can run it with the options you suggested (-snes_fd_color -snes_type newtonls -snes_monitor -snes_linesearch_monitor -ksp_monitor -pc_type lu).
> > > > > >
> > > > > > The output is 1800 lines long - how shall I share it?
> > > > >
> > > > >  Just email it, that is small enough for email.
> > > > >
> > > > >   Barry
> > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > On Sat, Aug 26, 2017 at 7:38 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> > > > > >
> > > > > > > On Aug 26, 2017, at 5:56 PM, zakaryah . <zakaryah at gmail.com> wrote:
> > > > > > >
> > > > > > > I'm using PETSc's SNES methods to solve PDEs which result from Euler-Lagrange equations for the strain energy of a 3D displacement field. There is an additional term in the Lagrangian which describes external forces which arise from various data sets, and that term contains nonlinearities (field terms higher than linear). The grid has about 1.6e6 elements, and the displacement field has 3 components at each grid element.
> > > > > > >
> > > > > > > I'm trying to solve a sequence of successively more complicated equations, and the latest equation is failing to converge on some data sets. In particular, the methods were successful for the infinitesimal bulk strain (compression) energy, as well as the full infinitesimal strain energy (bulk + shear), but I'm now trying to generalize to the finite strain, as certain data sets are known to result from displacement fields for which the infinitesimal strain is a poor approximation.
> > > > > > >
> > > > > > > I'm using a DMDA, closely following example 48, and my preferred solver is L-BFGS.
> > > > > >
> > > > > >    So you are using ?
> > > > > >
> > > > > > -snes_type qn -snes_qn_type lbfgs
> > > > > >
> > > > > > >
> > > > > > > I have read the FAQs "Why is Newton's method not converging?​" and "Why is my iterative linear solver not converging?​"​ which have raised a number of questions:
> > > > > >
> > > > > >   Quasi Newton methods either don't use Jacobians or use only the initial Jacobian (the idea behind quasi-Newton methods is to approximate Jacobian information from previous iterations without having the user compute a Jacobian at each iteration). With PETSc's qn it only uses the Jacobian if you use the option
> > > > > >
> > > > > >    -snes_qn_scale_type Jacobian
> > > > > >
> > > > > > otherwise the Jacobian is never computed or used
> > > > > >
> > > > > >
> > > > > > >
> > > > > > > Is there documentation for the DMDA/SNES methods somewhere?  I don't understand these very well.  For example, I am not allocating any matrix for the global Jacobian, and I believe this prevents me from changing the line search.  If I'm mistaken I would love to see an example of changing the line search type while using DMDA/SNES.
> > > > > >
> > > > > >    Whether you provide a Jacobian or not is orthogonal to the line search.
> > > > > >
> > > > > >     You should be able to change the line search with
> > > > > >
> > > > > > -snes_linesearch_type bt  or nleqerr or basic or  l2  or cp
> > > > > >
> > > > > > not all of them may work with qn
> > > > > >
> > > > > >
> > > > > > >
> > > > > > > I don't know how to interpret the linesearch monitor.  Even for problems which are converging properly, the linesearch monitor reports "lssucceed=0" on every iteration.  Is this a problem?
> > > > > >
> > > > > >   It returns a 0 if the line search does not believe it has achieved "sufficient decrease" in the function norm (or possibly some other measure of decrease) you should run -snes_linesearch_monitor also with the option -snes_monitor to see what is happening to the function norm
> > > > > >
> > > > > >    For qn you can add the option
> > > > > >
> > > > > >  -snes_qn_monitor
> > > > > >
> > > > > > to get more detailed monitoring
> > > > > >
> > > > > >
> > > > > > >
> > > > > > > I'm also having trouble understanding the methods for troubleshooting.  I suspect that I've made an error in the analytical Jacobian, which has a rather large number of non-zero elements, but I have no idea how to use -snes_type test -snes_test_display.  The FAQs mention that some troubleshooting tools are more useful for small test problems.  How small is small?
> > > > > >
> > > > > >    Tiny, at most a few dozen rows and columns in the Jacobin.
> > > > > >
> > > > > >    You should run without the -snes_test_display information, what does it say? Does it indicate the Jacobian or report there is likely a problem?
> > > > > >
> > > > > >     With DMDA you can also use -snes_fd_color to have PETSc compute the Jacobian for you instead of using your analytical form. If it works with this, but not your Jacobian then your Jacobian is wrong.
> > > > > >
> > > > > > >  When I try to run the program with -snes_type test -snes_test_display, I get errors like:
> > > > > > >
> > > > > > > [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: Local index 1076396032 too large 4979879 (max) at 0
> > > > > > >
> > > > > > > The second size is 1 less than the number of field elements, while the first number seems too large for any aspect of the problem - the Jacobian has at most 59 non-zero columns per row.
> > > > > > >
> > > > > > > Because I suspect a possible error in the Jacobian, I ran with -snes_mf_operator -pc_type ksp -ksp_ksp_rtol 1e-12 and observed very similar failure to converge (diverging residual) as with the explicit Jacobian.
> > > > > >
> > > > > >    What do you get with -ksp_monitor -ksp_ksp_monitor   it sounds like the true Jacobian is either very ill-conditioned or your function evaluation is wrong.
> > > > > >
> > > > > > > Do I need to set an SNES method which is somehow compatible with the "matrix-free" approach? If I instead use -snes_mf, the problem seems to converge, but horrendously slowly (true residual relative decrease by about 1e-5 per iteration).  I suppose this supports my suspicion that the Jacobian is incorrect but doesn't really suggest a solution.
> > > > > > >
> > > > > > > Is it possible that the analytical Jacobian is correct, but somehow pathological, which causes the SNES to diverge?
> > > > > >
> > > > > >    Yes
> > > > > >
> > > > > > >  Neither the Jacobian nor the function have singularities.
> > > > > > >
> > > > > > > Thanks for any help you can provide!
> > > > > >
> > > > > >     Try really hard to set up a small problem (like just use a very coarse grid) to experiment with as you try to get convergence. Using a big problem for debugging convergence is a recipe for pain.
> > > > > >
> > > > > >    Also since you have a Jacobian I would start with -snes_fd_color -snes_type ls -snes_monitor -snes_linesearch_monitor -ksp_monitor -pc_type lu (not on a huge problem), what happens? Send the output
> > > > > >
> > > > > >    Barry
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > >
> >
> >
>
>

```