bdsmatrix
Routines for Block Diagonal Symmetric Matrices