aboutsummaryrefslogtreecommitdiff
path: root/utils/dc
diff options
context:
space:
mode:
authortlatorre <tlatorre@uchicago.edu>2020-08-31 09:58:48 -0500
committertlatorre <tlatorre@uchicago.edu>2020-08-31 09:58:48 -0500
commitd0261d21c76a7d2fe740840d188775504195813e (patch)
tree53d88f96a607a0bfab552e03609a27c7c68b4417 /utils/dc
parentc2fea4840ee186eb87ede9f3ed91cfe60b118371 (diff)
downloadsddm-d0261d21c76a7d2fe740840d188775504195813e.tar.gz
sddm-d0261d21c76a7d2fe740840d188775504195813e.tar.bz2
sddm-d0261d21c76a7d2fe740840d188775504195813e.zip
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.
Diffstat (limited to 'utils/dc')
-rwxr-xr-xutils/dc2
1 files changed, 1 insertions, 1 deletions
diff --git a/utils/dc b/utils/dc
index a267090..d116c65 100755
--- a/utils/dc
+++ b/utils/dc
@@ -459,7 +459,7 @@ if __name__ == '__main__':
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)