Skip to content
Snippets Groups Projects
Commit 33f01d08 authored by Max Kahnt's avatar Max Kahnt
Browse files

Cleanup: Doc and asserts.

parent 0825792a
No related branches found
No related tags found
1 merge request!8Add generalized BlockGS steps with modular local solvers.
Pipeline #
...@@ -164,7 +164,8 @@ auto truncate(LinearSolver&& linearSolver) { ...@@ -164,7 +164,8 @@ auto truncate(LinearSolver&& linearSolver) {
auto mTruncated = m; auto mTruncated = m;
auto bTruncated = b; auto bTruncated = b;
assert(b.size() == m.N() && m.N() == m.M()); assert(b.size() == m.N());
assert(m.N() == m.M());
size_t blockSize = b.size(); size_t blockSize = b.size();
for (size_t j = 0; j < blockSize; ++j) { for (size_t j = 0; j < blockSize; ++j) {
if (not ignore[j]) if (not ignore[j])
...@@ -260,8 +261,7 @@ auto cg(size_t maxIter = LinearSolvers::defaultCgMaxIter, ...@@ -260,8 +261,7 @@ auto cg(size_t maxIter = LinearSolvers::defaultCgMaxIter,
/** /**
* \brief A Gauss--Seidel-type linear iteration step. * \brief A Gauss--Seidel-type linear iteration step.
* \param localSolver The solver how to solve the linear block correction * \param localSolver The solver for the linear block correction problems.
* problems.
*/ */
template <class Matrix, class Vector, class BitVector, class LocalSolver> template <class Matrix, class Vector, class BitVector, class LocalSolver>
struct BlockGSStep : public LinearIterationStep<Matrix, Vector, BitVector> { struct BlockGSStep : public LinearIterationStep<Matrix, Vector, BitVector> {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment