Implements direct sum (blkdiag) for sparse matrices.

- Provides new blkdiag export
- Provide simple test
- Documentation for blkdiag notes that it is only implemented for sparse
  matrices at the moment.