Is there a Haskell package that can compute the reduced row-echelon form of a matrix, with entries from either the rationals, or a finite field of prime size? I am looking for something that will be comparable in speed to computer algebra packages such as Sage.
I have looked at hmatrix and other modules but they appear to be focused on numerical analysis rather than linear algebra. The answer to this question suggests that the answer was a qualified 'no' six years ago.