Holdout cross-validation for large non-Gaussian covariance matrix estimation using Weingarten calculus