Using mUMFPACK to solve the sparse matrix