So I'm trying to invert a big (449x449) covariance matrix, which is thus symmetric and positive definite. (What I'm trying to do is to invert this matrix as part of a Gaussian Process fitting for the Mauna Loa CO2 dataset.)
This inversion is pretty long, so I wanted to use chol2inv instead of solve. But the chol2inv method gives me a very strange result : a matrix very close to 0 (sum of it is equal to 10^(-13)).
Why would chol2inv give me this?