diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-08-31 09:58:48 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-08-31 09:58:48 -0500 |
| commit | d0261d21c76a7d2fe740840d188775504195813e (patch) | |
| tree | 53d88f96a607a0bfab552e03609a27c7c68b4417 /utils/dc-closure-test | |
| parent | c2fea4840ee186eb87ede9f3ed91cfe60b118371 (diff) | |
| download | sddm-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-closure-test')
| -rwxr-xr-x | utils/dc-closure-test | 2 |
1 files changed, 1 insertions, 1 deletions
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) |