| /petsc/src/ksp/ksp/impls/cheby/ |
| H A D | cheby.c | 18 …rorCode KSPChebyshevComputeExtremeEigenvalues_Private(KSP kspest, PetscReal *emin, PetscReal *emax) in KSPChebyshevComputeExtremeEigenvalues_Private() argument 34 *emax = max; in KSPChebyshevComputeExtremeEigenvalues_Private() 40 static PetscErrorCode KSPChebyshevGetEigenvalues_Chebyshev(KSP ksp, PetscReal *emax, PetscReal *emi… in KSPChebyshevGetEigenvalues_Chebyshev() argument 45 *emax = 0; in KSPChebyshevGetEigenvalues_Chebyshev() 47 if (cheb->emax != 0.) { in KSPChebyshevGetEigenvalues_Chebyshev() 48 *emax = cheb->emax; in KSPChebyshevGetEigenvalues_Chebyshev() 50 *emax = cheb->tform[2] * cheb->emin_computed + cheb->tform[3] * cheb->emax_computed; in KSPChebyshevGetEigenvalues_Chebyshev() 52 *emax = cheb->tform[2] * cheb->emin_provided + cheb->tform[3] * cheb->emax_provided; in KSPChebyshevGetEigenvalues_Chebyshev() 64 static PetscErrorCode KSPChebyshevSetEigenvalues_Chebyshev(KSP ksp, PetscReal emax, PetscReal emin) in KSPChebyshevSetEigenvalues_Chebyshev() argument 69 …emax > emin || (emax == 0 && emin == 0) || (emax == -1 && emin == -1), PetscObjectComm((PetscObjec… in KSPChebyshevSetEigenvalues_Chebyshev() [all …]
|
| H A D | chebyshevimpl.h | 10 …PetscReal emin, emax; /* store user provided estimates of extreme eigenvalues o… member
|
| /petsc/src/snes/interface/noise/ |
| H A D | snesdnest.c | 22 static double emin, emax; in SNESNoise_dnest_() local 216 emax = PetscMax(d__1, eps[6]); in SNESNoise_dnest_() 218 if (emax <= emin * 4 && dnoise) { in SNESNoise_dnest_() 236 emax = PetscMax(d__1, eps[5]); in SNESNoise_dnest_() 238 if (emax <= emin * 4 && dnoise) { in SNESNoise_dnest_()
|
| /petsc/src/ksp/ksp/interface/saws/ |
| H A D | kspsaws.c | 81 PetscReal emax, emin; in KSPMonitorSAWs() local 86 PetscCall(KSPComputeExtremeSingularValues(ksp, &emax, &emin)); in KSPMonitorSAWs() 104 …tscCall(PetscInfo(ksp, "KSP extreme singular values min=%g max=%g\n", (double)emin, (double)emax)); in KSPMonitorSAWs()
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cgeig.c | 42 PetscErrorCode KSPComputeExtremeSingularValues_CG(KSP ksp, PetscReal *emax, PetscReal *emin) in KSPComputeExtremeSingularValues_CG() argument 52 *emax = *emin = 1.0; in KSPComputeExtremeSingularValues_CG() 72 *emax = dd[n - 1]; in KSPComputeExtremeSingularValues_CG()
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex56_seqaijmkl.out | 4 [0]main |b-Ax|/|b|=9.264005e-05, |b|=5.391826e+00, emax=9.987464e-01
|
| H A D | ex56_nns.out | 4 [0]main |b-Ax|/|b|=4.543854e-05, |b|=5.391826e+00, emax=9.949156e-01
|
| H A D | ex56_2.out | 4 [0]main |b-Ax|/|b|=9.023568e-05, |b|=4.969822e+00, emax=1.016793e+00
|
| H A D | ex56_hypre.out | 4 [0]main |b-Ax|/|b|=9.722297e-05, |b|=4.351790e+00, emax=9.991454e-01
|
| H A D | ex56_gamg.out | 4 [0]main |b-Ax|/|b|=4.993291e-05, |b|=4.351790e+00, emax=9.986640e-01
|
| H A D | ex56_1.out | 4 [0]main |b-Ax|/|b|=9.081303e-05, |b|=4.969822e+00, emax=9.982432e-01
|
| H A D | ex56_latebs-2.out | 285 [0]main |b-Ax|/|b|=2.425235e-04, |b|=5.391826e+00, emax=9.946388e-01
|
| H A D | ex56_latebs.out | 496 [0]main |b-Ax|/|b|=1.165490e-04, |b|=5.391826e+00, emax=9.990809e-01
|
| /petsc/lib/petsc/bin/maint/ |
| H A D | toclapack.sh | 623 static real emin, prec, emax; 1652 /* = 'L' or 'l', DLAMCH := emax */ 1665 /* emax = largest exponent before overflow */ 1666 /* rmax = overflow threshold - (base**emax)*(1-eps) */ 1699 emax = (doublereal) imax; 1728 rmach = emax; 1952 doublereal *eps, integer *emin, doublereal *rmin, integer *emax, 2252 *emax = lemax; 2405 logical *ieee, integer *emax, doublereal *rmax) 2508 *emax = expsum + *emin - 1; [all …]
|
| /petsc/src/ksp/pc/impls/gamg/ |
| H A D | gamg.c | 40 pc_gamg->emax = 0; in PCReset_GAMG() 931 PetscReal emax, emin; in PCSetUp_GAMG() local 934 emax = mg->max_eigen_DinvA[level]; in PCSetUp_GAMG() 935 …emin = %g\n", ((PetscObject)pc)->prefix, level, Aarr[level]->rmap->N, (double)emax, (double)emin)); in PCSetUp_GAMG() 937 cheb->emax_provided = emax; in PCSetUp_GAMG() 1222 PetscErrorCode PCGAMGSetEigenvalues(PC pc, PetscReal emax, PetscReal emin) in PCGAMGSetEigenvalues() argument 1226 PetscTryMethod(pc, "PCGAMGSetEigenvalues_C", (PC, PetscReal, PetscReal), (pc, emax, emin)); in PCGAMGSetEigenvalues() 1230 static PetscErrorCode PCGAMGSetEigenvalues_GAMG(PC pc, PetscReal emax, PetscReal emin) in PCGAMGSetEigenvalues_GAMG() argument 1236 …ck(emax > emin, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_INCOMP, "Maximum eigenvalue must b… in PCGAMGSetEigenvalues_GAMG() 1237 …ck(emax * emin > 0.0, PetscObjectComm((PetscObject)pc), PETSC_ERR_ARG_INCOMP, "Both eigenvalues mu… in PCGAMGSetEigenvalues_GAMG() [all …]
|
| H A D | agg.c | 1434 PetscReal alpha, emax, emin; in PCGAMGOptimizeProlongator_AGG() local 1443 if (pc_gamg->emax > 0) { in PCGAMGOptimizeProlongator_AGG() 1445 emax = pc_gamg->emax; in PCGAMGOptimizeProlongator_AGG() 1480 PetscCall(KSPComputeExtremeSingularValues(eksp, &emax, &emin)); in PCGAMGOptimizeProlongator_AGG() 1481 …ooth P0: max eigen=%e min=%e PC=%s\n", ((PetscObject)pc)->prefix, (double)emax, (double)emin, PCJA… in PCGAMGOptimizeProlongator_AGG() 1488 mg->max_eigen_DinvA[pc_gamg->current_level] = emax; in PCGAMGOptimizeProlongator_AGG() 1489 … spectra %g %g\n", ((PetscObject)pc)->prefix, pc_gamg->current_level, (double)emin, (double)emax)); in PCGAMGOptimizeProlongator_AGG() 1504 …PetscCheck(emax != 0.0, PetscObjectComm((PetscObject)pc), PETSC_ERR_PLIB, "Computed maximum singul… in PCGAMGOptimizeProlongator_AGG() 1526 alpha = -1.4 / emax; in PCGAMGOptimizeProlongator_AGG()
|
| /petsc/include/petsc/private/ |
| H A D | pcgamgimpl.h | 52 PetscReal emin, emax; member
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmreig.c | 4 PetscErrorCode KSPComputeExtremeSingularValues_GMRES(KSP ksp, PetscReal *emax, PetscReal *emin) in KSPComputeExtremeSingularValues_GMRES() argument 18 *emax = *emin = 1.0; in KSPComputeExtremeSingularValues_GMRES() 38 *emax = realpart[0]; in KSPComputeExtremeSingularValues_GMRES()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex56.c | 374 PetscReal emax, emin; in main() local 388 PetscCall(KSPComputeExtremeSingularValues(ksp, &emax, &emin)); in main() 410 … |b|=%e, emax=%e\n", 0, PETSC_FUNCTION_NAME, (double)(norm / norm2), (double)norm2, (double)emax)); in main()
|
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 661 PetscReal emax, emin; in SNESSolve_NEWTONTR() local 662 PetscCall(KSPComputeExtremeSingularValues(snes->ksp, &emax, &emin)); in SNESSolve_NEWTONTR() 663 if (emax > 0.0) beta_k = emax + 1; in SNESSolve_NEWTONTR()
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | iterativ.c | 1026 PetscReal emin, emax; in KSPMonitorSingularValue() local 1041 PetscCall(KSPComputeExtremeSingularValues(ksp, &emax, &emin)); in KSPMonitorSingularValue() 1042 …e min %14.12e max/min %14.12e\n", n, (double)rnorm, (double)emax, (double)emin, (double)(emax / em… in KSPMonitorSingularValue()
|
| H A D | itfunc.c | 55 PetscErrorCode KSPComputeExtremeSingularValues(KSP ksp, PetscReal *emax, PetscReal *emin) in KSPComputeExtremeSingularValues() argument 59 PetscAssertPointer(emax, 2); in KSPComputeExtremeSingularValues() 63 …s->computeextremesingularvalues) PetscUseTypeMethod(ksp, computeextremesingularvalues, emax, emin); in KSPComputeExtremeSingularValues() 66 *emax = -1.0; in KSPComputeExtremeSingularValues()
|
| /petsc/src/ksp/ksp/impls/minres/ |
| H A D | minres.c | 905 PetscErrorCode KSPComputeExtremeSingularValues_MINRES(KSP ksp, PetscReal *emax, PetscReal *emin) in KSPComputeExtremeSingularValues_MINRES() argument 915 *emax = *emin = 1.0; in KSPComputeExtremeSingularValues_MINRES() 937 *emax = dd[n - 1]; in KSPComputeExtremeSingularValues_MINRES()
|
| /petsc/src/ksp/ksp/impls/fetidp/ |
| H A D | fetidp.c | 248 static PetscErrorCode KSPComputeExtremeSingularValues_FETIDP(KSP ksp, PetscReal *emax, PetscReal *e… in KSPComputeExtremeSingularValues_FETIDP() argument 253 PetscCall(KSPComputeExtremeSingularValues(fetidp->innerksp, emax, emin)); in KSPComputeExtremeSingularValues_FETIDP()
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 190 KSPChebyshevSetEigenvalues(KSP ksp,PetscReal emax,PetscReal emin); 195 `scale=1.0, emax=0.01, emin=100.0`, and `max_steps=30`. The 1185 (smoothers), with `-pc_gamg_eigenvalues <emin,emax>`. GAMG sets the MG
|