Skip to content

Commit

Permalink
make cgs iterations match stopping checks
Browse files Browse the repository at this point in the history
  • Loading branch information
upsj committed Feb 5, 2021
1 parent 20a3487 commit d524823
Showing 1 changed file with 0 additions and 6 deletions.
6 changes: 0 additions & 6 deletions core/solver/cgs.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -139,9 +139,6 @@ void Cgs<ValueType>::apply_impl(const LinOp *b, LinOp *x) const

int iter = 0;
/* Memory movement summary:
* Per iteration:
* 14n * values + matrix/preconditioner storage
* Two iterations:
* 28n * values + 2 * matrix/preconditioner storage
* 2x SpMV: 4n * values + 2 * storage
* 2x Preconditioner: 4n * values + 2 * storage
Expand Down Expand Up @@ -169,9 +166,6 @@ void Cgs<ValueType>::apply_impl(const LinOp *b, LinOp *x) const
alpha.get(), rho.get(), gamma.get(),
&stop_status));

++iter;
this->template log<log::Logger::iteration_complete>(this, iter, r.get(),
dense_x);
get_preconditioner()->apply(t.get(), u_hat.get());
system_matrix_->apply(u_hat.get(), t.get());
// r = r - alpha * t
Expand Down

0 comments on commit d524823

Please sign in to comment.