added function for Multi Variable Block Diagonal Covariance Approximation for packed data measured on multidimensional grid