Matrix classes for matrices that are block-tridiagonal and sparse, and simply "block sparse". These talk together, and furthermore containts an algorithm for inversion of the block-tridiagonal version ...