Commit b046420f authored by Dave Moxey's avatar Dave Moxey
Browse files

Fixing iterative solver to avoid overwriting Dirichlet degrees of freedom when

solution is zero on input.

git-svn-id: 305cdda6-5ce1-45b3-a98d-dfc68c8b3305
parent 09d72266
......@@ -148,7 +148,7 @@ namespace Nektar
// If input vector is zero, set zero output and skip solve.
if (vExchange[0] < NekConstants::kNekZeroTol)
Vmath::Zero(nGlobal, pOutput, 1);
Vmath::Zero(nNonDir, tmp = pOutput+nDir, 1);
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