Skip to content
Snippets Groups Projects
Commit b605bfa3 authored by Elias Pipping's avatar Elias Pipping
Browse files

Make compile for blocksize > 1

parent c1d21118
No related branches found
No related tags found
No related merge requests found
Pipeline #
......@@ -99,7 +99,8 @@ static constexpr double defaultCgTol = 1e-15;
template <class MBlock, class VBlock>
VBlock cg(const MBlock& m, VBlock r, size_t maxIter = defaultCgMaxIter,
double tol = defaultCgTol) {
VBlock x = 0;
VBlock x;
x = 0;
VBlock p = r;
auto q = r;
auto rOldSquared = r * r;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment