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

Updated test results with new consistent transport-corrected total MGXS

Browse source
This commit is contained in:
Will Boyd 2017-03-07 09:23:12 -05:00
parent 3ff03261ca
commit 0352849b1a
14 changed files with 1858 additions and 1644 deletions
Download patch file Download diff file Expand all files Collapse all files

43
openmc/mgxs/mgxs.py
View file

@ -2715,7 +2715,10 @@ class TransportXS(MGXS):
super(TransportXS, self).__init__(domain, domain_type,
groups, by_nuclide, name, num_polar,
num_azimuthal)
self._estimator = 'analog'
# Use tracklength estimators for the total MGXS term, and
# analog estimators for the transport correction term
self._estimator = ['tracklength', 'tracklength', 'analog', 'analog']
self._valid_estimators = ['analog']
self.nu = nu
@ -2727,40 +2730,54 @@ class TransportXS(MGXS):
@property
def scores(self):
if not self.nu:
return ['flux', 'total', 'scatter-1']
return ['flux', 'total', 'flux', 'scatter-1']
else:
return ['flux', 'total', 'nu-scatter-1']
return ['flux', 'total', 'flux', 'nu-scatter-1']
@property
def tally_keys(self):
if not self.nu:
return super(TransportXS, self).tally_keys
else:
return ['flux', 'total', 'scatter-1']
return ['flux (tracklength)', 'total', 'flux (analog)', 'scatter-1']
@property
def filters(self):
group_edges = self.energy_groups.group_edges
energy_filter = openmc.EnergyFilter(group_edges)
energyout_filter = openmc.EnergyoutFilter(group_edges)
filters = [[energy_filter], [energy_filter], [energyout_filter]]
filters = [[energy_filter], [energy_filter],
[energy_filter], [energyout_filter]]
return self._add_angle_filters(filters)
@property
def rxn_rate_tally(self):
if self._rxn_rate_tally is None:
raise NotImplementedError('The reaction rate tally is poorly defined' \
' for the transport cross section')
@property
def xs_tally(self):
if self._xs_tally is None:
if self.tallies is None:
msg = 'Unable to get xs_tally since tallies have ' \
'not been loaded from a statepoint'
raise ValueError(msg)
# Switch EnergyoutFilter to EnergyFilter.
old_filt = self.tallies['scatter-1'].filters[-1]
new_filt = openmc.EnergyFilter(old_filt.bins)
new_filt.stride = old_filt.stride
self.tallies['scatter-1'].filters[-1] = new_filt
self._rxn_rate_tally = \
self.tallies['total'] - self.tallies['scatter-1']
self._rxn_rate_tally.sparse = self.sparse
# Compute total cross section
total_xs = self.tallies['total'] / self.tallies['flux (tracklength)']
return self._rxn_rate_tally
# Compute transport correction term
trans_corr = self.tallies['scatter-1'] / self.tallies['flux (analog)']
# Compute the transport-corrected total cross section
self._xs_tally = total_xs - trans_corr
self._compute_xs()
return self._xs_tally
@property
def nu(self):

1
src/xml/fox Submodule

@ -0,0 +1 @@
Subproject commit bdc852f4f43d969fb1b179cba79295c1e095a455

772
tests/test_mgxs_library_condense/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

16
tests/test_mgxs_library_condense/results_true.dat
View file

@ -1,9 +1,9 @@
material group in nuclide mean std. dev.
0 10000 1 total 0.453624 0.021053
material group in nuclide mean std. dev.
0 10000 1 total 0.400852 0.022858
material group in nuclide mean std. dev.
0 10000 1 total 0.400852 0.022858
material group in nuclide mean std. dev.
0 10000 1 total 0.4074 0.021863
material group in nuclide mean std. dev.
0 10000 1 total 0.4074 0.021863
material group in nuclide mean std. dev.
0 10000 1 total 0.064903 0.004313
material group in nuclide mean std. dev.
@ -92,9 +92,9 @@
material group in nuclide mean std. dev.
0 10001 1 total 0.311594 0.013793
material group in nuclide mean std. dev.
0 10001 1 total 0.279255 0.029189
0 10001 1 total 0.280977 0.015683
material group in nuclide mean std. dev.
0 10001 1 total 0.279255 0.029189
0 10001 1 total 0.280977 0.015683
material group in nuclide mean std. dev.
0 10001 1 total 0.00221 0.000286
material group in nuclide mean std. dev.
@ -183,9 +183,9 @@
material group in nuclide mean std. dev.
0 10002 1 total 0.904999 0.043964
material group in nuclide mean std. dev.
0 10002 1 total 0.499184 0.040914
0 10002 1 total 0.494581 0.046763
material group in nuclide mean std. dev.
0 10002 1 total 0.499184 0.040914
0 10002 1 total 0.494581 0.046763
material group in nuclide mean std. dev.
0 10002 1 total 0.00606 0.000555
material group in nuclide mean std. dev.

