Skip to content
Snippets Groups Projects

Abort CGSolver when the defect is invalid

Merged Abort CGSolver when the defect is invalid
feature/abort-on-invalid-defect into master
Merged Jö Fahlke requested to merge feature/abort-on-invalid-defect into master
  1. Apr 11, 2016
    • Jö Fahlke's avatar
      [test][SolverAbort] Check that the solver aborts when given an unsolvable system to solve. · c024db84
      Jö Fahlke authored 
      Here we try to solve A*x=b for x, with
```
      / 1  1 \       / 1 \
  A = |      |,  b = |   |,
      \ 1  1 /       \ 2 /
```
which has no solution.  In the CGSolver, the this leads to a NaN defect after
46 iterations.  In BiCGSTABSolver this leads to "abs(h) < EPSILON" after 1.5
iterations.

BiCGSTABSolver and RestartedGMResSolver also throw SolverAbort when they
detect "breakdown", checking for that is left as a todo.
      c024db84
  3. Mar 17, 2016