cg/stcg/stcg.c

18   PetscReal   norm_r, norm_d, norm_dp1, norm_p, dMp;
154     PetscCall(VecNorm(z, NORM_2, &norm_r)); /* norm_r = |z|      */
158     norm_r = PetscSqrtReal(rr); /* norm_r = |r|      */
162     norm_r = PetscSqrtReal(rz); /* norm_r = |r|_M    */
166     norm_r = 0.0;
170   PetscCall(KSPLogResidualHistory(ksp, norm_r));
171   PetscCall(KSPMonitor(ksp, ksp->its, norm_r));
172   ksp->rnorm = norm_r;
174   PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP));
400       PetscCall(VecNorm(z, NORM_2, &norm_r)); /* norm_r = |z|      */
404       PetscCall(VecNorm(r, NORM_2, &norm_r)); /* norm_r = |r|      */
408       norm_r = PetscSqrtReal(rz); /* norm_r = |r|_M    */
412       norm_r = 0.0;
416     PetscCall(KSPLogResidualHistory(ksp, norm_r));
417     PetscCall(KSPMonitor(ksp, ksp->its, norm_r));
418     ksp->rnorm = norm_r;
420     PetscCall((*ksp->converged)(ksp, ksp->its, norm_r, &ksp->reason, ksp->cnvP));
426       PetscCall(PetscInfo(ksp, "KSPCGSolve_STCG: truncated step: rnorm=%g, radius=%g\n", (double)norm_r, (double)cg->radius));