| /petsc/src/vec/vec/tests/ |
| H A D | ex47.c | 10 PetscReal norm, dnorm; in main() local
49 PetscCall(VecNorm(x, NORM_2, &dnorm)); in main()
50 …ck(norm / dnorm <= 1.e-6, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Vec read in 'noBsDim' does not match … in main()
57 PetscCall(VecNorm(x, NORM_2, &dnorm)); in main()
58 …eck(norm / dnorm <= 1.e-6, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Vec read in 'bsDim' does not match v… in main()
|
| /petsc/src/tao/bound/impls/bncg/ |
| H A D | bncg.c | 436 …l gnorm2, PetscReal step, PetscReal fold, PetscReal gnorm2_old, PetscReal dnorm, PetscBool pcgd_fa… in TaoBNCGStepDirectionUpdate() argument
452 if (step * dnorm < PETSC_MACHINE_EPSILON || step * dk_yk < PETSC_MACHINE_EPSILON) { in TaoBNCGStepDirectionUpdate()
502 cg->sts = step * step * dnorm * dnorm; in TaoBNCGStepDirectionUpdate()
518 cg->sts = step * step * dnorm * dnorm; in TaoBNCGStepDirectionUpdate()
538 …PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, step * step * dnorm * dnorm, &tau_k, c… in TaoBNCGStepDirectionUpdate()
551 snorm = step * dnorm; in TaoBNCGStepDirectionUpdate()
616 snorm = dnorm * step; in TaoBNCGStepDirectionUpdate()
629 …x(PetscMax(beta, cg->hz_eta * tau_k * gd_old / (dnorm * dnorm)), cg->dk_eta * tau_k * gd / (dnorm … in TaoBNCGStepDirectionUpdate()
651 …ta = PetscMax(PetscMax(beta, cg->hz_eta * gd_old / (dnorm * dnorm)), cg->dk_eta * gd / (dnorm * dn… in TaoBNCGStepDirectionUpdate()
664 snorm = step * dnorm; in TaoBNCGStepDirectionUpdate()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex204.c | 51 PetscReal dnorm; in main() local
86 PetscCall(VecNorm(d_vcl, NORM_INFINITY, &dnorm)); in main()
87 …PetscCheck(dnorm <= tol, PETSC_COMM_WORLD, PETSC_ERR_ARG_INCOMP, "Sequential CPU and MPI ViennaCL … in main()
104 PetscReal dnorm; in main() local
139 PetscCall(VecNorm(d_vcl, NORM_INFINITY, &dnorm)); in main()
140 …PetscCheck(dnorm <= tol, PETSC_COMM_WORLD, PETSC_ERR_ARG_INCOMP, "MPI CPU and MPI ViennaCL Vector … in main()
156 PetscReal dnorm; in main() local
190 PetscCall(VecNorm(d_vcl, NORM_INFINITY, &dnorm)); in main()
191 …PetscCheck(dnorm <= tol, PETSC_COMM_WORLD, PETSC_ERR_ARG_INCOMP, "MPI CPU and MPI ViennaCL Vector … in main()
208 PetscReal dnorm; in main() local
[all …]
|
| H A D | ex235.c | 65 PetscReal Aijnorm, Aijdiagnorm, Bnorm, dnorm; in main() local
87 PetscCall(VecNorm(Adiag, NORM_2, &dnorm)); in main()
89 …dnorm <= 100.0 * PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "Altered matrices diagon… in main()
|
| /petsc/src/tao/bound/impls/bnk/ |
| H A D | bnk.c | 451 PetscCall(KSPCGGetNormD(tao->ksp, &bnk->dnorm)); in TaoBNKComputeStep()
455 if (bnk->dnorm > 0.0) { in TaoBNKComputeStep()
456 tao->trust = bnk->dnorm; in TaoBNKComputeStep()
475 PetscCall(KSPCGGetNormD(tao->ksp, &bnk->dnorm)); in TaoBNKComputeStep()
477 …PetscCheck(bnk->dnorm != 0.0, PetscObjectComm((PetscObject)tao), PETSC_ERR_PLIB, "Initial directio… in TaoBNKComputeStep()
809 tao->trust = bnk->omega1 * PetscMin(bnk->dnorm, tao->trust); in TaoBNKUpdateTrustRadius()
812 tao->trust = bnk->omega2 * PetscMin(bnk->dnorm, tao->trust); in TaoBNKUpdateTrustRadius()
816 tao->trust = bnk->omega3 * PetscMin(bnk->dnorm, tao->trust); in TaoBNKUpdateTrustRadius()
818 tao->trust = PetscMax(bnk->omega3 * bnk->dnorm, tao->trust); in TaoBNKUpdateTrustRadius()
822 tao->trust = PetscMax(bnk->omega4 * bnk->dnorm, tao->trust); in TaoBNKUpdateTrustRadius()
[all …]
|
| H A D | bnk.h | 33 PetscReal fold, f, gnorm, dnorm; member
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | ngmresfunc.c | 115 …tscInt l, Vec X, Vec F, Vec XM, Vec FM, Vec XA, Vec FA, Vec D, PetscReal *dnorm, PetscReal *dminno… in SNESNGMRESNorms_Private() argument
135 if (dnorm) { in SNESNGMRESNorms_Private()
137 PetscCall(VecNormBegin(D, NORM_2, dnorm)); in SNESNGMRESNorms_Private()
151 if (dnorm) PetscCall(VecNormEnd(D, NORM_2, dnorm)); in SNESNGMRESNorms_Private()
163 …eal xAnorm, PetscReal fAnorm, PetscReal yAnorm, PetscReal objA, PetscReal dnorm, PetscReal objmin,… in SNESNGMRESSelect_Private() argument
221 …)(/* A */ (objA < ngmres->gammaA * objmin) && /* B */ (ngmres->epsilonB * dnorm < dminnorm || objA… in SNESNGMRESSelect_Private()
250 …nes, PetscInt l, PetscReal obj, PetscReal objM, PetscReal objA, PetscReal dnorm, PetscReal objmin,… in SNESNGMRESSelectRestart_Private() argument
263 if (ngmres->epsilonB * dnorm > dminnorm && objA > ngmres->deltaB * objmin && l > 0) { in SNESNGMRESSelectRestart_Private()
264 …es->monitor, "difference restart: %e > %e\n", (double)(ngmres->epsilonB * dnorm), (double)dminnorm… in SNESNGMRESSelectRestart_Private()
|
| H A D | anderson.c | 32 PetscReal dnorm, dminnorm = 0.0, fminnorm; in SNESSolve_Anderson() local
122 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, &dnorm, &dminnorm, NULL, NULL,… in SNESSolve_Anderson()
123 …PetscCall(SNESNGMRESSelectRestart_Private(snes, l, snes->norm, fMnorm, fnorm, dnorm, fminnorm, dmi… in SNESSolve_Anderson()
|
| H A D | snesngmres.c | 148 PetscReal dnorm = 0.0, dminnorm = 0.0; in SNESSolve_NGMRES() local
248 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, &dnorm, &dminnorm, &xMnorm, NU… in SNESSolve_NGMRES()
259 …, FM, xMnorm, fMnorm, yMnorm, objM, XA, FA, xAnorm, fAnorm, yAnorm, objA, dnorm, objmin, dminnorm,… in SNESSolve_NGMRES()
268 …PetscCall(SNESNGMRESSelectRestart_Private(snes, l, obj, objM, objA, dnorm, objmin, dminnorm, &sele… in SNESSolve_NGMRES()
|