Looks fine to me. Thanks. Barry Smith <[email protected]> writes:
> A proposed fix > https://bitbucket.org/petsc/petsc/pull-requests/645/do-not-assume-that-all-ksp-methods-support > > Needs Jed's approval. > > Barry > > > >> On Mar 8, 2017, at 10:33 AM, Barry Smith <[email protected]> wrote: >> >> >> Please tell us how you got this output. >> >> PETSc CG doesn't even implement right preconditioning. If you ask for it it >> should error out. CG supports no norm computation with left preconditioning. >> >> Barry >> >>> On Mar 8, 2017, at 10:26 AM, Kong, Fande <[email protected]> wrote: >>> >>> Hi All, >>> >>> The NONE norm type is supported only when CG is used with a right >>> preconditioner. Any reason for this? >>> >>> >>> >>> 0 Nonlinear |R| = 1.732051e+00 >>> 0 Linear |R| = 0.000000e+00 >>> 1 Linear |R| = 0.000000e+00 >>> 2 Linear |R| = 0.000000e+00 >>> 3 Linear |R| = 0.000000e+00 >>> 4 Linear |R| = 0.000000e+00 >>> 5 Linear |R| = 0.000000e+00 >>> 6 Linear |R| = 0.000000e+00 >>> 1 Nonlinear |R| = 1.769225e-08 >>> 0 Linear |R| = 0.000000e+00 >>> 1 Linear |R| = 0.000000e+00 >>> 2 Linear |R| = 0.000000e+00 >>> 3 Linear |R| = 0.000000e+00 >>> 4 Linear |R| = 0.000000e+00 >>> 5 Linear |R| = 0.000000e+00 >>> 6 Linear |R| = 0.000000e+00 >>> 7 Linear |R| = 0.000000e+00 >>> 8 Linear |R| = 0.000000e+00 >>> 9 Linear |R| = 0.000000e+00 >>> 10 Linear |R| = 0.000000e+00 >>> 2 Nonlinear |R| = 0.000000e+00 >>> SNES Object: 1 MPI processes >>> type: newtonls >>> maximum iterations=50, maximum function evaluations=10000 >>> tolerances: relative=1e-08, absolute=1e-50, solution=1e-50 >>> total number of linear solver iterations=18 >>> total number of function evaluations=23 >>> norm schedule ALWAYS >>> SNESLineSearch Object: 1 MPI processes >>> type: bt >>> interpolation: cubic >>> alpha=1.000000e-04 >>> maxstep=1.000000e+08, minlambda=1.000000e-12 >>> tolerances: relative=1.000000e-08, absolute=1.000000e-15, >>> lambda=1.000000e-08 >>> maximum iterations=40 >>> KSP Object: 1 MPI processes >>> type: cg >>> maximum iterations=10000, initial guess is zero >>> tolerances: relative=1e-05, absolute=1e-50, divergence=10000. >>> right preconditioning >>> using NONE norm type for convergence test >>> PC Object: 1 MPI processes >>> type: hypre >>> HYPRE BoomerAMG preconditioning >>> HYPRE BoomerAMG: Cycle type V >>> HYPRE BoomerAMG: Maximum number of levels 25 >>> HYPRE BoomerAMG: Maximum number of iterations PER hypre call 1 >>> HYPRE BoomerAMG: Convergence tolerance PER hypre call 0. >>> HYPRE BoomerAMG: Threshold for strong coupling 0.25 >>> HYPRE BoomerAMG: Interpolation truncation factor 0. >>> HYPRE BoomerAMG: Interpolation: max elements per row 0 >>> HYPRE BoomerAMG: Number of levels of aggressive coarsening 0 >>> HYPRE BoomerAMG: Number of paths for aggressive coarsening 1 >>> HYPRE BoomerAMG: Maximum row sums 0.9 >>> HYPRE BoomerAMG: Sweeps down 1 >>> HYPRE BoomerAMG: Sweeps up 1 >>> HYPRE BoomerAMG: Sweeps on coarse 1 >>> HYPRE BoomerAMG: Relax down symmetric-SOR/Jacobi >>> HYPRE BoomerAMG: Relax up symmetric-SOR/Jacobi >>> HYPRE BoomerAMG: Relax on coarse Gaussian-elimination >>> HYPRE BoomerAMG: Relax weight (all) 1. >>> HYPRE BoomerAMG: Outer relax weight (all) 1. >>> HYPRE BoomerAMG: Using CF-relaxation >>> HYPRE BoomerAMG: Not using more complex smoothers. >>> HYPRE BoomerAMG: Measure type local >>> HYPRE BoomerAMG: Coarsen type Falgout >>> HYPRE BoomerAMG: Interpolation type classical >>> linear system matrix followed by preconditioner matrix: >>> Mat Object: 1 MPI processes >>> type: mffd >>> rows=9, cols=9 >>> Matrix-free approximation: >>> err=1.49012e-08 (relative error in function evaluation) >>> Using wp compute h routine >>> Does not compute normU >>> Mat Object: () 1 MPI processes >>> type: seqaij >>> rows=9, cols=9 >>> total: nonzeros=49, allocated nonzeros=49 >>> total number of mallocs used during MatSetValues calls =0 >>> not using I-node routines >>> >>> Fande, >>> >>
signature.asc
Description: PGP signature
