#include <Eigen/src/SparseCholesky/SimplicialCholesky.h>
keeps off-diagonal entries; drops diagonal entries