From d0261d21c76a7d2fe740840d188775504195813e Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 31 Aug 2020 09:58:48 -0500 Subject: add estimate_errors to chi2 analysis This commit updates the estimate_errors() function so that it works without a list of constraints and uses arrays of low and high bounds passed in instead of hardcoded constraints. I can now call this function from the chi2 analysis to get the stepsizes before running the MCMC. --- utils/dc-closure-test | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'utils/dc-closure-test') diff --git a/utils/dc-closure-test b/utils/dc-closure-test index 72ae95f..bf26662 100755 --- a/utils/dc-closure-test +++ b/utils/dc-closure-test @@ -376,7 +376,7 @@ def fit(data, sacrifice, steps): print("nll(xopt) = ", nll(xopt)) #print("n = ", opt.get_numevals()) - stepsizes = estimate_errors(nll,xopt,constraints) + stepsizes = estimate_errors(nll,xopt,low,high,constraints) with printoptions(precision=3, suppress=True): print("Errors: ", stepsizes) -- cgit