compute_variance_block_diagonal.BM function

Tips Variances for the BM