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

[Barry Smith]

  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
> >
> >
> >
> 
>