This solver is of particular interest to us on cuML because of 1) its ability find the eigenvalues when given a sparse input and 2) its ability to accept a custom linear operation that we can use to perform mean centering at the vector level. It's important that we test sparse inputs here as well. 