Before changing blkmat
Fix blkdiags.

Implement csr_matrix instead of bsr_matrix.
Implement sparse.interleave.
Fix bugs in sparse.
Going back in time to fix bug
Fix shadow trajectory integration bug. Ignore manifolds not expanding/contracting.