Approximate a covariance matrix by a block diagonal matrix with blocks of approximately equal size using Ward's method for hierarchical clustering
divide.sdp(Sigma, max.size)