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