Work around to non positive-definite covariance matrix

This seems to be broken for the baseline but using the full covariance data fails on Cholensky Decomposition for the covariance data