Lines Matching refs:dp
93 PetscReal dp, dpold, w, dpest, tau, psi, cm, r0;
166 r0 = dp = PetscSqrtReal(PetscRealPart(dpi));
169 if (monitor) Kokkos::single(Kokkos::PerTeam(team), [=]() { printf("%3d KSP Residual norm %14.12e\n", 0, (double)dp); });
171 if (dp < atol) {
188 tau = dp;
189 dpold = dp;
229 dp = PetscSqrtReal(PetscRealPart(dpi));
231 if (!m) w = PetscSqrtReal(dp * dpold);
232 else w = dp;
305 dpold = dp;
335 PetscReal dp, r0;
395 /* dp <- r'*r */
398 r0 = dp = PetscSqrtReal(PetscRealPart(dpi));
400 if (monitor) Kokkos::single(Kokkos::PerTeam(team), [=]() { printf("%3d KSP Residual norm %14.12e\n", 0, (double)dp); });
402 if (dp < atol) {
466 /* dp <- r'*r */
469 dp = PetscSqrtReal(PetscRealPart(dpi));
471 if (monitor) Kokkos::single(Kokkos::PerTeam(team), [=]() { printf("%3d KSP Residual norm %14.12e\n", it + 1, (double)dp); });
473 if (dp < atol) {
477 if (dp / r0 < rtol) {
482 if (dp / r0 > dtol) {
484 Kokkos::single(Kokkos::PerTeam(team), [=]() { printf("ERROR block %d diverged: %d it, res=%e, r_0=%e (BICG does this)\n", team.league_rank(), it, dp, r0); });
488 if (dp / r0 > dtol) {
496 Kokkos::single(Kokkos::PerTeam(team), [=]() { printf("ERROR block %d diverged: BICG %d it, res=%e, r_0=%e r_res=%e\n", team.league_rank(), it, dp, r0, dp / r0); });