Got rid of the inversion of the Covariance matrix, since testing proved it to be wrong. Instead, I am solving "Ax = b" directly for use in the covariant chisq variables. However, that still has a nonsensical result, since the coChisqs end up being large and negative.
Probably, I will just discard the idea of the covariant chisq after this commit.