Routines for Block Diagonal Symmetric Matrices
