diff --git a/Common/src/linear_algebra/CSysSolve.cpp b/Common/src/linear_algebra/CSysSolve.cpp
index c1a356dfbcc..ae4c7fa72eb 100644
--- a/Common/src/linear_algebra/CSysSolve.cpp
+++ b/Common/src/linear_algebra/CSysSolve.cpp
@@ -1137,7 +1137,6 @@ unsigned long CSysSolve<ScalarType>::BCGSTAB_LinSolver(const CSysVector<ScalarTy
 
   ScalarType alpha = 1.0, omega = 1.0, rho = 1.0, rho_prime = 1.0;
   ScalarType norm_best = norm_r;
-  if (config->GetComm_Level() == COMM_FULL) x_best = x;
   p = ScalarType(0.0);
   v = ScalarType(0.0);
   r_0 = r;
@@ -1198,7 +1197,7 @@ unsigned long CSysSolve<ScalarType>::BCGSTAB_LinSolver(const CSysVector<ScalarTy
       /*--- Check if solution has converged, else output the relative residual if necessary ---*/
 
       norm_r = r.norm();
-      if (norm_r < norm_best) {
+      if ((i == 0) || (norm_r < norm_best)) {
         norm_best = norm_r;
         x_best = x;
       }
@@ -1211,7 +1210,8 @@ unsigned long CSysSolve<ScalarType>::BCGSTAB_LinSolver(const CSysVector<ScalarTy
     }
   }
 
-  /*--- BCGSTAB does not reduce the residual monotonically, return the best iterate. ---*/
+  /*--- BCGSTAB does not reduce the residual monotonically, return the best iterate.
+   * The initial guess is not a candidate, returning it would mean no update at all. ---*/
 
   if ((config->GetComm_Level() == COMM_FULL) && (norm_r > norm_best)) {
     x = x_best;
diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py
index 30824edc691..ae6c79c1664 100755
--- a/TestCases/parallel_regression.py
+++ b/TestCases/parallel_regression.py
@@ -957,7 +957,7 @@ def main():
     turbmod_sa_neg_rae2822.cfg_dir   = "turbulence_models/sa/rae2822"
     turbmod_sa_neg_rae2822.cfg_file  = "turb_SA_NEG_RAE2822.cfg"
     turbmod_sa_neg_rae2822.test_iter = 10
-    turbmod_sa_neg_rae2822.test_vals = [1.275994, 1.062588, -2.067822, 1.300662, 0.480881, 0.000000]
+    turbmod_sa_neg_rae2822.test_vals = [1.316535, 1.043996, 4.236921, 1.270979, 0.463385, 0.000000]
     turbmod_sa_neg_rae2822.test_vals_aarch64 = [-1.345593, 1.448310, 1.208721, -0.846597, 1.248410, 0.489117, 0.000000]
     test_list.append(turbmod_sa_neg_rae2822)