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

Regold MG survival biasing test for the 2-component fix.

Browse source
This commit is contained in:
nuclearkevin 2026-07-13 16:06:57 -05:00
parent a913f1ac09
commit 2fc96d0f10
2 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/eigenvalue.cpp
View file

@ -585,7 +585,7 @@ int openmc_get_keff(double* k_combined)
// block matrices of the three estimator equations reduces to scalars here
// See LA-12658-MS, Eqs. 36 and 40.
// We can use variance/covariances for the mean as the division by (n - 1)
// We can use the \Sigma matrix for the mean as the division by (n - 1)
// cancels.
const double f = cov(i, i) + cov(j, j) - 2.0 * cov(i, j);
const double w_1 = (cov(j, j) - cov(i, j)) / f;
@ -593,7 +593,7 @@ int openmc_get_keff(double* k_combined)
k_combined[0] = kv[i] * w_1 + kv[j] * w_2;
// We must use the sums of squares for the variance as the division by
// We must use the S matrix for the variance as the division by
// (n - 1) does not cancel.
const double s_ii = cov(i, i) * (n - 1);
const double s_jj = cov(j, j) * (n - 1);

2
tests/regression_tests/mg_survival_biasing/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
9.889968E-01 9.144186E-03
9.889968E-01 1.823442E-02