208
tests/test_mgxs_library_distribcell/inputs_true.dat
View file

@ -96,56 +96,56 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10003">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>total</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10004">
<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="10005">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter-1</scores>
<estimator>analog</estimator>
</tally>
<tally id="10005">
<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="10006">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>total</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10007">
<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>nu-scatter-1</scores>
<estimator>analog</estimator>
<scores>total</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10008">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10009">
<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>absorption</scores>
<estimator>tracklength</estimator>
<scores>nu-scatter-1</scores>
<estimator>analog</estimator>
</tally>
<tally id="10010">
<filter bins="10000" type="distribcell" />
@ -165,14 +165,14 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>fission</scores>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10013">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<scores>absorption</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10014">
@ -193,7 +193,7 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10017">
@ -207,7 +207,7 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>kappa-fission</scores>
<scores>nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10019">
@ -221,7 +221,7 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>kappa-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10021">
@ -229,14 +229,14 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10022">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10023">
<filter bins="10000" type="distribcell" />
@ -248,9 +248,8 @@
<tally id="10024">
<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-P3</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10025">
@ -265,15 +264,14 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter-P3</scores>
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10027">
<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>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10028">
@ -281,14 +279,15 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>nu-scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10029">
<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>flux</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10030">
@ -296,54 +295,47 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10031">
<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>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10032">
<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>flux</scores>
<estimator>tracklength</estimator>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10033">
<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>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10034">
<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-P3</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10035">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10036">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10037">
<tally id="10036">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -351,7 +343,29 @@
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10037">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10038">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10039">
<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-P3</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" />
@ -359,7 +373,7 @@
<scores>nu-scatter-0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10039">
<tally id="10041">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -367,47 +381,33 @@
<scores>scatter-0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10040">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10041">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10042">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10043">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10044">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10045">
<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>inverse-velocity</scores>
<estimator>tracklength</estimator>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10046">
<filter bins="10000" type="distribcell" />
@ -420,7 +420,7 @@
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>inverse-velocity</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10048">
@ -428,9 +428,23 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10049">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10050">
<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="10051">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -438,14 +452,14 @@
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10050">
<tally id="10052">
<filter bins="10000" type="distribcell" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10051">
<tally id="10053">
<filter bins="10000" type="distribcell" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
@ -453,7 +467,7 @@
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10052">
<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" />
@ -461,7 +475,7 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10053">
<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="energyout" />
@ -469,45 +483,45 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10054">
<tally id="10056">
<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="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>delayed-nu-fission</scores>
<estimator>tracklength</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>decay-rate</scores>
<scores>delayed-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>decay-rate</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10060">
<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="10059">
<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" />

4
tests/test_mgxs_library_distribcell/results_true.dat
View file

