[petsc-users] large number of iteration.

Josh L ysjosh.lo at gmail.com
Wed Aug 15 19:37:57 CDT 2018


Hi Matt,

Attach is the matrix and rhs in binary format.
I really appreciate your help.


Thanks,
Josh








2018-08-13 10:28 GMT-05:00 Matthew Knepley <knepley at gmail.com>:

> On Mon, Aug 13, 2018 at 11:25 AM Josh L <ysjosh.lo at gmail.com> wrote:
>
>> Hi Matt,
>>
>> Yes, the diagonal term is set to 1.0
>>
>
> Can you send the matrix and rhs? I think you can output them using
>
>  -ksp_view_mat binary:system.bin -ksp_view_rhs binary:system.bin::append
>
> and then we can run exactly what you do in KSP ex10.
>
>   Thanks,
>
>     Matt
>
>
>> Thanks,
>> Josh
>>
>>
>> 2018-08-13 9:13 GMT-05:00 Matthew Knepley <knepley at gmail.com>:
>>
>>> On Sun, Aug 12, 2018 at 11:02 PM Smith, Barry F. <bsmith at mcs.anl.gov>
>>> wrote:
>>>
>>>>
>>>>   Sending this to PETSc users because someone else (Mark?) will know
>>>> why the Eigen estimator fails at iteration 0.
>>>>
>>>
>>> Are you sure its the estimator? That should be GMRES, but the error is
>>> coming from Chebyshev.
>>>
>>> You are remembering to put something on the diagonal of the zeroed rows,
>>> right?
>>>
>>>    Matt
>>>
>>>
>>>>    I notice you are using a real old version of PETSc, for GAMG which
>>>> is a fast moving component of PETSc this may fail. You really should work
>>>> with the latest version of PETSc. Why aren't you?
>>>>
>>>>    Barry
>>>>
>>>>
>>>> > On Aug 12, 2018, at 9:25 PM, Josh L <ysjosh.lo at gmail.com> wrote:
>>>> >
>>>> > Hi Barry,
>>>> >
>>>> > I set the near null space following the example, but I get error when
>>>> solving it with -pc_type gamg(no error if solving with default PC).
>>>> > Moreover, there will be no error if I don't call
>>>> VecSetBlockSize(vec_coord,2) (it is a 2-D problem), but the near null space
>>>> will not affect iteration number.
>>>> >
>>>> > Is it because the way I set Dirichlet boundary conditions?
>>>> > I just zero out the corresponding row and modify right hand size.
>>>> > There is no Neumann bc.
>>>> >
>>>> >
>>>> > Thanks,
>>>> >
>>>> > The following is the error I get
>>>> >
>>>> > [0]PETSC ERROR: --------------------- Error Message
>>>> --------------------------------------------------------------
>>>> > [0]PETSC ERROR: Petsc has generated inconsistent data
>>>> > [0]PETSC ERROR: Eigen estimator failed: DIVERGED_NANORINF at
>>>> iteration 0
>>>> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/
>>>> documentation/faq.html for trouble shooting.
>>>> > [0]PETSC ERROR: Petsc Release Version 3.7.5, Jan, 01, 2017
>>>> > [0]PETSC ERROR: /work/03691/yslo/lonestar/inex/crack on a haswell
>>>> named nid00009 by yslo Sun Aug 12 21:14:17 2018
>>>> > [0]PETSC ERROR: Configure options --with-x=0 -with-pic --with-np=8
>>>> --with-external-packages-dir=/opt/apps/intel16/cray_mpich_7_3/petsc/3.7/externalpackages
>>>> --with-mpi-compilers=1 --with-mpi-dir=/opt/apps/intel16/cray_mpich/7.3.0
>>>> --with-scalar-type=real --with-shared-libraries=1 --with-precision=double
>>>> --with-chaco=1 --download-chaco --with-hypre=1 --download-hypre --with-ml=1
>>>> --download-ml --with-ml=1 --download-ml --with-plapack=1 --download-plapack
>>>> --with-spai=1 --download-spai --with-sundials=1 --download-sundials
>>>> --with-hdf5=1 --with-hdf5-dir=/opt/apps/intel16/cray_mpich_7_2/phdf5/1.8.16/x86_64
>>>> --with-mumps=1 --download-mumps --with-parmetis=1 --download-parmetis
>>>> --with-metis=1 --download-metis --with-scalapack=1 --download-scalapack
>>>> --with-blacs=1 --download-blacs --with-spooles=1 --download-spooles
>>>> --with-suitesparse=1 --download-suitesparse --with-superlu_dist=1
>>>> --download-superlu_dist --with-superlu=1 --download-superlu
>>>> --with-parmetis=1 --download-parmetis --with-metis=1 --download-metis
>>>> --with-debugging=no --with-blas-lapack-dir=/opt/apps/intel/16.0.1.150/
>>>> compilers_and_libraries_2016.1.150/linux/mkl --CFLAGS="-xhost -O2 -g"
>>>> --FFLAGS="-xhost -O2 -g" --CXXFLAGS="-xhost -O2 -g"
>>>> > [0]PETSC ERROR: #1 KSPSolve_Chebyshev() line 440 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/ksp/
>>>> impls/cheby/cheby.c
>>>> > [0]PETSC ERROR: #2 KSPSolve() line 656 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/ksp/
>>>> interface/itfunc.c
>>>> > [0]PETSC ERROR: #3 PCMGMCycle_Private() line 21 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/pc/impls/mg/mg.c
>>>> > [0]PETSC ERROR: #4 PCMGMCycle_Private() line 55 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/pc/impls/mg/mg.c
>>>> > [0]PETSC ERROR: #5 PCMGMCycle_Private() line 55 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/pc/impls/mg/mg.c
>>>> > [0]PETSC ERROR: #6 PCMGMCycle_Private() line 55 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/pc/impls/mg/mg.c
>>>> > [0]PETSC ERROR: #7 PCApply_MG() line 343 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/pc/impls/mg/mg.c
>>>> > [0]PETSC ERROR: #8 PCApply() line 482 in /opt/apps/intel16/cray_mpich_
>>>> 7_3/petsc/3.7/src/ksp/pc/interface/precon.c
>>>> > [0]PETSC ERROR: #9 KSP_PCApply() line 244 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/include/petsc/
>>>> private/kspimpl.h
>>>> > [0]PETSC ERROR: #10 KSPInitialResidual() line 69 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/ksp/
>>>> interface/itres.c
>>>> > [0]PETSC ERROR: #11 KSPSolve_GMRES() line 239 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/ksp/
>>>> impls/gmres/gmres.c
>>>> > [0]PETSC ERROR: #12 KSPSolve() line 656 in
>>>> /opt/apps/intel16/cray_mpich_7_3/petsc/3.7/src/ksp/ksp/
>>>> interface/itfunc.c
>>>> >
>>>> > 2018-08-11 14:27 GMT-05:00 Smith, Barry F. <bsmith at mcs.anl.gov>:
>>>> >
>>>> >
>>>> > > On Aug 9, 2018, at 9:11 PM, Josh L <ysjosh.lo at gmail.com> wrote:
>>>> > >
>>>> > >
>>>> > >
>>>> > > 2018-08-09 18:47 GMT-05:00 Smith, Barry F. <bsmith at mcs.anl.gov>:
>>>> > >
>>>> > >
>>>> > > > On Aug 9, 2018, at 3:03 PM, Josh L <ysjosh.lo at gmail.com> wrote:
>>>> > > >
>>>> > > > Hi Barry,
>>>> > > >
>>>> > > > I use the example ex49 of KSP(linear elasticity problem), and
>>>> it's output of ksp monitor is the same as my reduced problem, so I think my
>>>> code is correct , at least for the reduced problem.
>>>> > >
>>>> > >    Are you using -pc_type gamg ?
>>>> > >
>>>> > >       No, just with GMRES and default PC.
>>>> > >       strangely, with gamg, the final solution is still the same
>>>> but ksp output becomes different between my code and petsc example.
>>>> > >       with -pc_type gamg, my code take 35 iteration and petsc
>>>> example takes only 7.
>>>> >
>>>> >    You should use GAMG for elasticity type problems. Does you code
>>>> set the near null space like the example does?
>>>> >
>>>> > >
>>>> > >      Thanks.
>>>> > >
>>>> > >
>>>> > > >
>>>> > > > Within a Newton's method, several linear system(Ax=b) is solved,
>>>> and b is the residual vector whose norm becomes smaller and smaller.
>>>> > > > Should I change rtol for KSP ? since ||b|| is reducing, and
>>>> converge requires ||b-Ax|| / ||b||  < rtol ( || b - Ax || < atol=10^-50 is
>>>> rarely satisfied )
>>>> > > > I observe that for the first linear solve in Newton's method,
>>>> iterative solver converges with okay number of iteration
>>>> > > > , but for the following linear solve, iterative solver takes more
>>>> and more iterations.
>>>> > >
>>>> > >    You can probably get away with -ksp_rtol 1.e-3 or even try 1.e-2
>>>> > >
>>>> > >     Try also -ksp_pc_side right
>>>> > >
>>>> > >
>>>> > >    Barry
>>>> > >
>>>> > >
>>>> > > >
>>>> > > >
>>>> > > > Thanks,
>>>> > > > Josh
>>>> > > >
>>>> > > >
>>>> > > >
>>>> > > >
>>>> > > > 2018-08-09 10:26 GMT-05:00 Smith, Barry F. <bsmith at mcs.anl.gov>:
>>>> > > >
>>>> > > >    Ok, PCGAMG should work well for linear elasticity so on your
>>>> reduced problem use -pc_type gamg then the group (Mark) will have specific
>>>> suggestions for options. We want to see about 20 iterations for
>>>> "well-behaved" linear elasticity problems.
>>>> > > >
>>>> > > >    Barry
>>>> > > >
>>>> > > >
>>>> > > > > On Aug 8, 2018, at 5:48 PM, Josh L <ysjosh.lo at gmail.com> wrote:
>>>> > > > >
>>>> > > > > Hi Barry,
>>>> > > > >
>>>> > > > > Thanks for your help.
>>>> > > > >
>>>> > > > > I removed the phase field from my formulation by setting all
>>>> corresponding degrees of freedom as Dirichlet boundary condition, so
>>>> > > > > my formulation becomes simply linear elastic.
>>>> > > > > Then I use CG for KSP and algebraic multigrid for
>>>> precondiitoner. For a 150K dofs problem, the linear solver converges with
>>>> 200s iteration(still too much I think), but with GMRES it again takes 2000
>>>> iteration to converge. with GMRES and default PC, it barely converges
>>>> within 10K iteration.
>>>> > > > >
>>>> > > > > I will double check my hessian. below is the  ksp monitor
>>>> output from GMRES and gamg PC for first 500 iteration.
>>>> > > > >
>>>> > > > >  0 KSP preconditioned resid norm 9.175598437886e-07 true resid
>>>> norm 4.423520976130e-07 ||r(i)||/||b|| 1.000000000000e+00
>>>> > > > >   1 KSP preconditioned resid norm 9.171756840415e-07 true resid
>>>> norm 4.185847348836e-07 ||r(i)||/||b|| 9.462704871128e-01
>>>> > > > >   2 KSP preconditioned resid norm 9.171436192115e-07 true resid
>>>> norm 4.373193748329e-07 ||r(i)||/||b|| 9.886228124446e-01
>>>> > > > >   3 KSP preconditioned resid norm 9.047636402661e-07 true resid
>>>> norm 2.655845586443e-07 ||r(i)||/||b|| 6.003917695371e-01
>>>> > > > >   4 KSP preconditioned resid norm 9.032722605019e-07 true resid
>>>> norm 3.045614764490e-07 ||r(i)||/||b|| 6.885046506899e-01
>>>> > > > >   5 KSP preconditioned resid norm 8.926003608855e-07 true resid
>>>> norm 3.498775710352e-07 ||r(i)||/||b|| 7.909481449805e-01
>>>> > > > >   6 KSP preconditioned resid norm 8.788553543841e-07 true resid
>>>> norm 2.052341235884e-07 ||r(i)||/||b|| 4.639610045841e-01
>>>> > > > >   7 KSP preconditioned resid norm 8.725201948474e-07 true resid
>>>> norm 3.315801151567e-07 ||r(i)||/||b|| 7.495841365870e-01
>>>> > > > >   8 KSP preconditioned resid norm 8.704702843041e-07 true resid
>>>> norm 2.453167687304e-07 ||r(i)||/||b|| 5.545735400696e-01
>>>> > > > >   9 KSP preconditioned resid norm 8.673295560357e-07 true resid
>>>> norm 3.480581533019e-07 ||r(i)||/||b|| 7.868350917291e-01
>>>> > > > >  10 KSP preconditioned resid norm 8.650826882165e-07 true resid
>>>> norm 3.331076551593e-07 ||r(i)||/||b|| 7.530373586036e-01
>>>> > > > >  11 KSP preconditioned resid norm 8.632977966037e-07 true resid
>>>> norm 2.943531620635e-07 ||r(i)||/||b|| 6.654273002250e-01
>>>> > > > >  12 KSP preconditioned resid norm 8.622407761375e-07 true resid
>>>> norm 3.413473413778e-07 ||r(i)||/||b|| 7.716643443532e-01
>>>> > > > >  13 KSP preconditioned resid norm 8.612343318250e-07 true resid
>>>> norm 2.962692292083e-07 ||r(i)||/||b|| 6.697588432541e-01
>>>> > > > >  14 KSP preconditioned resid norm 8.579354639983e-07 true resid
>>>> norm 3.298514619407e-07 ||r(i)||/||b|| 7.456762694710e-01
>>>> > > > >  15 KSP preconditioned resid norm 8.578914052098e-07 true resid
>>>> norm 3.287293954137e-07 ||r(i)||/||b|| 7.431396780699e-01
>>>> > > > >  16 KSP preconditioned resid norm 8.564349944906e-07 true resid
>>>> norm 3.231662390796e-07 ||r(i)||/||b|| 7.305633698212e-01
>>>> > > > >  17 KSP preconditioned resid norm 8.544772830576e-07 true resid
>>>> norm 3.522175745975e-07 ||r(i)||/||b|| 7.962380567383e-01
>>>> > > > >  18 KSP preconditioned resid norm 8.518664697660e-07 true resid
>>>> norm 2.957457436154e-07 ||r(i)||/||b|| 6.685754294176e-01
>>>> > > > >  19 KSP preconditioned resid norm 8.502444318835e-07 true resid
>>>> norm 3.231647326483e-07 ||r(i)||/||b|| 7.305599643183e-01
>>>> > > > >  20 KSP preconditioned resid norm 8.488626061229e-07 true resid
>>>> norm 3.034773190650e-07 ||r(i)||/||b|| 6.860537583129e-01
>>>> > > > >  21 KSP preconditioned resid norm 8.472536047506e-07 true resid
>>>> norm 3.142174821592e-07 ||r(i)||/||b|| 7.103334286301e-01
>>>> > > > >  22 KSP preconditioned resid norm 8.465622251926e-07 true resid
>>>> norm 3.109996301173e-07 ||r(i)||/||b|| 7.030590151951e-01
>>>> > > > >  23 KSP preconditioned resid norm 8.455410030133e-07 true resid
>>>> norm 3.026169375479e-07 ||r(i)||/||b|| 6.841087431955e-01
>>>> > > > >  24 KSP preconditioned resid norm 8.441686432481e-07 true resid
>>>> norm 3.247156730398e-07 ||r(i)||/||b|| 7.340660862513e-01
>>>> > > > >  25 KSP preconditioned resid norm 8.422635872670e-07 true resid
>>>> norm 2.949977317008e-07 ||r(i)||/||b|| 6.668844418115e-01
>>>> > > > >  26 KSP preconditioned resid norm 8.391013229654e-07 true resid
>>>> norm 3.503429354312e-07 ||r(i)||/||b|| 7.920001675626e-01
>>>> > > > >  27 KSP preconditioned resid norm 8.369135179081e-07 true resid
>>>> norm 2.932796516614e-07 ||r(i)||/||b|| 6.630004768691e-01
>>>> > > > >  28 KSP preconditioned resid norm 8.344557411519e-07 true resid
>>>> norm 3.290576497205e-07 ||r(i)||/||b|| 7.438817437425e-01
>>>> > > > >  29 KSP preconditioned resid norm 8.255611852636e-07 true resid
>>>> norm 3.228082527538e-07 ||r(i)||/||b|| 7.297540906796e-01
>>>> > > > >  30 KSP preconditioned resid norm 8.191145448615e-07 true resid
>>>> norm 2.695709281706e-07 ||r(i)||/||b|| 6.094035263430e-01
>>>> > > > >  31 KSP preconditioned resid norm 8.190911648438e-07 true resid
>>>> norm 2.661586473845e-07 ||r(i)||/||b|| 6.016895790044e-01
>>>> > > > >  32 KSP preconditioned resid norm 8.188934260789e-07 true resid
>>>> norm 3.000926727757e-07 ||r(i)||/||b|| 6.784022826953e-01
>>>> > > > >  33 KSP preconditioned resid norm 8.166449659697e-07 true resid
>>>> norm 2.081220810731e-07 ||r(i)||/||b|| 4.704896443267e-01
>>>> > > > >  34 KSP preconditioned resid norm 8.130668560598e-07 true resid
>>>> norm 3.175998062677e-07 ||r(i)||/||b|| 7.179796546271e-01
>>>> > > > >  35 KSP preconditioned resid norm 8.103250854289e-07 true resid
>>>> norm 2.675762536241e-07 ||r(i)||/||b|| 6.048942800723e-01
>>>> > > > >  36 KSP preconditioned resid norm 8.079447876546e-07 true resid
>>>> norm 2.624343353237e-07 ||r(i)||/||b|| 5.932702404709e-01
>>>> > > > >  37 KSP preconditioned resid norm 8.064590891986e-07 true resid
>>>> norm 2.891085928184e-07 ||r(i)||/||b|| 6.535712035243e-01
>>>> > > > >  38 KSP preconditioned resid norm 8.039981484367e-07 true resid
>>>> norm 2.435760434885e-07 ||r(i)||/||b|| 5.506383824173e-01
>>>> > > > >  39 KSP preconditioned resid norm 8.030132361660e-07 true resid
>>>> norm 2.791572408348e-07 ||r(i)||/||b|| 6.310747532139e-01
>>>> > > > >  40 KSP preconditioned resid norm 8.025869823713e-07 true resid
>>>> norm 2.673367609434e-07 ||r(i)||/||b|| 6.043528727137e-01
>>>> > > > >  41 KSP preconditioned resid norm 8.011717744772e-07 true resid
>>>> norm 2.626259645676e-07 ||r(i)||/||b|| 5.937034457049e-01
>>>> > > > >  42 KSP preconditioned resid norm 7.994420264672e-07 true resid
>>>> norm 2.931234480129e-07 ||r(i)||/||b|| 6.626473562454e-01
>>>> > > > >  43 KSP preconditioned resid norm 7.976276065412e-07 true resid
>>>> norm 2.565425485040e-07 ||r(i)||/||b|| 5.799510161438e-01
>>>> > > > >  44 KSP preconditioned resid norm 7.963529254232e-07 true resid
>>>> norm 2.683594511633e-07 ||r(i)||/||b|| 6.066648098007e-01
>>>> > > > >  45 KSP preconditioned resid norm 7.951914923472e-07 true resid
>>>> norm 2.726317559874e-07 ||r(i)||/||b|| 6.163229641241e-01
>>>> > > > >  46 KSP preconditioned resid norm 7.943011197327e-07 true resid
>>>> norm 2.561291446355e-07 ||r(i)||/||b|| 5.790164577440e-01
>>>> > > > >  47 KSP preconditioned resid norm 7.926385476568e-07 true resid
>>>> norm 2.978344787798e-07 ||r(i)||/||b|| 6.732973131290e-01
>>>> > > > >  48 KSP preconditioned resid norm 7.909780991424e-07 true resid
>>>> norm 2.498059973427e-07 ||r(i)||/||b|| 5.647220815515e-01
>>>> > > > >  49 KSP preconditioned resid norm 7.863002375447e-07 true resid
>>>> norm 3.283898230799e-07 ||r(i)||/||b|| 7.423720263834e-01
>>>> > > > >  50 KSP preconditioned resid norm 7.778093966862e-07 true resid
>>>> norm 2.908211438411e-07 ||r(i)||/||b|| 6.574426693360e-01
>>>> > > > >  51 KSP preconditioned resid norm 7.727705792892e-07 true resid
>>>> norm 2.662262655905e-07 ||r(i)||/||b|| 6.018424396020e-01
>>>> > > > >  52 KSP preconditioned resid norm 7.580949386500e-07 true resid
>>>> norm 4.077725304210e-07 ||r(i)||/||b|| 9.218279570084e-01
>>>> > > > >  53 KSP preconditioned resid norm 7.223569984502e-07 true resid
>>>> norm 2.806806877450e-07 ||r(i)||/||b|| 6.345187222115e-01
>>>> > > > >  54 KSP preconditioned resid norm 6.833104315428e-07 true resid
>>>> norm 3.167664065828e-07 ||r(i)||/||b|| 7.160956357890e-01
>>>> > > > >  55 KSP preconditioned resid norm 6.480813314286e-07 true resid
>>>> norm 3.449389182846e-07 ||r(i)||/||b|| 7.797836161419e-01
>>>> > > > >  56 KSP preconditioned resid norm 5.628534804207e-07 true resid
>>>> norm 5.008589310744e-07 ||r(i)||/||b|| 1.132263040635e+00
>>>> > > > >  57 KSP preconditioned resid norm 4.676624755957e-07 true resid
>>>> norm 2.747589656053e-07 ||r(i)||/||b|| 6.211318248245e-01
>>>> > > > >  58 KSP preconditioned resid norm 4.266847724428e-07 true resid
>>>> norm 2.051356033363e-07 ||r(i)||/||b|| 4.637382855043e-01
>>>> > > > >  59 KSP preconditioned resid norm 3.963440719602e-07 true resid
>>>> norm 1.462894974865e-07 ||r(i)||/||b|| 3.307082712524e-01
>>>> > > > >  60 KSP preconditioned resid norm 3.710732500804e-07 true resid
>>>> norm 2.006153077396e-07 ||r(i)||/||b|| 4.535195126737e-01
>>>> > > > >  61 KSP preconditioned resid norm 3.700288732490e-07 true resid
>>>> norm 1.637656172799e-07 ||r(i)||/||b|| 3.702155322956e-01
>>>> > > > >  62 KSP preconditioned resid norm 3.681230892718e-07 true resid
>>>> norm 2.598403207310e-07 ||r(i)||/||b|| 5.874061005545e-01
>>>> > > > >  63 KSP preconditioned resid norm 3.447874366524e-07 true resid
>>>> norm 9.834277028700e-08 ||r(i)||/||b|| 2.223178567880e-01
>>>> > > > >  64 KSP preconditioned resid norm 3.308954457293e-07 true resid
>>>> norm 8.462173093522e-08 ||r(i)||/||b|| 1.912994905910e-01
>>>> > > > >  65 KSP preconditioned resid norm 3.238105866460e-07 true resid
>>>> norm 1.508442004566e-07 ||r(i)||/||b|| 3.410048268575e-01
>>>> > > > >  66 KSP preconditioned resid norm 3.149261250545e-07 true resid
>>>> norm 6.146735546031e-08 ||r(i)||/||b|| 1.389557228100e-01
>>>> > > > >  67 KSP preconditioned resid norm 3.080685419183e-07 true resid
>>>> norm 1.469742604456e-07 ||r(i)||/||b|| 3.322562755748e-01
>>>> > > > >  68 KSP preconditioned resid norm 3.017036477950e-07 true resid
>>>> norm 8.246002370676e-08 ||r(i)||/||b|| 1.864126431224e-01
>>>> > > > >  69 KSP preconditioned resid norm 2.951373085123e-07 true resid
>>>> norm 9.600472070625e-08 ||r(i)||/||b|| 2.170323622840e-01
>>>> > > > >  70 KSP preconditioned resid norm 2.854303567501e-07 true resid
>>>> norm 1.353606365274e-07 ||r(i)||/||b|| 3.060020225015e-01
>>>> > > > >  71 KSP preconditioned resid norm 2.720605973489e-07 true resid
>>>> norm 7.105192119695e-08 ||r(i)||/||b|| 1.606230005020e-01
>>>> > > > >  72 KSP preconditioned resid norm 2.679746464067e-07 true resid
>>>> norm 1.093157769657e-07 ||r(i)||/||b|| 2.471239032336e-01
>>>> > > > >  73 KSP preconditioned resid norm 2.560436669192e-07 true resid
>>>> norm 1.046442873472e-07 ||r(i)||/||b|| 2.365633347550e-01
>>>> > > > >  74 KSP preconditioned resid norm 2.486272155786e-07 true resid
>>>> norm 6.788331443968e-08 ||r(i)||/||b|| 1.534599130557e-01
>>>> > > > >  75 KSP preconditioned resid norm 2.450613653917e-07 true resid
>>>> norm 8.385393991346e-08 ||r(i)||/||b|| 1.895637894925e-01
>>>> > > > >  76 KSP preconditioned resid norm 2.381548789947e-07 true resid
>>>> norm 9.639184405623e-08 ||r(i)||/||b|| 2.179075098239e-01
>>>> > > > >  77 KSP preconditioned resid norm 2.318084665470e-07 true resid
>>>> norm 5.817254053914e-08 ||r(i)||/||b|| 1.315073238107e-01
>>>> > > > >  78 KSP preconditioned resid norm 2.301457113250e-07 true resid
>>>> norm 7.214348778091e-08 ||r(i)||/||b|| 1.630906424322e-01
>>>> > > > >  79 KSP preconditioned resid norm 2.288069605204e-07 true resid
>>>> norm 6.520316797685e-08 ||r(i)||/||b|| 1.474010597637e-01
>>>> > > > >  80 KSP preconditioned resid norm 2.272182407377e-07 true resid
>>>> norm 6.553558340584e-08 ||r(i)||/||b|| 1.481525322463e-01
>>>> > > > >  81 KSP preconditioned resid norm 2.261791202938e-07 true resid
>>>> norm 6.410731366977e-08 ||r(i)||/||b|| 1.449237248240e-01
>>>> > > > >  82 KSP preconditioned resid norm 2.250209367027e-07 true resid
>>>> norm 6.108112632184e-08 ||r(i)||/||b|| 1.380825967627e-01
>>>> > > > >  83 KSP preconditioned resid norm 2.237071893420e-07 true resid
>>>> norm 6.644497913123e-08 ||r(i)||/||b|| 1.502083509715e-01
>>>> > > > >  84 KSP preconditioned resid norm 2.234300550868e-07 true resid
>>>> norm 5.793844170260e-08 ||r(i)||/||b|| 1.309781100061e-01
>>>> > > > >  85 KSP preconditioned resid norm 2.229005171876e-07 true resid
>>>> norm 6.926007074308e-08 ||r(i)||/||b|| 1.565722670172e-01
>>>> > > > >  86 KSP preconditioned resid norm 2.225791427950e-07 true resid
>>>> norm 6.289120214757e-08 ||r(i)||/||b|| 1.421745312997e-01
>>>> > > > >  87 KSP preconditioned resid norm 2.224479658210e-07 true resid
>>>> norm 6.242310048343e-08 ||r(i)||/||b|| 1.411163207325e-01
>>>> > > > >  88 KSP preconditioned resid norm 2.217191418500e-07 true resid
>>>> norm 7.465963600280e-08 ||r(i)||/||b|| 1.687787543128e-01
>>>> > > > >  89 KSP preconditioned resid norm 2.208527637102e-07 true resid
>>>> norm 5.960066095409e-08 ||r(i)||/||b|| 1.347357936714e-01
>>>> > > > >  90 KSP preconditioned resid norm 2.204546869923e-07 true resid
>>>> norm 6.334569716267e-08 ||r(i)||/||b|| 1.432019820964e-01
>>>> > > > >  91 KSP preconditioned resid norm 2.204546270123e-07 true resid
>>>> norm 6.327475947831e-08 ||r(i)||/||b|| 1.430416173445e-01
>>>> > > > >  92 KSP preconditioned resid norm 2.204329606031e-07 true resid
>>>> norm 6.775206679705e-08 ||r(i)||/||b|| 1.531632090424e-01
>>>> > > > >  93 KSP preconditioned resid norm 2.203318409249e-07 true resid
>>>> norm 6.059983893344e-08 ||r(i)||/||b|| 1.369945779854e-01
>>>> > > > >  94 KSP preconditioned resid norm 2.202099673600e-07 true resid
>>>> norm 7.148771248833e-08 ||r(i)||/||b|| 1.616081688639e-01
>>>> > > > >  95 KSP preconditioned resid norm 2.200347850341e-07 true resid
>>>> norm 6.555792528003e-08 ||r(i)||/||b|| 1.482030392391e-01
>>>> > > > >  96 KSP preconditioned resid norm 2.196138742262e-07 true resid
>>>> norm 6.418387078605e-08 ||r(i)||/||b|| 1.450967931030e-01
>>>> > > > >  97 KSP preconditioned resid norm 2.193396877997e-07 true resid
>>>> norm 7.242613574407e-08 ||r(i)||/||b|| 1.637296084610e-01
>>>> > > > >  98 KSP preconditioned resid norm 2.186895327006e-07 true resid
>>>> norm 5.717301192712e-08 ||r(i)||/||b|| 1.292477468416e-01
>>>> > > > >  99 KSP preconditioned resid norm 2.183391464984e-07 true resid
>>>> norm 6.921630451950e-08 ||r(i)||/||b|| 1.564733272273e-01
>>>> > > > > 100 KSP preconditioned resid norm 2.182400010423e-07 true resid
>>>> norm 6.548314285564e-08 ||r(i)||/||b|| 1.480339829041e-01
>>>> > > > > 101 KSP preconditioned resid norm 2.177799012029e-07 true resid
>>>> norm 6.272281990396e-08 ||r(i)||/||b|| 1.417938792252e-01
>>>> > > > > 102 KSP preconditioned resid norm 2.169682429176e-07 true resid
>>>> norm 7.907535208604e-08 ||r(i)||/||b|| 1.787611102394e-01
>>>> > > > > 103 KSP preconditioned resid norm 2.164150347952e-07 true resid
>>>> norm 5.949404986863e-08 ||r(i)||/||b|| 1.344947841090e-01
>>>> > > > > 104 KSP preconditioned resid norm 2.163293952130e-07 true resid
>>>> norm 6.398348014122e-08 ||r(i)||/||b|| 1.446437814729e-01
>>>> > > > > 105 KSP preconditioned resid norm 2.160399869141e-07 true resid
>>>> norm 6.102742478406e-08 ||r(i)||/||b|| 1.379611967782e-01
>>>> > > > > 106 KSP preconditioned resid norm 2.152042538469e-07 true resid
>>>> norm 5.769020450098e-08 ||r(i)||/||b|| 1.304169344110e-01
>>>> > > > > 107 KSP preconditioned resid norm 2.150064693173e-07 true resid
>>>> norm 6.321432589105e-08 ||r(i)||/||b|| 1.429049986022e-01
>>>> > > > > 108 KSP preconditioned resid norm 2.144289854705e-07 true resid
>>>> norm 5.801860121360e-08 ||r(i)||/||b|| 1.311593220122e-01
>>>> > > > > 109 KSP preconditioned resid norm 2.131016574445e-07 true resid
>>>> norm 5.968888809818e-08 ||r(i)||/||b|| 1.349352437126e-01
>>>> > > > > 110 KSP preconditioned resid norm 2.131016340852e-07 true resid
>>>> norm 5.978688769843e-08 ||r(i)||/||b|| 1.351567857846e-01
>>>> > > > > 111 KSP preconditioned resid norm 2.128679373780e-07 true resid
>>>> norm 7.339671126222e-08 ||r(i)||/||b|| 1.659237328325e-01
>>>> > > > > 112 KSP preconditioned resid norm 2.124763157021e-07 true resid
>>>> norm 6.095583088857e-08 ||r(i)||/||b|| 1.377993485676e-01
>>>> > > > > 113 KSP preconditioned resid norm 2.114936950060e-07 true resid
>>>> norm 7.680678889067e-08 ||r(i)||/||b|| 1.736326996190e-01
>>>> > > > > 114 KSP preconditioned resid norm 2.111409526125e-07 true resid
>>>> norm 6.861513382847e-08 ||r(i)||/||b|| 1.551142951480e-01
>>>> > > > > 115 KSP preconditioned resid norm 2.104578242179e-07 true resid
>>>> norm 6.666552557721e-08 ||r(i)||/||b|| 1.507069276645e-01
>>>> > > > > 116 KSP preconditioned resid norm 2.095794193794e-07 true resid
>>>> norm 8.409379058143e-08 ||r(i)||/||b|| 1.901060061323e-01
>>>> > > > > 117 KSP preconditioned resid norm 2.088313732014e-07 true resid
>>>> norm 6.461621062155e-08 ||r(i)||/||b|| 1.460741589567e-01
>>>> > > > > 118 KSP preconditioned resid norm 2.076780694200e-07 true resid
>>>> norm 8.889797338722e-08 ||r(i)||/||b|| 2.009665464840e-01
>>>> > > > > 119 KSP preconditioned resid norm 2.059381548473e-07 true resid
>>>> norm 7.343872609510e-08 ||r(i)||/||b|| 1.660187133539e-01
>>>> > > > > 120 KSP preconditioned resid norm 2.036694513953e-07 true resid
>>>> norm 7.626153726963e-08 ||r(i)||/||b|| 1.724000805719e-01
>>>> > > > > 121 KSP preconditioned resid norm 2.036437947430e-07 true resid
>>>> norm 7.403000216323e-08 ||r(i)||/||b|| 1.673553772271e-01
>>>> > > > > 122 KSP preconditioned resid norm 2.034858116471e-07 true resid
>>>> norm 9.140544194375e-08 ||r(i)||/||b|| 2.066350367433e-01
>>>> > > > > 123 KSP preconditioned resid norm 2.013655754252e-07 true resid
>>>> norm 4.519978376641e-08 ||r(i)||/||b|| 1.021805570954e-01
>>>> > > > > 124 KSP preconditioned resid norm 1.998009157253e-07 true resid
>>>> norm 7.924252610751e-08 ||r(i)||/||b|| 1.791390309555e-01
>>>> > > > > 125 KSP preconditioned resid norm 1.978761477601e-07 true resid
>>>> norm 6.193706028697e-08 ||r(i)||/||b|| 1.400175575547e-01
>>>> > > > > 126 KSP preconditioned resid norm 1.966003019512e-07 true resid
>>>> norm 5.417132474463e-08 ||r(i)||/||b|| 1.224620048983e-01
>>>> > > > > 127 KSP preconditioned resid norm 1.957579266801e-07 true resid
>>>> norm 7.621152479630e-08 ||r(i)||/||b|| 1.722870202437e-01
>>>> > > > > 128 KSP preconditioned resid norm 1.950451407619e-07 true resid
>>>> norm 5.429984244536e-08 ||r(i)||/||b|| 1.227525374885e-01
>>>> > > > > 129 KSP preconditioned resid norm 1.948052418686e-07 true resid
>>>> norm 6.458532086853e-08 ||r(i)||/||b|| 1.460043282648e-01
>>>> > > > > 130 KSP preconditioned resid norm 1.946142511609e-07 true resid
>>>> norm 5.988337322193e-08 ||r(i)||/||b|| 1.353749050701e-01
>>>> > > > > 131 KSP preconditioned resid norm 1.942556458837e-07 true resid
>>>> norm 5.860812144886e-08 ||r(i)||/||b|| 1.324920165748e-01
>>>> > > > > 132 KSP preconditioned resid norm 1.935241717066e-07 true resid
>>>> norm 8.129956448465e-08 ||r(i)||/||b|| 1.837892595590e-01
>>>> > > > > 133 KSP preconditioned resid norm 1.924599958511e-07 true resid
>>>> norm 5.461025246990e-08 ||r(i)||/||b|| 1.234542636162e-01
>>>> > > > > 134 KSP preconditioned resid norm 1.920525995377e-07 true resid
>>>> norm 6.091841929726e-08 ||r(i)||/||b|| 1.377147743302e-01
>>>> > > > > 135 KSP preconditioned resid norm 1.912474104817e-07 true resid
>>>> norm 6.146257286777e-08 ||r(i)||/||b|| 1.389449110775e-01
>>>> > > > > 136 KSP preconditioned resid norm 1.899112302838e-07 true resid
>>>> norm 5.363752248243e-08 ||r(i)||/||b|| 1.212552687596e-01
>>>> > > > > 137 KSP preconditioned resid norm 1.871603374495e-07 true resid
>>>> norm 1.033541760028e-07 ||r(i)||/||b|| 2.336468540796e-01
>>>> > > > > 138 KSP preconditioned resid norm 1.848601023477e-07 true resid
>>>> norm 4.885699993060e-08 ||r(i)||/||b|| 1.104482157861e-01
>>>> > > > > 139 KSP preconditioned resid norm 1.827750189949e-07 true resid
>>>> norm 8.191132589117e-08 ||r(i)||/||b|| 1.851722334610e-01
>>>> > > > > 140 KSP preconditioned resid norm 1.793562029798e-07 true resid
>>>> norm 6.003412567942e-08 ||r(i)||/||b|| 1.357157024989e-01
>>>> > > > > 141 KSP preconditioned resid norm 1.762198590085e-07 true resid
>>>> norm 5.955815320828e-08 ||r(i)||/||b|| 1.346396988500e-01
>>>> > > > > 142 KSP preconditioned resid norm 1.735812332311e-07 true resid
>>>> norm 7.966870349397e-08 ||r(i)||/||b|| 1.801024657142e-01
>>>> > > > > 143 KSP preconditioned resid norm 1.688412385715e-07 true resid
>>>> norm 5.688299429377e-08 ||r(i)||/||b|| 1.285921206223e-01
>>>> > > > > 144 KSP preconditioned resid norm 1.661476005287e-07 true resid
>>>> norm 7.082505578352e-08 ||r(i)||/||b|| 1.601101388819e-01
>>>> > > > > 145 KSP preconditioned resid norm 1.610383377332e-07 true resid
>>>> norm 5.920400760597e-08 ||r(i)||/||b|| 1.338391022117e-01
>>>> > > > > 146 KSP preconditioned resid norm 1.544803228895e-07 true resid
>>>> norm 6.998127061136e-08 ||r(i)||/||b|| 1.582026421690e-01
>>>> > > > > 147 KSP preconditioned resid norm 1.526491664454e-07 true resid
>>>> norm 5.645428190677e-08 ||r(i)||/||b|| 1.276229551333e-01
>>>> > > > > 148 KSP preconditioned resid norm 1.514873753547e-07 true resid
>>>> norm 5.864109468824e-08 ||r(i)||/||b|| 1.325665572848e-01
>>>> > > > > 149 KSP preconditioned resid norm 1.504821676672e-07 true resid
>>>> norm 6.213039985019e-08 ||r(i)||/||b|| 1.404546292093e-01
>>>> > > > > 150 KSP preconditioned resid norm 1.486089733443e-07 true resid
>>>> norm 5.614398416742e-08 ||r(i)||/||b|| 1.269214828423e-01
>>>> > > > > 151 KSP preconditioned resid norm 1.485963531298e-07 true resid
>>>> norm 5.477576327252e-08 ||r(i)||/||b|| 1.238284243888e-01
>>>> > > > > 152 KSP preconditioned resid norm 1.485376571004e-07 true resid
>>>> norm 6.371560914841e-08 ||r(i)||/||b|| 1.440382208929e-01
>>>> > > > > 153 KSP preconditioned resid norm 1.474539380136e-07 true resid
>>>> norm 3.640242097893e-08 ||r(i)||/||b|| 8.229286393207e-02
>>>> > > > > 154 KSP preconditioned resid norm 1.468139621410e-07 true resid
>>>> norm 4.869307908912e-08 ||r(i)||/||b|| 1.100776493474e-01
>>>> > > > > 155 KSP preconditioned resid norm 1.460083083948e-07 true resid
>>>> norm 4.630532787515e-08 ||r(i)||/||b|| 1.046797972136e-01
>>>> > > > > 156 KSP preconditioned resid norm 1.443569894841e-07 true resid
>>>> norm 3.603261718441e-08 ||r(i)||/||b|| 8.145686971725e-02
>>>> > > > > 157 KSP preconditioned resid norm 1.423163427629e-07 true resid
>>>> norm 6.442826410227e-08 ||r(i)||/||b|| 1.456492790470e-01
>>>> > > > > 158 KSP preconditioned resid norm 1.409955526676e-07 true resid
>>>> norm 3.683864351150e-08 ||r(i)||/||b|| 8.327900717616e-02
>>>> > > > > 159 KSP preconditioned resid norm 1.392132231651e-07 true resid
>>>> norm 6.633950516791e-08 ||r(i)||/||b|| 1.499699120359e-01
>>>> > > > > 160 KSP preconditioned resid norm 1.373585912794e-07 true resid
>>>> norm 4.602835576127e-08 ||r(i)||/||b|| 1.040536622515e-01
>>>> > > > > 161 KSP preconditioned resid norm 1.358554686788e-07 true resid
>>>> norm 4.404395475469e-08 ||r(i)||/||b|| 9.956764078290e-02
>>>> > > > > 162 KSP preconditioned resid norm 1.356117729059e-07 true resid
>>>> norm 5.220262213455e-08 ||r(i)||/||b|| 1.180114718936e-01
>>>> > > > > 163 KSP preconditioned resid norm 1.345828776136e-07 true resid
>>>> norm 4.364444677578e-08 ||r(i)||/||b|| 9.866449602318e-02
>>>> > > > > 164 KSP preconditioned resid norm 1.335743316358e-07 true resid
>>>> norm 4.365564289163e-08 ||r(i)||/||b|| 9.868980643971e-02
>>>> > > > > 165 KSP preconditioned resid norm 1.335181864567e-07 true resid
>>>> norm 4.195338231508e-08 ||r(i)||/||b|| 9.484160364891e-02
>>>> > > > > 166 KSP preconditioned resid norm 1.334081155239e-07 true resid
>>>> norm 4.320705771627e-08 ||r(i)||/||b|| 9.767571567858e-02
>>>> > > > > 167 KSP preconditioned resid norm 1.330652717803e-07 true resid
>>>> norm 4.671646276077e-08 ||r(i)||/||b|| 1.056092262541e-01
>>>> > > > > 168 KSP preconditioned resid norm 1.327117008156e-07 true resid
>>>> norm 4.076640307294e-08 ||r(i)||/||b|| 9.215826779826e-02
>>>> > > > > 169 KSP preconditioned resid norm 1.323823831437e-07 true resid
>>>> norm 4.609403230274e-08 ||r(i)||/||b|| 1.042021334396e-01
>>>> > > > > 170 KSP preconditioned resid norm 1.321738715249e-07 true resid
>>>> norm 4.115763769210e-08 ||r(i)||/||b|| 9.304270944840e-02
>>>> > > > > 171 KSP preconditioned resid norm 1.317312129565e-07 true resid
>>>> norm 4.587805680503e-08 ||r(i)||/||b|| 1.037138900270e-01
>>>> > > > > 172 KSP preconditioned resid norm 1.315622581976e-07 true resid
>>>> norm 4.346795666887e-08 ||r(i)||/||b|| 9.826551496744e-02
>>>> > > > > 173 KSP preconditioned resid norm 1.313903001504e-07 true resid
>>>> norm 4.268100404068e-08 ||r(i)||/||b|| 9.648649632497e-02
>>>> > > > > 174 KSP preconditioned resid norm 1.310657837357e-07 true resid
>>>> norm 4.888365697755e-08 ||r(i)||/||b|| 1.105084778423e-01
>>>> > > > > 175 KSP preconditioned resid norm 1.308286004287e-07 true resid
>>>> norm 4.218058360777e-08 ||r(i)||/||b|| 9.535522457196e-02
>>>> > > > > 176 KSP preconditioned resid norm 1.305438769953e-07 true resid
>>>> norm 5.026215795117e-08 ||r(i)||/||b|| 1.136247758797e-01
>>>> > > > > 177 KSP preconditioned resid norm 1.303521375300e-07 true resid
>>>> norm 4.518428840914e-08 ||r(i)||/||b|| 1.021455276305e-01
>>>> > > > > 178 KSP preconditioned resid norm 1.302430349652e-07 true resid
>>>> norm 4.699834591242e-08 ||r(i)||/||b|| 1.062464633174e-01
>>>> > > > > 179 KSP preconditioned resid norm 1.299950133828e-07 true resid
>>>> norm 4.780011567673e-08 ||r(i)||/||b|| 1.080589782091e-01
>>>> > > > > 180 KSP preconditioned resid norm 1.297967984970e-07 true resid
>>>> norm 4.439369158202e-08 ||r(i)||/||b|| 1.003582707567e-01
>>>> > > > > 181 KSP preconditioned resid norm 1.297963590960e-07 true resid
>>>> norm 4.419939374367e-08 ||r(i)||/||b|| 9.991903278447e-02
>>>> > > > > 182 KSP preconditioned resid norm 1.297926195884e-07 true resid
>>>> norm 4.616612604888e-08 ||r(i)||/||b|| 1.043651116339e-01
>>>> > > > > 183 KSP preconditioned resid norm 1.297319003442e-07 true resid
>>>> norm 3.981403448576e-08 ||r(i)||/||b|| 9.000530279069e-02
>>>> > > > > 184 KSP preconditioned resid norm 1.296224619385e-07 true resid
>>>> norm 5.101974490162e-08 ||r(i)||/||b|| 1.153374092198e-01
>>>> > > > > 185 KSP preconditioned resid norm 1.294198893701e-07 true resid
>>>> norm 4.272250450111e-08 ||r(i)||/||b|| 9.658031403412e-02
>>>> > > > > 186 KSP preconditioned resid norm 1.292654432794e-07 true resid
>>>> norm 4.341411289813e-08 ||r(i)||/||b|| 9.814379344509e-02
>>>> > > > > 187 KSP preconditioned resid norm 1.291166048974e-07 true resid
>>>> norm 4.903770559588e-08 ||r(i)||/||b|| 1.108567267127e-01
>>>> > > > > 188 KSP preconditioned resid norm 1.287267461311e-07 true resid
>>>> norm 3.838768901036e-08 ||r(i)||/||b|| 8.678084543401e-02
>>>> > > > > 189 KSP preconditioned resid norm 1.285403517589e-07 true resid
>>>> norm 4.836219780019e-08 ||r(i)||/||b|| 1.093296450071e-01
>>>> > > > > 190 KSP preconditioned resid norm 1.283960766544e-07 true resid
>>>> norm 4.373642319734e-08 ||r(i)||/||b|| 9.887242184076e-02
>>>> > > > > 191 KSP preconditioned resid norm 1.281846985469e-07 true resid
>>>> norm 4.196631979171e-08 ||r(i)||/||b|| 9.487085065985e-02
>>>> > > > > 192 KSP preconditioned resid norm 1.279291087436e-07 true resid
>>>> norm 5.017898447715e-08 ||r(i)||/||b|| 1.134367503804e-01
>>>> > > > > 193 KSP preconditioned resid norm 1.277583114588e-07 true resid
>>>> norm 4.217096345089e-08 ||r(i)||/||b|| 9.533347683541e-02
>>>> > > > > 194 KSP preconditioned resid norm 1.275868702813e-07 true resid
>>>> norm 4.784228026600e-08 ||r(i)||/||b|| 1.081542972762e-01
>>>> > > > > 195 KSP preconditioned resid norm 1.273992295764e-07 true resid
>>>> norm 4.379242155010e-08 ||r(i)||/||b|| 9.899901410304e-02
>>>> > > > > 196 KSP preconditioned resid norm 1.272302392639e-07 true resid
>>>> norm 4.236252438739e-08 ||r(i)||/||b|| 9.576652765067e-02
>>>> > > > > 197 KSP preconditioned resid norm 1.269756464572e-07 true resid
>>>> norm 5.265147070813e-08 ||r(i)||/||b|| 1.190261581040e-01
>>>> > > > > 198 KSP preconditioned resid norm 1.266399406059e-07 true resid
>>>> norm 3.658633711709e-08 ||r(i)||/||b|| 8.270863259045e-02
>>>> > > > > 199 KSP preconditioned resid norm 1.262216614405e-07 true resid
>>>> norm 4.823243517156e-08 ||r(i)||/||b|| 1.090362980798e-01
>>>> > > > > 200 KSP preconditioned resid norm 1.256037114182e-07 true resid
>>>> norm 4.118749777942e-08 ||r(i)||/||b|| 9.311021243410e-02
>>>> > > > > 201 KSP preconditioned resid norm 1.252357182659e-07 true resid
>>>> norm 4.297385591425e-08 ||r(i)||/||b|| 9.714852974845e-02
>>>> > > > > 202 KSP preconditioned resid norm 1.242188985573e-07 true resid
>>>> norm 4.735199275670e-08 ||r(i)||/||b|| 1.070459324421e-01
>>>> > > > > 203 KSP preconditioned resid norm 1.231408970960e-07 true resid
>>>> norm 3.893514112125e-08 ||r(i)||/||b|| 8.801843900221e-02
>>>> > > > > 204 KSP preconditioned resid norm 1.220276184666e-07 true resid
>>>> norm 4.828090901382e-08 ||r(i)||/||b|| 1.091458801130e-01
>>>> > > > > 205 KSP preconditioned resid norm 1.200057581901e-07 true resid
>>>> norm 3.901824500784e-08 ||r(i)||/||b|| 8.820630718920e-02
>>>> > > > > 206 KSP preconditioned resid norm 1.190463271301e-07 true resid
>>>> norm 3.812404409291e-08 ||r(i)||/||b|| 8.618483850000e-02
>>>> > > > > 207 KSP preconditioned resid norm 1.176408472739e-07 true resid
>>>> norm 4.100938913911e-08 ||r(i)||/||b|| 9.270757245280e-02
>>>> > > > > 208 KSP preconditioned resid norm 1.166353857162e-07 true resid
>>>> norm 3.610482888020e-08 ||r(i)||/||b|| 8.162011455361e-02
>>>> > > > > 209 KSP preconditioned resid norm 1.159620049487e-07 true resid
>>>> norm 3.566307135323e-08 ||r(i)||/||b|| 8.062145866534e-02
>>>> > > > > 210 KSP preconditioned resid norm 1.146020921987e-07 true resid
>>>> norm 3.592344723620e-08 ||r(i)||/||b|| 8.121007548070e-02
>>>> > > > > 211 KSP preconditioned resid norm 1.145847448331e-07 true resid
>>>> norm 3.473565983922e-08 ||r(i)||/||b|| 7.852491268077e-02
>>>> > > > > 212 KSP preconditioned resid norm 1.145177079586e-07 true resid
>>>> norm 4.122532839361e-08 ||r(i)||/||b|| 9.319573393247e-02
>>>> > > > > 213 KSP preconditioned resid norm 1.132872800274e-07 true resid
>>>> norm 2.141692237998e-08 ||r(i)||/||b|| 4.841600728366e-02
>>>> > > > > 214 KSP preconditioned resid norm 1.123613362087e-07 true resid
>>>> norm 3.246863698820e-08 ||r(i)||/||b|| 7.339998422842e-02
>>>> > > > > 215 KSP preconditioned resid norm 1.111053137881e-07 true resid
>>>> norm 2.903254997529e-08 ||r(i)||/||b|| 6.563221951914e-02
>>>> > > > > 216 KSP preconditioned resid norm 1.099813982910e-07 true resid
>>>> norm 2.173780265025e-08 ||r(i)||/||b|| 4.914140289501e-02
>>>> > > > > 217 KSP preconditioned resid norm 1.089711486676e-07 true resid
>>>> norm 4.266914180936e-08 ||r(i)||/||b|| 9.645968005942e-02
>>>> > > > > 218 KSP preconditioned resid norm 1.077041044884e-07 true resid
>>>> norm 2.174207697189e-08 ||r(i)||/||b|| 4.915106560862e-02
>>>> > > > > 219 KSP preconditioned resid norm 1.068087285500e-07 true resid
>>>> norm 3.726239717827e-08 ||r(i)||/||b|| 8.423696277093e-02
>>>> > > > > 220 KSP preconditioned resid norm 1.060944039406e-07 true resid
>>>> norm 2.612138847397e-08 ||r(i)||/||b|| 5.905112378786e-02
>>>> > > > > 221 KSP preconditioned resid norm 1.052145789908e-07 true resid
>>>> norm 2.848574746885e-08 ||r(i)||/||b|| 6.439609447443e-02
>>>> > > > > 222 KSP preconditioned resid norm 1.034013372786e-07 true resid
>>>> norm 4.496176952774e-08 ||r(i)||/||b|| 1.016424919659e-01
>>>> > > > > 223 KSP preconditioned resid norm 1.016447882738e-07 true resid
>>>> norm 2.629128290689e-08 ||r(i)||/||b|| 5.943519438194e-02
>>>> > > > > 224 KSP preconditioned resid norm 1.002724751991e-07 true resid
>>>> norm 3.297488178025e-08 ||r(i)||/||b|| 7.454442277585e-02
>>>> > > > > 225 KSP preconditioned resid norm 9.969279422176e-08 true resid
>>>> norm 3.069952269998e-08 ||r(i)||/||b|| 6.940064908846e-02
>>>> > > > > 226 KSP preconditioned resid norm 9.811249548568e-08 true resid
>>>> norm 3.219822477037e-08 ||r(i)||/||b|| 7.278867884683e-02
>>>> > > > > 227 KSP preconditioned resid norm 9.637808562080e-08 true resid
>>>> norm 4.091298315961e-08 ||r(i)||/||b|| 9.248963298779e-02
>>>> > > > > 228 KSP preconditioned resid norm 9.526012021860e-08 true resid
>>>> norm 2.720466192827e-08 ||r(i)||/||b|| 6.150001791574e-02
>>>> > > > > 229 KSP preconditioned resid norm 9.392714618074e-08 true resid
>>>> norm 4.018517470597e-08 ||r(i)||/||b|| 9.084431818639e-02
>>>> > > > > 230 KSP preconditioned resid norm 9.288540879331e-08 true resid
>>>> norm 2.742429822661e-08 ||r(i)||/||b|| 6.199653709024e-02
>>>> > > > > 231 KSP preconditioned resid norm 9.182683719708e-08 true resid
>>>> norm 3.418196050419e-08 ||r(i)||/||b|| 7.727319637148e-02
>>>> > > > > 232 KSP preconditioned resid norm 9.152618596319e-08 true resid
>>>> norm 3.061702948039e-08 ||r(i)||/||b|| 6.921416140130e-02
>>>> > > > > 233 KSP preconditioned resid norm 9.028012081487e-08 true resid
>>>> norm 3.309940751576e-08 ||r(i)||/||b|| 7.482593095946e-02
>>>> > > > > 234 KSP preconditioned resid norm 8.964512632411e-08 true resid
>>>> norm 3.657895201743e-08 ||r(i)||/||b|| 8.269193751948e-02
>>>> > > > > 235 KSP preconditioned resid norm 8.927591956786e-08 true resid
>>>> norm 2.996402299113e-08 ||r(i)||/||b|| 6.773794710778e-02
>>>> > > > > 236 KSP preconditioned resid norm 8.863790942206e-08 true resid
>>>> norm 3.770401488950e-08 ||r(i)||/||b|| 8.523530258577e-02
>>>> > > > > 237 KSP preconditioned resid norm 8.842416675824e-08 true resid
>>>> norm 3.294247881070e-08 ||r(i)||/||b|| 7.447117124223e-02
>>>> > > > > 238 KSP preconditioned resid norm 8.793063593991e-08 true resid
>>>> norm 3.836231748404e-08 ||r(i)||/||b|| 8.672348948056e-02
>>>> > > > > 239 KSP preconditioned resid norm 8.778266322588e-08 true resid
>>>> norm 3.572179823322e-08 ||r(i)||/||b|| 8.075421915252e-02
>>>> > > > > 240 KSP preconditioned resid norm 8.738590572964e-08 true resid
>>>> norm 3.574931934089e-08 ||r(i)||/||b|| 8.081643454117e-02
>>>> > > > > 241 KSP preconditioned resid norm 8.738404146064e-08 true resid
>>>> norm 3.532576961772e-08 ||r(i)||/||b|| 7.985893998998e-02
>>>> > > > > 242 KSP preconditioned resid norm 8.737492315747e-08 true resid
>>>> norm 3.850858638399e-08 ||r(i)||/||b|| 8.705415118814e-02
>>>> > > > > 243 KSP preconditioned resid norm 8.704894133677e-08 true resid
>>>> norm 2.725999899728e-08 ||r(i)||/||b|| 6.162511525180e-02
>>>> > > > > 244 KSP preconditioned resid norm 8.691957267756e-08 true resid
>>>> norm 3.445507385918e-08 ||r(i)||/||b|| 7.789060805883e-02
>>>> > > > > 245 KSP preconditioned resid norm 8.664547872292e-08 true resid
>>>> norm 3.233777274784e-08 ||r(i)||/||b|| 7.310414695067e-02
>>>> > > > > 246 KSP preconditioned resid norm 8.647910251595e-08 true resid
>>>> norm 2.869258593620e-08 ||r(i)||/||b|| 6.486368232688e-02
>>>> > > > > 247 KSP preconditioned resid norm 8.602870326813e-08 true resid
>>>> norm 3.830342792446e-08 ||r(i)||/||b|| 8.659036123294e-02
>>>> > > > > 248 KSP preconditioned resid norm 8.584186792773e-08 true resid
>>>> norm 2.914852049114e-08 ||r(i)||/||b|| 6.589438740865e-02
>>>> > > > > 249 KSP preconditioned resid norm 8.568202493555e-08 true resid
>>>> norm 3.726679821110e-08 ||r(i)||/||b|| 8.424691193326e-02
>>>> > > > > 250 KSP preconditioned resid norm 8.535775976392e-08 true resid
>>>> norm 3.374879488937e-08 ||r(i)||/||b|| 7.629396372591e-02
>>>> > > > > 251 KSP preconditioned resid norm 8.505569181999e-08 true resid
>>>> norm 3.019219829750e-08 ||r(i)||/||b|| 6.825376992767e-02
>>>> > > > > 252 KSP preconditioned resid norm 8.489565939881e-08 true resid
>>>> norm 3.532505625863e-08 ||r(i)||/||b|| 7.985732734003e-02
>>>> > > > > 253 KSP preconditioned resid norm 8.454497959931e-08 true resid
>>>> norm 2.833135878246e-08 ||r(i)||/||b|| 6.404707683166e-02
>>>> > > > > 254 KSP preconditioned resid norm 8.448142425824e-08 true resid
>>>> norm 3.023601426215e-08 ||r(i)||/||b|| 6.835282216432e-02
>>>> > > > > 255 KSP preconditioned resid norm 8.447314709431e-08 true resid
>>>> norm 2.955565650261e-08 ||r(i)||/||b|| 6.681477642380e-02
>>>> > > > > 256 KSP preconditioned resid norm 8.434191340398e-08 true resid
>>>> norm 3.011718192322e-08 ||r(i)||/||b|| 6.808418471561e-02
>>>> > > > > 257 KSP preconditioned resid norm 8.420519688672e-08 true resid
>>>> norm 3.194610164662e-08 ||r(i)||/||b|| 7.221871857057e-02
>>>> > > > > 258 KSP preconditioned resid norm 8.414920999588e-08 true resid
>>>> norm 2.927199782498e-08 ||r(i)||/||b|| 6.617352553076e-02
>>>> > > > > 259 KSP preconditioned resid norm 8.411142889741e-08 true resid
>>>> norm 3.116811073402e-08 ||r(i)||/||b|| 7.045995916421e-02
>>>> > > > > 260 KSP preconditioned resid norm 8.393401673483e-08 true resid
>>>> norm 2.913633916610e-08 ||r(i)||/||b|| 6.586684978622e-02
>>>> > > > > 261 KSP preconditioned resid norm 8.386932870420e-08 true resid
>>>> norm 2.931516814567e-08 ||r(i)||/||b|| 6.627111819715e-02
>>>> > > > > 262 KSP preconditioned resid norm 8.383102906154e-08 true resid
>>>> norm 2.947670084156e-08 ||r(i)||/||b|| 6.663628589221e-02
>>>> > > > > 263 KSP preconditioned resid norm 8.371299202072e-08 true resid
>>>> norm 2.900846382950e-08 ||r(i)||/||b|| 6.557776935168e-02
>>>> > > > > 264 KSP preconditioned resid norm 8.354612585259e-08 true resid
>>>> norm 3.153902284033e-08 ||r(i)||/||b|| 7.129845887590e-02
>>>> > > > > 265 KSP preconditioned resid norm 8.348749379055e-08 true resid
>>>> norm 2.912311287048e-08 ||r(i)||/||b|| 6.583694985881e-02
>>>> > > > > 266 KSP preconditioned resid norm 8.328479663117e-08 true resid
>>>> norm 3.201390360666e-08 ||r(i)||/||b|| 7.237199457042e-02
>>>> > > > > 267 KSP preconditioned resid norm 8.304748250656e-08 true resid
>>>> norm 3.103708076427e-08 ||r(i)||/||b|| 7.016374723157e-02
>>>> > > > > 268 KSP preconditioned resid norm 8.294564373794e-08 true resid
>>>> norm 2.975578688471e-08 ||r(i)||/||b|| 6.726719969290e-02
>>>> > > > > 269 KSP preconditioned resid norm 8.286083099532e-08 true resid
>>>> norm 3.122272154943e-08 ||r(i)||/||b|| 7.058341470044e-02
>>>> > > > > 270 KSP preconditioned resid norm 8.249515525349e-08 true resid
>>>> norm 2.958445476571e-08 ||r(i)||/||b|| 6.687987900443e-02
>>>> > > > > 271 KSP preconditioned resid norm 8.249243482197e-08 true resid
>>>> norm 2.915759140675e-08 ||r(i)||/||b|| 6.591489350698e-02
>>>> > > > > 272 KSP preconditioned resid norm 8.247967811026e-08 true resid
>>>> norm 3.226477388325e-08 ||r(i)||/||b|| 7.293912260697e-02
>>>> > > > > 273 KSP preconditioned resid norm 8.225540420138e-08 true resid
>>>> norm 2.335784346654e-08 ||r(i)||/||b|| 5.280373619246e-02
>>>> > > > > 274 KSP preconditioned resid norm 8.208649428128e-08 true resid
>>>> norm 3.103095297226e-08 ||r(i)||/||b|| 7.014989448385e-02
>>>> > > > > 275 KSP preconditioned resid norm 8.177662092059e-08 true resid
>>>> norm 2.780330970969e-08 ||r(i)||/||b|| 6.285334659815e-02
>>>> > > > > 276 KSP preconditioned resid norm 8.156719698267e-08 true resid
>>>> norm 2.529092246235e-08 ||r(i)||/||b|| 5.717373693676e-02
>>>> > > > > 277 KSP preconditioned resid norm 8.129665601936e-08 true resid
>>>> norm 3.216359520748e-08 ||r(i)||/||b|| 7.271039378141e-02
>>>> > > > > 278 KSP preconditioned resid norm 8.111942990953e-08 true resid
>>>> norm 2.578429842977e-08 ||r(i)||/||b|| 5.828908367091e-02
>>>> > > > > 279 KSP preconditioned resid norm 8.103355449576e-08 true resid
>>>> norm 3.114448259432e-08 ||r(i)||/||b|| 7.040654438485e-02
>>>> > > > > 280 KSP preconditioned resid norm 8.096927476093e-08 true resid
>>>> norm 2.905799246198e-08 ||r(i)||/||b|| 6.568973588864e-02
>>>> > > > > 281 KSP preconditioned resid norm 8.069018681340e-08 true resid
>>>> norm 2.806750078967e-08 ||r(i)||/||b|| 6.345058821044e-02
>>>> > > > > 282 KSP preconditioned resid norm 8.047055052412e-08 true resid
>>>> norm 3.356282539492e-08 ||r(i)||/||b|| 7.587355316280e-02
>>>> > > > > 283 KSP preconditioned resid norm 8.023864541552e-08 true resid
>>>> norm 2.697346140135e-08 ||r(i)||/||b|| 6.097735615340e-02
>>>> > > > > 284 KSP preconditioned resid norm 8.014006658762e-08 true resid
>>>> norm 3.033393715647e-08 ||r(i)||/||b|| 6.857419083160e-02
>>>> > > > > 285 KSP preconditioned resid norm 8.000866139550e-08 true resid
>>>> norm 2.744969989335e-08 ||r(i)||/||b|| 6.205396118040e-02
>>>> > > > > 286 KSP preconditioned resid norm 7.963315124505e-08 true resid
>>>> norm 2.853888338671e-08 ||r(i)||/||b|| 6.451621579442e-02
>>>> > > > > 287 KSP preconditioned resid norm 7.901906690217e-08 true resid
>>>> norm 3.213929486191e-08 ||r(i)||/||b|| 7.265545938481e-02
>>>> > > > > 288 KSP preconditioned resid norm 7.847487693433e-08 true resid
>>>> norm 2.531887675170e-08 ||r(i)||/||b|| 5.723693159436e-02
>>>> > > > > 289 KSP preconditioned resid norm 7.776109704697e-08 true resid
>>>> norm 3.151424977919e-08 ||r(i)||/||b|| 7.124245583834e-02
>>>> > > > > 290 KSP preconditioned resid norm 7.697887176101e-08 true resid
>>>> norm 2.452083040958e-08 ||r(i)||/||b|| 5.543283402950e-02
>>>> > > > > 291 KSP preconditioned resid norm 7.630186207736e-08 true resid
>>>> norm 2.646550673478e-08 ||r(i)||/||b|| 5.982905219075e-02
>>>> > > > > 292 KSP preconditioned resid norm 7.539569101676e-08 true resid
>>>> norm 2.639665425804e-08 ||r(i)||/||b|| 5.967340134812e-02
>>>> > > > > 293 KSP preconditioned resid norm 7.387540581378e-08 true resid
>>>> norm 2.766591414879e-08 ||r(i)||/||b|| 6.254274433891e-02
>>>> > > > > 294 KSP preconditioned resid norm 7.148854762848e-08 true resid
>>>> norm 4.506538736910e-08 ||r(i)||/||b|| 1.018767348732e-01
>>>> > > > > 295 KSP preconditioned resid norm 6.875818761377e-08 true resid
>>>> norm 1.988285545247e-08 ||r(i)||/||b|| 4.494803022244e-02
>>>> > > > > 296 KSP preconditioned resid norm 6.720148343751e-08 true resid
>>>> norm 2.185917455369e-08 ||r(i)||/||b|| 4.941578139143e-02
>>>> > > > > 297 KSP preconditioned resid norm 6.559210493364e-08 true resid
>>>> norm 2.579380002545e-08 ||r(i)||/||b|| 5.831056338297e-02
>>>> > > > > 298 KSP preconditioned resid norm 6.296172453909e-08 true resid
>>>> norm 2.690729261887e-08 ||r(i)||/||b|| 6.082777218435e-02
>>>> > > > > 299 KSP preconditioned resid norm 6.126664666716e-08 true resid
>>>> norm 2.523269850506e-08 ||r(i)||/||b|| 5.704211337804e-02
>>>> > > > > 300 KSP preconditioned resid norm 5.996078539496e-08 true resid
>>>> norm 1.432165407581e-08 ||r(i)||/||b|| 3.237614143370e-02
>>>> > > > > 301 KSP preconditioned resid norm 5.994360410331e-08 true resid
>>>> norm 1.354971586959e-08 ||r(i)||/||b|| 3.063106503328e-02
>>>> > > > > 302 KSP preconditioned resid norm 5.968927298570e-08 true resid
>>>> norm 2.119292714347e-08 ||r(i)||/||b|| 4.790963410783e-02
>>>> > > > > 303 KSP preconditioned resid norm 5.871687993468e-08 true resid
>>>> norm 8.809042798461e-09 ||r(i)||/||b|| 1.991409749382e-02
>>>> > > > > 304 KSP preconditioned resid norm 5.799650512507e-08 true resid
>>>> norm 1.745679956812e-08 ||r(i)||/||b|| 3.946358491871e-02
>>>> > > > > 305 KSP preconditioned resid norm 5.745928074883e-08 true resid
>>>> norm 1.274997712124e-08 ||r(i)||/||b|| 2.882314154277e-02
>>>> > > > > 306 KSP preconditioned resid norm 5.632393272566e-08 true resid
>>>> norm 1.202295750268e-08 ||r(i)||/||b|| 2.717961001553e-02
>>>> > > > > 307 KSP preconditioned resid norm 5.396643228430e-08 true resid
>>>> norm 2.502414535309e-08 ||r(i)||/||b|| 5.657064923649e-02
>>>> > > > > 308 KSP preconditioned resid norm 5.160090231555e-08 true resid
>>>> norm 1.262817440249e-08 ||r(i)||/||b|| 2.854778912689e-02
>>>> > > > > 309 KSP preconditioned resid norm 4.907022431548e-08 true resid
>>>> norm 2.171207321606e-08 ||r(i)||/||b|| 4.908323783977e-02
>>>> > > > > 310 KSP preconditioned resid norm 4.776335659862e-08 true resid
>>>> norm 1.151540470517e-08 ||r(i)||/||b|| 2.603221453522e-02
>>>> > > > > 311 KSP preconditioned resid norm 4.697621190386e-08 true resid
>>>> norm 1.508650195148e-08 ||r(i)||/||b|| 3.410518913076e-02
>>>> > > > > 312 KSP preconditioned resid norm 4.595176479120e-08 true resid
>>>> norm 1.933671761706e-08 ||r(i)||/||b|| 4.371340776138e-02
>>>> > > > > 313 KSP preconditioned resid norm 4.509046368009e-08 true resid
>>>> norm 1.186794873424e-08 ||r(i)||/||b|| 2.682919058886e-02
>>>> > > > > 314 KSP preconditioned resid norm 4.450909170016e-08 true resid
>>>> norm 1.579009874879e-08 ||r(i)||/||b|| 3.569577003024e-02
>>>> > > > > 315 KSP preconditioned resid norm 4.337869810365e-08 true resid
>>>> norm 1.647588479037e-08 ||r(i)||/||b|| 3.724608717643e-02
>>>> > > > > 316 KSP preconditioned resid norm 4.223810751226e-08 true resid
>>>> norm 1.498137614524e-08 ||r(i)||/||b|| 3.386753725389e-02
>>>> > > > > 317 KSP preconditioned resid norm 4.156026722749e-08 true resid
>>>> norm 1.333299778266e-08 ||r(i)||/||b|| 3.014114289185e-02
>>>> > > > > 318 KSP preconditioned resid norm 4.099199515842e-08 true resid
>>>> norm 1.310678580056e-08 ||r(i)||/||b|| 2.962975844646e-02
>>>> > > > > 319 KSP preconditioned resid norm 4.066386009421e-08 true resid
>>>> norm 1.506644976724e-08 ||r(i)||/||b|| 3.405985830867e-02
>>>> > > > > 320 KSP preconditioned resid norm 4.023087548652e-08 true resid
>>>> norm 1.240117063260e-08 ||r(i)||/||b|| 2.803461473228e-02
>>>> > > > > 321 KSP preconditioned resid norm 3.991268177919e-08 true resid
>>>> norm 1.312932289209e-08 ||r(i)||/||b|| 2.968070675585e-02
>>>> > > > > 322 KSP preconditioned resid norm 3.956889976017e-08 true resid
>>>> norm 1.421972275585e-08 ||r(i)||/||b|| 3.214571114861e-02
>>>> > > > > 323 KSP preconditioned resid norm 3.900185975518e-08 true resid
>>>> norm 1.347489145169e-08 ||r(i)||/||b|| 3.046191376598e-02
>>>> > > > > 324 KSP preconditioned resid norm 3.877682049038e-08 true resid
>>>> norm 1.279539901604e-08 ||r(i)||/||b|| 2.892582421352e-02
>>>> > > > > 325 KSP preconditioned resid norm 3.869150965109e-08 true resid
>>>> norm 1.224801722014e-08 ||r(i)||/||b|| 2.768838960238e-02
>>>> > > > > 326 KSP preconditioned resid norm 3.851624335552e-08 true resid
>>>> norm 1.309261686557e-08 ||r(i)||/||b|| 2.959772754831e-02
>>>> > > > > 327 KSP preconditioned resid norm 3.836778709819e-08 true resid
>>>> norm 1.256963552910e-08 ||r(i)||/||b|| 2.841545365542e-02
>>>> > > > > 328 KSP preconditioned resid norm 3.810190892142e-08 true resid
>>>> norm 1.532295260746e-08 ||r(i)||/||b|| 3.463971955858e-02
>>>> > > > > 329 KSP preconditioned resid norm 3.779873511822e-08 true resid
>>>> norm 1.357026597966e-08 ||r(i)||/||b|| 3.067752148771e-02
>>>> > > > > 330 KSP preconditioned resid norm 3.758209190896e-08 true resid
>>>> norm 1.326727470061e-08 ||r(i)||/||b|| 2.999256649217e-02
>>>> > > > > 331 KSP preconditioned resid norm 3.758046617049e-08 true resid
>>>> norm 1.304204653926e-08 ||r(i)||/||b|| 2.948340611389e-02
>>>> > > > > 332 KSP preconditioned resid norm 3.756576997625e-08 true resid
>>>> norm 1.547320887283e-08 ||r(i)||/||b|| 3.497939527431e-02
>>>> > > > > 333 KSP preconditioned resid norm 3.735750946435e-08 true resid
>>>> norm 9.401475437163e-09 ||r(i)||/||b|| 2.125337595977e-02
>>>> > > > > 334 KSP preconditioned resid norm 3.726816411327e-08 true resid
>>>> norm 1.293737018744e-08 ||r(i)||/||b|| 2.924677029283e-02
>>>> > > > > 335 KSP preconditioned resid norm 3.715836120670e-08 true resid
>>>> norm 1.137076112403e-08 ||r(i)||/||b|| 2.570522709260e-02
>>>> > > > > 336 KSP preconditioned resid norm 3.711352115906e-08 true resid
>>>> norm 1.100868545692e-08 ||r(i)||/||b|| 2.488670341188e-02
>>>> > > > > 337 KSP preconditioned resid norm 3.697154209112e-08 true resid
>>>> norm 1.292178506561e-08 ||r(i)||/||b|| 2.921153790237e-02
>>>> > > > > 338 KSP preconditioned resid norm 3.687250442553e-08 true resid
>>>> norm 1.038256633339e-08 ||r(i)||/||b|| 2.347127184298e-02
>>>> > > > > 339 KSP preconditioned resid norm 3.681163672322e-08 true resid
>>>> norm 1.217819776625e-08 ||r(i)||/||b|| 2.753055276999e-02
>>>> > > > > 340 KSP preconditioned resid norm 3.676709918842e-08 true resid
>>>> norm 1.163124964184e-08 ||r(i)||/||b|| 2.629409853509e-02
>>>> > > > > 341 KSP preconditioned resid norm 3.664698023565e-08 true resid
>>>> norm 1.141439938468e-08 ||r(i)||/||b|| 2.580387760401e-02
>>>> > > > > 342 KSP preconditioned resid norm 3.656190729093e-08 true resid
>>>> norm 1.232477473873e-08 ||r(i)||/||b|| 2.786191091946e-02
>>>> > > > > 343 KSP preconditioned resid norm 3.650243259584e-08 true resid
>>>> norm 1.132768943490e-08 ||r(i)||/||b|| 2.560785739692e-02
>>>> > > > > 344 KSP preconditioned resid norm 3.650130795451e-08 true resid
>>>> norm 1.141084292815e-08 ||r(i)||/||b|| 2.579583772683e-02
>>>> > > > > 345 KSP preconditioned resid norm 3.647229125147e-08 true resid
>>>> norm 1.138306478808e-08 ||r(i)||/||b|| 2.573304127980e-02
>>>> > > > > 346 KSP preconditioned resid norm 3.643856164392e-08 true resid
>>>> norm 1.116515940173e-08 ||r(i)||/||b|| 2.524043507871e-02
>>>> > > > > 347 KSP preconditioned resid norm 3.639892567434e-08 true resid
>>>> norm 1.174359887576e-08 ||r(i)||/||b|| 2.654808000036e-02
>>>> > > > > 348 KSP preconditioned resid norm 3.632894120883e-08 true resid
>>>> norm 1.064588448828e-08 ||r(i)||/||b|| 2.406654008362e-02
>>>> > > > > 349 KSP preconditioned resid norm 3.632109543324e-08 true resid
>>>> norm 1.070526312693e-08 ||r(i)||/||b|| 2.420077396422e-02
>>>> > > > > 350 KSP preconditioned resid norm 3.628173432975e-08 true resid
>>>> norm 1.071879321292e-08 ||r(i)||/||b|| 2.423136065311e-02
>>>> > > > > 351 KSP preconditioned resid norm 3.621714360676e-08 true resid
>>>> norm 1.066953807038e-08 ||r(i)||/||b|| 2.412001237918e-02
>>>> > > > > 352 KSP preconditioned resid norm 3.614598721323e-08 true resid
>>>> norm 1.067226627219e-08 ||r(i)||/||b|| 2.412617986843e-02
>>>> > > > > 353 KSP preconditioned resid norm 3.608734513646e-08 true resid
>>>> norm 1.072773727635e-08 ||r(i)||/||b|| 2.425157998401e-02
>>>> > > > > 354 KSP preconditioned resid norm 3.593202477726e-08 true resid
>>>> norm 1.166649044369e-08 ||r(i)||/||b|| 2.637376539332e-02
>>>> > > > > 355 KSP preconditioned resid norm 3.574816645908e-08 true resid
>>>> norm 1.093827518487e-08 ||r(i)||/||b|| 2.472753095079e-02
>>>> > > > > 356 KSP preconditioned resid norm 3.551731486762e-08 true resid
>>>> norm 1.076438469683e-08 ||r(i)||/||b|| 2.433442670423e-02
>>>> > > > > 357 KSP preconditioned resid norm 3.526513581484e-08 true resid
>>>> norm 1.110249798589e-08 ||r(i)||/||b|| 2.509878001212e-02
>>>> > > > > 358 KSP preconditioned resid norm 3.465142744209e-08 true resid
>>>> norm 1.363275954561e-08 ||r(i)||/||b|| 3.081879710568e-02
>>>> > > > > 359 KSP preconditioned resid norm 3.445104570595e-08 true resid
>>>> norm 1.002161340638e-08 ||r(i)||/||b|| 2.265528627639e-02
>>>> > > > > 360 KSP preconditioned resid norm 3.412299596043e-08 true resid
>>>> norm 1.132489003434e-08 ||r(i)||/||b|| 2.560152895272e-02
>>>> > > > > 361 KSP preconditioned resid norm 3.411529696237e-08 true resid
>>>> norm 1.086147186780e-08 ||r(i)||/||b|| 2.455390609972e-02
>>>> > > > > 362 KSP preconditioned resid norm 3.406577468016e-08 true resid
>>>> norm 1.456732293497e-08 ||r(i)||/||b|| 3.293151092439e-02
>>>> > > > > 363 KSP preconditioned resid norm 3.365555787056e-08 true resid
>>>> norm 7.711588827294e-09 ||r(i)||/||b|| 1.743314628529e-02
>>>> > > > > 364 KSP preconditioned resid norm 3.345269937851e-08 true resid
>>>> norm 1.136806437719e-08 ||r(i)||/||b|| 2.569913071178e-02
>>>> > > > > 365 KSP preconditioned resid norm 3.314065779377e-08 true resid
>>>> norm 9.636397417285e-09 ||r(i)||/||b|| 2.178445059780e-02
>>>> > > > > 366 KSP preconditioned resid norm 3.276476852537e-08 true resid
>>>> norm 8.198503499564e-09 ||r(i)||/||b|| 1.853388634033e-02
>>>> > > > > 367 KSP preconditioned resid norm 3.259579551722e-08 true resid
>>>> norm 1.047105927610e-08 ||r(i)||/||b|| 2.367132275987e-02
>>>> > > > > 368 KSP preconditioned resid norm 3.250921252208e-08 true resid
>>>> norm 8.145570438537e-09 ||r(i)||/||b|| 1.841422360715e-02
>>>> > > > > 369 KSP preconditioned resid norm 3.241361225044e-08 true resid
>>>> norm 1.107380056199e-08 ||r(i)||/||b|| 2.503390539290e-02
>>>> > > > > 370 KSP preconditioned resid norm 3.226640160009e-08 true resid
>>>> norm 8.914996078156e-09 ||r(i)||/||b|| 2.015361999245e-02
>>>> > > > > 371 KSP preconditioned resid norm 3.219199485113e-08 true resid
>>>> norm 9.028612338383e-09 ||r(i)||/||b|| 2.041046575138e-02
>>>> > > > > 372 KSP preconditioned resid norm 3.215039853884e-08 true resid
>>>> norm 9.957541741887e-09 ||r(i)||/||b|| 2.251044314161e-02
>>>> > > > > 373 KSP preconditioned resid norm 3.191438320482e-08 true resid
>>>> norm 8.943736857091e-09 ||r(i)||/||b|| 2.021859262192e-02
>>>> > > > > 374 KSP preconditioned resid norm 3.186200739666e-08 true resid
>>>> norm 8.758295038024e-09 ||r(i)||/||b|| 1.979937494427e-02
>>>> > > > > 375 KSP preconditioned resid norm 3.183870059099e-08 true resid
>>>> norm 8.943256269651e-09 ||r(i)||/||b|| 2.021750618548e-02
>>>> > > > > 376 KSP preconditioned resid norm 3.179122668815e-08 true resid
>>>> norm 8.590862092304e-09 ||r(i)||/||b|| 1.942086889304e-02
>>>> > > > > 377 KSP preconditioned resid norm 3.167956879702e-08 true resid
>>>> norm 9.582658100383e-09 ||r(i)||/||b|| 2.166296520824e-02
>>>> > > > > 378 KSP preconditioned resid norm 3.136955166729e-08 true resid
>>>> norm 9.977269021243e-09 ||r(i)||/||b|| 2.255503946987e-02
>>>> > > > > 379 KSP preconditioned resid norm 3.130691641774e-08 true resid
>>>> norm 8.509125545196e-09 ||r(i)||/||b|| 1.923609177194e-02
>>>> > > > > 380 KSP preconditioned resid norm 3.117879107848e-08 true resid
>>>> norm 9.053973136302e-09 ||r(i)||/||b|| 2.046779745176e-02
>>>> > > > > 381 KSP preconditioned resid norm 3.088492490292e-08 true resid
>>>> norm 1.048308325237e-08 ||r(i)||/||b|| 2.369850467295e-02
>>>> > > > > 382 KSP preconditioned resid norm 3.050052142818e-08 true resid
>>>> norm 9.220813554635e-09 ||r(i)||/||b|| 2.084496401033e-02
>>>> > > > > 383 KSP preconditioned resid norm 2.985444875009e-08 true resid
>>>> norm 1.234239992416e-08 ||r(i)||/||b|| 2.790175516463e-02
>>>> > > > > 384 KSP preconditioned resid norm 2.952644478108e-08 true resid
>>>> norm 8.796660599472e-09 ||r(i)||/||b|| 1.988610576719e-02
>>>> > > > > 385 KSP preconditioned resid norm 2.923012320168e-08 true resid
>>>> norm 9.047640736274e-09 ||r(i)||/||b|| 2.045348215843e-02
>>>> > > > > 386 KSP preconditioned resid norm 2.885571723790e-08 true resid
>>>> norm 1.000194856842e-08 ||r(i)||/||b|| 2.261083110579e-02
>>>> > > > > 387 KSP preconditioned resid norm 2.854406928323e-08 true resid
>>>> norm 7.861779671019e-09 ||r(i)||/||b|| 1.777267410608e-02
>>>> > > > > 388 KSP preconditioned resid norm 2.797557513811e-08 true resid
>>>> norm 1.546960768219e-08 ||r(i)||/||b|| 3.497125426932e-02
>>>> > > > > 389 KSP preconditioned resid norm 2.745355913117e-08 true resid
>>>> norm 7.687966698265e-09 ||r(i)||/||b|| 1.737974509390e-02
>>>> > > > > 390 KSP preconditioned resid norm 2.661217744985e-08 true resid
>>>> norm 1.052014137922e-08 ||r(i)||/||b|| 2.378227985352e-02
>>>> > > > > 391 KSP preconditioned resid norm 2.658120474029e-08 true resid
>>>> norm 9.465616738864e-09 ||r(i)||/||b|| 2.139837651939e-02
>>>> > > > > 392 KSP preconditioned resid norm 2.623787329422e-08 true resid
>>>> norm 1.909951824281e-08 ||r(i)||/||b|| 4.317718474915e-02
>>>> > > > > 393 KSP preconditioned resid norm 2.543325628802e-08 true resid
>>>> norm 5.948456163537e-09 ||r(i)||/||b|| 1.344733346046e-02
>>>> > > > > 394 KSP preconditioned resid norm 2.517323260201e-08 true resid
>>>> norm 1.034609461157e-08 ||r(i)||/||b|| 2.338882231462e-02
>>>> > > > > 395 KSP preconditioned resid norm 2.490676685335e-08 true resid
>>>> norm 8.641670445923e-09 ||r(i)||/||b|| 1.953572842212e-02
>>>> > > > > 396 KSP preconditioned resid norm 2.439513980463e-08 true resid
>>>> norm 7.519572419040e-09 ||r(i)||/||b|| 1.699906581119e-02
>>>> > > > > 397 KSP preconditioned resid norm 2.414249722638e-08 true resid
>>>> norm 1.073069104260e-08 ||r(i)||/||b|| 2.425825739384e-02
>>>> > > > > 398 KSP preconditioned resid norm 2.392817846978e-08 true resid
>>>> norm 7.116187521104e-09 ||r(i)||/||b|| 1.608715672313e-02
>>>> > > > > 399 KSP preconditioned resid norm 2.371648249227e-08 true resid
>>>> norm 1.098047911448e-08 ||r(i)||/||b|| 2.482293895233e-02
>>>> > > > > 400 KSP preconditioned resid norm 2.345879128128e-08 true resid
>>>> norm 8.016833131810e-09 ||r(i)||/||b|| 1.812319456621e-02
>>>> > > > > 401 KSP preconditioned resid norm 2.331728390350e-08 true resid
>>>> norm 7.808056352388e-09 ||r(i)||/||b|| 1.765122488290e-02
>>>> > > > > 402 KSP preconditioned resid norm 2.325356933312e-08 true resid
>>>> norm 9.927418776284e-09 ||r(i)||/||b|| 2.244234588206e-02
>>>> > > > > 403 KSP preconditioned resid norm 2.314519860272e-08 true resid
>>>> norm 7.444871113908e-09 ||r(i)||/||b|| 1.683019285787e-02
>>>> > > > > 404 KSP preconditioned resid norm 2.308308693430e-08 true resid
>>>> norm 7.801054524067e-09 ||r(i)||/||b|| 1.763539625145e-02
>>>> > > > > 405 KSP preconditioned resid norm 2.306358661113e-08 true resid
>>>> norm 7.805058505977e-09 ||r(i)||/||b|| 1.764444782357e-02
>>>> > > > > 406 KSP preconditioned resid norm 2.305322796876e-08 true resid
>>>> norm 7.571079190824e-09 ||r(i)||/||b|| 1.711550421413e-02
>>>> > > > > 407 KSP preconditioned resid norm 2.298488633980e-08 true resid
>>>> norm 8.830501095637e-09 ||r(i)||/||b|| 1.996260703473e-02
>>>> > > > > 408 KSP preconditioned resid norm 2.287294808434e-08 true resid
>>>> norm 7.013326899093e-09 ||r(i)||/||b|| 1.585462561823e-02
>>>> > > > > 409 KSP preconditioned resid norm 2.286549605315e-08 true resid
>>>> norm 7.082627748739e-09 ||r(i)||/||b|| 1.601129007177e-02
>>>> > > > > 410 KSP preconditioned resid norm 2.286053015844e-08 true resid
>>>> norm 7.033545134087e-09 ||r(i)||/||b|| 1.590033182173e-02
>>>> > > > > 411 KSP preconditioned resid norm 2.280043141482e-08 true resid
>>>> norm 7.319703904833e-09 ||r(i)||/||b|| 1.654723453179e-02
>>>> > > > > 412 KSP preconditioned resid norm 2.278834021186e-08 true resid
>>>> norm 7.137620865339e-09 ||r(i)||/||b|| 1.613560985436e-02
>>>> > > > > 413 KSP preconditioned resid norm 2.274834174042e-08 true resid
>>>> norm 7.536215034436e-09 ||r(i)||/||b|| 1.703668881668e-02
>>>> > > > > 414 KSP preconditioned resid norm 2.271641465026e-08 true resid
>>>> norm 7.111204145035e-09 ||r(i)||/||b|| 1.607589109085e-02
>>>> > > > > 415 KSP preconditioned resid norm 2.271597283066e-08 true resid
>>>> norm 7.105850673247e-09 ||r(i)||/||b|| 1.606378880442e-02
>>>> > > > > 416 KSP preconditioned resid norm 2.264562496921e-08 true resid
>>>> norm 7.590206017147e-09 ||r(i)||/||b|| 1.715874313269e-02
>>>> > > > > 417 KSP preconditioned resid norm 2.264130303752e-08 true resid
>>>> norm 7.322331610853e-09 ||r(i)||/||b|| 1.655317483599e-02
>>>> > > > > 418 KSP preconditioned resid norm 2.257869016612e-08 true resid
>>>> norm 8.689110661905e-09 ||r(i)||/||b|| 1.964297379574e-02
>>>> > > > > 419 KSP preconditioned resid norm 2.252448770837e-08 true resid
>>>> norm 7.589690206564e-09 ||r(i)||/||b|| 1.715757706931e-02
>>>> > > > > 420 KSP preconditioned resid norm 2.251854717911e-08 true resid
>>>> norm 7.606351482008e-09 ||r(i)||/||b|| 1.719524225849e-02
>>>> > > > > 421 KSP preconditioned resid norm 2.251852680166e-08 true resid
>>>> norm 7.587363564342e-09 ||r(i)||/||b|| 1.715231736276e-02
>>>> > > > > 422 KSP preconditioned resid norm 2.251740722171e-08 true resid
>>>> norm 8.042777274478e-09 ||r(i)||/||b|| 1.818184500057e-02
>>>> > > > > 423 KSP preconditioned resid norm 2.251066763399e-08 true resid
>>>> norm 7.161199860439e-09 ||r(i)||/||b|| 1.618891353535e-02
>>>> > > > > 424 KSP preconditioned resid norm 2.250468324723e-08 true resid
>>>> norm 7.892076688728e-09 ||r(i)||/||b|| 1.784116483524e-02
>>>> > > > > 425 KSP preconditioned resid norm 2.248157980181e-08 true resid
>>>> norm 7.381796123970e-09 ||r(i)||/||b|| 1.668760284806e-02
>>>> > > > > 426 KSP preconditioned resid norm 2.245901138193e-08 true resid
>>>> norm 7.147464780092e-09 ||r(i)||/||b|| 1.615786342748e-02
>>>> > > > > 427 KSP preconditioned resid norm 2.241049755791e-08 true resid
>>>> norm 8.351728267326e-09 ||r(i)||/||b|| 1.888027277907e-02
>>>> > > > > 428 KSP preconditioned resid norm 2.234407003838e-08 true resid
>>>> norm 6.648970113909e-09 ||r(i)||/||b|| 1.503094514480e-02
>>>> > > > > 429 KSP preconditioned resid norm 2.225463790450e-08 true resid
>>>> norm 9.144039069449e-09 ||r(i)||/||b|| 2.067140433784e-02
>>>> > > > > 430 KSP preconditioned resid norm 2.219515360063e-08 true resid
>>>> norm 7.964312082788e-09 ||r(i)||/||b|| 1.800446324492e-02
>>>> > > > > 431 KSP preconditioned resid norm 2.207622538143e-08 true resid
>>>> norm 7.744159741884e-09 ||r(i)||/||b|| 1.750677748263e-02
>>>> > > > > 432 KSP preconditioned resid norm 2.189830543256e-08 true resid
>>>> norm 1.088914113364e-08 ||r(i)||/||b|| 2.461645642103e-02
>>>> > > > > 433 KSP preconditioned resid norm 2.136231265645e-08 true resid
>>>> norm 6.923381943847e-09 ||r(i)||/||b|| 1.565129221995e-02
>>>> > > > > 434 KSP preconditioned resid norm 2.037470395914e-08 true resid
>>>> norm 1.030114953089e-08 ||r(i)||/||b|| 2.328721755017e-02
>>>> > > > > 435 KSP preconditioned resid norm 1.961890283205e-08 true resid
>>>> norm 7.295267301635e-09 ||r(i)||/||b|| 1.649199210539e-02
>>>> > > > > 436 KSP preconditioned resid norm 1.913560642295e-08 true resid
>>>> norm 6.039280397468e-09 ||r(i)||/||b|| 1.365265459361e-02
>>>> > > > > 437 KSP preconditioned resid norm 1.877609800000e-08 true resid
>>>> norm 8.984892222755e-09 ||r(i)||/||b|| 2.031163019513e-02
>>>> > > > > 438 KSP preconditioned resid norm 1.811081731215e-08 true resid
>>>> norm 6.564934723165e-09 ||r(i)||/||b|| 1.484097115983e-02
>>>> > > > > 439 KSP preconditioned resid norm 1.785755067179e-08 true resid
>>>> norm 5.336927073942e-09 ||r(i)||/||b|| 1.206488474394e-02
>>>> > > > > 440 KSP preconditioned resid norm 1.772971743003e-08 true resid
>>>> norm 4.970154886573e-09 ||r(i)||/||b|| 1.123574390942e-02
>>>> > > > > 441 KSP preconditioned resid norm 1.735301571457e-08 true resid
>>>> norm 5.652698869454e-09 ||r(i)||/||b|| 1.277873191957e-02
>>>> > > > > 442 KSP preconditioned resid norm 1.631139023345e-08 true resid
>>>> norm 7.160283137945e-09 ||r(i)||/||b|| 1.618684115342e-02
>>>> > > > > 443 KSP preconditioned resid norm 1.483956291991e-08 true resid
>>>> norm 7.988466840254e-09 ||r(i)||/||b|| 1.805906851886e-02
>>>> > > > > 444 KSP preconditioned resid norm 1.410013314968e-08 true resid
>>>> norm 7.482143551338e-09 ||r(i)||/||b|| 1.691445251806e-02
>>>> > > > > 445 KSP preconditioned resid norm 1.356330060794e-08 true resid
>>>> norm 3.083065574605e-09 ||r(i)||/||b|| 6.969709403983e-03
>>>> > > > > 446 KSP preconditioned resid norm 1.314115350470e-08 true resid
>>>> norm 4.164976230172e-09 ||r(i)||/||b|| 9.415522730979e-03
>>>> > > > > 447 KSP preconditioned resid norm 1.236636402318e-08 true resid
>>>> norm 5.761202097878e-09 ||r(i)||/||b|| 1.302401894094e-02
>>>> > > > > 448 KSP preconditioned resid norm 1.180595146541e-08 true resid
>>>> norm 4.038199940007e-09 ||r(i)||/||b|| 9.128926847635e-03
>>>> > > > > 449 KSP preconditioned resid norm 1.131424900987e-08 true resid
>>>> norm 5.090038236142e-09 ||r(i)||/||b|| 1.150675731755e-02
>>>> > > > > 450 KSP preconditioned resid norm 1.063388583397e-08 true resid
>>>> norm 3.010230442125e-09 ||r(i)||/||b|| 6.805055200074e-03
>>>> > > > > 451 KSP preconditioned resid norm 1.051458985758e-08 true resid
>>>> norm 1.649000916776e-09 ||r(i)||/||b|| 3.727801734579e-03
>>>> > > > > 452 KSP preconditioned resid norm 1.012147631812e-08 true resid
>>>> norm 5.884734670038e-09 ||r(i)||/||b|| 1.330328193716e-02
>>>> > > > > 453 KSP preconditioned resid norm 9.390240831193e-09 true resid
>>>> norm 2.031511507263e-09 ||r(i)||/||b|| 4.592521473788e-03
>>>> > > > > 454 KSP preconditioned resid norm 9.098339938653e-09 true resid
>>>> norm 1.769428017679e-09 ||r(i)||/||b|| 4.000044370146e-03
>>>> > > > > 455 KSP preconditioned resid norm 8.890748233216e-09 true resid
>>>> norm 2.050221736547e-09 ||r(i)||/||b|| 4.634818615330e-03
>>>> > > > > 456 KSP preconditioned resid norm 8.527203662681e-09 true resid
>>>> norm 2.051638846865e-09 ||r(i)||/||b|| 4.638022195295e-03
>>>> > > > > 457 KSP preconditioned resid norm 8.268992410262e-09 true resid
>>>> norm 2.382052265392e-09 ||r(i)||/||b|| 5.384968847771e-03
>>>> > > > > 458 KSP preconditioned resid norm 7.908973261474e-09 true resid
>>>> norm 3.438027649085e-09 ||r(i)||/||b|| 7.772151794095e-03
>>>> > > > > 459 KSP preconditioned resid norm 7.427882465376e-09 true resid
>>>> norm 2.945117076322e-09 ||r(i)||/||b|| 6.657857150931e-03
>>>> > > > > 460 KSP preconditioned resid norm 7.015919500900e-09 true resid
>>>> norm 1.769207642661e-09 ||r(i)||/||b|| 3.999546180990e-03
>>>> > > > > 461 KSP preconditioned resid norm 6.803968849346e-09 true resid
>>>> norm 1.305695673120e-09 ||r(i)||/||b|| 2.951711272911e-03
>>>> > > > > 462 KSP preconditioned resid norm 6.710918799980e-09 true resid
>>>> norm 1.795032965660e-09 ||r(i)||/||b|| 4.057928006550e-03
>>>> > > > > 463 KSP preconditioned resid norm 6.577483627377e-09 true resid
>>>> norm 1.782001044516e-09 ||r(i)||/||b|| 4.028467490336e-03
>>>> > > > > 464 KSP preconditioned resid norm 6.417912605066e-09 true resid
>>>> norm 1.290960823243e-09 ||r(i)||/||b|| 2.918401043442e-03
>>>> > > > > 465 KSP preconditioned resid norm 6.365751616365e-09 true resid
>>>> norm 1.263290494523e-09 ||r(i)||/||b|| 2.855848319336e-03
>>>> > > > > 466 KSP preconditioned resid norm 6.293400160765e-09 true resid
>>>> norm 1.425147805149e-09 ||r(i)||/||b|| 3.221749852300e-03
>>>> > > > > 467 KSP preconditioned resid norm 5.891998607109e-09 true resid
>>>> norm 2.318264750150e-09 ||r(i)||/||b|| 5.240768072899e-03
>>>> > > > > 468 KSP preconditioned resid norm 5.803460471252e-09 true resid
>>>> norm 1.081820742140e-09 ||r(i)||/||b|| 2.445610064874e-03
>>>> > > > > 469 KSP preconditioned resid norm 5.754085062170e-09 true resid
>>>> norm 1.656644764768e-09 ||r(i)||/||b|| 3.745081743044e-03
>>>> > > > > 470 KSP preconditioned resid norm 5.525556632876e-09 true resid
>>>> norm 1.979760762585e-09 ||r(i)||/||b|| 4.475531535326e-03
>>>> > > > > 471 KSP preconditioned resid norm 5.443965422137e-09 true resid
>>>> norm 9.955797379882e-10 ||r(i)||/||b|| 2.250649976253e-03
>>>> > > > > 472 KSP preconditioned resid norm 5.414797083501e-09 true resid
>>>> norm 1.143450718543e-09 ||r(i)||/||b|| 2.584933415516e-03
>>>> > > > > 473 KSP preconditioned resid norm 5.408283291827e-09 true resid
>>>> norm 1.029764590400e-09 ||r(i)||/||b|| 2.327929710194e-03
>>>> > > > > 474 KSP preconditioned resid norm 5.355045752288e-09 true resid
>>>> norm 1.886388210190e-09 ||r(i)||/||b|| 4.264449564882e-03
>>>> > > > > 475 KSP preconditioned resid norm 5.241533530016e-09 true resid
>>>> norm 1.425455398363e-09 ||r(i)||/||b|| 3.222445210625e-03
>>>> > > > > 476 KSP preconditioned resid norm 5.125189682735e-09 true resid
>>>> norm 1.354368288121e-09 ||r(i)||/||b|| 3.061742660270e-03
>>>> > > > > 477 KSP preconditioned resid norm 5.057977457429e-09 true resid
>>>> norm 1.740562818023e-09 ||r(i)||/||b|| 3.934790469888e-03
>>>> > > > > 478 KSP preconditioned resid norm 4.979286088604e-09 true resid
>>>> norm 1.417735808423e-09 ||r(i)||/||b|| 3.204993976684e-03
>>>> > > > > 479 KSP preconditioned resid norm 4.894533209216e-09 true resid
>>>> norm 1.326984715102e-09 ||r(i)||/||b|| 2.999838188318e-03
>>>> > > > > 480 KSP preconditioned resid norm 4.777786489681e-09 true resid
>>>> norm 1.456212644054e-09 ||r(i)||/||b|| 3.291976350766e-03
>>>> > > > > 481 KSP preconditioned resid norm 4.775308813328e-09 true resid
>>>> norm 1.364283691096e-09 ||r(i)||/||b|| 3.084157842718e-03
>>>> > > > > 482 KSP preconditioned resid norm 4.760397485605e-09 true resid
>>>> norm 2.049838358389e-09 ||r(i)||/||b|| 4.633951934330e-03
>>>> > > > > 483 KSP preconditioned resid norm 4.670996943034e-09 true resid
>>>> norm 8.889456596556e-10 ||r(i)||/||b|| 2.009588435214e-03
>>>> > > > > 484 KSP preconditioned resid norm 4.597559289814e-09 true resid
>>>> norm 1.386989122943e-09 ||r(i)||/||b|| 3.135486709406e-03
>>>> > > > > 485 KSP preconditioned resid norm 4.570207496163e-09 true resid
>>>> norm 1.099843879320e-09 ||r(i)||/||b|| 2.486353936728e-03
>>>> > > > > 486 KSP preconditioned resid norm 4.542158022055e-09 true resid
>>>> norm 1.011236566485e-09 ||r(i)||/||b|| 2.286044469873e-03
>>>> > > > > 487 KSP preconditioned resid norm 4.526766510181e-09 true resid
>>>> norm 1.399400517345e-09 ||r(i)||/||b|| 3.163544436426e-03
>>>> > > > > 488 KSP preconditioned resid norm 4.505705850329e-09 true resid
>>>> norm 9.117728324580e-10 ||r(i)||/||b|| 2.061192514692e-03
>>>> > > > > 489 KSP preconditioned resid norm 4.487150083935e-09 true resid
>>>> norm 1.381092775475e-09 ||r(i)||/||b|| 3.122157175082e-03
>>>> > > > > 490 KSP preconditioned resid norm 4.473919835893e-09 true resid
>>>> norm 9.652867482479e-10 ||r(i)||/||b|| 2.182168352895e-03
>>>> > > > > 491 KSP preconditioned resid norm 4.452517186816e-09 true resid
>>>> norm 1.131636320873e-09 ||r(i)||/||b|| 2.558225284743e-03
>>>> > > > > 492 KSP preconditioned resid norm 4.400952072877e-09 true resid
>>>> norm 1.786271233792e-09 ||r(i)||/||b|| 4.038120862161e-03
>>>> > > > > 493 KSP preconditioned resid norm 4.366719039095e-09 true resid
>>>> norm 1.028221901066e-09 ||r(i)||/||b|| 2.324442240954e-03
>>>> > > > > 494 KSP preconditioned resid norm 4.331818466820e-09 true resid
>>>> norm 1.573868226869e-09 ||r(i)||/||b|| 3.557953574452e-03
>>>> > > > > 495 KSP preconditioned resid norm 4.302438670088e-09 true resid
>>>> norm 1.043943800033e-09 ||r(i)||/||b|| 2.359983835651e-03
>>>> > > > > 496 KSP preconditioned resid norm 4.254710406786e-09 true resid
>>>> norm 1.390625661617e-09 ||r(i)||/||b|| 3.143707623679e-03
>>>> > > > > 497 KSP preconditioned resid norm 4.221208729086e-09 true resid
>>>> norm 1.442390433810e-09 ||r(i)||/||b|| 3.260729273340e-03
>>>> > > > > 498 KSP preconditioned resid norm 4.194726771303e-09 true resid
>>>> norm 1.122208863419e-09 ||r(i)||/||b|| 2.536913172729e-03
>>>> > > > > 499 KSP preconditioned resid norm 4.161552261100e-09 true resid
>>>> norm 1.685040004550e-09 ||r(i)||/||b|| 3.809273231987e-03
>>>> > > > > 500 KSP preconditioned resid norm 4.157891682720e-09 true resid
>>>> norm 1.444023641507e-09 ||r(i)||/||b|| 3.264421372248e-03
>>>> > > > >
>>>> > > > >
>>>> > > > >
>>>> > > > >
>>>> > > > > 2018-08-06 15:40 GMT-05:00 Smith, Barry F. <bsmith at mcs.anl.gov
>>>> >:
>>>> > > > >
>>>> > > > >     Ok, I think you want to use PCFIELDSPLIT with two fields u
>>>> and c. Start by using PC_COMPOSITE_MULTIPLICATIVE with direct solvers (LU)
>>>> for each field. Then try replacing the direct solvers LU with hypre
>>>> BoomerAMG first one at a time.
>>>> > > > >
>>>> > > > >     Let us know how it goes and we'll make additional
>>>> suggestions,
>>>> > > > >
>>>> > > > >     Barry
>>>> > > > >
>>>> > > > >
>>>> > > > > > On Aug 6, 2018, at 12:58 PM, Josh L <ysjosh.lo at gmail.com>
>>>> wrote:
>>>> > > > > >
>>>> > > > > > Hi Barry,
>>>> > > > > >
>>>> > > > > > (1)
>>>> > > > > > I am basically solving the equilibrium coupled with phase
>>>> field equation:
>>>> > > > > >
>>>> > > > > > 1-D example
>>>> > > > > > u is displacement and c is phase-field parameter.
>>>> > > > > >
>>>> > > > > > equilibrium: (c^2)*E* (\Delta u) = 0
>>>> > > > > >
>>>> > > > > > phase-field equation:  2*G_c*L*(\Delta
>>>> c)-E*(\epsilon)^2*c+G_c/2L(1-c)=0
>>>> > > > > >
>>>> > > > > > I solve the above 2 equations coupled.
>>>> > > > > > it is the formulation for quasi-static fracture problem,
>>>> where G_c is the critical energy release rate and L is the length scale in
>>>> phase-field method (both constant).
>>>> > > > > >
>>>> > > > > > The free energy I am using now is
>>>> > > > > > 1-D:
>>>> > > > > >
>>>> > > > > > \psi=c^2*0.5*E*\epsilon^2+G_c/4L*(1-c)^2+G_c*L*grad c * grad
>>>> c
>>>> > > > > >
>>>> > > > > >
>>>> > > > > > it is not convex w.r.t u and c simultaneously.
>>>> > > > > >
>>>> > > > > > (2)
>>>> > > > > >
>>>> > > > > > with an appropriate step size. Newton's method usually
>>>> converges with 5 iteration, and I can see the quadratic convergence of
>>>> Newton's method.
>>>> > > > > >
>>>> > > > > >
>>>> > > > > >
>>>> > > > > >
>>>> > > > > > Thanks,
>>>> > > > > > Josh
>>>> > > > > >
>>>> > > > > >
>>>> > > > > >
>>>> > > > > > 2018-08-05 22:23 GMT-05:00 Smith, Barry F. <
>>>> bsmith at mcs.anl.gov>:
>>>> > > > > >
>>>> > > > > >
>>>> > > > > > > On Aug 5, 2018, at 9:38 PM, Josh L <ysjosh.lo at gmail.com>
>>>> wrote:
>>>> > > > > > >
>>>> > > > > > > Hi,
>>>> > > > > > >
>>>> > > > > > > I am using GRMES to solve the linear systems within Newton
>>>> Raphson.
>>>> > > > > > > I am using finite element method to solve fracture problem
>>>> with phase field method.(the formulation gives symmetric but indefinite
>>>> hessian so can't use CG)
>>>> > > > > >
>>>> > > > > >     What exact equations are you using for the phase field
>>>> method?
>>>> > > > > >
>>>> > > > > >      With the direct LU solver roughly how many iterations
>>>> does your Newton's method take?
>>>> > > > > >
>>>> > > > > >
>>>> > > > > >    Barry
>>>> > > > > >
>>>> > > > > > >
>>>> > > > > > > The total iteration number to solve one linear system is
>>>> often very large( > 2000) even for a small sized problem(120K degree of
>>>> freedom). I've tried different preconditioner, but it doesn't help a lot.
>>>> I have checked my hessian and right hand side, and sure they are correct.
>>>> > > > > > >
>>>> > > > > > > I think it might stem from the great difference in element
>>>> size(max/min=100) and damage stiffness around crack.
>>>> > > > > > >
>>>> > > > > > > Anyone has idea to reduce the number of iteration?
>>>> > > > > > >
>>>> > > > > > > Now, I get away with it by using LU, which is not ideal for
>>>> large size problem (I will have to run it in the future)
>>>> > > > > > >
>>>> > > > > > >
>>>> > > > > > > Thanks,
>>>> > > > > > > Josh
>>>> > > > > > >
>>>> > > > > >
>>>> > > > > >
>>>> > > > >
>>>> > > > >
>>>> > > >
>>>> > > >
>>>> >
>>>> >
>>>>
>>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which their
>>> experiments lead.
>>> -- Norbert Wiener
>>>
>>> https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
>>>
>>
>>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180815/a6e59140/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mat
Type: application/octet-stream
Size: 17555128 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180815/a6e59140/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rhs
Type: application/octet-stream
Size: 646424 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180815/a6e59140/attachment-0003.obj>


More information about the petsc-users mailing list