DOC: improvements.
Much faster implementation of set_submatrix
Use linalg to solve, speed up considerable.
Allow solving with sparse matrices. (speed up)