@ -1,9 +1,9 @@
sum(distribcell) group in nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.457353 0.010474
sum(distribcell) group in nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.405649 0.015784
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.410174 0.011573
sum(distribcell) group in nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.405641 0.015787
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.410166 0.011577
sum(distribcell) group in nuclide mean std. dev.
0 ((0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13... 1 total 0.066556 0.00251
sum(distribcell) group in nuclide mean std. dev.

772
tests/test_mgxs_library_hdf5/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

24
tests/test_mgxs_library_hdf5/results_true.dat
View file

@ -2,11 +2,11 @@ domain=10000 type=total
[4.14825464e-01 6.60169863e-01]
[2.27929105e-02 4.75188999e-02]
domain=10000 type=transport
[3.56859612e-01 6.47647614e-01]
[2.54935937e-02 2.37037335e-02]
[3.63092031e-01 6.44850709e-01]
[2.38384843e-02 4.76746408e-02]
domain=10000 type=nu-transport
[3.56859612e-01 6.47647614e-01]
[2.54935937e-02 2.37037335e-02]
[3.63092031e-01 6.44850709e-01]
[2.38384843e-02 4.76746408e-02]
domain=10000 type=absorption
[2.74078431e-02 2.64510714e-01]
[2.69249666e-03 2.33670618e-02]
@ -195,11 +195,11 @@ domain=10001 type=total
[3.13737666e-01 3.00821380e-01]
[1.55819223e-02 2.80524816e-02]
domain=10001 type=transport
[2.73227852e-01 3.12374814e-01]
[3.31153641e-02 4.96058281e-02]
[2.75508079e-01 3.12035015e-01]
[1.77418855e-02 3.23843473e-02]
domain=10001 type=nu-transport
[2.73227852e-01 3.12374814e-01]
[3.31153641e-02 4.96058281e-02]
[2.75508079e-01 3.12035015e-01]
[1.77418855e-02 3.23843473e-02]
domain=10001 type=absorption
[1.57499139e-03 5.40037826e-03]
[3.22547919e-04 6.18139027e-04]
@ -388,11 +388,11 @@ domain=10002 type=total
[6.64572195e-01 2.05238389e+00]
[3.12147473e-02 2.24342891e-01]
domain=10002 type=transport
[2.90565237e-01 1.51643790e+00]
[2.38518529e-02 2.35197252e-01]
[2.83322749e-01 1.49973953e+00]
[3.52061127e-02 2.30902118e-01]
domain=10002 type=nu-transport
[2.90565237e-01 1.51643790e+00]
[2.38518529e-02 2.35197252e-01]
[2.83322749e-01 1.49973953e+00]
[3.52061127e-02 2.30902118e-01]
domain=10002 type=absorption
[6.90399495e-04 3.16872549e-02]
[4.41475663e-05 3.74655831e-03]

208
tests/test_mgxs_library_mesh/inputs_true.dat
View file

@ -333,56 +333,56 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10003">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>total</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10004">
<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="10005">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter-1</scores>
<estimator>analog</estimator>
</tally>
<tally id="10005">
<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="10006">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>total</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10007">
<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>nu-scatter-1</scores>
<estimator>analog</estimator>
<scores>total</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10008">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10009">
<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>absorption</scores>
<estimator>tracklength</estimator>
<scores>nu-scatter-1</scores>
<estimator>analog</estimator>
</tally>
<tally id="10010">
<filter bins="1" type="mesh" />
@ -402,14 +402,14 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>fission</scores>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10013">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<scores>absorption</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10014">
@ -430,7 +430,7 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10017">
@ -444,7 +444,7 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>kappa-fission</scores>
<scores>nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10019">
@ -458,7 +458,7 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>kappa-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10021">
@ -466,14 +466,14 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10022">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10023">
<filter bins="1" type="mesh" />
@ -485,9 +485,8 @@
<tally id="10024">
<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-P3</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10025">
@ -502,15 +501,14 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-scatter-P3</scores>
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10027">
<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>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10028">
@ -518,14 +516,15 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<scores>nu-scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10029">
<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>flux</scores>
<scores>nu-scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10030">
@ -533,54 +532,47 @@
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<scores>scatter</scores>
<estimator>analog</estimator>
</tally>
<tally id="10031">
<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>flux</scores>
<estimator>analog</estimator>
</tally>
<tally id="10032">
<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>flux</scores>
<estimator>tracklength</estimator>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10033">
<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>
<estimator>tracklength</estimator>
<estimator>analog</estimator>
</tally>
<tally id="10034">
<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-P3</scores>
<estimator>analog</estimator>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10035">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10036">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10037">
<tally id="10036">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -588,7 +580,29 @@
<scores>scatter-P3</scores>
<estimator>analog</estimator>
</tally>
<tally id="10037">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10038">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>scatter</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10039">
<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-P3</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" />
@ -596,7 +610,7 @@
<scores>nu-scatter-0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10039">
<tally id="10041">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -604,47 +618,33 @@
<scores>scatter-0</scores>
<estimator>analog</estimator>
</tally>
<tally id="10040">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10041">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10042">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10043">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energyout" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10044">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10045">
<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>inverse-velocity</scores>
<estimator>tracklength</estimator>
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10046">
<filter bins="1" type="mesh" />
@ -657,7 +657,7 @@
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<scores>inverse-velocity</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10048">
@ -665,9 +665,23 @@
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>analog</estimator>
<estimator>tracklength</estimator>
</tally>
<tally id="10049">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>prompt-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10050">
<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="10051">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<filter bins="0.0 20000000.0" type="energyout" />
@ -675,14 +689,14 @@
<scores>prompt-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10050">
<tally id="10052">
<filter bins="1" type="mesh" />
<filter bins="0.0 20000000.0" type="energy" />
<nuclides>total</nuclides>
<scores>flux</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10051">
<tally id="10053">
<filter bins="1" type="mesh" />
<filter bins="1 2 3 4 5 6" type="delayedgroup" />
<filter bins="0.0 20000000.0" type="energy" />
@ -690,7 +704,7 @@
<scores>delayed-nu-fission</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10052">
<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" />
@ -698,7 +712,7 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10053">
<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="energyout" />
@ -706,45 +720,45 @@
<scores>delayed-nu-fission</scores>
<estimator>analog</estimator>
</tally>
<tally id="10054">
<tally id="10056">
<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="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>delayed-nu-fission</scores>
<estimator>tracklength</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>decay-rate</scores>
<scores>delayed-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>decay-rate</scores>
<estimator>tracklength</estimator>
</tally>
<tally id="10060">
<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="10059">
<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" />

16
tests/test_mgxs_library_mesh/results_true.dat
View file

@ -6,16 +6,16 @@
3 2 2 1 1 total 0.641095 0.091519
mesh 1 group in nuclide mean std. dev.
x y z
0 1 1 1 1 total 0.413423 0.087250
1 1 2 1 1 total 0.392074 0.244272
2 2 1 1 1 total 0.458841 0.087921
3 2 2 1 1 total 0.403898 0.074343
0 1 1 1 1 total 0.407867 0.104648
1 1 2 1 1 total 0.417805 0.300173
2 2 1 1 1 total 0.451699 0.087229
3 2 2 1 1 total 0.396449 0.095884
mesh 1 group in nuclide mean std. dev.
x y z
0 1 1 1 1 total 0.413423 0.087250
1 1 2 1 1 total 0.392074 0.244272
2 2 1 1 1 total 0.458841 0.087921
3 2 2 1 1 total 0.403898 0.074343
0 1 1 1 1 total 0.407867 0.104648
1 1 2 1 1 total 0.417805 0.300173
2 2 1 1 1 total 0.451699 0.087229
3 2 2 1 1 total 0.396449 0.095884
mesh 1 group in nuclide mean std. dev.
x y z
0 1 1 1 1 total 0.021476 0.004248

772
tests/test_mgxs_library_no_nuclides/inputs_true.dat
View file

File diff suppressed because it is too large Load diff

24
tests/test_mgxs_library_no_nuclides/results_true.dat
View file

@ -2,11 +2,11 @@
1 10000 1 total 0.414825 0.022793
0 10000 2 total 0.660170 0.047519
material group in nuclide mean std. dev.
1 10000 1 total 0.356860 0.025494
0 10000 2 total 0.647648 0.023704
1 10000 1 total 0.363092 0.023838
0 10000 2 total 0.644851 0.047675
material group in nuclide mean std. dev.
1 10000 1 total 0.356860 0.025494
0 10000 2 total 0.647648 0.023704
1 10000 1 total 0.363092 0.023838
0 10000 2 total 0.644851 0.047675
material group in nuclide mean std. dev.
1 10000 1 total 0.027408 0.002692
0 10000 2 total 0.264511 0.023367
@ -209,11 +209,11 @@
1 10001 1 total 0.313738 0.015582
0 10001 2 total 0.300821 0.028052
material group in nuclide mean std. dev.
1 10001 1 total 0.273228 0.033115
0 10001 2 total 0.312375 0.049606
1 10001 1 total 0.275508 0.017742
0 10001 2 total 0.312035 0.032384
material group in nuclide mean std. dev.
1 10001 1 total 0.273228 0.033115
0 10001 2 total 0.312375 0.049606
1 10001 1 total 0.275508 0.017742
0 10001 2 total 0.312035 0.032384
material group in nuclide mean std. dev.
1 10001 1 total 0.001575 0.000323
0 10001 2 total 0.005400 0.000618
@ -416,11 +416,11 @@
1 10002 1 total 0.664572 0.031215
0 10002 2 total 2.052384 0.224343
material group in nuclide mean std. dev.
1 10002 1 total 0.290565 0.023852
0 10002 2 total 1.516438 0.235197
1 10002 1 total 0.283323 0.035206
0 10002 2 total 1.499740 0.230902
material group in nuclide mean std. dev.
1 10002 1 total 0.290565 0.023852
0 10002 2 total 1.516438 0.235197
1 10002 1 total 0.283323 0.035206
0 10002 2 total 1.499740 0.230902
material group in nuclide mean std. dev.
1 10002 1 total 0.000690 0.000044
0 10002 2 total 0.031687 0.003747

640
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 @@
5df54054ab37cd69b7096b9f785be57a83a8c3872d3e065ff605aa6de5688640a5551f99313a3e60a9ab8f902c1d4927630e1838db26a2a09d41e29338634b43
107576b21fa8ed72f71ac65866e4ad1100cc62df527f6f4e6bb8a6318f3694614b6ba1c72634b2e97474a75ea1b3112dcb9b4b36e58f0417e5b04aee5d3c7570