Skip to content
Snippets Groups Projects
Commit 538d99fa authored by Jö Fahlke's avatar Jö Fahlke
Browse files

Merge branch 'feature/abort-on-invalid-defect' into 'master'

Abort CGSolver when the defect is invalid

The first commit introduces the exception `SolverAbort`, derived from `ISTLError`.

The second commit makes `CGSolver` abort by throwing `SolverAbort` when the defect becomes invalid (infinite or NaN).  The check is done in every iteration before the convergence check.

These two commits fix #12.

The third and fourth commit make `BiCGStab` and `GMRes` throw `SolverAbort` when they detect a breakdown.  `SolverAbort` seemed more specific when the `ISTLError` they were throwing before.

The fifth commit documents when `SolverAbort` is thrown at the moment.

The sixth commit adds a unit tests that tries to trigger `SolverAbort` and makes sure that it is actually thrown.  It does this for the NaN-check in `CGSolver` and for the "abs(h) < EPSILON"-check in `BiCGSTABSolver` -- there a more conditions in `BiCGSTABSolver` and `RestartedGMResSolver` that throw, but I don't know how to trigger them.



See merge request !36
parents dcade4ee c024db84
No related branches found
No related tags found
Loading
Loading
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment