Commit 697541fc authored by Chris Cantwell's avatar Chris Cantwell

Disabled the updating of m_rhs_magnitude on each solve in SC iter CG.

parent 46013265
......@@ -513,7 +513,7 @@ namespace Nektar
rho_new = vExchange[0];
mu = vExchange[1];
eps = vExchange[2];
cout << "CG Iteration " << k << ", " << eps << ", " << m_rhs_magnitude << endl;
m_totalIterations++;
// test if norm is within tolerance
if (eps < m_tolerance * m_tolerance * m_rhs_magnitude)
......@@ -546,7 +546,9 @@ namespace Nektar
m_expList.lock()->GetComm()->GetRowComm()->AllReduce(
vExchange, Nektar::LibUtilities::ReduceSum);
cout << "Rank " << m_expList.lock()->GetComm()->GetRank() << ", vExchange[0] = " << vExchange[0] << endl;
int x;
cin >> x;
// To ensure that very different rhs values are not being
// used in subsequent solvers such as the velocit solve in
// INC NS. If this works we then need to work out a better
......
......@@ -540,7 +540,9 @@ namespace Nektar
m_precon->BuildPreconditioner();
}
if (m_rhs_magnitude == NekConstants::kNekUnsetDouble) {
Set_Rhs_Magnitude(F_GlobBnd);
}
return m_S1Blk;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment