explicitly set beta = 0 on first iter

parthenon-hpc-lab · Sep 18, 2024 · 38543b0 · 38543b0
1 parent 507cbd3
commit 38543b0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/solvers/cg_solver.hpp b/src/solvers/cg_solver.hpp
@@ -181,7 +181,7 @@ class CGSolver {
     auto correct_p = itl.AddTask(
         get_ru, "p <- u + beta p",
         [](CGSolver *solver, std::shared_ptr<MeshData<Real>> &md) {
-          Real beta = solver->ru.val / solver->ru_old;
+          Real beta = solver->iter_counter > 0 ? solver->ru.val / solver->ru_old : 0.0;
           return AddFieldsAndStore<u, p, p>(md, 1.0, beta);
         },
         this, md);