Searched refs:norm_p (Results 1 – 4 of 4) sorted by relevance
| /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()
231 norm_p = rz; in KSPCGSolve_STCG()
235 PetscCall(VecDot(p, p, &norm_p)); in KSPCGSolve_STCG()
253 if (cg->radius != 0.0 && norm_p > 0.0) { in KSPCGSolve_STCG()
259 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_STCG()
318 norm_dp1 = norm_d + alpha * (2.0 * dMp + alpha * norm_p); in KSPCGSolve_STCG()
328 if (norm_p > 0.0) { in KSPCGSolve_STCG()
334 step = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPCGSolve_STCG()
463 dMp = beta * (dMp + alpha * norm_p); in KSPCGSolve_STCG()
464 norm_p = beta * (rzm1 + beta * norm_p); in KSPCGSolve_STCG()
[all …]
|
| /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()
234 norm_p = rz; in KSPCGSolve_NASH()
238 PetscCall(VecDot(p, p, &norm_p)); in KSPCGSolve_NASH()
256 if (cg->radius && norm_p > 0.0) { in KSPCGSolve_NASH()
262 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_NASH()
320 norm_dp1 = norm_d + alpha * (2.0 * dMp + alpha * norm_p); in KSPCGSolve_NASH()
330 if (norm_p > 0.0) { in KSPCGSolve_NASH()
335 step = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPCGSolve_NASH()
464 dMp = beta * (dMp + alpha * norm_p); in KSPCGSolve_NASH()
465 norm_p = beta * (rzm1 + beta * norm_p); in KSPCGSolve_NASH()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/ |
| H A D | cg.c | 124 PetscReal r2, norm_p, norm_d, dMp; in KSPSolve_CG() local
237 PetscCall(VecNorm(P, NORM_2, &norm_p)); in KSPSolve_CG()
238 norm_p *= norm_p; in KSPSolve_CG()
252 PetscCall(VecNorm(P, NORM_2, &norm_p)); in KSPSolve_CG()
253 norm_p *= norm_p; in KSPSolve_CG()
266 if (norm_p > 0.0) { in KSPSolve_CG()
267 a = PetscSqrtReal(r2 / norm_p); in KSPSolve_CG()
272 } else if (norm_p > 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()
[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()
271 norm_p = rz; in KSPCGSolve_GLTR()
275 PetscCall(VecDot(p, p, &norm_p)); in KSPCGSolve_GLTR()
295 if (cg->radius && norm_p > 0.0) { in KSPCGSolve_GLTR()
298 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_GLTR()
351 if (cg->radius && norm_p > 0.0) { in KSPCGSolve_GLTR()
353 step = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPCGSolve_GLTR()
368 norm_dp1 = norm_d + alpha * (2.0 * dMp + alpha * norm_p); in KSPCGSolve_GLTR()
375 if (norm_p > 0.0) { in KSPCGSolve_GLTR()
378 step = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPCGSolve_GLTR()
[all …]
|