Hi, a few minutes ago I did a pull request for the addition of the new solver mkl-pardiso. This solver supports sequential aij matrices. So any question or suggestion is welcome. Thanks