It has a whole collection of sparse matrix solvers. The one I used was CHOLMOD, which is the one used in MatLab IIRC. Which solver to choose depends on the nature of your matrix.

In this paper a new algorithm for reducing an arbitrary real square matrix to tri-diagonal form using real similarity transformations is described. The method is essentially a generalization of a ...