crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Changing check for Tabular to use std. dev.

Browse source
This commit is contained in:
Patrick Shriwise 2022-06-25 14:38:19 +02:00
parent f92ea393f2
commit 3e61076715
1 changed files with 14 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

16
tests/unit_tests/test_stats.py
View file

@ -165,14 +165,26 @@ def test_tabular():
n_samples = 100_000
samples = d.sample(n_samples, seed=100)
assert samples.mean() == pytest.approx(d.mean(), rel=1e-03)
diff = np.abs(samples - d.mean())
# within_1_sigma = np.count_nonzero(diff < samples.std())
# assert within_1_sigma / n_samples >= 0.68
within_2_sigma = np.count_nonzero(diff < 2*samples.std())
assert within_2_sigma / n_samples >= 0.95
within_3_sigma = np.count_nonzero(diff < 3*samples.std())
assert within_3_sigma / n_samples >= 0.99
# test histogram sampling
d = openmc.stats.Tabular(x, p, interpolation='histogram')
d.normalize()
samples = d.sample(n_samples, seed=100)
assert samples.mean() == pytest.approx(d.mean(), rel=1e-03)
diff = np.abs(samples - d.mean())
# within_1_sigma = np.count_nonzero(diff < samples.std())
# assert within_1_sigma / n_samples >= 0.68
within_2_sigma = np.count_nonzero(diff < 2*samples.std())
assert within_2_sigma / n_samples >= 0.95
within_3_sigma = np.count_nonzero(diff < 3*samples.std())
assert within_3_sigma / n_samples >= 0.99
def test_legendre():