| /petsc/src/ksp/ksp/impls/cg/ ! |
| H A D | cg.c | 122 PetscScalar dpi = 0.0, a = 1.0, beta, betaold = 1.0, b = 0, *e = NULL, *d = NULL, dpiold; in KSPSolve_CG() local
256 dpiold = dpi; in KSPSolve_CG()
258 PetscCall(VecXDot(P, W, &dpi)); /* dpi <- p'w */ in KSPSolve_CG()
259 KSPCheckDot(ksp, dpi); in KSPSolve_CG()
262 …if ((dpi == 0.0) || ((i > 0) && ((PetscSign(PetscRealPart(dpi)) * PetscSign(PetscRealPart(dpiold))… in KSPSolve_CG()
276 cg->obj += PetscRealPart(a * (0.5 * a * dpi - betaold)); in KSPSolve_CG()
280 …tscCall(PetscInfo(ksp, "converged due to negative curvature: %g\n", (double)(PetscRealPart(dpi)))); in KSPSolve_CG()
283 …verged due to indefinite matrix, dpi %g, dpiold %g", (double)PetscRealPart(dpi), (double)PetscReal… in KSPSolve_CG()
289 a = beta / dpi; /* a = beta/p'w */ in KSPSolve_CG()
299 cg->obj += PetscRealPart(dp * (0.5 * dp * dpi - beta)); in KSPSolve_CG()
[all …]
|
| /petsc/src/ksp/ksp/impls/bicg/ ! |
| H A D | bicg.c | 14 PetscScalar dpi, a = 1.0, beta, betaold = 1.0, b, ma; in KSPSolve_BiCG() local
79 PetscCall(VecDot(Zr, Pl, &dpi)); /* dpi <- z'p */ in KSPSolve_BiCG()
80 KSPCheckDot(ksp, dpi); in KSPSolve_BiCG()
81 a = beta / dpi; /* a = beta/p'z */ in KSPSolve_BiCG()
|
| /petsc/src/ksp/ksp/impls/fcg/ ! |
| H A D | fcg.c | 76 PetscScalar alpha = 0.0, beta = 0.0, dpi = 0.0, dpiold, s; in KSPSolve_FCG() local
205 dpiold = dpi; in KSPSolve_FCG()
206 PetscCall(VecXDot(Pcurr, Ccurr, &dpi)); /* dpi <- pi'*w */ in KSPSolve_FCG()
207 …if ((dpi == 0.0) || ((i > 0) && ((PetscSign(PetscRealPart(dpi)) * PetscSign(PetscRealPart(dpiold))… in KSPSolve_FCG()
208 …verged due to indefinite matrix, dpi %g, dpiold %g", (double)PetscRealPart(dpi), (double)PetscReal… in KSPSolve_FCG()
214 alpha = beta / dpi; /* alpha <- beta/dpi */ in KSPSolve_FCG()
263 PetscCall(VecScale(Ccurr, 1.0 / dpi)); /* w <- ci/dpi */ in KSPSolve_FCG()
|
| /petsc/src/ksp/ksp/impls/cg/cgne/ ! |
| H A D | cgne.c | 44 PetscScalar dpi, a = 1.0, beta, betaold = 1.0, b = 0, *e = NULL, *d = NULL; in KSPSolve_CGNE() local
136 PetscCall(VecXDot(P, Z, &dpi)); /* dpi <- z'p */ in KSPSolve_CGNE()
137 KSPCheckDot(ksp, dpi); in KSPSolve_CGNE()
138 a = beta / dpi; /* a = beta/p'z */ in KSPSolve_CGNE()
|
| /petsc/src/ksp/pc/impls/bjacobi/bjkokkos/ ! |
| H A D | bjkokkos.kokkos.cxx | 95 PetscScalar *ptr = work_space_shared, rho, rhoold, a, s, b, eta, etaold, psiold, cf, dpi; in BJSolve_TFQMR() local
164 …e(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += R[idx] * PetscConj(R[idx]); }, dpi); in BJSolve_TFQMR()
166 r0 = dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_TFQMR()
227 …e(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += R[idx] * PetscConj(R[idx]); }, dpi); in BJSolve_TFQMR()
229 dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_TFQMR()
337 PetscScalar *ptr = work_space_shared, dpi, a = 1.0, beta, betaold = 1.0, t1, t2; in BJSolve_BICG() local
396 …team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += Rr[idx] * PetscConj(Rr[idx]); }, dpi); in BJSolve_BICG()
398 r0 = dp = PetscSqrtReal(PetscRealPart(dpi)); in BJSolve_BICG()
448 …team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += Zr[idx] * PetscConj(Pl[idx]); }, dpi); in BJSolve_BICG()
450 if (dpi == 0) { in BJSolve_BICG()
[all …]
|