Covariance Calculation methods seem to break

[arthurmloureiro]

Methods like calling TJPCov and using the internal gaussian seem to be quite unstable leading to non-positive definite matrices and breaking the cholesky decomposition.

I haven't yet found a way to get a new covariance without causing the problem above.

[bruno-moraes]

Can confirm. I am having issues with 'gaus_internal' covariance generation in a LCDM + neutrinos run. I can provide a config_test.yml file if needed, but it should be easily reproducible.

[arthurmloureiro]

Tagging @felipeaoli here in case he has some insights on what is going on :D

[felipeaoli]

please, send us @bruno-moraes, share the config_test.yml. Thank you!