BUG: matrix.approximate: case where alpha is finite but alpha**2 is infinity is handled now correctly
BUG: matrix.approximate: case where alpha is finite but alpha**2 is infinity is handled now correctly

jor-/matrix-decomposition