crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Updated test results with Legendre moments for consistent scattering matrices

Browse source
This commit is contained in:
Will Boyd 2017-03-02 09:49:54 -05:00
parent 1dc199e2f4
commit 6d092292bc
13 changed files with 1538 additions and 1569 deletions
Download patch file Download diff file Expand all files Collapse all files

23
openmc/mgxs/mgxs.py
View file

@ -4574,10 +4574,11 @@ class ScatterProbabilityMatrix(MatrixMGXS):
.. math::
\langle \sigma_{s,g'\rightarrow g} \phi \rangle &= \int_{r \in V} dr
\langle \sigma_{s,\ell,g'\rightarrow g} \phi \rangle &= \int_{r \in V} dr
\int_{4\pi} d\Omega' \int_{E_{g'}}^{E_{g'-1}} dE' \int_{4\pi} d\Omega
\int_{E_g}^{E_{g-1}} dE \; \sigma_s (r, E'
\rightarrow E, \Omega' \cdot \Omega) \psi(r, E', \Omega')\\
\int_{E_g}^{E_{g-1}} dE \; P_\ell (\Omega \cdot \Omega')
\sigma_{s} (r, E' \rightarrow E, \Omega' \cdot \Omega) \psi(r, E',
\Omega')\\
\langle \sigma_{s,g'} \phi \rangle &= \int_{r \in V} dr
\int_{4\pi} d\Omega' \int_{E_{g'}}^{E_{g'-1}} dE' \int_{4\pi} d\Omega
\int_{0}^{\infty} dE \; \sigma_s (r, E'
@ -4586,11 +4587,19 @@ class ScatterProbabilityMatrix(MatrixMGXS):
\sigma_{s,g'\rightarrow g} \phi \rangle}{\langle
\sigma_{s,g'} \phi \rangle}
.. math::
This class can also be used to gather a prompt-nu-fission cross section
(which only includes the contributions from prompt neutrons). This is
accomplished by setting the :attr:`NuFissionMatrixXS.prompt` attribute to
`True`.
\langle \sigma_{s,\ell,g'\rightarrow g} \phi \rangle &= \int_{r \in V} dr
\int_{4\pi} d\Omega' \int_{E_{g'}}^{E_{g'-1}} dE' \int_{4\pi} d\Omega
\int_{E_g}^{E_{g-1}} dE \; P_\ell (\Omega \cdot \Omega') \sigma_s (r, E'
\rightarrow E, \Omega' \cdot \Omega) \psi(r, E', \Omega')\\
\langle \phi \rangle &= \int_{r \in V} dr \int_{4\pi} d\Omega
\int_{E_g}^{E_{g-1}} dE \; \psi (r, E, \Omega) \\
\sigma_{s,\ell,g'\rightarrow g} &= \frac{\langle
\sigma_{s,\ell,g'\rightarrow g} \phi \rangle}{\langle \phi \rangle}
To incorporate the effect of neutron multiplication from (n,xn) reactions
in the above probaility matrix, the `nu` parameter can be set to `True`.
Parameters
----------

676
tests/test_mgxs_library_condense/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

42
tests/test_mgxs_library_condense/results_true.dat
View file

@ -36,10 +36,16 @@
0 10000 1 1 total 1.0 0.066111
material group in group out nuclide mean std. dev.
0 10000 1 1 total 1.0 0.066111
material group in group out nuclide mean std. dev.
0 10000 1 1 total 0.341361 0.031842
material group in group out nuclide mean std. dev.
0 10000 1 1 total 0.343079 0.043498
material group in group out nuclide moment mean std. dev.
0 10000 1 1 total P0 0.388721 0.031279
1 10000 1 1 total P1 0.046155 0.006407
2 10000 1 1 total P2 0.017957 0.003039
3 10000 1 1 total P3 0.006618 0.002480
material group in group out nuclide moment mean std. dev.
0 10000 1 1 total P0 0.389304 0.043096
1 10000 1 1 total P1 0.046224 0.007316
2 10000 1 1 total P2 0.017984 0.003336
3 10000 1 1 total P3 0.006628 0.002534
material group out nuclide mean std. dev.
0 10000 1 total 1.0 0.046071
material group out nuclide mean std. dev.
@ -123,10 +129,16 @@
0 10001 1 1 total 1.0 0.095039
material group in group out nuclide mean std. dev.
0 10001 1 1 total 1.0 0.095039
material group in group out nuclide mean std. dev.
0 10001 1 1 total 0.276654 0.033287
material group in group out nuclide mean std. dev.
0 10001 1 1 total 0.27737 0.051266
material group in group out nuclide moment mean std. dev.
0 10001 1 1 total P0 0.309384 0.032376
1 10001 1 1 total P1 0.030756 0.007617
2 10001 1 1 total P2 0.018997 0.004420
3 10001 1 1 total P3 0.006263 0.003364
material group in group out nuclide moment mean std. dev.
0 10001 1 1 total P0 0.307987 0.050720
1 10001 1 1 total P1 0.030617 0.008524
2 10001 1 1 total P2 0.018911 0.005015
3 10001 1 1 total P3 0.006235 0.003442
material group out nuclide mean std. dev.
0 10001 1 total 0.0 0.0
material group out nuclide mean std. dev.
@ -210,10 +222,16 @@
0 10002 1 1 total 1.0 0.056867
material group in group out nuclide mean std. dev.
0 10002 1 1 total 1.0 0.056867
material group in group out nuclide mean std. dev.
0 10002 1 1 total 0.489924 0.069059
material group in group out nuclide mean std. dev.
0 10002 1 1 total 0.492997 0.086399
material group in group out nuclide moment mean std. dev.
0 10002 1 1 total P0 0.898938 0.067118
1 10002 1 1 total P1 0.408384 0.028127
2 10002 1 1 total P2 0.142591 0.010824
3 10002 1 1 total P3 0.008696 0.003588
material group in group out nuclide moment mean std. dev.
0 10002 1 1 total P0 0.903415 0.084918
1 10002 1 1 total P1 0.410417 0.036718
2 10002 1 1 total P2 0.143301 0.013612
3 10002 1 1 total P3 0.008739 0.003640
material group out nuclide mean std. dev.
0 10002 1 total 0.0 0.0
material group out nuclide mean std. dev.

110
tests/test_mgxs_library_distribcell/inputs_true.dat
View file

@ -304,7 +304,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>scatter-P0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10032">
@ -312,7 +312,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>nu-scatter-P0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10033">
@ -334,7 +334,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10036">
@ -355,23 +355,24 @@
</tally>
<tally id="10038">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter-1</scores>
<scores>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10039">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10040">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>nu-scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10041">
@ -387,51 +388,50 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10043">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10044">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter-1</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10045">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10046">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10047">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10048">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>inverse-velocity</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10049">
<filter bins="10000" type="distribcell" />
@ -444,7 +444,7 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>inverse-velocity</scores>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10051">
@ -452,54 +452,40 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10052">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10053">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10054">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10055">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10056">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10057">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10058">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -507,30 +493,30 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10059">
<tally id="10057">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10058">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10059">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10060">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10061">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10062">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
@ -538,14 +524,14 @@
<scores>decay-rate</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10063">
<tally id="10061">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10064">
<tally id="10062">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />

14
tests/test_mgxs_library_distribcell/results_true.dat
View file

@ -36,10 +36,16 @@
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total 1.0 0.037213
sum(distribcell) group in group out nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total 1.0 0.037208
sum(distribcell) group in group out nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total 0.342636 0.017656
sum(distribcell) group in group out nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total 0.340468 0.025369
sum(distribcell) group in group out nuclide moment mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P0 0.390797 0.016955
1 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P1 0.047641 0.005091
2 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P2 0.015866 0.003708
3 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P3 0.005430 0.003170
sum(distribcell) group in group out nuclide moment mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P0 0.387654 0.024885
1 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P1 0.047226 0.005527
2 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P2 0.015740 0.003750
3 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 1 total P3 0.005392 0.003157
sum(distribcell) group out nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 1.0 0.080455
sum(distribcell) group out nuclide mean std. dev.

676
tests/test_mgxs_library_hdf5/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

84
tests/test_mgxs_library_hdf5/results_true.dat
View file

@ -71,15 +71,27 @@ domain=10000 type=nu-scatter probability matrix
[[7.82243018e-02 1.25560869e-03]
[2.24310192e-03 4.10531468e-02]]
domain=10000 type=consistent scatter matrix
[[3.32987194e-01 9.94653712e-04]
[8.87128136e-04 3.79756741e-01]]
[[3.73090623e-02 4.89318749e-04]
[8.89289900e-04 3.01139605e-02]]
[[[3.86422967e-01 5.21704775e-02 2.01849914e-02 9.53256688e-03]
[9.94653712e-04 -2.08433942e-04 -1.03964400e-04 2.35646553e-04]]
[[8.87128136e-04 -7.36559899e-04 4.73756321e-04 -1.64539748e-04]
[3.94772020e-01 1.58130798e-02 6.11300510e-03 -1.00731826e-02]]]
[[[3.66286904e-02 7.76748968e-03 3.13767806e-03 2.32683668e-03]
[4.89318749e-04 1.50458419e-04 1.85500930e-04 1.29783344e-04]]
[[8.89289900e-04 7.38354757e-04 4.74910776e-04 1.64940700e-04]
[2.98710064e-02 4.44330993e-03 1.01307463e-02 1.00367467e-02]]]
domain=10000 type=consistent nu-scatter matrix
[[3.32465997e-01 9.88930322e-04]
[9.24639842e-04 3.96145576e-01]]
[[5.08819059e-02 8.36973463e-04]
[1.60212263e-03 2.87181296e-02]]
[[[3.84199430e-01 5.18702806e-02 2.00688439e-02 9.47771503e-03]
[9.88930322e-04 -2.07234582e-04 -1.03366173e-04 2.34290606e-04]]
[[9.24639842e-04 -7.67704913e-04 4.93788836e-04 -1.71497217e-04]
[4.11464730e-01 1.64817268e-02 6.37149004e-03 -1.04991213e-02]]]
[[[5.04005137e-02 9.04259889e-03 3.61169756e-03 2.46795126e-03]
[8.36973463e-04 2.06752354e-04 1.97694759e-04 2.06602832e-04]]
[[1.60212263e-03 1.33020162e-03 8.55587476e-04 2.97153075e-04]
[2.84588375e-02 4.60349295e-03 1.05573088e-02 1.04561824e-02]]]
domain=10000 type=chi
[1.00000000e+00 0.00000000e+00]
[4.60705493e-02 0.00000000e+00]
@ -257,15 +269,27 @@ domain=10001 type=nu-scatter probability matrix
[[1.08778697e-01 0.00000000e+00]
[0.00000000e+00 1.42427173e-01]]
domain=10001 type=consistent scatter matrix
[[2.70706971e-01 0.00000000e+00]
[0.00000000e+00 3.06542088e-01]]
[[3.82655052e-02 0.00000000e+00]
[0.00000000e+00 5.27278829e-02]]
[[[3.12162675e-01 3.84813069e-02 2.08815337e-02 8.01673640e-03]
[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
[[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[2.95421002e-01 -1.11817183e-02 8.81141444e-03 -3.26075959e-03]]]
[[[3.72534018e-02 8.74305413e-03 4.83468235e-03 3.77642682e-03]
[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
[[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[5.02358893e-02 1.61616720e-02 1.14951123e-02 7.31312479e-03]]]
domain=10001 type=consistent nu-scatter matrix
[[2.71891125e-01 0.00000000e+00]
[0.00000000e+00 3.07477884e-01]]
[[5.90732661e-02 0.00000000e+00]
[0.00000000e+00 7.57666234e-02]]
[[[3.10120713e-01 3.82295876e-02 2.07449405e-02 7.96429620e-03]
[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
[[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[2.96264249e-01 -1.12136353e-02 8.83656566e-03 -3.27006707e-03]]]
[[[5.84608636e-02 1.03230576e-02 5.67743887e-03 3.92760936e-03]
[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
[[0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[7.40187611e-02 1.63372534e-02 1.16408402e-02 7.35838382e-03]]]
domain=10001 type=chi
[0.00000000e+00 0.00000000e+00]
[0.00000000e+00 0.00000000e+00]
@ -443,15 +467,27 @@ domain=10002 type=nu-scatter probability matrix
[[3.60184962e-02 2.54736726e-03]
[2.18820864e-04 1.35884974e-01]]
domain=10002 type=consistent scatter matrix
[[2.51505773e-01 3.10225138e-02]
[4.40143026e-04 1.47923453e+00]]
[[4.18067487e-02 2.23201039e-03]
[4.44773843e-04 3.55391574e-01]]
[[[6.32859281e-01 3.76972649e-01 1.50714804e-01 9.04734705e-03]
[3.10225138e-02 8.66134326e-03 -2.53964096e-03 -3.74315061e-03]]
[[4.40143026e-04 3.97073448e-04 3.17256062e-04 2.12303394e-04]
[2.02025649e+00 5.06259696e-01 1.10372136e-01 2.48080660e-02]]]
[[[3.81421848e-02 2.37145043e-02 1.06635009e-02 3.86848985e-03]
[2.23201039e-03 9.99377011e-04 1.00968851e-03 8.26439590e-04]]
[[4.44773843e-04 4.01251123e-04 3.20593966e-04 2.14537073e-04]
[3.52193929e-01 7.91402819e-02 1.84875925e-02 8.77085752e-03]]]
domain=10002 type=consistent nu-scatter matrix
[[2.58651993e-01 3.13677176e-02]
[4.43343102e-04 1.48230049e+00]]
[[4.56422796e-02 2.98774961e-03]
[7.71125112e-04 4.71762821e-01]]
[[[6.39901439e-01 3.81167422e-01 1.52391887e-01 9.14802163e-03]
[3.13677176e-02 8.75772258e-03 -2.56790088e-03 -3.78480261e-03]]
[[4.43343102e-04 3.99960386e-04 3.19562684e-04 2.13846954e-04]
[2.03494484e+00 5.09940477e-01 1.11174601e-01 2.49884339e-02]]]
[[[4.26392749e-02 2.63117825e-02 1.16194647e-02 3.92016710e-03]
[2.98774961e-03 1.14887360e-03 1.03342895e-03 8.68386243e-04]]
[[7.71125112e-04 6.95667747e-04 5.55828678e-04 3.71952908e-04]
[4.68587349e-01 1.10634806e-01 2.50303502e-02 9.60120658e-03]]]
domain=10002 type=chi
[0.00000000e+00 0.00000000e+00]
[0.00000000e+00 0.00000000e+00]

110
tests/test_mgxs_library_mesh/inputs_true.dat
View file

@ -541,7 +541,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>scatter-P0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10032">
@ -549,7 +549,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>nu-scatter-P0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10033">
@ -571,7 +571,7 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10036">
@ -592,23 +592,24 @@
</tally>
<tally id="10038">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter-1</scores>
<scores>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10039">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10040">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>nu-scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10041">
@ -624,51 +625,50 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<scores>scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10043">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10044">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter-1</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10045">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10046">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10047">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10048">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>inverse-velocity</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10049">
<filter bins="1" type="mesh" />
@ -681,7 +681,7 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>inverse-velocity</scores>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10051">
@ -689,54 +689,40 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10052">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10053">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10054">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10055">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10056">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10057">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10058">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -744,30 +730,30 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10059">
<tally id="10057">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10058">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10059">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10060">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10061">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10062">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
@ -775,14 +761,14 @@
<scores>decay-rate</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10063">
<tally id="10061">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10064">
<tally id="10062">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />

48
tests/test_mgxs_library_mesh/results_true.dat
View file

@ -118,18 +118,42 @@
1 1 2 1 1 1 total 1.0 0.454973
2 2 1 1 1 1 total 1.0 0.146747
3 2 2 1 1 1 total 1.0 0.141824
mesh 1 group in group out nuclide mean std. dev.
x y z
0 1 1 1 1 1 total 0.385281 0.140163
1 1 2 1 1 1 total 0.381031 0.407858
2 2 1 1 1 1 total 0.429729 0.130705
3 2 2 1 1 1 total 0.373069 0.127649
mesh 1 group in group out nuclide mean std. dev.
x y z
0 1 1 1 1 1 total 0.391249 0.172830
1 1 2 1 1 1 total 0.366825 0.473935
2 2 1 1 1 1 total 0.436538 0.175458
3 2 2 1 1 1 total 0.380998 0.149007
mesh 1 group in group out nuclide moment mean std. dev.
x y z
0 1 1 1 1 1 total P0 0.633490 0.135514
1 1 1 1 1 1 total P1 0.245219 0.051009
2 1 1 1 1 1 total P2 0.091493 0.017479
3 1 1 1 1 1 total P3 0.013124 0.004906
4 1 2 1 1 1 total P0 0.618705 0.392783
5 1 2 1 1 1 total P1 0.232972 0.149703
6 1 2 1 1 1 total P2 0.073396 0.049002
7 1 2 1 1 1 total P3 -0.003195 0.010229
8 2 1 1 1 1 total P0 0.686150 0.126578
9 2 1 1 1 1 total P1 0.257252 0.047890
10 2 1 1 1 1 total P2 0.094522 0.018315
11 2 1 1 1 1 total P3 0.016676 0.005187
12 2 2 1 1 1 total P0 0.619269 0.124057
13 2 2 1 1 1 total P1 0.242153 0.048064
14 2 2 1 1 1 total P2 0.087677 0.018646
15 2 2 1 1 1 total P3 0.019785 0.014168
mesh 1 group in group out nuclide moment mean std. dev.
x y z
0 1 1 1 1 1 total P0 0.638348 0.169129
1 1 1 1 1 1 total P1 0.247099 0.064297
2 1 1 1 1 1 total P2 0.092195 0.022758
3 1 1 1 1 1 total P3 0.013224 0.005359
4 1 2 1 1 1 total P0 0.588378 0.462704
5 1 2 1 1 1 total P1 0.221552 0.175616
6 1 2 1 1 1 total P2 0.069798 0.056753
7 1 2 1 1 1 total P3 -0.003039 0.009829
8 2 1 1 1 1 total P0 0.698373 0.171880
9 2 1 1 1 1 total P1 0.261835 0.064773
10 2 1 1 1 1 total P2 0.096206 0.024355
11 2 1 1 1 1 total P3 0.016973 0.005960
12 2 2 1 1 1 total P0 0.625643 0.146236
13 2 2 1 1 1 total P1 0.244646 0.056797
14 2 2 1 1 1 total P2 0.088580 0.021648
15 2 2 1 1 1 total P3 0.019989 0.014515
mesh 1 group out nuclide mean std. dev.
x y z
0 1 1 1 1 total 1.0 0.135958

676
tests/test_mgxs_library_no_nuclides/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

132
tests/test_mgxs_library_no_nuclides/results_true.dat
View file

@ -82,16 +82,40 @@
2 10000 1 2 total 0.002567 0.001256
1 10000 2 1 total 0.002242 0.002243
0 10000 2 2 total 0.997758 0.041053
material group in group out nuclide mean std. dev.
3 10000 1 1 total 0.332987 0.037309
2 10000 1 2 total 0.000995 0.000489
1 10000 2 1 total 0.000887 0.000889
0 10000 2 2 total 0.379757 0.030114
material group in group out nuclide mean std. dev.
3 10000 1 1 total 0.332466 0.050882
2 10000 1 2 total 0.000989 0.000837
1 10000 2 1 total 0.000925 0.001602
0 10000 2 2 total 0.396146 0.028718
material group in group out nuclide moment mean std. dev.
12 10000 1 1 total P0 0.386423 0.036629
13 10000 1 1 total P1 0.052170 0.007767
14 10000 1 1 total P2 0.020185 0.003138
15 10000 1 1 total P3 0.009533 0.002327
8 10000 1 2 total P0 0.000995 0.000489
9 10000 1 2 total P1 -0.000208 0.000150
10 10000 1 2 total P2 -0.000104 0.000186
11 10000 1 2 total P3 0.000236 0.000130
4 10000 2 1 total P0 0.000887 0.000889
5 10000 2 1 total P1 -0.000737 0.000738
6 10000 2 1 total P2 0.000474 0.000475
7 10000 2 1 total P3 -0.000165 0.000165
0 10000 2 2 total P0 0.394772 0.029871
1 10000 2 2 total P1 0.015813 0.004443
2 10000 2 2 total P2 0.006113 0.010131
3 10000 2 2 total P3 -0.010073 0.010037
material group in group out nuclide moment mean std. dev.
12 10000 1 1 total P0 0.384199 0.050401
13 10000 1 1 total P1 0.051870 0.009043
14 10000 1 1 total P2 0.020069 0.003612
15 10000 1 1 total P3 0.009478 0.002468
8 10000 1 2 total P0 0.000989 0.000837
9 10000 1 2 total P1 -0.000207 0.000207
10 10000 1 2 total P2 -0.000103 0.000198
11 10000 1 2 total P3 0.000234 0.000207
4 10000 2 1 total P0 0.000925 0.001602
5 10000 2 1 total P1 -0.000768 0.001330
6 10000 2 1 total P2 0.000494 0.000856
7 10000 2 1 total P3 -0.000171 0.000297
0 10000 2 2 total P0 0.411465 0.028459
1 10000 2 2 total P1 0.016482 0.004603
2 10000 2 2 total P2 0.006371 0.010557
3 10000 2 2 total P3 -0.010499 0.010456
material group out nuclide mean std. dev.
1 10000 1 total 1.0 0.046071
0 10000 2 total 0.0 0.000000
@ -270,16 +294,40 @@
2 10001 1 2 total 0.0 0.000000
1 10001 2 1 total 0.0 0.000000
0 10001 2 2 total 1.0 0.142427
material group in group out nuclide mean std. dev.
3 10001 1 1 total 0.270707 0.038266
2 10001 1 2 total 0.000000 0.000000
1 10001 2 1 total 0.000000 0.000000
0 10001 2 2 total 0.306542 0.052728
material group in group out nuclide mean std. dev.
3 10001 1 1 total 0.271891 0.059073
2 10001 1 2 total 0.000000 0.000000
1 10001 2 1 total 0.000000 0.000000
0 10001 2 2 total 0.307478 0.075767
material group in group out nuclide moment mean std. dev.
12 10001 1 1 total P0 0.312163 0.037253
13 10001 1 1 total P1 0.038481 0.008743
14 10001 1 1 total P2 0.020882 0.004835
15 10001 1 1 total P3 0.008017 0.003776
8 10001 1 2 total P0 0.000000 0.000000
9 10001 1 2 total P1 0.000000 0.000000
10 10001 1 2 total P2 0.000000 0.000000
11 10001 1 2 total P3 0.000000 0.000000
4 10001 2 1 total P0 0.000000 0.000000
5 10001 2 1 total P1 0.000000 0.000000
6 10001 2 1 total P2 0.000000 0.000000
7 10001 2 1 total P3 0.000000 0.000000
0 10001 2 2 total P0 0.295421 0.050236
1 10001 2 2 total P1 -0.011182 0.016162
2 10001 2 2 total P2 0.008811 0.011495
3 10001 2 2 total P3 -0.003261 0.007313
material group in group out nuclide moment mean std. dev.
12 10001 1 1 total P0 0.310121 0.058461
13 10001 1 1 total P1 0.038230 0.010323
14 10001 1 1 total P2 0.020745 0.005677
15 10001 1 1 total P3 0.007964 0.003928
8 10001 1 2 total P0 0.000000 0.000000
9 10001 1 2 total P1 0.000000 0.000000
10 10001 1 2 total P2 0.000000 0.000000
11 10001 1 2 total P3 0.000000 0.000000
4 10001 2 1 total P0 0.000000 0.000000
5 10001 2 1 total P1 0.000000 0.000000
6 10001 2 1 total P2 0.000000 0.000000
7 10001 2 1 total P3 0.000000 0.000000
0 10001 2 2 total P0 0.296264 0.074019
1 10001 2 2 total P1 -0.011214 0.016337
2 10001 2 2 total P2 0.008837 0.011641
3 10001 2 2 total P3 -0.003270 0.007358
material group out nuclide mean std. dev.
1 10001 1 total 0.0 0.0
0 10001 2 total 0.0 0.0
@ -458,16 +506,40 @@
2 10002 1 2 total 0.046729 0.002547
1 10002 2 1 total 0.000218 0.000219
0 10002 2 2 total 0.999782 0.135885
material group in group out nuclide mean std. dev.
3 10002 1 1 total 0.251506 0.041807
2 10002 1 2 total 0.031023 0.002232
1 10002 2 1 total 0.000440 0.000445
0 10002 2 2 total 1.479235 0.355392
material group in group out nuclide mean std. dev.
3 10002 1 1 total 0.258652 0.045642
2 10002 1 2 total 0.031368 0.002988
1 10002 2 1 total 0.000443 0.000771
0 10002 2 2 total 1.482300 0.471763
material group in group out nuclide moment mean std. dev.
12 10002 1 1 total P0 0.632859 0.038142
13 10002 1 1 total P1 0.376973 0.023715
14 10002 1 1 total P2 0.150715 0.010664
15 10002 1 1 total P3 0.009047 0.003868
8 10002 1 2 total P0 0.031023 0.002232
9 10002 1 2 total P1 0.008661 0.000999
10 10002 1 2 total P2 -0.002540 0.001010
11 10002 1 2 total P3 -0.003743 0.000826
4 10002 2 1 total P0 0.000440 0.000445
5 10002 2 1 total P1 0.000397 0.000401
6 10002 2 1 total P2 0.000317 0.000321
7 10002 2 1 total P3 0.000212 0.000215
0 10002 2 2 total P0 2.020256 0.352194
1 10002 2 2 total P1 0.506260 0.079140
2 10002 2 2 total P2 0.110372 0.018488
3 10002 2 2 total P3 0.024808 0.008771
material group in group out nuclide moment mean std. dev.
12 10002 1 1 total P0 0.639901 0.042639
13 10002 1 1 total P1 0.381167 0.026312
14 10002 1 1 total P2 0.152392 0.011619
15 10002 1 1 total P3 0.009148 0.003920
8 10002 1 2 total P0 0.031368 0.002988
9 10002 1 2 total P1 0.008758 0.001149
10 10002 1 2 total P2 -0.002568 0.001033
11 10002 1 2 total P3 -0.003785 0.000868
4 10002 2 1 total P0 0.000443 0.000771
5 10002 2 1 total P1 0.000400 0.000696
6 10002 2 1 total P2 0.000320 0.000556
7 10002 2 1 total P3 0.000214 0.000372
0 10002 2 2 total P0 2.034945 0.468587
1 10002 2 2 total P1 0.509940 0.110635
2 10002 2 2 total P2 0.111175 0.025030
3 10002 2 2 total P3 0.024988 0.009601
material group out nuclide mean std. dev.
1 10002 1 total 0.0 0.0
0 10002 2 total 0.0 0.0

514
tests/test_mgxs_library_nuclides/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

2
tests/test_mgxs_library_nuclides/results_true.dat
View file

@ -1 +1 @@
5c57e011264efb348fa333166e42da34f564feb2c2507f89c72fadb1056ed871e251911800fe3976cac9c66fa456edb63ee3bc0c4a40d8d82abe2c58c6cc3128
98b6de0cd546baf457c878d9695806ed53d1af1f380e164d1803041947346a54830aef4a30eab1d261abce700a4f45e2f81d6e5f78c124110c0f60601161cf9b