Skip to content
Snippets Groups Projects
Commit 19660295 authored by oliver.sander_at_tu-dresden.de's avatar oliver.sander_at_tu-dresden.de
Browse files

Fix bug in method lowerTriangularSolve

This makes triangularsolvetest.cc pass.
parent d34a21fe
No related branches found
No related tags found
1 merge request!12Implement addToDiagonal for scalars (interpreted as 1x1 matrices)
Pipeline #33209 failed
...@@ -37,7 +37,7 @@ namespace MatrixVector { ...@@ -37,7 +37,7 @@ namespace MatrixVector {
// Note: We could drop the check for ignore nodes here bcs. b[j] will // Note: We could drop the check for ignore nodes here bcs. b[j] will
// be ignored anyway due to the check above. // be ignored anyway due to the check above.
if (ignore == nullptr or (*ignore)[j].none()) if (ignore == nullptr or (*ignore)[j].none())
b[j] -= x[i] * *cIt; cIt->mmv(x[i], b[j]);
} }
} }
} else { } else {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment