-
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.
Jö Fahlke authoredHere 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.