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

zakaryah . zakaryah at gmail.com
Sun Aug 27 18:51:59 CDT 2017


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.



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:
> > > Sure - it was this:
> > >
> > > 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
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > >
> >
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170827/24628550/attachment-0001.html>


More information about the petsc-users mailing list