123   PetscReal   dp = 0.0;  in KSPSolve_CG()  local
171     PetscCall(VecNorm(Z, NORM_2, &dp)); /*    dp <- z'*z = e'*A'*B'*B*A*e       */  in KSPSolve_CG()
172     KSPCheckNorm(ksp, dp);  in KSPSolve_CG()
175     PetscCall(VecNorm(R, NORM_2, &dp)); /*    dp <- r'*r = e'*A'*A*e            */  in KSPSolve_CG()
176     KSPCheckNorm(ksp, dp);  in KSPSolve_CG()
182     dp = PetscSqrtReal(PetscAbsScalar(beta)); /*    dp <- r'*z = r'*B*r = e'*A'*B*A*e */  in KSPSolve_CG()
185     dp = 0.0;  in KSPSolve_CG()
200   PetscCall(KSPLogResidualHistory(ksp, dp));  in KSPSolve_CG()
201   PetscCall(KSPMonitor(ksp, ksp->its, dp));  in KSPSolve_CG()
202   ksp->rnorm = dp;  in KSPSolve_CG()
204 …PetscCall((*ksp->converged)(ksp, ksp->its, dp, &ksp->reason, ksp->cnvP)); /* test for convergence …  in KSPSolve_CG()
269             PetscCall(VecNorm(R, NORM_2, &dp));  in KSPSolve_CG()
270             a = cg->radius > dp ? 1.0 : cg->radius / dp;  in KSPSolve_CG()
297           dp = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p;  in KSPSolve_CG()
298           PetscCall(VecAXPY(X, dp, P)); /*     x <- x + ap                      */  in KSPSolve_CG()
299           cg->obj += PetscRealPart(dp * (0.5 * dp * dpi - beta));  in KSPSolve_CG()
309       PetscCall(VecNorm(Z, NORM_2, &dp)); /*     dp <- z'*z                       */  in KSPSolve_CG()
310       KSPCheckNorm(ksp, dp);  in KSPSolve_CG()
312       PetscCall(VecNorm(R, NORM_2, &dp)); /*     dp <- r'*r                       */  in KSPSolve_CG()
313       KSPCheckNorm(ksp, dp);  in KSPSolve_CG()
318       dp = PetscSqrtReal(PetscAbsScalar(beta));  in KSPSolve_CG()
320       dp = 0.0;  in KSPSolve_CG()
325     ksp->rnorm = dp;  in KSPSolve_CG()
326     PetscCall(KSPLogResidualHistory(ksp, dp));  in KSPSolve_CG()
327     PetscCall(KSPMonitor(ksp, i + 1, dp));  in KSPSolve_CG()
328     PetscCall((*ksp->converged)(ksp, i + 1, dp, &ksp->reason, ksp->cnvP));  in KSPSolve_CG()
368   PetscReal   dp = 0.0;  in KSPSolve_CG_SingleReduction()  local
407     PetscCall(VecNorm(Z, NORM_2, &dp)); /*    dp <- z'*z = e'*A'*B'*B*A'*e'     */  in KSPSolve_CG_SingleReduction()
408     KSPCheckNorm(ksp, dp);  in KSPSolve_CG_SingleReduction()
411     PetscCall(VecNorm(R, NORM_2, &dp)); /*    dp <- r'*r = e'*A'*A*e            */  in KSPSolve_CG_SingleReduction()
412     KSPCheckNorm(ksp, dp);  in KSPSolve_CG_SingleReduction()
420     dp = PetscSqrtReal(PetscAbsScalar(beta)); /*    dp <- r'*z = r'*B*r = e'*A'*B*A*e */  in KSPSolve_CG_SingleReduction()
423     dp = 0.0;  in KSPSolve_CG_SingleReduction()
428   PetscCall(KSPLogResidualHistory(ksp, dp));  in KSPSolve_CG_SingleReduction()
429   PetscCall(KSPMonitor(ksp, 0, dp));  in KSPSolve_CG_SingleReduction()
430   ksp->rnorm = dp;  in KSPSolve_CG_SingleReduction()
432   PetscCall((*ksp->converged)(ksp, 0, dp, &ksp->reason, ksp->cnvP)); /* test for convergence */  in KSPSolve_CG_SingleReduction()
493       PetscCall(VecNorm(Z, NORM_2, &dp)); /*    dp <- z'*z                        */  in KSPSolve_CG_SingleReduction()
494       KSPCheckNorm(ksp, dp);  in KSPSolve_CG_SingleReduction()
496       PetscCall(VecNorm(R, NORM_2, &dp)); /*    dp <- r'*r                        */  in KSPSolve_CG_SingleReduction()
497       KSPCheckNorm(ksp, dp);  in KSPSolve_CG_SingleReduction()
507       dp = PetscSqrtReal(PetscAbsScalar(beta)); /*    dp <- r'*z = r'*B*r = e'*A'*B*A*e */  in KSPSolve_CG_SingleReduction()
509       dp = 0.0;  in KSPSolve_CG_SingleReduction()
511     ksp->rnorm = dp;  in KSPSolve_CG_SingleReduction()
512     PetscCall(KSPLogResidualHistory(ksp, dp));  in KSPSolve_CG_SingleReduction()
513     PetscCall(KSPMonitor(ksp, i + 1, dp));  in KSPSolve_CG_SingleReduction()
514     PetscCall((*ksp->converged)(ksp, i + 1, dp, &ksp->reason, ksp->cnvP));  in KSPSolve_CG_SingleReduction()