Try with SparseLU solver in solver.h.
Build option, solver position and sanity check for solver.
Change build option from Release to Debug.
Move solver_ from Solver to Derived solver classes in case
different solvers are needed (e.g., the matrix is not PSD).
Sanity check the norm of the residual from solver.