[WIP] Reintroducing cholesky
Created by: tgiani
Hi, here as discussed in the last pc I m trying to do the following:
- reading the theory cov matrix from file at the level of
ComputeCovMat_basic
inchisquare.cc
, using the functionread_total_covmat
(which btw should now be calledread_theory_covmat
) - construct the total cov matrix inside
ComputeCovMat_basic
, adding to the additive errors also the part coming from the theory cov mat, if present - performing the replica sampling doing a gaussian sampling using just the additive part of the total cov mat, and adding the mjultiplicative uncertainties by hand, like in the old implementation. The gaussian sampling is performed using cholesky decomposition
- use the same function
ComputeCovMat_basic
to build the fit covariance matrix inExperiment::LoadFitCovMat