Fix MG mode delayed nu-fission indexing bug

This commit is contained in:
Sterling Harper 2019-02-07 13:13:43 -05:00
parent 4878a4601d
commit 3ff2410dfa
2 changed files with 2 additions and 2 deletions

View file

@ -77,7 +77,7 @@ class XsData {
// [angle][delayed group]
xt::xtensor<double, 2> decay_rate;
// delayed_nu_fission has the following dimensions:
// [angle][incoming group][delayed group]
// [angle][delayed group][incoming group]
xt::xtensor<double, 3> delayed_nu_fission;
// chi_prompt has the following dimensions:
// [angle][incoming group][outgoing group]

View file

@ -467,7 +467,7 @@ Mgxs::get_xs(int xstype, int gin, int* gout, double* mu, int* dg)
val = xs_t->delayed_nu_fission(a, *dg, gin);
} else {
val = 0.;
for (int d = 0; d < xs_t->delayed_nu_fission.shape()[2]; d++) {
for (int d = 0; d < xs_t->delayed_nu_fission.shape()[1]; d++) {
val += xs_t->delayed_nu_fission(a, d, gin);
}
}