Routines for block diagonal symmetric matrices
