Fix local implementation of CG with Null preconditioner

Issue/feature addressed

Currently, there is no assembly performed with the Null preconditioner when solving in local space. This MR fix this problem.

Proposed solution

Implementation

Tests

Suggested reviewers

Please suggest any people who would be appropriate to review your code.

Notes

Please add any other information that could be useful for reviewers.

Checklist

• [ ] Functions and classes, or changes to them, are documented.
• [ ] User guide/documentation is updated.
• Changelog is updated.
• Suitable tests added for new functionality.
• Contributed code is correctly formatted. (See the contributing guidelines).
• [ ] License added to any new files.
• No extraneous files have been added (e.g. compiler output or test data files).

library/LibUtilities/LinearAlgebra/NekLinSysIterCG.cpp

@@ -157,11 +157,9 @@ void NekLinSysIterCG::DoConjugateGradient(
     }
 
     m_operator.DoNekSysPrecon(r_A, tmp = w_A + nDir);
-
     m_operator.DoNekSysLhsEval(w_A, s_A);
 
     vExchange[0] = Vmath::Dot2(nNonDir, r_A, w_A + nDir, m_map + nDir);
-
     vExchange[1] = Vmath::Dot2(nNonDir, s_A + nDir, w_A + nDir, m_map + nDir);
 
     m_rowComm->AllReduce(vExchange, Nektar::LibUtilities::ReduceSum);