Routines for Block Diagonal Symmetric matrices