| /petsc/src/ksp/ksp/impls/cg/stcg/ |
| H A D | stcg.c | 18 PetscReal norm_r, norm_d, norm_dp1, norm_p, dMp; in KSPCGSolve_STCG()
57 cg->norm_d = 0.0; in KSPCGSolve_STCG()
89 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
92 cg->norm_d = cg->radius; in KSPCGSolve_STCG()
124 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
127 cg->norm_d = cg->radius; in KSPCGSolve_STCG()
202 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
205 cg->norm_d = cg->radius; in KSPCGSolve_STCG()
228 norm_d = 0.0; in KSPCGSolve_STCG()
260 cg->norm_d = cg->radius; in KSPCGSolve_STCG()
[all …]
|
| H A D | stcgimpl.h | 20 PetscReal norm_d; member
|
| /petsc/src/ksp/ksp/impls/cg/nash/ |
| H A D | nash.c | 19 PetscReal norm_r, norm_d, norm_dp1, norm_p, dMp; in KSPCGSolve_NASH()
60 cg->norm_d = 0.0; in KSPCGSolve_NASH()
92 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
95 cg->norm_d = cg->radius; in KSPCGSolve_NASH()
127 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
130 cg->norm_d = cg->radius; in KSPCGSolve_NASH()
205 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
208 cg->norm_d = cg->radius; in KSPCGSolve_NASH()
231 norm_d = 0.0; in KSPCGSolve_NASH()
263 cg->norm_d = cg->radius; in KSPCGSolve_NASH()
[all …]
|
| H A D | nashimpl.h | 21 PetscReal norm_d; member
|
| /petsc/src/tao/unconstrained/impls/ntr/ |
| H A D | ntr.c | 54 PetscReal norm_d; in TaoSolve_NTR() local
249 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NTR()
253 if (norm_d > 0.0) { in TaoSolve_NTR()
254 tao->trust = norm_d; in TaoSolve_NTR()
273 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NTR()
275 …PetscCheck(norm_d != 0.0, PetscObjectComm((PetscObject)tao), PETSC_ERR_PLIB, "Initial direction ze… in TaoSolve_NTR()
294 tao->trust = tr->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTR()
302 tao->trust = tr->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTR()
316 tao->trust = tr->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTR()
321 tao->trust = tr->alpha2 * PetscMin(tao->trust, norm_d); in TaoSolve_NTR()
[all …]
|
| /petsc/src/tao/unconstrained/impls/nls/ |
| H A D | nls.c | 51 PetscReal norm_d = 0.0, e_min; in TaoSolve_NLS() local
282 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NLS()
286 if (norm_d > 0.0) { in TaoSolve_NLS()
287 tao->trust = norm_d; in TaoSolve_NLS()
306 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NLS()
308 PetscCheck(norm_d != 0.0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Initial direction zero"); in TaoSolve_NLS()
535 tao->trust = nlsP->omega1 * PetscMin(norm_d, tao->trust); in TaoSolve_NLS()
538 tao->trust = nlsP->omega2 * PetscMin(norm_d, tao->trust); in TaoSolve_NLS()
542 tao->trust = nlsP->omega3 * PetscMin(norm_d, tao->trust); in TaoSolve_NLS()
544 tao->trust = PetscMax(nlsP->omega3 * norm_d, tao->trust); in TaoSolve_NLS()
[all …]
|
| /petsc/src/tao/unconstrained/impls/ntl/ |
| H A D | ntl.c | 42 PetscReal norm_d = 0.0; in TaoSolve_NTL() local
248 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NTL()
252 if (norm_d > 0.0) { in TaoSolve_NTL()
253 tao->trust = norm_d; in TaoSolve_NTL()
272 PetscCall(KSPCGGetNormD(tao->ksp, &norm_d)); in TaoSolve_NTL()
274 …PetscCheck(norm_d != 0.0, PetscObjectComm((PetscObject)tao), PETSC_ERR_PLIB, "Initial direction ze… in TaoSolve_NTL()
297 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
306 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
321 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
327 tao->trust = tl->alpha2 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/gltr/ |
| H A D | gltr.c | 69 PetscReal norm_r, norm_d, norm_dp1, norm_p, dMp; in KSPCGSolve_GLTR()
109 cg->norm_d = 0.0; in KSPCGSolve_GLTR()
143 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
146 cg->norm_d = cg->radius; in KSPCGSolve_GLTR()
174 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
177 cg->norm_d = cg->radius; in KSPCGSolve_GLTR()
243 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
246 cg->norm_d = cg->radius; in KSPCGSolve_GLTR()
268 norm_d = 0.0; in KSPCGSolve_GLTR()
299 cg->norm_d = cg->radius; in KSPCGSolve_GLTR()
[all …]
|
| H A D | gltrimpl.h | 33 PetscReal norm_d; member
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cgtype.c | 134 PetscErrorCode KSPCGGetNormD(KSP ksp, PetscReal *norm_d) in KSPCGGetNormD() argument
138 PetscAssertPointer(norm_d, 2); in KSPCGGetNormD()
139 PetscUseMethod(ksp, "KSPCGGetNormD_C", (KSP, PetscReal *), (ksp, norm_d)); in KSPCGGetNormD()
|
| H A D | cg.c | 124 PetscReal r2, norm_p, norm_d, dMp; in KSPSolve_CG() local
158 PetscCall(VecNorm(X, NORM_2, &norm_d)); in KSPSolve_CG()
159 norm_d *= norm_d; in KSPSolve_CG()
163 norm_d = 0.0; in KSPSolve_CG()
273 a = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
292 PetscReal norm_dp1 = norm_d + PetscRealPart(a) * (2.0 * dMp + PetscRealPart(a) * norm_p); in KSPSolve_CG()
297 dp = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
338 PetscCall(VecNorm(X, NORM_2, &norm_d)); in KSPSolve_CG()
339 norm_d *= norm_d; in KSPSolve_CG()
|
| /petsc/src/snes/impls/vi/ss/ |
| H A D | vissimpl.h | 21 PetscScalar norm_d; /* two norm of the descent direction */ member
|