I've been using cholmod from SuiteSparse for over a year now to solve one of my linear systems that is difficult to solve with an iterative method. Are there other sparse direct solvers, such as clique, that I should try sometime? My linear system is SPD and I am solving it on a single node. My code has not yet been parallelized using MPI - it is only a 2D code. Thanks, Dave