| /petsc/src/ksp/ksp/utils/lmvm/rescale/ |
| H A D | symbrdnrescale.c | 25 for (PetscInt i = 0; i < next - oldest; ++i) signew += ldb->yts[i] / ldb->yty[i]; in SymBroydenRescaleUpdateScalar()
30 for (PetscInt i = 0; i < next - oldest; ++i) signew += ldb->sts[i] / ldb->yts[i]; in SymBroydenRescaleUpdateScalar()
36 b += ldb->yts[i]; in SymBroydenRescaleUpdateScalar()
65 …nRescale ldb, Vec D, Vec s, Vec y, Vec V, Vec W, Vec BFGS, Vec DFP, PetscReal theta, PetscReal yts) in DiagonalUpdate() argument
77 PetscCall(VecAXPY(D, 1.0 / yts, ldb->V)); in DiagonalUpdate()
108 PetscCall(VecAXPBY(DFP, stDs / yts, -1.0, V)); in DiagonalUpdate()
114 PetscCall(VecAXPY(D, 1.0 / yts, DFP)); in DiagonalUpdate()
117 PetscCall(VecAXPBYPCZ(D, 1.0 - theta, theta / yts, 1.0, BFGS, DFP)); in DiagonalUpdate()
128 PetscScalar yts; in SymBroydenRescaleUpdateDiagonal() local
134 PetscCall(MatLMVMProductsGetDiagonalValue(B, LMBASIS_Y, LMBASIS_S, next - 1, &yts)); in SymBroydenRescaleUpdateDiagonal()
[all …]
|
| H A D | symbrdnrescale.h | 12 …PetscReal *yty, *sts, *yts; /* scalar arrays for recycling dot product… member
|
| /petsc/src/tao/bound/impls/bncg/ |
| H A D | bncg.c | 257 PetscErrorCode TaoBNCGComputeScalarScaling(PetscReal yty, PetscReal yts, PetscReal sts, PetscReal *… in TaoBNCGComputeScalarScaling() argument
263 if (1.0 == alpha) *scale = yts / yty; in TaoBNCGComputeScalarScaling()
264 else if (0.0 == alpha) *scale = sts / yts; in TaoBNCGComputeScalarScaling()
268 b = yts; in TaoBNCGComputeScalarScaling()
594 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, step * dk_yk, cg->yts, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
617 cg->yts = step * dk_yk; in TaoBNCGStepDirectionUpdate()
624 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
665 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
668 PetscCall(TaoBNCGComputeScalarScaling(ynorm2, cg->yts, snorm * snorm, &tau_k, cg->alpha)); in TaoBNCGStepDirectionUpdate()
707 cg->yts = dk_yk * step; in TaoBNCGStepDirectionUpdate()
[all …]
|
| H A D | bncg.h | 28 PetscReal as_step, as_tol, yts, yty, sts; member
|
| /petsc/src/ksp/ksp/utils/lmvm/symbrdn/ |
| H A D | symbrdn.c | 6 static inline PetscScalar PhiFromPsi(PetscScalar psi, PetscScalar yts, PetscScalar stBs, PetscScala… in PhiFromPsi() argument
8 PetscScalar numer = (1.0 - psi) * PetscRealPart(PetscConj(yts) * yts); in PhiFromPsi()
260 PetscScalar yts; in SymBroydenCompactDenseUpdateArrays() local
277 yts = YtS[i + i * ldays]; in SymBroydenCompactDenseUpdateArrays()
279 phi = PhiFromPsi(Psi[i], yts, sBis, YtHkY[i]); in SymBroydenCompactDenseUpdateArrays()
284 beta = -phi / yts; in SymBroydenCompactDenseUpdateArrays()
285 delta = (1.0 + phi * sBis / yts) / yts; in SymBroydenCompactDenseUpdateArrays()
310 yts = YtS[i + i * ldays]; in SymBroydenCompactDenseUpdateArrays()
312 phi = PhiFromPsi(Psi[i], yts, sBis, YtHkY[i]); in SymBroydenCompactDenseUpdateArrays()
318 beta = -phi / yts; in SymBroydenCompactDenseUpdateArrays()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/dense/ |
| H A D | denseqn.h | 35 PetscReal *ytq, *stp, *yts; member
|
| H A D | denseqn.c | 81 PetscCall(PetscFree(lqn->yts)); in MatLMVMDQNResetDestructive()
196 PetscCall(PetscMalloc1(lmvm->m, &lqn->yts)); in MatAllocate_LMVMDQN_Internal()
334 if (lqn->use_recursive && (is_dbfgs || is_ddfp)) lqn->yts[idx] = PetscRealPart(curvature); in MatUpdate_LMVMDQN()
503 mlqn->yts[i] = blqn->yts[i]; in MatCopy_LMVMDQN()
908 …PetscCall(VecAXPBYPCZ(lbfgs->PQ[idx], -sjtpi / lbfgs->stp[idx_j], yjtsi / lbfgs->yts[idx_j], 1.0, … in MatMult_LMVMDBFGS()
934 …PetscCall(VecAXPBYPCZ(Z, -stz / lbfgs->stp[idx], ytx / lbfgs->yts[idx], 1.0, lbfgs->PQ[idx], lbfgs… in MatMult_LMVMDBFGS()
1273 …PetscCall(VecAXPBYPCZ(ldfp->PQ[idx], -yjtqi / ldfp->ytq[idx_j], sjtyi / ldfp->yts[idx_j], 1.0, ldf… in MatSolve_LMVMDDFP()
1299 …PetscCall(VecAXPBYPCZ(dX, -ytx / ldfp->ytq[idx], stf / ldfp->yts[idx], 1.0, ldfp->PQ[idx], ldfp->c… in MatSolve_LMVMDDFP()
|
| /petsc/src/ksp/ksp/utils/lmvm/tests/ |
| H A D | ex1.c | 234 PetscScalar stBs, yts; in SymmetricBroydenUpdate_Explicit() local
240 PetscCall(VecDot(s, y, &yts)); in SymmetricBroydenUpdate_Explicit()
241 PetscCall(RankOneAXPY(B, (yts + phi * stBs) / (yts * yts), y, y)); in SymmetricBroydenUpdate_Explicit()
242 PetscCall(RankOneAXPY(B, -phi / yts, y, Bs)); in SymmetricBroydenUpdate_Explicit()
243 PetscCall(RankOneAXPY(B, -phi / yts, Bs, y)); in SymmetricBroydenUpdate_Explicit()
|