There are a number of different libraries out there that solve a sparse linear system of equations, however I'm finding it difficult to figure out what the differences are.
As far as I can tell there are three major packages: Trilinos, PETSc, and Intel MKL. They can all do sparse matrix solves, they are all fast (as far as I can tell, I haven't been able to find solid benchmarks on any of them), and they are all parallelizable. What I can't find is the differences.
So, what are the differences between the different sparse linear system solvers out there?