check_full_blockdiag function

Replace a full with a block-diagonal variance covariance matrix