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

Merged with develop

Browse source
This commit is contained in:
Will Boyd 2015-09-05 17:59:12 -04:00
commit bcc9ded29b
138 changed files with 12158 additions and 8403 deletions
Download patch file Download diff file Expand all files Collapse all files

3
.gitignore vendored
View file

@ -60,3 +60,6 @@ data/nndc
# PyCharm project configuration files
.idea
.idea/*
# IPython notebook checkpoints
.ipynb_checkpoints

22
CMakeLists.txt
View file

@ -287,6 +287,22 @@ foreach(test ${TESTS})
WORKING_DIRECTORY ${TEST_PATH}
COMMAND $<TARGET_FILE:openmc> -p ${TEST_PATH})
elseif(${test} MATCHES "test_filter_distribcell")
# Add each case for distribcell tests
add_test(NAME ${TEST_NAME}_case-1
WORKING_DIRECTORY ${TEST_PATH}/case-1
COMMAND $<TARGET_FILE:openmc> ${TEST_PATH}/case-1)
add_test(NAME ${TEST_NAME}_case-2
WORKING_DIRECTORY ${TEST_PATH}/case-2
COMMAND $<TARGET_FILE:openmc> ${TEST_PATH}/case-2)
add_test(NAME ${TEST_NAME}_case-3
WORKING_DIRECTORY ${TEST_PATH}/case-3
COMMAND $<TARGET_FILE:openmc> ${TEST_PATH}/case-3)
add_test(NAME ${TEST_NAME}_case-4
WORKING_DIRECTORY ${TEST_PATH}/case-4
COMMAND $<TARGET_FILE:openmc> ${TEST_PATH}/case-4)
# If a restart test is encounted, need to run with -r and restart file(s)
elseif(${test} MATCHES "restart")
@ -299,9 +315,9 @@ foreach(test ${TESTS})
elseif(${test} MATCHES "test_sourcepoint_restart")
set(RESTART_FILE statepoint.07.h5 source.07.h5)
elseif(${test} MATCHES "test_particle_restart_eigval")
set(RESTART_FILE particle_12_616.h5)
set(RESTART_FILE particle_9_555.h5)
elseif(${test} MATCHES "test_particle_restart_fixed")
set(RESTART_FILE particle_7_6144.h5)
set(RESTART_FILE particle_7_928.h5)
else(${test} MATCHES "test_statepoint_restart")
message(FATAL_ERROR "Restart test ${test} not recognized")
endif(${test} MATCHES "test_statepoint_restart")
@ -314,7 +330,7 @@ foreach(test ${TESTS})
elseif(${test} MATCHES "test_sourcepoint_restart")
set(RESTART_FILE statepoint.07.binary source.07.binary)
elseif(${test} MATCHES "test_particle_restart_eigval")
set(RESTART_FILE particle_12_616.binary)
set(RESTART_FILE particle_9_555.binary)
elseif(${test} MATCHES "test_particle_restart_fixed")
set(RESTART_FILE particle_7_6144.binary)
else(${test} MATCHES "test_statepoint_restart")

7
docs/source/conf.py
View file

@ -28,7 +28,8 @@ extensions = ['sphinx.ext.autodoc',
'sphinx.ext.napoleon',
'sphinx.ext.pngmath',
'sphinxcontrib.tikz',
'sphinx_numfig']
'sphinx_numfig',
'notebook_sphinxext']
# Add any paths that contain templates here, relative to this directory.
templates_path = ['_templates']
@ -51,9 +52,9 @@ copyright = u'2011-2015, Massachusetts Institute of Technology'
# built documents.
#
# The short X.Y version.
version = "0.6"
version = "0.7"
# The full version, including alpha/beta/rc tags.
release = "0.6.2"
release = "0.7.0"
# The language for content autogenerated by Sphinx. Refer to documentation
# for a list of supported languages.

2
docs/source/developers.rst
View file

@ -11,7 +11,7 @@ Active development of the OpenMC Monte Carlo code is currently led by:
* `Nick Horelik <https://github.com/nhorelik>`_
* `Adam Nelson <https://github.com/nelsonag>`_
* `Jon Walsh <https://github.com/walshjon>`_
* `Sterling Harper <https://github.com/walshjon>`_
* `Sterling Harper <https://github.com/smharper>`_
* `Will Boyd <https://github.com/wbinventor>`_
* `Benoit Forget <http://web.mit.edu/nse/people/faculty/forget.html>`_
* `Kord Smith <http://web.mit.edu/nse/people/faculty/smith.html>`_

2
docs/source/index.rst
View file

@ -28,7 +28,7 @@ free to send a message to the User's Group `mailing list`_.
:maxdepth: 1
quickinstall
releasenotes/index
releasenotes
methods/index
usersguide/index
devguide/index

10
docs/source/methods/physics.rst
View file

@ -187,11 +187,11 @@ secondary photons from nuclear de-excitation are tracked in OpenMC.
------------------------
These types of reactions are just treated as inelastic scattering and as such
are subject to the same procedure as described in
:ref:`inelastic-scatter`. Rather than tracking multiple secondary neutrons, the
weight of the outgoing neutron is multiplied by the number of secondary
neutrons, e.g. for :math:`(n,2n)`, only one outgoing neutron is tracked but its
weight is doubled.
are subject to the same procedure as described in :ref:`inelastic-scatter`. For
reactions with integral multiplicity, e.g., :math:`(n,2n)`, an appropriate
number of secondary neutrons are created. For reactions that have a multiplicity
given as a function of the incoming neutron energy (which occasionally occurs
for MT=5), the weight of the outgoing neutron is multiplied by the multiplcity.
.. _fission:

2467
docs/source/pythonapi/examples/pandas-dataframes.ipynb Normal file
View file

File diff suppressed because one or more lines are too long

11
docs/source/pythonapi/examples/pandas-dataframes.rst Normal file
View file

@ -0,0 +1,11 @@
=================
Pandas Dataframes
=================
.. only:: html
.. notebook:: pandas-dataframes.ipynb
.. only:: latex
IPython notebooks must be viewed in the online HTML documentation.

1657
docs/source/pythonapi/examples/tally-arithmetic.ipynb Normal file
View file

File diff suppressed because it is too large Load diff

11
docs/source/pythonapi/examples/tally-arithmetic.rst Normal file
View file

@ -0,0 +1,11 @@
================
Tally Arithmetic
================
.. only:: html
.. notebook:: tally-arithmetic.ipynb
.. only:: latex
IPython notebooks must be viewed in the online HTML documentation.

53
docs/source/pythonapi/index.rst
View file

@ -4,29 +4,68 @@
Python API
==========
--------
Contents
--------
OpenMC includes a rich Python API that enables programmatic pre- and
post-processing. The easiest way to begin using the API is to take a look at the
example Jupyter_ notebooks provided. However, this assumes that you are already
familiar with Python and common third-party packages such as NumPy_. If you have
never programmed in Python before, there are many good tutorials available
online. We recommend going through the modules from Codecademy_ and/or the
`Scipy lectures`_. The full API documentation serves to provide more information
on a given module or class.
**Handling nuclear data:**
.. toctree::
:maxdepth: 1
ace
**Creating input files:**
.. toctree::
:maxdepth: 1
cmfd
element
executor
filter
geometry
material
mesh
nuclide
opencg_compatible
particle_restart
plots
settings
statepoint
summary
surface
tallies
trigger
universe
**Running OpenMC:**
.. toctree::
:maxdepth: 1
executor
**Post-processing:**
.. toctree::
:maxdepth: 1
particle_restart
statepoint
summary
tallies
**Example Jupyter Notebooks:**
.. toctree::
:maxdepth: 1
examples/pandas-dataframes
examples/tally-arithmetic
.. _Jupyter: https://jupyter.org/
.. _NumPy: http://www.numpy.org/
.. _Codecademy: https://www.codecademy.com/tracks/python
.. _Scipy lectures: https://scipy-lectures.github.io/

67
docs/source/releasenotes.rst Normal file
View file

@ -0,0 +1,67 @@
.. _releasenotes:
==============================
Release Notes for OpenMC 0.7.0
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Complete Python API
- Python 3 compatability for all scripts
- All scripts consistently named openmc-* and installed together
- New 'distribcell' tally filter for repeated cells
- Ability to specify outer lattice universe
- XML input validation utility (openmc-validate-xml)
- Support for hexagonal lattices
- Material union energy grid method
- Tally triggers
- Remove dependence on PETSc
- Significant OpenMP performance improvements
- Support for Fortran 2008 MPI interface
- Use of Travis CI for continuous integration
- Simplifications and improvements to test suite
---------
Bug Fixes
---------
- b5f712_: Fix bug in spherical harmonics tallies
- e6675b_: Ensure all constants are double precision
- 04e2c1_: Fix potential bug in sample_nuclide routine
- 6121d9_: Fix bugs related to particle track files
- 2f0e89_: Fixes for nuclide specification in tallies
.. _b5f712: https://github.com/mit-crpg/openmc/commit/b5f712
.. _e6675b: https://github.com/mit-crpg/openmc/commit/e6675b
.. _04e2c1: https://github.com/mit-crpg/openmc/commit/04e2c1
.. _6121d9: https://github.com/mit-crpg/openmc/commit/6121d9
.. _2f0e89: https://github.com/mit-crpg/openmc/commit/2f0e89
------------
Contributors
------------
This release contains new contributions from the following people:
- `Will Boyd <wbinventor@gmail.com>`_
- `Matt Ellis <mellis13@mit.edu>`_
- `Sterling Harper <sterlingmharper@mit.edu>`_
- `Bryan Herman <bherman@mit.edu>`_
- `Nicholas Horelik <nicholas.horelik@gmail.com>`_
- `Colin Josey <cjosey@mit.edu>`_
- `William Lyu <PaleNeutron@users.noreply.github.com>`_
- `Adam Nelson <nelsonag@umich.edu>`_
- `Paul Romano <paul.k.romano@gmail.com>`_
- `Anthony Scopatz <scopatz@gmail.com>`_
- `Jon Walsh <walshjon@mit.edu>`_

25
docs/source/releasenotes/index.rst
View file

@ -1,25 +0,0 @@
.. _releasenotes:
=============
Release Notes
=============
The release notes for OpenMC give a list of system requirements, new features,
bugs fixed, and known issues for each successive release.
.. toctree::
:maxdepth: 1
notes_0.6.2
notes_0.6.1
notes_0.6.0
notes_0.5.4
notes_0.5.3
notes_0.5.2
notes_0.5.1
notes_0.5.0
notes_0.4.4
notes_0.4.3
notes_0.4.2
notes_0.4.1
notes_0.4.0

36
docs/source/releasenotes/notes_0.4.0.rst
View file

@ -1,36 +0,0 @@
.. _notes_0.4.0:
==============================
Release Notes for OpenMC 0.4.0
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions as well as
Mac OS X. However, it has not been tested yet on any releases of Microsoft
Windows. Memory requirements will vary depending on the size of the problem at
hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- The probability table method for treatment of energy self-shielding in the
unresolved resonance range has been implemented and is now turned on by
default.
- Calculation of Shannon entropy for assessing convergence of the fission source
distribution.
- Ability to compile with the PGI Fortran compiler.
- Ability to run on IBM BlueGene/P machines.
- Completely rewrote how nested universes are handled. Geometry is now much more
robust.
---------
Bug Fixes
---------
- Many geometry errors have been fixed. The Monte Carlo performance benchmark
can now be successfully run in OpenMC.

55
docs/source/releasenotes/notes_0.4.1.rst
View file

@ -1,55 +0,0 @@
.. _notes_0.4.1:
==============================
Release Notes for OpenMC 0.4.1
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions as well as
Mac OS X. However, it has not been tested yet on any releases of Microsoft
Windows. Memory requirements will vary depending on the size of the problem at
hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- A batching method has been implemented so that statistics can be calculated
based on multiple generations instead of a single generation. This can help to
overcome problems with underpredicted variance in problems where there is
correlation between successive fission source iterations.
- Users now have the option to select a non-unionized energy grid for problems
with many nuclides where the use of a unionized grid is not feasible.
- Improved plotting capability (Nick Horelik). The plotting input is now in
``plots.xml`` instead of ``plot.xml``.
- Added multiple estimators for k-effective and added a global tally for
leakage.
- Moved cross section-related output into cross_sections.out.
- Improved timing capabilities.
- Can now use more than 2**31 - 1 particles per generation.
- Improved fission bank synchronization method. This also necessitated changing
the source bank to be of type Bank rather than of type Particle.
- Added HDF5 output (not complete yet).
- Major changes to tally implementation.
---------
Bug Fixes
---------
- `b206a8`_: Fixed subtle error in the sampling of energy distributions.
- `800742`_: Fixed error in sampling of angle and rotating angles.
- `a07c08`_: Fixed bug in linear-linear interpolation during sampling energy.
- `a75283`_: Fixed energy and energyout tally filters to support many bins.
- `95cfac`_: Fixed error in cell neighbor searches.
- `83a803`_: Fixed bug related to probability tables.
.. _b206a8: https://github.com/mit-crpg/openmc/commit/b206a8
.. _800742: https://github.com/mit-crpg/openmc/commit/800742
.. _a07c08: https://github.com/mit-crpg/openmc/commit/a07c08
.. _a75283: https://github.com/mit-crpg/openmc/commit/a75283
.. _95cfac: https://github.com/mit-crpg/openmc/commit/95cfac
.. _83a803: https://github.com/mit-crpg/openmc/commit/83a803

56
docs/source/releasenotes/notes_0.4.2.rst
View file

@ -1,56 +0,0 @@
.. _notes_0.4.2:
==============================
Release Notes for OpenMC 0.4.2
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Ability to specify void materials.
- Option to not reduce tallies across processors at end of each batch.
- Uniform fission site method for reducing variance on local tallies.
- Reading/writing binary source files.
- Added more messages for <trace> or high verbosity.
- Estimator for diffusion coefficient.
- Ability to specify 'point' source type.
- Ability to change random number seed.
- Users can now specify units='sum' on a <density> tag. This tells the code that
the total material density is the sum of the atom fractions listed for each
nuclide on the material.
---------
Bug Fixes
---------
- a27f8f_: Fixed runtime error bug when using Intel compiler with DEBUG on.
- afe121_: Fixed minor bug in fission bank algorithms.
- e0968e_: Force re-evaluation of cross-sections when each particle is born.
- 298db8_: Fixed bug in surface currents when using energy-in filter.
- 2f3bbe_: Fixed subtle bug in S(a,b) cross section calculation.
- 671f30_: Fixed surface currents on mesh not encompassing geometry.
- b2c40e_: Fixed bug in incoming energy filter for track-length tallies.
- 5524fd_: Mesh filter now works with track-length tallies.
- d050c7_: Added Bessel's correction to make estimate of variance unbiased.
- 2a5b9c_: Fixed regression in plotting.
.. _a27f8f: https://github.com/mit-crpg/openmc/commit/a27f8f
.. _afe121: https://github.com/mit-crpg/openmc/commit/afe121
.. _e0968e: https://github.com/mit-crpg/openmc/commit/e0968e
.. _298db8: https://github.com/mit-crpg/openmc/commit/298db8
.. _2f3bbe: https://github.com/mit-crpg/openmc/commit/2f3bbe
.. _671f30: https://github.com/mit-crpg/openmc/commit/671f30
.. _b2c40e: https://github.com/mit-crpg/openmc/commit/b2c40e
.. _5524fd: https://github.com/mit-crpg/openmc/commit/5524fd
.. _d050c7: https://github.com/mit-crpg/openmc/commit/d050c7
.. _2a5b9c: https://github.com/mit-crpg/openmc/commit/2a5b9c

53
docs/source/releasenotes/notes_0.4.3.rst
View file

@ -1,53 +0,0 @@
.. _notes_0.4.3:
==============================
Release Notes for OpenMC 0.4.3
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Option to report confidence intervals for tally results.
- Rotation and translation for filled cells.
- Ability to explicitly specify <estimator> for tallies.
- Ability to store state points and use them to restart runs.
- Fixed source calculations (no subcritical multiplication however).
- Expanded options for external source distribution.
- Ability to tally reaction rates for individual nuclides within a material.
- Reduced memory usage by removing redundant storage or some cross-sections.
- 3bd35b_: Log-log interpolation for URR probability tables.
- Support to specify labels on tallies (nelsonag_).
---------
Bug Fixes
---------
- 33f29a_: Handle negative values in probability table.
- 1c472d_: Fixed survival biasing with probability tables.
- 3c6e80_: Fixed writing tallies with no filters.
- 460ef1_: Invalid results for duplicate tallies.
- 0069d5_: Fixed bug with 0 inactive batches.
- 7af2cf_: Fixed bug in score_analog_tallies.
- 85a60e_: Pick closest angular distribution for law 61.
- 3212f5_: Fixed issue with blank line at beginning of XML files.
.. _nelsonag: https://github.com/nelsonag
.. _33f29a: https://github.com/mit-crpg/openmc/commit/33f29a
.. _1c472d: https://github.com/mit-crpg/openmc/commit/1c472d
.. _3c6e80: https://github.com/mit-crpg/openmc/commit/3c6e80
.. _3bd35b: https://github.com/mit-crpg/openmc/commit/3bd35b
.. _0069d5: https://github.com/mit-crpg/openmc/commit/0069d5
.. _7af2cf: https://github.com/mit-crpg/openmc/commit/7af2cf
.. _460ef1: https://github.com/mit-crpg/openmc/commit/460ef1
.. _85a60e: https://github.com/mit-crpg/openmc/commit/85a60e
.. _3212f5: https://github.com/mit-crpg/openmc/commit/3212f5

45
docs/source/releasenotes/notes_0.4.4.rst
View file

@ -1,45 +0,0 @@
.. _notes_0.4.4:
==============================
Release Notes for OpenMC 0.4.4
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Ability to write state points when using <no_reduce>.
- Real-time XML validation in GNU Emacs with RELAX NG schemata.
- Writing state points every n batches with <state_point interval="..." />
- Suppress creation of summary.out and cross_sections.out by default with option
to turn them on with <output> tag in settings.xml file.
- Ability to create HDF5 state points.
- Binary source file is now part of state point file by default.
- Enhanced state point usage and added state point Python scripts.
- Turning confidence intervals on affects k-effective.
- Option to specify <upper_right> for tally meshes.
---------
Bug Fixes
---------
- 4654ee_: Fixed plotting with void cells.
- 7ee461_: Fixed bug with multi-line input using type='word'.
- 792eb3_: Fixed degrees of freedom for confidence intervals.
- 7fd617_: Fixed bug with restart runs in parallel.
- dc4a8f_: Fixed bug with fixed source restart runs.
.. _4654ee: https://github.com/mit-crpg/openmc/commit/4654ee
.. _7ee461: https://github.com/mit-crpg/openmc/commit/7ee461
.. _792eb3: https://github.com/mit-crpg/openmc/commit/792eb3
.. _7fd617: https://github.com/mit-crpg/openmc/commit/7fd617
.. _dc4a8f: https://github.com/mit-crpg/openmc/commit/dc4a8f

52
docs/source/releasenotes/notes_0.5.0.rst
View file

@ -1,52 +0,0 @@
.. _notes_0.5.0:
==============================
Release Notes for OpenMC 0.5.0
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- All user input options that formerly accepted "off" or "on" should now be
"false" or "true" (the proper XML schema datatype).
- The <criticality> element is deprecated and was replaced with <eigenvalue>.
- Added 'events' score that returns number of events that scored to a tally.
- Restructured tally filter implementation and user input.
- Source convergence acceleration via CMFD (implemented with PETSc).
- Ability to read source files in parallel when number of particles is greater
than that number of source sites.
- Cone surface types.
---------
Bug Fixes
---------
- 737b90_: Coincident surfaces from separate universes / particle traveling
tangent to a surface.
- a819b4_: Output of surface neighbors in summary.out file.
- b11696_: Reading long attribute lists in XML input.
- 2bd46a_: Search for tallying nuclides when no default_xs specified.
- 7a1f08_: Fix word wrapping when writing messages.
- c0e3ec_: Prevent underflow when compiling with MPI=yes and DEBUG=yes.
- 6f8d9d_: Set default tally labels.
- 6a3a5e_: Fix problem with corner-crossing in lattices.
.. _737b90: https://github.com/mit-crpg/openmc/commit/737b90
.. _a819b4: https://github.com/mit-crpg/openmc/commit/a819b4
.. _b11696: https://github.com/mit-crpg/openmc/commit/b11696
.. _2bd46a: https://github.com/mit-crpg/openmc/commit/2bd46a
.. _7a1f08: https://github.com/mit-crpg/openmc/commit/7a1f08
.. _c0e3ec: https://github.com/mit-crpg/openmc/commit/c0e3ec
.. _6f8d9d: https://github.com/mit-crpg/openmc/commit/6f8d9d
.. _6a3a5e: https://github.com/mit-crpg/openmc/commit/6a3a5e

45
docs/source/releasenotes/notes_0.5.1.rst
View file

@ -1,45 +0,0 @@
.. _notes_0.5.1:
==============================
Release Notes for OpenMC 0.5.1
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Absorption and combined estimators for k-effective.
- Natural elements can now be specified in materials using <element> rather than
<nuclide>.
- Support for multiple S(a,b) tables in a single material (e.g. BeO).
- Test suite using Python nosetests.
- Proper install capability with 'make install'.
- Lattices can now be 2 or 3 dimensions.
- New scatter-PN score type.
- New kappa-fission score type.
- Ability to tally any reaction by specifying MT.
---------
Bug Fixes
---------
- 94103e_: Two checks for outgoing energy filters.
- e77059_: Fix reaction name for MT=849.
- b0fe88_: Fix distance to surface for cones.
- 63bfd2_: Fix tracklength tallies with cell filter and universes.
- 88daf7_: Fix analog tallies with survival biasing.
.. _94103e: https://github.com/mit-crpg/openmc/commit/94103e
.. _e77059: https://github.com/mit-crpg/openmc/commit/e77059
.. _b0fe88: https://github.com/mit-crpg/openmc/commit/b0fe88
.. _63bfd2: https://github.com/mit-crpg/openmc/commit/63bfd2
.. _88daf7: https://github.com/mit-crpg/openmc/commit/88daf7

57
docs/source/releasenotes/notes_0.5.2.rst
View file

@ -1,57 +0,0 @@
.. _notes_0.5.2:
==============================
Release Notes for OpenMC 0.5.2
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Python script for mesh tally plotting
- Isotopic abundances based on IUPAC 2009 when using <element>
- Particle restart files for debugging
- Code will abort after certain number of lost particles (defaults to 10)
- Region outside lattice can be filled with material (void by default)
- 3D voxel plots
- Full HDF5/PHDF5 support (including support in statepoint.py)
- Cell overlap checking with -g command line flag (or when plotting)
---------
Bug Fixes
---------
- 7632f3_: Fixed bug in statepoint.py for multiple generations per batch.
- f85ac4_: Fix infinite loop bug in error module.
- 49c36b_: Don't convert surface ids if surface filter is for current tallies.
- 5ccc78_: Fix bug in reassignment of bins for mesh filter.
- b1f52f_: Fixed bug in plot color specification.
- eae7e5_: Fixed many memory leaks.
- 10c1cc_: Minor CMFD fixes.
- afdb50_: Add compatibility for XML comments without whitespace.
- a3c593_: Fixed bug in use of free gas scattering for H-1.
- 3a66e3_: Fixed bug in 2D mesh tally implementation.
- ab0793_: Corrected PETSC_NULL references to their correct types.
- 182ebd_: Use analog estimator with energyout filter.
.. _7632f3: https://github.com/mit-crpg/openmc/commit/7632f3
.. _f85ac4: https://github.com/mit-crpg/openmc/commit/f85ac4
.. _49c36b: https://github.com/mit-crpg/openmc/commit/49c36b
.. _5ccc78: https://github.com/mit-crpg/openmc/commit/5ccc78
.. _b1f52f: https://github.com/mit-crpg/openmc/commit/b1f52f
.. _eae7e5: https://github.com/mit-crpg/openmc/commit/eae7e5
.. _10c1cc: https://github.com/mit-crpg/openmc/commit/10c1cc
.. _afdb50: https://github.com/mit-crpg/openmc/commit/afdb50
.. _a3c593: https://github.com/mit-crpg/openmc/commit/a3c593
.. _3a66e3: https://github.com/mit-crpg/openmc/commit/3a66e3
.. _ab0793: https://github.com/mit-crpg/openmc/commit/ab0793
.. _182ebd: https://github.com/mit-crpg/openmc/commit/182ebd

49
docs/source/releasenotes/notes_0.5.3.rst
View file

@ -1,49 +0,0 @@
.. _notes_0.5.3:
==============================
Release Notes for OpenMC 0.5.3
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Output interface enhanced to allow multiple files handles to be opened
- Particle restart file linked to output interface
- Particle restarts and state point restarts are both identified with the -r
command line flag.
- Particle instance no longer global, passed to all physics routines
- Physics routines refactored to rely less on global memory, more arguments
passed in
- CMFD routines refactored and now can compute dominance ratio on the fly
- PETSc 3.4.2 or higher must be used and compiled with fortran datatype support
- Memory leaks fixed except for ones from xml-fortran package
- Test suite enhanced to test output with different compiler options
- Description of OpenMC development workflow added
- OpenMP shared-memory parallelism added
- Special run mode --tallies removed.
---------
Bug Fixes
---------
- 2b1e8a_: Normalize direction vector after reflecting particle.
- 5853d2_: Set blank default for cross section listing alias.
- e178c7_: Fix infinite loop with words greater than 80 characters in write_message.
- c18a6e_: Check for valid secondary mode on S(a,b) tables.
- 82c456_: Fix bug where last process could have zero particles.
.. _2b1e8a: https://github.com/mit-crpg/openmc/commit/2b1e8a
.. _5853d2: https://github.com/mit-crpg/openmc/commit/5853d2
.. _e178c7: https://github.com/mit-crpg/openmc/commit/e178c7
.. _c18a6e: https://github.com/mit-crpg/openmc/commit/c18a6e
.. _82c456: https://github.com/mit-crpg/openmc/commit/82c456

64
docs/source/releasenotes/notes_0.5.4.rst
View file

@ -1,64 +0,0 @@
.. _notes_0.5.4:
==============================
Release Notes for OpenMC 0.5.4
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Source sites outside geometry are resampled
- XML-Fortran backend replaced by FoX XML
- Ability to write particle track files
- Handle lost particles more gracefully (via particle track files)
- Multiple random number generator streams
- Mesh tally plotting utility converted to use Tkinter rather than PyQt
- Script added to download ACE data from NNDC
- Mixed ASCII/binary cross_sections.xml now allowed
- Expanded options for writing source bank
- Re-enabled ability to use source file as starting source
- S(a,b) recalculation avoided when same nuclide and S(a,b) table are accessed
---------
Bug Fixes
---------
- 32c03c_: Check for valid data in cross_sections.xml
- c71ef5_: Fix bug in statepoint.py
- 8884fb_: Check for all ZAIDs for S(a,b) tables
- b38af0_: Fix XML reading on multiple levels of input
- d28750_: Fix bug in convert_xsdir.py
- cf567c_: ENDF/B-VI data checked for compatibility
- 6b9461_: Fix p_valid sampling inside of sample_energy
.. _32c03c: https://github.com/mit-crpg/openmc/commit/32c03c
.. _c71ef5: https://github.com/mit-crpg/openmc/commit/c71ef5
.. _8884fb: https://github.com/mit-crpg/openmc/commit/8884fb
.. _b38af0: https://github.com/mit-crpg/openmc/commit/b38af0
.. _d28750: https://github.com/mit-crpg/openmc/commit/d28750
.. _cf567c: https://github.com/mit-crpg/openmc/commit/cf567c
.. _6b9461: https://github.com/mit-crpg/openmc/commit/6b9461
------------
Contributors
------------
This release contains new contributions from the following people:
- `Sterling Harper <smharper@mit.edu>`_
- `Bryan Herman <bherman@mit.edu>`_
- `Nick Horelik <nhorelik@mit.edu>`_
- `Adam Nelson <nelsonag@umich.edu>`_
- `Paul Romano <paul.k.romano@gmail.com>`_
- `Tuomas Viitanen <tuomas.viitanen@vtt.fi>`_
- `Jon Walsh <walshjon@mit.edu>`_

59
docs/source/releasenotes/notes_0.6.0.rst
View file

@ -1,59 +0,0 @@
.. _notes_0.6.0:
==============================
Release Notes for OpenMC 0.6.0
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Legendre and spherical harmonic expansion tally scores
- CMake is now default build system
- Regression test suite based on CTests and NNDC cross sections
- FoX is now a git submodule
- Support for older cross sections (e.g. MCNP 66c)
- Progress bar for plots
- Expanded support for natural elements via <natural_elements> in settings.xml
---------
Bug Fixes
---------
- 41f7ca_: Fixed erroneous results from survival biasing
- 038736_: Fix tallies over void materials
- 46f9e8_: Check for negative values in probability tables
- d1ca35_: Fixed sampling of angular distribution
- 0291c0_: Fixed indexing error in plotting
- d7a7d0_: Fix bug with <element> specifying xs attribute
- 85b3cb_: Fix out-of-bounds error with OpenMP threading
.. _41f7ca: https://github.com/mit-crpg/openmc/commit/41f7ca
.. _038736: https://github.com/mit-crpg/openmc/commit/038736
.. _46f9e8: https://github.com/mit-crpg/openmc/commit/46f9e8
.. _d1ca35: https://github.com/mit-crpg/openmc/commit/d1ca35
.. _0291c0: https://github.com/mit-crpg/openmc/commit/0291c0
.. _d7a7d0: https://github.com/mit-crpg/openmc/commit/d7a7d0
.. _85b3cb: https://github.com/mit-crpg/openmc/commit/85b3cb
------------
Contributors
------------
This release contains new contributions from the following people:
- `Sterling Harper <smharper@mit.edu>`_
- `Bryan Herman <bherman@mit.edu>`_
- `Nick Horelik <nhorelik@mit.edu>`_
- `Adam Nelson <nelsonag@umich.edu>`_
- `Paul Romano <paul.k.romano@gmail.com>`_
- `Jon Walsh <walshjon@mit.edu>`_

65
docs/source/releasenotes/notes_0.6.1.rst
View file

@ -1,65 +0,0 @@
.. _notes_0.6.1:
==============================
Release Notes for OpenMC 0.6.1
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Coarse mesh finite difference (CMFD) acceleration no longer requires PETSc
- Statepoint file numbering is now zero-padded
- Python scripts now compatible with Python 2 or 3
- Ability to run particle restarts in fixed source calculations
- Capability to filter box source by fissionable materials
- Nuclide/element names are now case insensitive in input files
- Improved treatment of resonance scattering for heavy nuclides
---------
Bug Fixes
---------
- 03e890_: Check for energy-dependent multiplicities in ACE files
- 4439de_: Fix distance-to-surface calculation for general plane surface
- 5808ed_: Account for differences in URR band probabilities at different energies
- 2e60c0_: Allow zero atom/weight percents in materials
- 3e0870_: Don't use PWD environment variable when setting path to input files
- dc4776_: Handle probability table resampling correctly
- 01178b_: Fix metastables nuclides in NNDC cross_sections.xml file
- 62ec43_: Don't read tallies.xml when OpenMC is run in plotting mode
- 2a95ef_: Prevent segmentation fault on "current" score without mesh filter
- 93e482_: Check for negative values in probability tables
.. _03e890: https://github.com/mit-crpg/openmc/commit/03e890
.. _4439de: https://github.com/mit-crpg/openmc/commit/4439de
.. _5808ed: https://github.com/mit-crpg/openmc/commit/5808ed
.. _2e60c0: https://github.com/mit-crpg/openmc/commit/2e60c0
.. _3e0870: https://github.com/mit-crpg/openmc/commit/3e0870
.. _dc4776: https://github.com/mit-crpg/openmc/commit/dc4776
.. _01178b: https://github.com/mit-crpg/openmc/commit/01178b
.. _62ec43: https://github.com/mit-crpg/openmc/commit/62ec43
.. _2a95ef: https://github.com/mit-crpg/openmc/commit/2a95ef
.. _93e482: https://github.com/mit-crpg/openmc/commit/93e482
------------
Contributors
------------
This release contains new contributions from the following people:
- `Sterling Harper <smharper@mit.edu>`_
- `Bryan Herman <bherman@mit.edu>`_
- `Adam Nelson <nelsonag@umich.edu>`_
- `Paul Romano <paul.k.romano@gmail.com>`_
- `Jon Walsh <walshjon@mit.edu>`_
- `Will Boyd <wbinventor@gmail.com>`_

58
docs/source/releasenotes/notes_0.6.2.rst
View file

@ -1,58 +0,0 @@
.. _notes_0.6.2:
==============================
Release Notes for OpenMC 0.6.2
==============================
-------------------
System Requirements
-------------------
There are no special requirements for running the OpenMC code. As of this
release, OpenMC has been tested on a variety of Linux distributions, Mac OS X,
and Microsoft Windows 7. Memory requirements will vary depending on the size of
the problem at hand (mostly on the number of nuclides in the problem).
------------
New Features
------------
- Meshline plotting capability
- Support for plotting cells/materials on middle universe levels
- Ability to model cells with no surfaces
- Compatibility with PETSc 3.5
- Compatability with OpenMPI 1.7/1.8
- Improved overall performance via logarithmic-mapped energy grid search
- Improved multi-threaded performance with atomic operations
- Support for fixed source problems with fissionable materials
---------
Bug Fixes
---------
- 26fb93_: Fix problem with -t, --track command-line flag
- 2f07c0_: Improved evaporation spectrum algorithm
- e6abb9_: Fix segfault when tallying in a void material
- 291b45_: Handle metastable nuclides in NNDC data and multiplicities in MT=5 data
.. _26fb93: https://github.com/mit-crpg/openmc/commit/26fb93
.. _2f07c0: https://github.com/mit-crpg/openmc/commit/2f07c0
.. _e6abb9: https://github.com/mit-crpg/openmc/commit/e6abb9
.. _291b45: https://github.com/mit-crpg/openmc/commit/291b45
------------
Contributors
------------
This release contains new contributions from the following people:
- `Will Boyd <wbinventor@gmail.com>`_
- `Matt Ellis <mellis13@mit.edu>`_
- `Sterling Harper <smharper@mit.edu>`_
- `Bryan Herman <bherman@mit.edu>`_
- `Nicholas Horelik <nicholas.horelik@gmail.com>`_
- `Anton Leontiev <bunder@t-25.ru>`_
- `Adam Nelson <nelsonag@umich.edu>`_
- `Paul Romano <paul.k.romano@gmail.com>`_
- `Jon Walsh <walshjon@mit.edu>`_
- `John Xia <john.danger.xia@gmail.com>`_

68
docs/sphinxext/LICENSE Normal file
View file

@ -0,0 +1,68 @@
The file notebook_sphinxext.py was derived from code in PyNE and yt.
PyNE has the following license:
-------------------------------------------------------------------------------
Copyright 2011-2015, the PyNE Development Team. All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are
permitted provided that the following conditions are met:
1. Redistributions of source code must retain the above copyright notice, this list of
conditions and the following disclaimer.
2. Redistributions in binary form must reproduce the above copyright notice, this list
of conditions and the following disclaimer in the documentation and/or other materials
provided with the distribution.
THIS SOFTWARE IS PROVIDED BY THE PYNE DEVELOPMENT TEAM ``AS IS'' AND ANY EXPRESS OR IMPLIED
WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND
FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL <COPYRIGHT HOLDER> OR
CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
The views and conclusions contained in the software and documentation are those of the
authors and should not be interpreted as representing official policies, either expressed
or implied, of the stakeholders of the PyNE project or the employers of PyNE developers.
-------------------------------------------------------------------------------
yt has the following license:
-------------------------------------------------------------------------------
yt is licensed under the terms of the Modified BSD License (also known as New
or Revised BSD), as follows:
Copyright (c) 2013-, yt Development Team
Copyright (c) 2006-2013, Matthew Turk <matthewturk@gmail.com>
All rights reserved.
Redistribution and use in source and binary forms, with or without
modification, are permitted provided that the following conditions are met:
Redistributions of source code must retain the above copyright notice, this
list of conditions and the following disclaimer.
Redistributions in binary form must reproduce the above copyright notice, this
list of conditions and the following disclaimer in the documentation and/or
other materials provided with the distribution.
Neither the name of the yt Development Team nor the names of its
contributors may be used to endorse or promote products derived from this
software without specific prior written permission.
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE
FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-------------------------------------------------------------------------------

120
docs/sphinxext/notebook_sphinxext.py Normal file
View file

@ -0,0 +1,120 @@
import sys
import os.path
import re
import time
from docutils import io, nodes, statemachine, utils
try:
from docutils.utils.error_reporting import ErrorString # the new way
except ImportError:
from docutils.error_reporting import ErrorString # the old way
from docutils.parsers.rst import Directive, convert_directive_function
from docutils.parsers.rst import directives, roles, states
from docutils.parsers.rst.roles import set_classes
from docutils.transforms import misc
try:
from IPython.nbconver.exporters import html
except ImportError:
from IPython.nbconvert import html
class Notebook(Directive):
"""Use nbconvert to insert a notebook into the environment.
This is based on the Raw directive in docutils
"""
required_arguments = 1
optional_arguments = 0
final_argument_whitespace = True
option_spec = {}
has_content = False
def run(self):
# check if raw html is supported
if not self.state.document.settings.raw_enabled:
raise self.warning('"%s" directive disabled.' % self.name)
# set up encoding
attributes = {'format': 'html'}
encoding = self.options.get(
'encoding', self.state.document.settings.input_encoding)
e_handler = self.state.document.settings.input_encoding_error_handler
# get path to notebook
source_dir = os.path.dirname(
os.path.abspath(self.state.document.current_source))
nb_path = os.path.normpath(os.path.join(source_dir,
self.arguments[0]))
nb_path = utils.relative_path(None, nb_path)
# convert notebook to html
exporter = html.HTMLExporter(template_file='full')
output, resources = exporter.from_filename(nb_path)
header = output.split('<head>', 1)[1].split('</head>',1)[0]
body = output.split('<body>', 1)[1].split('</body>',1)[0]
# add HTML5 scoped attribute to header style tags
header = header.replace('<style', '<style scoped="scoped"')
header = header.replace('body {\n overflow: visible;\n padding: 8px;\n}\n',
'')
header = header.replace("code,pre{", "code{")
# Filter out styles that conflict with the sphinx theme.
filter_strings = [
'navbar',
'body{',
'alert{',
'uneditable-input{',
'collapse{',
]
filter_strings.extend(['h%s{' % (i+1) for i in range(6)])
line_begin = [
'pre{',
'p{margin'
]
filterfunc = lambda x: not any([s in x for s in filter_strings])
header_lines = filter(filterfunc, header.split('\n'))
filterfunc = lambda x: not any([x.startswith(s) for s in line_begin])
header_lines = filter(filterfunc, header_lines)
header = '\n'.join(header_lines)
# concatenate raw html lines
lines = ['<div class="ipynotebook">']
lines.append(header)
lines.append(body)
lines.append('</div>')
text = '\n'.join(lines)
# add dependency
self.state.document.settings.record_dependencies.add(nb_path)
attributes['source'] = nb_path
# create notebook node
nb_node = notebook('', text, **attributes)
(nb_node.source, nb_node.line) = \
self.state_machine.get_source_and_line(self.lineno)
return [nb_node]
class notebook(nodes.raw):
pass
def visit_notebook_node(self, node):
self.visit_raw(node)
def depart_notebook_node(self, node):
self.depart_raw(node)
def setup(app):
app.add_node(notebook,
html=(visit_notebook_node, depart_notebook_node))
app.add_directive('notebook', Notebook)

8
openmc/clean_xml.py
View file

@ -82,11 +82,11 @@ def clean_xml_indentation(element, level=0):
if not element.tail or not element.tail.strip():
element.tail = i
for element in element:
clean_xml_indentation(element, level+1)
for sub_element in element:
clean_xml_indentation(sub_element, level+1)
if not element.tail or not element.tail.strip():
element.tail = i
if not sub_element.tail or not sub_element.tail.strip():
sub_element.tail = i
else:
if level and (not element.tail or not element.tail.strip()):

30
openmc/tallies.py
View file

@ -823,10 +823,10 @@ class Tally(object):
A list of filter type strings
(e.g., ['mesh', 'energy']; default is [])
filter_bins : list
filter_bins : list of Iterables
A list of the filter bins corresponding to the filter_types
parameter (e.g., [1, (0., 0.625e-6)]; default is []). Each bin in
the list is the integer ID for 'material', 'surface', 'cell',
parameter (e.g., [(1,), (0., 0.625e-6)]; default is []). Each bin
in the list is the integer ID for 'material', 'surface', 'cell',
'cellborn', and 'universe' Filters. Each bin is an integer for the
cell instance ID for 'distribcell Filters. Each bin is a 2-tuple of
floats for 'energy' and 'energyout' filters corresponding to the
@ -914,7 +914,8 @@ class Tally(object):
filter_indices[i][j] = filter_index
# Account for stride in each of the previous filters
filter_indices[:i] *= filter.num_bins
for indices in filter_indices[:i]:
indices *= filter.num_bins
# Apply outer product sum between all filter bin indices
filter_indices = list(map(sum, itertools.product(*filter_indices)))
@ -1456,7 +1457,7 @@ class Tally(object):
Raises
------
ValueError
When this method is called before the other tally is populated
When this method is called before the other tally is populated
with data by the StatePoint.read_results() method.
"""
@ -1470,6 +1471,7 @@ class Tally(object):
new_name = '({0} {1} {2})'.format(self.name, binary_op, other.name)
new_tally = Tally(name=new_name)
new_tally.with_batch_statistics = True
new_tally._derived = True
data = self._align_tally_data(other)
@ -1497,7 +1499,7 @@ class Tally(object):
new_tally._std_dev = np.abs(new_tally.mean) * \
np.sqrt(self_rel_err**2 + other_rel_err**2)
elif binary_op == '^':
data = self._align_tally_data(power)
data = self._align_tally_data(other)
mean_ratio = data['other']['mean'] / data['self']['mean']
first_term = mean_ratio * data['self']['std. dev.']
second_term = \
@ -1512,6 +1514,7 @@ class Tally(object):
new_tally.with_summary = self.with_summary
if self.num_realizations == other.num_realizations:
new_tally.num_realizations = self.num_realizations
new_tally.num_score_bins = self.num_score_bins * other.num_score_bins
# Generate filter "outer products"
if self.filters == other.filters:
@ -1671,6 +1674,7 @@ class Tally(object):
elif isinstance(other, Real):
new_tally = Tally(name='derived')
new_tally._derived = True
new_tally.with_batch_statistics = True
new_tally.name = self.name
new_tally._mean = self._mean + other
@ -1739,6 +1743,8 @@ class Tally(object):
new_tally = self._outer_product(other, binary_op='-')
elif isinstance(other, Real):
new_tally = Tally(name='derived')
new_tally._derived = True
new_tally.name = self.name
new_tally._mean = self._mean - other
new_tally._std_dev = self._std_dev
@ -1807,6 +1813,8 @@ class Tally(object):
new_tally = self._outer_product(other, binary_op='*')
elif isinstance(other, Real):
new_tally = Tally(name='derived')
new_tally._derived = True
new_tally.name = self.name
new_tally._mean = self._mean * other
new_tally._std_dev = self._std_dev * np.abs(other)
@ -1875,6 +1883,8 @@ class Tally(object):
new_tally = self._outer_product(other, binary_op='/')
elif isinstance(other, Real):
new_tally = Tally(name='derived')
new_tally._derived = True
new_tally.name = self.name
new_tally._mean = self._mean / other
new_tally._std_dev = self._std_dev * np.abs(1. / other)
@ -1943,6 +1953,8 @@ class Tally(object):
new_tally = self._outer_product(power, binary_op='^')
elif isinstance(power, Real):
new_tally = Tally(name='derived')
new_tally._derived = True
new_tally.name = self.name
new_tally._mean = self._mean ** power
self_rel_err = self.std_dev / self.mean
@ -2087,10 +2099,10 @@ class Tally(object):
A list of filter type strings
(e.g., ['mesh', 'energy']; default is [])
filter_bins : list
filter_bins : list of Iterables
A list of the filter bins corresponding to the filter_types
parameter (e.g., [1, (0., 0.625e-6)]; default is []). Each bin in
the list is the integer ID for 'material', 'surface', 'cell',
parameter (e.g., [(1,), (0., 0.625e-6)]; default is []). Each bin
in the list is the integer ID for 'material', 'surface', 'cell',
'cellborn', and 'universe' Filters. Each bin is an integer for the
cell instance ID for 'distribcell Filters. Each bin is a 2-tuple of
floats for 'energy' and 'energyout' filters corresponding to the

3
openmc/universe.py
View file

@ -714,7 +714,8 @@ class Lattice(object):
assert isinstance(u, Universe)
univs[u._id] = u
univs[self._outer._id] = self._outer
if self.outer is not None:
univs[self.outer._id] = self.outer
return univs

18
scripts/openmc-plot-mesh-tally
View file

@ -133,7 +133,11 @@ class MeshPlotter(tk.Frame):
self.mesh = selectedTally.filters_by_name['mesh'].mesh
# Get mesh dimensions
self.nx, self.ny, self.nz = self.mesh.dimension
if len(self.mesh.dimension) == 2:
self.nx, self.ny = self.mesh.dimension
self.nz = 1
else:
self.nx, self.ny, self.nz = self.mesh.dimension
# Repopulate comboboxes baesd on current basis selection
text = self.basisBox.get()
@ -210,7 +214,7 @@ class MeshPlotter(tk.Frame):
mesh_filter = f
continue
index = self.filterBoxes[f.type].current()
spec_list.append((f, index))
spec_list.append((f.type, (index,)))
text = self.basisBox.get()
if text == 'xy':
@ -229,11 +233,11 @@ class MeshPlotter(tk.Frame):
else:
meshtuple = (i + 1, axial_level, j + 1)
filters, filter_bins = zip(*spec_list + [
(mesh_filter, meshtuple)])
mean = selectedTally.get_value(
self.scoreBox.get(), filters, filter_bins)
stdev = selectedTally.get_value(
self.scoreBox.get(), filters, filter_bins,
(mesh_filter.type, (meshtuple,))])
mean = selectedTally.get_values(
[self.scoreBox.get()], filters, filter_bins)
stdev = selectedTally.get_values(
[self.scoreBox.get()], filters, filter_bins,
value='std_dev')
if mbvalue == 'Mean':
matrix[i, j] = mean

2
scripts/openmc-statepoint-3d
View file

@ -219,7 +219,7 @@ def main(file_, o):
for y in range(1,ny+1):
for z in range(1,nz+1):
filterspec[0][1] = (x,y,z)
val = sp.get_value(tally.id-1, filterspec, sid)[o.valerr]
val = sp.get_values(tally.id-1, filterspec, sid)[o.valerr]
if o.vtk:
# vtk cells go z, y, x, so we store it now and enter it later
i = (z-1)*nx*ny + (y-1)*nx + x-1

53
scripts/openmc-track-to-vtk
View file

@ -1,4 +1,4 @@
#!/usr/bin/env python2
#!/usr/bin/env python
"""Convert binary particle track to VTK poly data.
Usage information can be obtained by running 'track.py --help':
@ -64,26 +64,43 @@ def main():
for fname in args.input:
# Write coordinate values to points array.
if fname.endswith('.binary'):
track = open(fname, 'rb').read()
coords = [struct.unpack("ddd", track[24*i : 24*(i+1)])
for i in range(len(track)/24)]
n_points = len(coords)
for triplet in coords:
points.InsertNextPoint(triplet)
track = open(fname, 'rb')
# Determine number of particles and tracks/particle
n_particles = struct.unpack('i', track.read(4))[0]
n_coords = struct.unpack('i'*n_particles, track.read(4*n_particles))
coords = []
for i in range(n_particles):
# Read coordinates for each particle
coords.append([struct.unpack('ddd', track.read(24))
for j in range(n_coords[i])])
# Add coordinates to points data
for triplet in coords[i]:
points.InsertNextPoint(triplet)
else:
coords = h5py.File(fname).get('coordinates')
n_points = coords.shape[0]
for i in range(n_points):
points.InsertNextPoint(coords[i, :])
track = h5py.File(fname)
n_particles = track['n_particles'].value[0]
n_coords = track['n_coords']
coords = []
for i in range(n_particles):
coords.append(track['coordinates_' + str(i + 1)].value)
for j in range(n_coords[i]):
points.InsertNextPoint(coords[i][j,:])
# Create VTK line and assign points to line.
line = vtk.vtkPolyLine()
line.GetPointIds().SetNumberOfIds(n_points)
for i in range(n_points):
line.GetPointIds().SetId(i, point_offset+i)
for i in range(n_particles):
# Create VTK line and assign points to line.
line = vtk.vtkPolyLine()
line.GetPointIds().SetNumberOfIds(n_coords[i])
for j in range(n_coords[i]):
line.GetPointIds().SetId(j, point_offset + j)
# Add line to cell array
cells.InsertNextCell(line)
point_offset += n_coords[i]
cells.InsertNextCell(line)
point_offset += n_points
data = vtk.vtkPolyData()
data.SetPoints(points)
data.SetLines(cells)

2
setup.py
View file

@ -10,7 +10,7 @@ except ImportError:
have_setuptools = False
kwargs = {'name': 'openmc',
'version': '0.6.2',
'version': '0.7.0',
'packages': ['openmc', 'openmc.mgxs'],
'scripts': glob.glob('scripts/openmc-*'),

7
src/constants.F90
View file

@ -7,8 +7,8 @@ module constants
! OpenMC major, minor, and release numbers
integer, parameter :: VERSION_MAJOR = 0
integer, parameter :: VERSION_MINOR = 6
integer, parameter :: VERSION_RELEASE = 2
integer, parameter :: VERSION_MINOR = 7
integer, parameter :: VERSION_RELEASE = 0
! Revision numbers for binary files
integer, parameter :: REVISION_STATEPOINT = 13
@ -44,6 +44,9 @@ module constants
integer, parameter :: MAX_EVENTS = 10000
integer, parameter :: MAX_SAMPLE = 100000
! Maximum number of secondary particles created
integer, parameter :: MAX_SECONDARY = 1000
! Maximum number of words in a single line, length of line, and length of
! single word
integer, parameter :: MAX_WORDS = 500

288
src/eigenvalue.F90
View file

@ -4,273 +4,24 @@ module eigenvalue
use message_passing
#endif
use cmfd_execute, only: cmfd_init_batch, execute_cmfd
use constants, only: ZERO
use error, only: fatal_error, warning
use global
use math, only: t_percentile
use mesh, only: count_bank_sites
use mesh_header, only: StructuredMesh
use output, only: write_message, header, print_columns, &
print_batch_keff, print_generation
use particle_header, only: Particle
use random_lcg, only: prn, set_particle_seed, prn_skip
use search, only: binary_search
use source, only: get_source_particle
use state_point, only: write_state_point, write_source_point
use string, only: to_str
use tally, only: synchronize_tallies, setup_active_usertallies, &
reset_result
use trigger, only: check_triggers
use tracking, only: transport
implicit none
private
public :: run_eigenvalue
real(8) :: keff_generation ! Single-generation k on each
! processor
real(8) :: k_sum(2) = ZERO ! Used to reduce sum and sum_sq
real(8) :: keff_generation ! Single-generation k on each processor
real(8) :: k_sum(2) = ZERO ! Used to reduce sum and sum_sq
contains
!===============================================================================
! RUN_EIGENVALUE encompasses all the main logic where iterations are performed
! over the batches, generations, and histories in a k-eigenvalue calculation.
!===============================================================================
subroutine run_eigenvalue()
type(Particle) :: p
integer(8) :: i_work
if (master) call header("K EIGENVALUE SIMULATION", level=1)
! Display column titles
if(master) call print_columns()
! Turn on inactive timer
call time_inactive % start()
! ==========================================================================
! LOOP OVER BATCHES
BATCH_LOOP: do current_batch = 1, n_max_batches
call initialize_batch()
! Handle restart runs
if (restart_run .and. current_batch <= restart_batch) then
call replay_batch_history()
cycle BATCH_LOOP
end if
! =======================================================================
! LOOP OVER GENERATIONS
GENERATION_LOOP: do current_gen = 1, gen_per_batch
call initialize_generation()
! Start timer for transport
call time_transport % start()
! ====================================================================
! LOOP OVER PARTICLES
!$omp parallel do schedule(static) firstprivate(p)
PARTICLE_LOOP: do i_work = 1, work
current_work = i_work
! grab source particle from bank
call get_source_particle(p, current_work)
! transport particle
call transport(p)
end do PARTICLE_LOOP
!$omp end parallel do
! Accumulate time for transport
call time_transport % stop()
call finalize_generation()
end do GENERATION_LOOP
call finalize_batch()
if (satisfy_triggers) exit BATCH_LOOP
end do BATCH_LOOP
call time_active % stop()
! ==========================================================================
! END OF RUN WRAPUP
if (master) call header("SIMULATION FINISHED", level=1)
! Clear particle
call p % clear()
end subroutine run_eigenvalue
!===============================================================================
! INITIALIZE_BATCH
!===============================================================================
subroutine initialize_batch()
call write_message("Simulating batch " // trim(to_str(current_batch)) &
&// "...", 8)
! Reset total starting particle weight used for normalizing tallies
total_weight = ZERO
if (current_batch == n_inactive + 1) then
! Switch from inactive batch timer to active batch timer
call time_inactive % stop()
call time_active % start()
! Enable active batches (and tallies_on if it hasn't been enabled)
active_batches = .true.
tallies_on = .true.
! Add user tallies to active tallies list
!$omp parallel
call setup_active_usertallies()
!$omp end parallel
end if
! check CMFD initialize batch
if (cmfd_run) call cmfd_init_batch()
end subroutine initialize_batch
!===============================================================================
! INITIALIZE_GENERATION
!===============================================================================
subroutine initialize_generation()
! set overall generation number
overall_gen = gen_per_batch*(current_batch - 1) + current_gen
! Reset number of fission bank sites
n_bank = 0
! Count source sites if using uniform fission source weighting
if (ufs) call count_source_for_ufs()
! Store current value of tracklength k
keff_generation = global_tallies(K_TRACKLENGTH) % value
end subroutine initialize_generation
!===============================================================================
! FINALIZE_GENERATION
!===============================================================================
subroutine finalize_generation()
! Update global tallies with the omp private accumulation variables
!$omp parallel
!$omp critical
global_tallies(K_TRACKLENGTH) % value = &
global_tallies(K_TRACKLENGTH) % value + global_tally_tracklength
global_tallies(K_COLLISION) % value = &
global_tallies(K_COLLISION) % value + global_tally_collision
global_tallies(LEAKAGE) % value = &
global_tallies(LEAKAGE) % value + global_tally_leakage
global_tallies(K_ABSORPTION) % value = &
global_tallies(K_ABSORPTION) % value + global_tally_absorption
!$omp end critical
! reset private tallies
global_tally_tracklength = 0
global_tally_collision = 0
global_tally_leakage = 0
global_tally_absorption = 0
!$omp end parallel
#ifdef _OPENMP
! Join the fission bank from each thread into one global fission bank
call join_bank_from_threads()
#endif
! Distribute fission bank across processors evenly
call time_bank % start()
call synchronize_bank()
call time_bank % stop()
! Calculate shannon entropy
if (entropy_on) call shannon_entropy()
! Collect results and statistics
call calculate_generation_keff()
call calculate_average_keff()
! Write generation output
if (master .and. current_gen /= gen_per_batch) call print_generation()
end subroutine finalize_generation
!===============================================================================
! FINALIZE_BATCH handles synchronization and accumulation of tallies,
! calculation of Shannon entropy, getting single-batch estimate of keff, and
! turning on tallies when appropriate
!===============================================================================
subroutine finalize_batch()
! Collect tallies
call time_tallies % start()
call synchronize_tallies()
call time_tallies % stop()
! Reset global tally results
if (.not. active_batches) then
call reset_result(global_tallies)
n_realizations = 0
end if
! Perform CMFD calculation if on
if (cmfd_on) call execute_cmfd()
! Display output
if (master) call print_batch_keff()
! Calculate combined estimate of k-effective
if (master) call calculate_combined_keff()
! Check_triggers
if (master) call check_triggers()
#ifdef MPI
call MPI_BCAST(satisfy_triggers, 1, MPI_LOGICAL, 0, &
MPI_COMM_WORLD, mpi_err)
#endif
if (satisfy_triggers .or. &
(trigger_on .and. current_batch == n_max_batches)) then
call statepoint_batch % add(current_batch)
end if
! Write out state point if it's been specified for this batch
if (statepoint_batch % contains(current_batch)) then
call write_state_point()
end if
! Write out source point if it's been specified for this batch
if ((sourcepoint_batch % contains(current_batch) .or. source_latest) .and. &
source_write) then
call write_source_point()
end if
if (master .and. current_batch == n_max_batches) then
! Make sure combined estimate of k-effective is calculated at the last
! batch in case no state point is written
call calculate_combined_keff()
end if
end subroutine finalize_batch
!===============================================================================
! SYNCHRONIZE_BANK samples source sites from the fission sites that were
! accumulated during the generation. This routine is what allows this Monte
@ -830,40 +581,11 @@ contains
end subroutine count_source_for_ufs
!===============================================================================
! REPLAY_BATCH_HISTORY displays keff and entropy for each generation within a
! batch using data read from a state point file
!===============================================================================
subroutine replay_batch_history
! Write message at beginning
if (current_batch == 1) then
call write_message("Replaying history from state point...", 1)
end if
do current_gen = 1, gen_per_batch
overall_gen = overall_gen + 1
call calculate_average_keff()
! print out batch keff
if (current_gen < gen_per_batch) then
if (master) call print_generation()
else
if (master) call print_batch_keff()
end if
end do
! Write message at end
if (current_batch == restart_batch) then
call write_message("Resuming simulation...", 1)
end if
end subroutine replay_batch_history
#ifdef _OPENMP
!===============================================================================
! JOIN_BANK_FROM_THREADS
! JOIN_BANK_FROM_THREADS joins threadprivate fission banks into a single fission
! bank that can be sampled. Note that this operation is necessarily sequential
! to preserve the order of the bank when using varying numbers of threads.
!===============================================================================
subroutine join_bank_from_threads()

165
src/fixed_source.F90
View file

@ -1,165 +0,0 @@
module fixed_source
#ifdef MPI
use message_passing
#endif
use constants, only: ZERO, MAX_LINE_LEN
use global
use output, only: write_message, header
use particle_header, only: Particle
use random_lcg, only: set_particle_seed
use source, only: sample_external_source, copy_source_attributes
use state_point, only: write_state_point
use string, only: to_str
use tally, only: synchronize_tallies, setup_active_usertallies
use trigger, only: check_triggers
use tracking, only: transport
implicit none
contains
subroutine run_fixedsource()
integer(8) :: i ! index over histories in single cycle
type(Particle) :: p
if (master) call header("FIXED SOURCE TRANSPORT SIMULATION", level=1)
! Allocate particle and dummy source site
!$omp parallel
allocate(source_site)
!$omp end parallel
! Turn timer and tallies on
tallies_on = .true.
!$omp parallel
call setup_active_usertallies()
!$omp end parallel
call time_active % start()
! ==========================================================================
! LOOP OVER BATCHES
BATCH_LOOP: do current_batch = 1, n_max_batches
! In a restart run, skip any batches that have already been simulated
if (restart_run .and. current_batch <= restart_batch) then
if (current_batch > n_inactive) n_realizations = n_realizations + 1
cycle BATCH_LOOP
end if
call initialize_batch()
! Start timer for transport
call time_transport % start()
! =======================================================================
! LOOP OVER PARTICLES
!$omp parallel do schedule(static) firstprivate(p)
PARTICLE_LOOP: do i = 1, work
! Set unique particle ID
p % id = (current_batch - 1)*n_particles + work_index(rank) + i
! set particle trace
trace = .false.
if (current_batch == trace_batch .and. current_gen == trace_gen .and. &
work_index(rank) + i == trace_particle) trace = .true.
! set random number seed
call set_particle_seed(p % id)
! grab source particle from bank
call sample_source_particle(p)
! transport particle
call transport(p)
end do PARTICLE_LOOP
!$omp end parallel do
! Accumulate time for transport
call time_transport % stop()
call finalize_batch()
if (satisfy_triggers) exit BATCH_LOOP
end do BATCH_LOOP
call time_active % stop()
! ==========================================================================
! END OF RUN WRAPUP
if (master) call header("SIMULATION FINISHED", level=1)
end subroutine run_fixedsource
!===============================================================================
! INITIALIZE_BATCH
!===============================================================================
subroutine initialize_batch()
call write_message("Simulating batch " // trim(to_str(current_batch)) &
&// "...", 1)
! Reset total starting particle weight used for normalizing tallies
total_weight = ZERO
end subroutine initialize_batch
!===============================================================================
! FINALIZE_BATCH
!===============================================================================
subroutine finalize_batch()
! Collect and accumulate tallies
call time_tallies % start()
call synchronize_tallies()
call time_tallies % stop()
! Check_triggers
if (master) call check_triggers()
#ifdef MPI
call MPI_BCAST(satisfy_triggers, 1, MPI_LOGICAL, 0, &
MPI_COMM_WORLD, mpi_err)
#endif
if (satisfy_triggers .or. &
(trigger_on .and. current_batch == n_max_batches)) then
call statepoint_batch % add(current_batch)
end if
! Write out state point if it's been specified for this batch
if (statepoint_batch % contains(current_batch)) then
call write_state_point()
end if
end subroutine finalize_batch
!===============================================================================
! SAMPLE_SOURCE_PARTICLE
!===============================================================================
subroutine sample_source_particle(p)
type(Particle), intent(inout) :: p
! Set particle
call p % initialize()
! Sample the external source distribution
call sample_external_source(source_site)
! Copy source attributes to the particle
call copy_source_attributes(p, source_site)
! Determine whether to create track file
if (write_all_tracks) p % write_track = .true.
end subroutine sample_source_particle
end module fixed_source

4
src/global.F90
View file

@ -283,10 +283,6 @@ module global
! Mode to run in (fixed source, eigenvalue, plotting, etc)
integer :: run_mode = NONE
! Fixed source particle bank
type(Bank), pointer :: source_site => null()
!$omp threadprivate(source_site)
! Restart run
logical :: restart_run = .false.
integer :: restart_batch

45
src/initialize.F90
View file

@ -20,7 +20,6 @@ module initialize
write_message
use output_interface
use random_lcg, only: initialize_prng
use source, only: initialize_source
use state_point, only: load_state_point
use string, only: to_str, str_to_int, starts_with, ends_with
use tally_header, only: TallyObject, TallyResult, TallyFilter
@ -141,13 +140,9 @@ contains
! Determine how much work each processor should do
call calculate_work()
! Allocate banks and create source particles -- for a fixed source
! calculation, the external source distribution is sampled during the
! run, not at initialization
if (run_mode == MODE_EIGENVALUE) then
call allocate_banks()
if (.not. restart_run) call initialize_source()
end if
! Allocate source bank, and for eigenvalue simulations also allocate the
! fission bank
call allocate_banks()
! If this is a restart run, load the state point data and binary source
! file
@ -904,31 +899,33 @@ contains
call fatal_error("Failed to allocate source bank.")
end if
if (run_mode == MODE_EIGENVALUE) then
#ifdef _OPENMP
! If OpenMP is being used, each thread needs its own private fission
! bank. Since the private fission banks need to be combined at the end of a
! generation, there is also a 'master_fission_bank' that is used to collect
! the sites from each thread.
! If OpenMP is being used, each thread needs its own private fission
! bank. Since the private fission banks need to be combined at the end of
! a generation, there is also a 'master_fission_bank' that is used to
! collect the sites from each thread.
n_threads = omp_get_max_threads()
n_threads = omp_get_max_threads()
!$omp parallel
thread_id = omp_get_thread_num()
thread_id = omp_get_thread_num()
if (thread_id == 0) then
allocate(fission_bank(3*work))
else
allocate(fission_bank(3*work/n_threads))
end if
if (thread_id == 0) then
allocate(fission_bank(3*work))
else
allocate(fission_bank(3*work/n_threads))
end if
!$omp end parallel
allocate(master_fission_bank(3*work), STAT=alloc_err)
allocate(master_fission_bank(3*work), STAT=alloc_err)
#else
allocate(fission_bank(3*work), STAT=alloc_err)
allocate(fission_bank(3*work), STAT=alloc_err)
#endif
! Check for allocation errors
if (alloc_err /= 0) then
call fatal_error("Failed to allocate fission bank.")
! Check for allocation errors
if (alloc_err /= 0) then
call fatal_error("Failed to allocate fission bank.")
end if
end if
end subroutine allocate_banks

9
src/main.F90
View file

@ -1,13 +1,12 @@
program main
use constants
use eigenvalue, only: run_eigenvalue
use finalize, only: finalize_run
use fixed_source, only: run_fixedsource
use global
use initialize, only: initialize_run
use particle_restart, only: run_particle_restart
use plot, only: run_plot
use simulation, only: run_simulation
implicit none
@ -16,10 +15,8 @@ program main
! start problem based on mode
select case (run_mode)
case (MODE_FIXEDSOURCE)
call run_fixedsource()
case (MODE_EIGENVALUE)
call run_eigenvalue()
case (MODE_FIXEDSOURCE, MODE_EIGENVALUE)
call run_simulation()
case (MODE_PLOTTING)
call run_plot()
case (MODE_PARTICLE)

14
src/output.F90
View file

@ -7,7 +7,7 @@ module output
use endf, only: reaction_name
use error, only: fatal_error, warning
use geometry_header, only: Cell, Universe, Surface, Lattice, RectLattice, &
&HexLattice, BASE_UNIVERSE
HexLattice, BASE_UNIVERSE
use global
use math, only: t_percentile
use mesh_header, only: StructuredMesh
@ -1547,11 +1547,15 @@ contains
write(ou,100) "Total time in simulation", time_inactive % elapsed + &
time_active % elapsed
write(ou,100) " Time in transport only", time_transport % elapsed
write(ou,100) " Time in inactive batches", time_inactive % elapsed
if (run_mode == MODE_EIGENVALUE) then
write(ou,100) " Time in inactive batches", time_inactive % elapsed
end if
write(ou,100) " Time in active batches", time_active % elapsed
write(ou,100) " Time synchronizing fission bank", time_bank % elapsed
write(ou,100) " Sampling source sites", time_bank_sample % elapsed
write(ou,100) " SEND/RECV source sites", time_bank_sendrecv % elapsed
if (run_mode == MODE_EIGENVALUE) then
write(ou,100) " Time synchronizing fission bank", time_bank % elapsed
write(ou,100) " Sampling source sites", time_bank_sample % elapsed
write(ou,100) " SEND/RECV source sites", time_bank_sendrecv % elapsed
end if
write(ou,100) " Time accumulating tallies", time_tallies % elapsed
if (cmfd_run) write(ou,100) " Time in CMFD", time_cmfd % elapsed
if (cmfd_run) write(ou,100) " Building matrices", &

66
src/particle_header.F90
View file

@ -1,6 +1,8 @@
module particle_header
use constants, only: NEUTRON, ONE, NONE, ZERO
use bank_header, only: Bank
use constants, only: NEUTRON, ONE, NONE, ZERO, MAX_SECONDARY
use error, only: fatal_error
use geometry_header, only: BASE_UNIVERSE
implicit none
@ -79,9 +81,15 @@ module particle_header
! Track output
logical :: write_track = .false.
! Secondary particles created
integer :: n_secondary = 0
type(Bank) :: secondary_bank(MAX_SECONDARY)
contains
procedure :: initialize => initialize_particle
procedure :: clear => clear_particle
procedure :: initialize_from_source
procedure :: create_secondary
end type Particle
contains
@ -122,7 +130,7 @@ contains
end subroutine initialize_particle
!===============================================================================
! CLEAR_PARTICLE
! CLEAR_PARTICLE resets all coordinate levels for the particle
!===============================================================================
subroutine clear_particle(this)
@ -138,7 +146,7 @@ contains
end subroutine clear_particle
!===============================================================================
! RESET_COORD
! RESET_COORD clears data from a single coordinate level
!===============================================================================
elemental subroutine reset_coord(this)
@ -154,4 +162,56 @@ contains
end subroutine reset_coord
!===============================================================================
! INITIALIZE_FROM_SOURCE initializes a particle from data stored in a source
! site. The source site may have been produced from an external source, from
! fission, or simply as a secondary particle.
!===============================================================================
subroutine initialize_from_source(this, src)
class(Particle), intent(inout) :: this
type(Bank), intent(in) :: src
! set defaults
call this % initialize()
! copy attributes from source bank site
this % wgt = src % wgt
this % last_wgt = src % wgt
this % coord(1) % xyz = src % xyz
this % coord(1) % uvw = src % uvw
this % last_xyz = src % xyz
this % last_uvw = src % uvw
this % E = src % E
this % last_E = src % E
end subroutine initialize_from_source
!===============================================================================
! CREATE_SECONDARY stores the current phase space attributes of the particle in
! the secondary bank and increments the number of sites in the secondary bank.
!===============================================================================
subroutine create_secondary(this, uvw, type)
class(Particle), intent(inout) :: this
real(8), intent(in) :: uvw(3)
integer, intent(in) :: type
integer :: n
! Check to make sure that the hard-limit on secondary particles is not
! exceeded.
if (this % n_secondary == MAX_SECONDARY) then
call fatal_error("Too many secondary particles created.")
end if
n = this % n_secondary + 1
this % secondary_bank(n) % wgt = this % wgt
this % secondary_bank(n) % xyz(:) = this % coord(1) % xyz
this % secondary_bank(n) % uvw(:) = uvw
this % secondary_bank(n) % E = this % E
this % n_secondary = n
end subroutine create_secondary
end module particle_header

7
src/particle_restart_write.F90
View file

@ -43,12 +43,7 @@ contains
call pr % file_create(filename)
! Get information about source particle
select case (run_mode)
case (MODE_EIGENVALUE)
src => source_bank(current_work)
case (MODE_FIXEDSOURCE)
src => source_site
end select
src => source_bank(current_work)
! Write data to file
call pr % write_data(FILETYPE_PARTICLE_RESTART, 'filetype')

44
src/physics.F90
View file

@ -382,8 +382,7 @@ contains
end do
! Perform collision physics for inelastic scattering
call inelastic_scatter(nuc, rxn, p % E, p % coord(1) % uvw, &
p % mu, p % wgt)
call inelastic_scatter(nuc, rxn, p)
p % event_MT = rxn % MT
end if
@ -1268,24 +1267,23 @@ contains
! than fission), i.e. level scattering, (n,np), (n,na), etc.
!===============================================================================
subroutine inelastic_scatter(nuc, rxn, E, uvw, mu, wgt)
subroutine inelastic_scatter(nuc, rxn, p)
type(Nuclide), pointer :: nuc
type(Reaction), pointer :: rxn
type(Particle), intent(inout) :: p
type(Nuclide), pointer :: nuc
type(Reaction), pointer :: rxn
real(8), intent(inout) :: E ! energy in lab (incoming/outgoing)
real(8), intent(inout) :: uvw(3) ! directional cosines
real(8), intent(out) :: mu ! cosine of scattering angle in lab
real(8), intent(inout) :: wgt ! particle weight
integer :: law ! secondary energy distribution law
real(8) :: A ! atomic weight ratio of nuclide
real(8) :: E_in ! incoming energy
real(8) :: E_cm ! outgoing energy in center-of-mass
real(8) :: Q ! Q-value of reaction
real(8) :: yield ! neutron yield
integer :: i ! loop index
integer :: law ! secondary energy distribution law
real(8) :: E ! energy in lab (incoming/outgoing)
real(8) :: mu ! cosine of scattering angle in lab
real(8) :: A ! atomic weight ratio of nuclide
real(8) :: E_in ! incoming energy
real(8) :: E_cm ! outgoing energy in center-of-mass
real(8) :: Q ! Q-value of reaction
real(8) :: yield ! neutron yield
! copy energy of neutron
E_in = E
E_in = p % E
! determine A and Q
A = nuc % awr
@ -1319,16 +1317,22 @@ contains
mu = mu * sqrt(E_cm/E) + ONE/(A+ONE) * sqrt(E_in/E)
end if
! Set outgoing energy and scattering angle
p % E = E
p % mu = mu
! change direction of particle
uvw = rotate_angle(uvw, mu)
p % coord(1) % uvw = rotate_angle(p % coord(1) % uvw, mu)
! change weight of particle based on yield
if (rxn % multiplicity_with_E) then
yield = interpolate_tab1(rxn % multiplicity_E, E_in)
p % wgt = yield * p % wgt
else
yield = rxn % multiplicity
do i = 1, rxn % multiplicity - 1
call p % create_secondary(p % coord(1) % uvw, NEUTRON)
end do
end if
wgt = yield * wgt
end subroutine inelastic_scatter

371
src/simulation.F90 Normal file
View file

@ -0,0 +1,371 @@
module simulation
#ifdef MPI
use mpi
#endif
use cmfd_execute, only: cmfd_init_batch, execute_cmfd
use constants, only: ZERO
use eigenvalue, only: count_source_for_ufs, calculate_average_keff, &
calculate_combined_keff, calculate_generation_keff, &
shannon_entropy, synchronize_bank, keff_generation
#ifdef _OPENMP
use eigenvalue, only: join_bank_from_threads
#endif
use global
use output, only: write_message, header, print_columns, &
print_batch_keff, print_generation
use particle_header, only: Particle
use random_lcg, only: set_particle_seed
use source, only: initialize_source
use state_point, only: write_state_point, write_source_point
use string, only: to_str
use tally, only: synchronize_tallies, setup_active_usertallies, &
reset_result
use trigger, only: check_triggers
use tracking, only: transport
implicit none
private
public :: run_simulation
contains
!===============================================================================
! RUN_SIMULATION encompasses all the main logic where iterations are performed
! over the batches, generations, and histories in a fixed source or k-eigenvalue
! calculation.
!===============================================================================
subroutine run_simulation()
type(Particle) :: p
integer(8) :: i_work
if (.not. restart_run) call initialize_source()
! Display header
if (master) then
if (run_mode == MODE_FIXEDSOURCE) then
call header("FIXED SOURCE TRANSPORT SIMULATION", level=1)
elseif (run_mode == MODE_EIGENVALUE) then
call header("K EIGENVALUE SIMULATION", level=1)
call print_columns()
end if
end if
! Turn on inactive timer
call time_inactive % start()
! ==========================================================================
! LOOP OVER BATCHES
BATCH_LOOP: do current_batch = 1, n_max_batches
call initialize_batch()
! Handle restart runs
if (restart_run .and. current_batch <= restart_batch) then
call replay_batch_history()
cycle BATCH_LOOP
end if
! =======================================================================
! LOOP OVER GENERATIONS
GENERATION_LOOP: do current_gen = 1, gen_per_batch
call initialize_generation()
! Start timer for transport
call time_transport % start()
! ====================================================================
! LOOP OVER PARTICLES
!$omp parallel do schedule(static) firstprivate(p)
PARTICLE_LOOP: do i_work = 1, work
current_work = i_work
! grab source particle from bank
call initialize_history(p, current_work)
! transport particle
call transport(p)
end do PARTICLE_LOOP
!$omp end parallel do
! Accumulate time for transport
call time_transport % stop()
call finalize_generation()
end do GENERATION_LOOP
call finalize_batch()
if (satisfy_triggers) exit BATCH_LOOP
end do BATCH_LOOP
call time_active % stop()
! ==========================================================================
! END OF RUN WRAPUP
if (master) call header("SIMULATION FINISHED", level=1)
! Clear particle
call p % clear()
end subroutine run_simulation
!===============================================================================
! INITIALIZE_HISTORY
!===============================================================================
subroutine initialize_history(p, index_source)
type(Particle), intent(inout) :: p
integer(8), intent(in) :: index_source
integer(8) :: particle_seed ! unique index for particle
integer :: i
! set defaults
call p % initialize_from_source(source_bank(index_source))
! set identifier for particle
p % id = work_index(rank) + index_source
! set random number seed
particle_seed = (overall_gen - 1)*n_particles + p % id
call set_particle_seed(particle_seed)
! set particle trace
trace = .false.
if (current_batch == trace_batch .and. current_gen == trace_gen .and. &
p % id == trace_particle) trace = .true.
! Set particle track.
p % write_track = .false.
if (write_all_tracks) then
p % write_track = .true.
else if (allocated(track_identifiers)) then
do i=1, size(track_identifiers(1,:))
if (current_batch == track_identifiers(1,i) .and. &
&current_gen == track_identifiers(2,i) .and. &
&p % id == track_identifiers(3,i)) then
p % write_track = .true.
exit
end if
end do
end if
end subroutine initialize_history
!===============================================================================
! INITIALIZE_BATCH
!===============================================================================
subroutine initialize_batch()
if (run_mode == MODE_FIXEDSOURCE) then
call write_message("Simulating batch " // trim(to_str(current_batch)) &
// "...", 1)
end if
! Reset total starting particle weight used for normalizing tallies
total_weight = ZERO
if (current_batch == n_inactive + 1) then
! Switch from inactive batch timer to active batch timer
call time_inactive % stop()
call time_active % start()
! Enable active batches (and tallies_on if it hasn't been enabled)
active_batches = .true.
tallies_on = .true.
! Add user tallies to active tallies list
!$omp parallel
call setup_active_usertallies()
!$omp end parallel
end if
! check CMFD initialize batch
if (run_mode == MODE_EIGENVALUE) then
if (cmfd_run) call cmfd_init_batch()
end if
end subroutine initialize_batch
!===============================================================================
! INITIALIZE_GENERATION
!===============================================================================
subroutine initialize_generation()
! set overall generation number
overall_gen = gen_per_batch*(current_batch - 1) + current_gen
if (run_mode == MODE_EIGENVALUE) then
! Reset number of fission bank sites
n_bank = 0
! Count source sites if using uniform fission source weighting
if (ufs) call count_source_for_ufs()
! Store current value of tracklength k
keff_generation = global_tallies(K_TRACKLENGTH) % value
end if
end subroutine initialize_generation
!===============================================================================
! FINALIZE_GENERATION
!===============================================================================
subroutine finalize_generation()
! Update global tallies with the omp private accumulation variables
!$omp parallel
!$omp critical
if (run_mode == MODE_EIGENVALUE) then
global_tallies(K_COLLISION) % value = &
global_tallies(K_COLLISION) % value + global_tally_collision
global_tallies(K_ABSORPTION) % value = &
global_tallies(K_ABSORPTION) % value + global_tally_absorption
global_tallies(K_TRACKLENGTH) % value = &
global_tallies(K_TRACKLENGTH) % value + global_tally_tracklength
end if
global_tallies(LEAKAGE) % value = &
global_tallies(LEAKAGE) % value + global_tally_leakage
!$omp end critical
! reset private tallies
if (run_mode == MODE_EIGENVALUE) then
global_tally_collision = 0
global_tally_absorption = 0
global_tally_tracklength = 0
end if
global_tally_leakage = 0
!$omp end parallel
if (run_mode == MODE_EIGENVALUE) then
#ifdef _OPENMP
! Join the fission bank from each thread into one global fission bank
call join_bank_from_threads()
#endif
! Distribute fission bank across processors evenly
call time_bank % start()
call synchronize_bank()
call time_bank % stop()
! Calculate shannon entropy
if (entropy_on) call shannon_entropy()
! Collect results and statistics
call calculate_generation_keff()
call calculate_average_keff()
! Write generation output
if (master .and. current_gen /= gen_per_batch) call print_generation()
end if
end subroutine finalize_generation
!===============================================================================
! FINALIZE_BATCH handles synchronization and accumulation of tallies,
! calculation of Shannon entropy, getting single-batch estimate of keff, and
! turning on tallies when appropriate
!===============================================================================
subroutine finalize_batch()
! Collect tallies
call time_tallies % start()
call synchronize_tallies()
call time_tallies % stop()
! Reset global tally results
if (.not. active_batches) then
call reset_result(global_tallies)
n_realizations = 0
end if
if (run_mode == MODE_EIGENVALUE) then
! Perform CMFD calculation if on
if (cmfd_on) call execute_cmfd()
! Display output
if (master) call print_batch_keff()
! Calculate combined estimate of k-effective
if (master) call calculate_combined_keff()
end if
! Check_triggers
if (master) call check_triggers()
#ifdef MPI
call MPI_BCAST(satisfy_triggers, 1, MPI_LOGICAL, 0, &
MPI_COMM_WORLD, mpi_err)
#endif
if (satisfy_triggers .or. &
(trigger_on .and. current_batch == n_max_batches)) then
call statepoint_batch % add(current_batch)
end if
! Write out state point if it's been specified for this batch
if (statepoint_batch % contains(current_batch)) then
call write_state_point()
end if
! Write out source point if it's been specified for this batch
if ((sourcepoint_batch % contains(current_batch) .or. source_latest) .and. &
source_write) then
call write_source_point()
end if
if (master .and. current_batch == n_max_batches .and. &
run_mode == MODE_EIGENVALUE) then
! Make sure combined estimate of k-effective is calculated at the last
! batch in case no state point is written
call calculate_combined_keff()
end if
end subroutine finalize_batch
!===============================================================================
! REPLAY_BATCH_HISTORY displays keff and entropy for each generation within a
! batch using data read from a state point file
!===============================================================================
subroutine replay_batch_history
! Write message at beginning
if (current_batch == 1) then
call write_message("Replaying history from state point...", 1)
end if
if (run_mode == MODE_EIGENVALUE) then
do current_gen = 1, gen_per_batch
overall_gen = overall_gen + 1
call calculate_average_keff()
! print out batch keff
if (current_gen < gen_per_batch) then
if (master) call print_generation()
else
if (master) call print_batch_keff()
end if
end do
end if
! Write message at end
if (current_batch == restart_batch) then
call write_message("Resuming simulation...", 1)
end if
end subroutine replay_batch_history
end module simulation

75
src/source.F90
View file

@ -78,7 +78,7 @@ contains
! Write out initial source
if (write_initial_source) then
call write_message('Writing out initial source guess...', 1)
call write_message('Writing out initial source...', 1)
#ifdef HDF5
filename = trim(path_output) // 'initial_source.h5'
#else
@ -92,7 +92,8 @@ contains
end subroutine initialize_source
!===============================================================================
! SAMPLE_EXTERNAL_SOURCE
! SAMPLE_EXTERNAL_SOURCE samples the user-specified external source and stores
! the position, angle, and energy in a Bank type.
!===============================================================================
subroutine sample_external_source(site)
@ -245,74 +246,4 @@ contains
end subroutine sample_external_source
!===============================================================================
! GET_SOURCE_PARTICLE returns the next source particle
!===============================================================================
subroutine get_source_particle(p, index_source)
type(Particle), intent(inout) :: p
integer(8), intent(in) :: index_source
integer(8) :: particle_seed ! unique index for particle
integer :: i
type(Bank), pointer :: src
! set defaults
call p % initialize()
! Copy attributes from source to particle
src => source_bank(index_source)
call copy_source_attributes(p, src)
! set identifier for particle
p % id = work_index(rank) + index_source
! set random number seed
particle_seed = (overall_gen - 1)*n_particles + p % id
call set_particle_seed(particle_seed)
! set particle trace
trace = .false.
if (current_batch == trace_batch .and. current_gen == trace_gen .and. &
p % id == trace_particle) trace = .true.
! Set particle track.
p % write_track = .false.
if (write_all_tracks) then
p % write_track = .true.
else if (allocated(track_identifiers)) then
do i=1, size(track_identifiers(1,:))
if (current_batch == track_identifiers(1,i) .and. &
&current_gen == track_identifiers(2,i) .and. &
&p % id == track_identifiers(3,i)) then
p % write_track = .true.
exit
end if
end do
end if
end subroutine get_source_particle
!===============================================================================
! COPY_SOURCE_ATTRIBUTES
!===============================================================================
subroutine copy_source_attributes(p, src)
type(Particle), intent(inout) :: p
type(Bank), pointer :: src
! copy attributes from source bank site
p % wgt = src % wgt
p % last_wgt = src % wgt
p % coord(1) % xyz = src % xyz
p % coord(1) % uvw = src % uvw
p % last_xyz = src % xyz
p % last_uvw = src % uvw
p % E = src % E
p % last_E = src % E
end subroutine copy_source_attributes
end module source

80
src/track_output.F90
View file

@ -12,18 +12,22 @@ module track_output
implicit none
integer, private :: n_tracks ! total number of tracks
real(8), private, allocatable :: coords(:,:) ! track coordinates
!$omp threadprivate(n_tracks, coords)
type, private :: TrackCoordinates
real(8), allocatable :: coords(:,:)
end type TrackCoordinates
type(TrackCoordinates), private, allocatable :: tracks(:)
!$omp threadprivate(tracks)
contains
!===============================================================================
! INITIALIZE_PARTICLE_TRACK
! INITIALIZE_PARTICLE_TRACK allocates the array to store particle track
! information
!===============================================================================
subroutine initialize_particle_track()
n_tracks = 0
allocate(tracks(1))
end subroutine initialize_particle_track
!===============================================================================
@ -32,23 +36,50 @@ contains
subroutine write_particle_track(p)
type(Particle), intent(in) :: p
real(8), allocatable :: new_coords(:, :)
integer :: i
integer :: n_tracks
! Add another column to coords
n_tracks = n_tracks + 1
if (allocated(coords)) then
allocate(new_coords(3, n_tracks))
new_coords(:, 1:n_tracks-1) = coords
call move_alloc(FROM=new_coords, TO=coords)
i = size(tracks)
if (allocated(tracks(i) % coords)) then
n_tracks = size(tracks(i) % coords, 2)
allocate(new_coords(3, n_tracks + 1))
new_coords(:, 1:n_tracks) = tracks(i) % coords
call move_alloc(FROM=new_coords, TO=tracks(i) % coords)
else
allocate(coords(3,1))
n_tracks = 0
allocate(tracks(i) % coords(3, 1))
end if
! Write current coordinates into the newest column.
coords(:, n_tracks) = p % coord(1) % xyz
n_tracks = n_tracks + 1
tracks(i) % coords(:, n_tracks) = p % coord(1) % xyz
end subroutine write_particle_track
!===============================================================================
! ADD_PARTICLE_TRACK creates a new entry in the track coordinates for a
! secondary particle
!===============================================================================
subroutine add_particle_track()
type(TrackCoordinates), allocatable :: new_tracks(:)
integer :: i
! Determine current number of particle tracks
i = size(tracks)
! Create array one larger than current
allocate(new_tracks(i + 1))
! Copy memory and move allocation
new_tracks(1:i) = tracks(i)
call move_alloc(FROM=new_tracks, TO=tracks)
end subroutine add_particle_track
!===============================================================================
! FINALIZE_PARTICLE_TRACK writes the particle track array to disk.
!===============================================================================
@ -60,6 +91,10 @@ contains
character(MAX_FILE_LEN) :: fname
type(BinaryOutput) :: binout
integer :: i
integer, allocatable :: n_coords(:)
integer :: n_particle_tracks
#ifdef HDF5
fname = trim(path_output) // 'track_' // trim(to_str(current_batch)) &
// '_' // trim(to_str(current_gen)) // '_' // trim(to_str(p % id)) &
@ -69,13 +104,26 @@ contains
// '_' // trim(to_str(current_gen)) // '_' // trim(to_str(p % id)) &
// '.binary'
#endif
! Determine total number of particles and number of coordinates for each
n_particle_tracks = size(tracks)
allocate(n_coords(n_particle_tracks))
do i = 1, n_particle_tracks
n_coords(i) = size(tracks(i) % coords, 2)
end do
!$omp critical (FinalizeParticleTrack)
call binout % file_create(fname)
length = [3, n_tracks]
call binout % write_data(coords, 'coordinates', length=length)
call binout % write_data(n_particle_tracks, 'n_particles')
call binout % write_data(n_coords, 'n_coords', length=n_particle_tracks)
do i = 1, n_particle_tracks
length(:) = [3, n_coords(i)]
call binout % write_data(tracks(i) % coords, 'coordinates_' // &
trim(to_str(i)), length=length)
end do
call binout % file_close()
!$omp end critical (FinalizeParticleTrack)
deallocate(coords)
deallocate(tracks)
end subroutine finalize_particle_track
end module track_output

44
src/tracking.F90
View file

@ -15,7 +15,7 @@ module tracking
use tally, only: score_analog_tally, score_tracklength_tally, &
score_surface_current
use track_output, only: initialize_particle_track, write_particle_track, &
finalize_particle_track
add_particle_track, finalize_particle_track
implicit none
@ -45,20 +45,6 @@ contains
call write_message("Simulating Particle " // trim(to_str(p % id)))
end if
! If the cell hasn't been determined based on the particle's location,
! initiate a search for the current cell
if (p % coord(p % n_coord) % cell == NONE) then
call find_cell(p, found_cell)
! Particle couldn't be located
if (.not. found_cell) then
call fatal_error("Could not locate particle " // trim(to_str(p % id)))
end if
! set birth cell attribute
p % cell_born = p % coord(p % n_coord) % cell
end if
! Initialize number of events to zero
n_event = 0
@ -74,7 +60,19 @@ contains
call initialize_particle_track()
endif
do while (p % alive)
EVENT_LOOP: do
! If the cell hasn't been determined based on the particle's location,
! initiate a search for the current cell. This generally happens at the
! beginning of the history and again for any secondary particles
if (p % coord(p % n_coord) % cell == NONE) then
call find_cell(p, found_cell)
if (.not. found_cell) then
call fatal_error("Could not locate particle " // trim(to_str(p % id)))
end if
! set birth cell attribute
if (p % cell_born == NONE) p % cell_born = p % coord(p % n_coord) % cell
end if
! Write particle track.
if (p % write_track) call write_particle_track(p)
@ -197,7 +195,19 @@ contains
p % alive = .false.
end if
end do
! Check for secondary particles if this particle is dead
if (.not. p % alive) then
if (p % n_secondary > 0) then
call p % initialize_from_source(p % secondary_bank(p % n_secondary))
p % n_secondary = p % n_secondary - 1
! Enter new particle in particle track file
if (p % write_track) call add_particle_track()
else
exit EVENT_LOOP
end if
end if
end do EVENT_LOOP
! Finish particle track output.
if (p % write_track) then

2
tests/test_basic/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.011726E-01 1.841644E-03
3.021779E-01 3.813358E-03

500
tests/test_cmfd_feed/results_true.dat
View file

@ -1,128 +1,128 @@
k-combined:
1.177396E+00 4.883437E-03
1.172666E+00 8.502438E-03
tally 1:
1.107335E+01
1.235221E+01
2.002676E+01
4.017045E+01
2.844184E+01
8.111731E+01
3.428492E+01
1.177848E+02
3.735839E+01
1.397523E+02
3.894180E+01
1.517824E+02
3.528006E+01
1.246308E+02
2.863448E+01
8.222321E+01
2.125384E+01
4.535192E+01
1.112124E+01
1.241089E+01
1.170812E+01
1.376785E+01
2.179886E+01
4.765478E+01
2.945614E+01
8.709999E+01
3.527293E+01
1.245879E+02
3.829349E+01
1.470691E+02
3.709040E+01
1.379455E+02
3.380335E+01
1.145311E+02
2.801351E+01
7.871047E+01
2.029625E+01
4.131602E+01
1.084302E+01
1.180329E+01
tally 2:
2.243113E+01
2.535057E+01
1.557550E+01
1.222813E+01
2.164704E+00
2.389690E-01
4.049814E+01
8.218116E+01
2.854669E+01
4.085264E+01
3.934452E+00
7.833208E-01
5.706024E+01
1.633739E+02
4.049737E+01
8.234850E+01
5.231368E+00
1.386110E+00
6.798682E+01
2.320807E+02
4.819178E+01
1.166820E+02
6.247056E+00
1.968990E+00
7.449317E+01
2.782374E+02
5.315156E+01
1.417287E+02
6.667584E+00
2.250879E+00
7.530107E+01
2.849330E+02
5.378916E+01
1.454619E+02
7.071012E+00
2.522919E+00
6.820303E+01
2.335083E+02
4.850330E+01
1.181013E+02
6.241747E+00
1.973368E+00
5.831373E+01
1.704779E+02
4.149124E+01
8.633656E+01
5.385636E+00
1.468642E+00
4.184350E+01
8.792430E+01
2.971699E+01
4.435896E+01
3.784806E+00
7.343249E-01
2.202673E+01
2.444732E+01
1.531988E+01
1.183160E+01
2.073873E+00
2.272561E-01
2.270565E+01
2.599927E+01
1.590852E+01
1.276260E+01
2.252857E+00
2.614120E-01
4.313167E+01
9.326539E+01
3.044479E+01
4.648169E+01
4.023051E+00
8.172006E-01
5.859113E+01
1.725665E+02
4.171599E+01
8.755981E+01
5.512216E+00
1.531207E+00
6.892516E+01
2.383198E+02
4.904413E+01
1.207096E+02
6.542718E+00
2.155749E+00
7.421495E+01
2.764539E+02
5.288881E+01
1.405388E+02
6.811354E+00
2.358827E+00
7.278191E+01
2.661597E+02
5.169924E+01
1.343999E+02
6.516967E+00
2.148745E+00
6.655238E+01
2.222812E+02
4.729758E+01
1.123214E+02
6.102046E+00
1.890147E+00
5.708495E+01
1.636585E+02
4.068603E+01
8.317681E+01
5.394757E+00
1.465413E+00
4.136562E+01
8.598520E+01
2.958591E+01
4.402226E+01
3.765802E+00
7.200302E-01
2.275517E+01
2.614738E+01
1.589295E+01
1.276624E+01
2.232715E+00
2.558645E-01
tally 3:
1.500439E+01
1.135466E+01
9.942586E-01
5.051831E-02
2.744453E+01
3.776958E+01
1.781992E+00
1.612021E-01
3.901013E+01
7.642401E+01
2.472924E+00
3.077607E-01
4.642527E+01
1.082730E+02
2.924369E+00
4.335819E-01
5.109083E+01
1.309660E+02
3.325865E+00
5.592642E-01
5.182512E+01
1.350762E+02
3.259912E+00
5.350517E-01
4.672480E+01
1.095974E+02
3.097309E+00
4.854729E-01
3.997405E+01
8.014326E+01
2.589176E+00
3.374615E-01
2.861624E+01
4.114210E+01
1.806237E+00
1.656944E-01
1.474531E+01
1.096355E+01
9.807860E-01
4.934364E-02
1.529144E+01
1.179942E+01
1.023883E+00
5.386625E-02
2.936854E+01
4.326483E+01
1.881629E+00
1.788063E-01
4.015056E+01
8.114284E+01
2.594958E+00
3.407980E-01
4.720593E+01
1.118311E+02
3.161769E+00
5.053887E-01
5.095790E+01
1.304930E+02
3.308202E+00
5.528151E-01
4.979520E+01
1.246892E+02
3.163884E+00
5.062497E-01
4.554330E+01
1.041770E+02
3.019145E+00
4.618487E-01
3.921119E+01
7.727273E+01
2.472070E+00
3.099171E-01
2.843166E+01
4.067093E+01
1.823607E+00
1.688171E-01
1.530477E+01
1.184246E+01
1.047996E+00
5.549017E-02
tally 4:
0.000000E+00
0.000000E+00
@ -160,8 +160,8 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
3.050887E+00
4.698799E-01
3.111592E+00
4.883699E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -208,10 +208,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
5.388951E+00
1.456464E+00
2.664528E+00
3.577345E-01
5.536088E+00
1.540676E+00
2.727975E+00
3.757452E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -256,10 +256,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
7.165280E+00
2.573230E+00
4.930263E+00
1.218988E+00
7.518115E+00
2.840502E+00
5.271874E+00
1.398895E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -304,10 +304,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
8.550278E+00
3.668687E+00
6.985934E+00
2.450395E+00
8.764240E+00
3.855378E+00
7.176540E+00
2.591613E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -352,10 +352,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
9.236725E+00
4.281659E+00
8.429932E+00
3.565265E+00
9.381092E+00
4.414024E+00
8.597689E+00
3.710217E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -400,10 +400,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
9.396613E+00
4.431004E+00
9.307625E+00
4.351276E+00
9.158655E+00
4.215178E+00
9.188880E+00
4.244766E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -448,10 +448,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
8.721885E+00
3.821463E+00
9.438995E+00
4.479948E+00
8.362511E+00
3.509173E+00
9.159213E+00
4.209143E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -496,10 +496,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
7.096946E+00
2.525113E+00
8.605114E+00
3.710134E+00
7.029505E+00
2.479106E+00
8.613258E+00
3.719199E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -544,10 +544,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
5.222832E+00
1.373166E+00
7.480684E+00
2.804650E+00
5.119892E+00
1.320586E+00
7.401001E+00
2.749355E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -592,10 +592,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
2.726690E+00
3.773397E-01
5.470375E+00
1.504111E+00
2.765680E+00
3.903229E-01
5.461998E+00
1.501206E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -642,8 +642,8 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
3.026290E+00
4.597502E-01
3.044921E+00
4.656739E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -662,114 +662,114 @@ k cmfd
0.000000E+00
0.000000E+00
0.000000E+00
1.150583E+00
1.159578E+00
1.162798E+00
1.171003E+00
1.162732E+00
1.165591E+00
1.166032E+00
1.165387E+00
1.165451E+00
1.170726E+00
1.170820E+00
1.167945E+00
1.166050E+00
1.165757E+00
1.167631E+00
1.168366E+00
1.177990E+00
1.160010E+00
1.155990E+00
1.160167E+00
1.162166E+00
1.161566E+00
1.164454E+00
1.166269E+00
1.168529E+00
1.168622E+00
1.170296E+00
1.168644E+00
1.172975E+00
1.176543E+00
1.173389E+00
1.178422E+00
cmfd entropy
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.215201E+00
3.212850E+00
3.214094E+00
3.208054E+00
3.212891E+00
3.213349E+00
3.208290E+00
3.206895E+00
3.208786E+00
3.209630E+00
3.212321E+00
3.211750E+00
3.212453E+00
3.213370E+00
3.211791E+00
3.212685E+00
3.214145E+00
3.225292E+00
3.229509E+00
3.228530E+00
3.224203E+00
3.225547E+00
3.224720E+00
3.224546E+00
3.224527E+00
3.223579E+00
3.224380E+00
3.223483E+00
3.222819E+00
3.223067E+00
3.224007E+00
3.220616E+00
cmfd balance
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.348026E-03
5.102747E-03
4.043947E-03
3.952742E-03
2.533934E-03
2.215030E-03
2.945032E-03
2.597862E-03
2.300873E-03
1.991974E-03
1.679430E-03
1.636923E-03
1.581336E-03
1.431472E-03
1.419168E-03
1.256242E-03
4.801684E-03
2.802571E-03
1.828029E-03
2.220542E-03
1.709900E-03
2.008246E-03
2.578373E-03
2.000076E-03
1.645365E-03
1.462882E-03
1.208273E-03
1.146126E-03
1.214196E-03
1.082376E-03
8.967163E-04
1.154433E-03
cmfd dominance ratio
0.000E+00
0.000E+00
0.000E+00
0.000E+00
5.524E-01
5.476E-01
5.472E-01
5.521E-01
5.445E-01
5.527E-01
5.488E-01
5.078E-01
5.474E-01
5.419E-01
5.443E-01
5.448E-01
5.408E-01
5.391E-01
5.400E-01
5.396E-01
5.408E-01
5.398E-01
5.398E-01
5.399E-01
5.388E-01
5.389E-01
5.475E-01
5.473E-01
5.469E-01
5.461E-01
5.455E-01
5.454E-01
5.459E-01
5.460E-01
5.432E-01
cmfd openmc source comparison
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.388554E-03
7.868057E-03
6.387209E-03
6.979018E-03
5.634468E-03
5.579466E-03
5.647779E-03
5.289856E-03
4.550547E-03
4.373716E-03
4.042350E-03
3.813510E-03
3.749151E-03
3.358126E-03
3.562360E-03
3.904942E-03
9.186654E-03
6.033650E-03
3.920380E-03
4.218939E-03
4.591972E-03
4.042772E-03
4.100500E-03
3.664495E-03
3.266803E-03
3.164213E-03
3.310474E-03
3.165822E-03
3.849586E-03
2.718170E-03
2.431480E-03
3.322902E-03
cmfd source
4.142294E-02
7.484382E-02
1.050784E-01
1.256771E-01
1.444717E-01
1.420230E-01
1.350735E-01
1.118744E-01
7.755343E-02
4.198177E-02
4.296288E-02
7.964357E-02
1.107722E-01
1.359821E-01
1.425321E-01
1.356719E-01
1.285829E-01
1.040603E-01
7.630230E-02
4.348975E-02

502
tests/test_cmfd_nofeed/results_true.dat
View file

@ -1,128 +1,128 @@
k-combined:
1.170519E+00 8.422960E-03
1.171115E+00 6.173328E-03
tally 1:
1.078122E+01
1.170828E+01
1.999878E+01
4.018561E+01
2.812898E+01
7.936860E+01
3.362276E+01
1.133841E+02
3.714459E+01
1.382628E+02
3.828125E+01
1.470408E+02
3.599263E+01
1.299451E+02
3.090749E+01
9.596105E+01
2.144103E+01
4.630323E+01
1.143002E+01
1.314355E+01
1.151618E+01
1.331859E+01
2.120660E+01
4.514836E+01
2.759616E+01
7.639131E+01
3.216668E+01
1.036501E+02
3.664720E+01
1.345450E+02
3.771246E+01
1.424209E+02
3.523750E+01
1.245225E+02
2.973298E+01
8.860064E+01
2.152108E+01
4.647187E+01
1.169538E+01
1.375047E+01
tally 2:
2.247115E+01
2.548076E+01
1.574700E+01
1.251937E+01
2.119907E+00
2.284737E-01
4.160187E+01
8.733627E+01
2.945600E+01
4.383448E+01
4.027136E+00
8.279807E-01
5.722117E+01
1.643758E+02
4.057100E+01
8.267151E+01
5.388803E+00
1.465180E+00
6.769175E+01
2.300722E+02
4.800300E+01
1.157663E+02
6.213554E+00
1.946061E+00
7.376629E+01
2.729166E+02
5.253400E+01
1.385018E+02
6.438590E+00
2.103693E+00
7.454128E+01
2.790080E+02
5.311800E+01
1.417584E+02
6.784745E+00
2.328936E+00
6.907125E+01
2.397912E+02
4.912000E+01
1.213330E+02
6.405754E+00
2.075535E+00
6.012056E+01
1.815248E+02
4.280600E+01
9.207864E+01
5.534450E+00
1.555396E+00
4.229808E+01
8.999380E+01
3.003200E+01
4.541718E+01
3.844618E+00
7.537006E-01
2.272774E+01
2.597695E+01
1.577800E+01
1.252728E+01
2.221451E+00
2.528223E-01
2.274639E+01
2.606952E+01
1.588200E+01
1.270445E+01
2.140989E+00
2.357207E-01
4.205792E+01
8.880940E+01
2.970000E+01
4.427086E+01
3.919645E+00
7.773724E-01
5.560960E+01
1.559764E+02
3.947900E+01
7.872700E+01
5.238942E+00
1.400918E+00
6.492259E+01
2.117369E+02
4.612200E+01
1.069035E+02
5.989449E+00
1.813201E+00
7.217377E+01
2.608499E+02
5.148500E+01
1.327923E+02
6.607336E+00
2.205529E+00
7.305896E+01
2.681514E+02
5.187500E+01
1.352457E+02
6.722921E+00
2.290262E+00
6.884269E+01
2.380550E+02
4.904800E+01
1.208314E+02
6.177320E+00
1.927173E+00
5.902100E+01
1.748370E+02
4.201000E+01
8.858460E+01
5.542381E+00
1.549108E+00
4.268091E+01
9.151405E+01
3.029500E+01
4.614050E+01
3.822093E+00
7.420139E-01
2.362279E+01
2.812041E+01
1.653100E+01
1.377737E+01
2.336090E+00
2.851840E-01
tally 3:
1.517200E+01
1.162361E+01
9.172562E-01
4.342120E-02
2.835000E+01
4.061736E+01
1.868192E+00
1.761390E-01
3.911800E+01
7.686442E+01
2.478956E+00
3.088180E-01
4.617600E+01
1.071239E+02
2.980518E+00
4.488545E-01
5.059100E+01
1.284662E+02
3.257651E+00
5.341674E-01
5.115000E+01
1.314866E+02
3.365441E+00
5.729274E-01
4.732100E+01
1.126317E+02
3.051596E+00
4.705140E-01
4.120700E+01
8.535705E+01
2.704451E+00
3.705800E-01
2.900900E+01
4.238065E+01
1.872660E+00
1.780668E-01
1.520000E+01
1.162785E+01
1.007084E+00
5.171907E-02
1.523800E+01
1.170551E+01
1.071050E+00
5.839198E-02
2.862100E+01
4.111143E+01
1.892774E+00
1.812712E-01
3.804200E+01
7.314552E+01
2.423654E+00
2.968521E-01
4.433500E+01
9.878201E+01
2.823929E+00
4.033633E-01
4.954300E+01
1.229796E+02
3.226029E+00
5.265680E-01
4.999000E+01
1.256279E+02
3.232464E+00
5.286388E-01
4.723500E+01
1.120638E+02
3.015553E+00
4.606928E-01
4.050800E+01
8.237529E+01
2.592073E+00
3.412174E-01
2.911800E+01
4.263022E+01
1.875109E+00
1.785438E-01
1.592800E+01
1.279461E+01
1.038638E+00
5.538157E-02
tally 4:
0.000000E+00
0.000000E+00
@ -160,8 +160,8 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
3.017000E+00
4.575810E-01
3.065000E+00
4.742170E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -208,10 +208,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
5.447000E+00
1.491865E+00
2.697000E+00
3.709330E-01
5.420000E+00
1.474674E+00
2.693000E+00
3.667090E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -256,10 +256,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
7.403000E+00
2.751813E+00
5.151000E+00
1.334701E+00
7.243000E+00
2.637431E+00
5.092000E+00
1.305200E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -304,10 +304,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
8.635000E+00
3.741143E+00
7.048000E+00
2.495310E+00
8.280000E+00
3.445670E+00
6.765000E+00
2.307253E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -352,10 +352,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
9.381000E+00
4.415619E+00
8.456000E+00
3.587840E+00
8.980000E+00
4.046484E+00
8.108000E+00
3.299338E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -400,10 +400,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
9.297000E+00
4.334551E+00
9.198000E+00
4.240540E+00
9.016000E+00
4.079320E+00
8.962000E+00
4.034032E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -448,10 +448,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
8.737000E+00
3.843559E+00
9.508000E+00
4.553256E+00
8.465000E+00
3.595665E+00
9.296000E+00
4.340524E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -496,10 +496,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
7.338000E+00
2.709162E+00
8.888000E+00
3.969398E+00
7.247000E+00
2.638527E+00
8.865000E+00
3.946315E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -544,10 +544,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
5.354000E+00
1.442386E+00
7.584000E+00
2.887298E+00
5.179000E+00
1.353661E+00
7.492000E+00
2.817588E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -592,10 +592,10 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
2.707000E+00
3.709270E-01
5.507000E+00
1.522587E+00
2.821000E+00
4.067990E-01
5.617000E+00
1.587757E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -642,8 +642,8 @@ tally 4:
0.000000E+00
0.000000E+00
0.000000E+00
3.112000E+00
4.868600E-01
3.134000E+00
4.937920E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -662,114 +662,114 @@ k cmfd
0.000000E+00
0.000000E+00
0.000000E+00
1.150583E+00
1.160876E+00
1.160893E+00
1.157393E+00
1.157826E+00
1.163206E+00
1.169286E+00
1.169322E+00
1.169866E+00
1.177576E+00
1.183172E+00
1.184784E+00
1.186581E+00
1.183233E+00
1.181032E+00
1.180107E+00
1.177990E+00
1.160491E+00
1.145875E+00
1.148719E+00
1.140676E+00
1.141509E+00
1.143597E+00
1.141954E+00
1.150311E+00
1.155088E+00
1.155464E+00
1.152786E+00
1.156950E+00
1.159040E+00
1.160571E+00
1.161251E+00
cmfd entropy
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.215201E+00
3.214833E+00
3.211409E+00
3.218241E+00
3.220068E+00
3.217667E+00
3.214425E+00
3.215427E+00
3.214339E+00
3.212343E+00
3.211075E+00
3.211281E+00
3.209704E+00
3.210597E+00
3.213481E+00
3.213943E+00
3.214145E+00
3.222082E+00
3.225870E+00
3.230292E+00
3.228784E+00
3.228863E+00
3.228331E+00
3.230222E+00
3.231212E+00
3.230979E+00
3.229831E+00
3.229258E+00
3.228559E+00
3.227915E+00
3.227427E+00
3.229561E+00
cmfd balance
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.348026E-03
4.651411E-03
3.551519E-03
2.893286E-03
2.930836E-03
2.342081E-03
2.443793E-03
2.380658E-03
2.091259E-03
2.144887E-03
2.055334E-03
1.907704E-03
1.879350E-03
1.598308E-03
1.247988E-03
1.179680E-03
4.801684E-03
3.228380E-03
2.568997E-03
2.195796E-03
2.248884E-03
3.405416E-03
2.332198E-03
2.576061E-03
2.326651E-03
2.324425E-03
2.205364E-03
2.112702E-03
1.864656E-03
1.804877E-03
1.557106E-03
1.312058E-03
cmfd dominance ratio
0.000E+00
0.000E+00
0.000E+00
0.000E+00
5.524E-01
5.492E-01
5.459E-01
5.486E-01
5.475E-01
5.455E-01
5.433E-01
5.425E-01
5.412E-01
5.404E-01
5.399E-01
5.400E-01
5.401E-01
5.414E-01
5.433E-01
5.437E-01
5.472E-01
5.510E-01
5.519E-01
5.535E-01
5.535E-01
5.505E-01
5.488E-01
5.505E-01
5.510E-01
5.513E-01
5.510E-01
5.508E-01
5.487E-01
5.489E-01
5.481E-01
5.499E-01
cmfd openmc source comparison
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.388554E-03
6.525044E-03
5.810856E-03
3.651620E-03
3.448625E-03
3.805838E-03
4.064235E-03
3.244723E-03
3.963513E-03
4.186768E-03
3.785555E-03
2.732731E-03
2.308430E-03
2.074797E-03
1.581512E-03
1.729988E-03
9.186654E-03
5.964812E-03
4.465905E-03
4.119425E-03
4.973577E-03
4.092492E-03
4.063342E-03
2.804589E-03
3.632667E-03
5.005042E-03
3.428575E-03
3.007070E-03
3.091465E-03
3.030625E-03
2.751739E-03
1.762364E-03
cmfd source
3.816410E-02
7.860013E-02
1.050204E-01
1.270331E-01
1.389768E-01
1.438216E-01
1.304530E-01
1.155470E-01
7.975741E-02
4.262650E-02
4.538792E-02
8.103354E-02
1.045198E-01
1.221411E-01
1.398214E-01
1.401011E-01
1.305055E-01
1.120110E-01
8.032924E-02
4.414939E-02

2
tests/test_density_atombcm/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
1.760126E+00 1.038820E-02
1.752274E+00 4.032481E-02

2
tests/test_density_atomcm3/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
1.092203E+00 1.990176E-02
1.092376E+00 1.759788E-02

2
tests/test_density_kgm3/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
8.085745E-01 9.674599E-03
7.994522E-01 1.065745E-02

2
tests/test_density_sum/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.319139E-01 1.688777E-02
3.231215E-01 6.421320E-03

2
tests/test_eigenvalue_genperbatch/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.012381E-01 1.890480E-03
3.015627E-01 5.978844E-03

2
tests/test_eigenvalue_no_inactive/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.066374E-01 7.794575E-03
3.130246E-01 6.960311E-03

2
tests/test_energy_grid/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.215828E-01 2.966835E-03
3.155788E-01 7.559348E-03

16
tests/test_entropy/results_true.dat
View file

@ -1,13 +1,13 @@
k-combined:
3.011726E-01 1.841644E-03
3.021779E-01 3.813358E-03
entropy:
7.608094E+00
8.167702E+00
8.273634E+00
8.238974E+00
8.307173E+00
8.239618E+00
8.230443E+00
8.201657E+00
8.289158E+00
8.364683E+00
8.239452E+00
8.234598E+00
8.278421E+00
8.260773E+00
8.351860E+00
8.303719E+00
8.271058E+00

14
tests/test_filter_cell/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
1.517577E+01
4.747271E+01
3.151504E+00
2.051857E+00
4.536316E+01
4.258781E+02
1.423676E+01
4.330937E+01
2.914798E+00
1.831649E+00
4.088282E+01
3.662539E+02

6
tests/test_filter_cellborn/results_true.dat
View file

@ -1,10 +1,10 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
7.449502E+01
1.145793E+03
7.007584E+01
1.050827E+03
0.000000E+00
0.000000E+00
0.000000E+00

18
tests/test_filter_energy/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.687671E+01
1.475192E+02
4.148025E+01
3.443331E+02
5.223662E+01
5.466362E+02
1.050254E+01
2.218108E+01
2.770022E+01
1.559167E+02
4.285040E+01
3.679177E+02
5.288354E+01
5.619308E+02
1.134904E+01
2.579927E+01

18
tests/test_filter_energyout/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.740000E+01
1.516786E+02
4.163000E+01
3.469463E+02
5.041000E+01
5.097599E+02
6.990000E+00
9.850700E+00
2.835000E+01
1.631463E+02
4.252000E+01
3.622146E+02
5.192000E+01
5.401444E+02
7.390000E+00
1.097710E+01

66
tests/test_filter_group_transfer/results_true.dat
View file

@ -1,54 +1,54 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.470000E+01
1.233286E+02
2.571000E+01
1.344145E+02
0.000000E+00
0.000000E+00
7.000000E-02
2.100000E-03
2.000000E-02
2.000000E-04
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.590851E-01
1.552999E-01
8.909862E-01
1.660072E-01
0.000000E+00
0.000000E+00
2.474615E+00
1.227479E+00
2.700000E+00
1.469800E+00
2.269580E+00
1.046926E+00
2.640000E+00
1.397800E+00
0.000000E+00
0.000000E+00
3.712000E+01
2.758636E+02
3.799000E+01
2.892451E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.370712E-01
2.596756E-02
3.267964E-01
2.170374E-02
0.000000E+00
0.000000E+00
9.771334E-01
1.963983E-01
8.811378E-01
1.642094E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.440000E+00
3.947400E+00
4.510000E+00
4.068900E+00
0.000000E+00
0.000000E+00
4.688000E+01
4.409726E+02
2.867366E-02
2.783031E-04
4.842000E+01
4.700146E+02
6.079426E-02
8.532151E-04
0.000000E+00
0.000000E+00
1.080794E-01
2.580055E-03
1.182446E-01
3.979568E-03
0.000000E+00
0.000000E+00
0.000000E+00
@ -57,11 +57,11 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
3.530000E+00
2.494900E+00
1.572766E-01
5.137814E-03
6.990000E+00
9.850700E+00
3.134136E-01
2.138680E-02
3.500000E+00
2.455800E+00
1.158638E-01
2.993934E-03
7.390000E+00
1.097710E+01
3.126942E-01
2.185793E-02

18
tests/test_filter_material/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.819256E+01
1.591068E+02
6.599750E+00
8.721650E+00
5.383171E+01
5.950793E+02
4.140672E+01
3.649396E+02
2.923791E+01
1.711686E+02
6.753642E+00
9.145065E+00
4.921142E+01
5.186991E+02
4.731999E+01
4.905260E+02

590
tests/test_filter_mesh_2d/results_true.dat
View file

@ -1,5 +1,5 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
@ -45,8 +45,10 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
8.751430E-01
7.658753E-01
3.228098E-02
1.042062E-03
3.222708E-01
1.038585E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -73,14 +75,26 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
8.182335E-01
3.748630E-01
2.711997E-01
5.338821E-02
3.359680E-01
5.168399E-02
0.000000E+00
0.000000E+00
1.287286E-01
1.657106E-02
7.512292E-01
1.499955E-01
2.231206E+00
1.531548E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.706070E-01
1.373496E-01
0.000000E+00
0.000000E+00
0.000000E+00
@ -95,120 +109,106 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.855163E-02
4.699326E-03
6.967594E-02
4.854737E-03
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.993440E-02
9.986884E-03
1.978123E-01
3.912972E-02
8.758664E-02
7.671419E-03
4.204632E-01
7.722514E-02
2.486465E+00
1.923132E+00
1.697091E-01
1.631523E-02
0.000000E+00
0.000000E+00
3.931419E-01
9.002841E-02
1.092722E+00
3.733055E-01
2.384227E+00
1.926937E+00
9.101131E-01
3.634496E-01
3.284661E-01
1.078900E-01
0.000000E+00
0.000000E+00
6.885295E-02
4.740728E-03
0.000000E+00
0.000000E+00
6.175782E-01
2.022672E-01
1.204200E-02
8.137904E-05
2.136417E-01
1.678274E-02
2.932366E-01
4.285227E-02
6.387464E-01
1.526040E-01
2.419633E-02
5.854622E-04
9.286912E-02
7.247209E-03
5.629729E-01
9.281109E-02
7.345786E-01
1.755550E-01
1.219449E-01
1.487057E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.062751E+00
3.899247E-01
1.264401E+00
5.654864E-01
1.615537E+00
5.605431E-01
4.267750E-01
1.013728E-01
2.024264E+00
9.373843E-01
1.122195E-01
7.011342E-03
8.223839E-01
2.170405E-01
1.239203E+00
4.165595E-01
1.169828E+00
4.406334E-01
1.373464E+00
6.398480E-01
2.824398E+00
2.991665E+00
6.158518E-01
1.951901E-01
7.188149E-01
1.091668E-01
7.595527E-01
2.108412E-01
5.784902E-01
1.362812E-01
7.167515E-01
3.225769E-01
2.081496E-02
4.332627E-04
9.212853E-01
2.277306E-01
1.266470E+00
6.777774E-01
1.983303E-01
3.797884E-02
4.840974E-02
1.266547E-03
1.183485E+00
4.184814E-01
3.027254E-01
8.598449E-02
9.889868E-01
4.608531E-01
8.698103E-01
4.598559E-01
1.332831E+00
4.809984E-01
1.564949E+00
5.782651E-01
1.143572E+00
4.399439E-01
1.326651E+00
6.376565E-01
1.716813E+00
1.314280E+00
7.673229E-01
2.364966E-01
2.539284E+00
1.945563E+00
1.263219E+00
4.919339E-01
5.430042E-01
1.300266E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
7.553587E-01
1.738170E-01
2.048487E+00
1.239856E+00
3.761862E-01
8.452912E-02
1.350273E-02
1.823238E-04
3.348203E-01
6.603737E-02
2.382593E-01
2.884302E-02
1.548117E+00
7.256982E-01
9.635803E-01
4.452538E-01
1.030188E+00
3.496685E-01
7.053043E-01
2.814424E-01
1.260458E-01
8.365754E-03
1.484515E+00
7.183840E-01
1.887042E+00
1.143800E+00
3.798783E+00
4.487065E+00
2.274180E+00
1.991784E+00
1.903719E-01
3.237886E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.675232E-02
9.361012E-03
2.319594E-01
2.331698E-02
1.573495E+00
6.394722E-01
4.432570E-01
1.005943E-01
9.353148E-01
3.125416E-01
8.359366E-01
2.985072E-01
1.657665E+00
9.207020E-01
3.737550E+00
3.558505E+00
1.742376E+00
8.732217E-01
5.153816E+00
6.543973E+00
1.653035E+00
1.061068E+00
9.963191E-01
3.716347E-01
2.282805E-01
2.383414E-02
8.749983E-01
2.714666E-01
1.728411E-01
1.190218E-02
9.250054E-02
@ -219,28 +219,28 @@ tally 1:
1.711154E-02
3.218800E-01
7.906855E-02
9.721141E-01
3.463278E-01
8.364857E-01
2.612988E-01
6.317454E-01
1.837078E-01
2.134394E-01
1.710045E-02
3.743814E-02
1.401614E-03
1.314651E+00
7.700853E-01
7.389940E-01
2.322340E-01
6.147456E-01
1.001719E-01
3.157573E+00
2.794673E+00
6.671644E-01
1.996570E-01
0.000000E+00
0.000000E+00
1.057036E+00
3.778638E-01
9.231639E-01
2.991795E-01
3.375678E-01
1.041383E-01
1.181686E-01
8.023920E-03
4.912969E-01
2.073894E-01
9.395786E-01
4.653730E-01
6.998437E-01
2.917085E-01
3.074214E+00
2.819088E+00
2.570673E+00
1.358321E+00
1.108912E+00
3.843307E-01
4.950896E-02
2.451137E-03
0.000000E+00
0.000000E+00
0.000000E+00
@ -255,24 +255,24 @@ tally 1:
7.670183E-03
0.000000E+00
0.000000E+00
9.451595E-01
2.560533E-01
2.255990E+00
1.237092E+00
9.318430E-01
2.664794E-01
5.320317E-01
2.267827E-01
9.971891E-01
5.070904E-01
7.088786E-02
5.025089E-03
3.601108E-02
1.296798E-03
3.683352E+00
3.243246E+00
3.708141E+00
3.629996E+00
6.469411E-01
1.789549E-01
7.829878E-01
2.033989E-01
8.994770E-01
2.351438E-01
4.712797E-01
7.213659E-02
2.133532E+00
9.765422E-01
4.533607E-01
1.511497E-01
1.878729E+00
2.099266E+00
4.287190E+00
4.748691E+00
1.961229E+00
1.085868E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -285,28 +285,28 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.238471E-01
1.533812E-02
2.580606E+00
1.878196E+00
2.843921E+00
1.997749E+00
1.145124E+00
4.571704E-01
1.135664E+00
4.319756E-01
1.061692E-01
1.127190E-02
1.912282E-01
3.656822E-02
3.289827E-01
1.075283E-01
1.750908E+00
8.092384E-01
2.156426E+00
9.498067E-01
1.480596E+00
6.002164E-01
3.249216E-01
1.005106E-01
8.875810E-02
7.878000E-03
7.494206E-02
5.616313E-03
1.007353E+00
3.619658E-01
1.649444E+00
7.918743E-01
2.458176E-01
4.377921E-02
2.766784E+00
2.677426E+00
2.703501E+00
2.548083E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -319,28 +319,28 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.286351E-01
1.805028E-01
1.595058E+00
6.351621E-01
1.568755E+00
6.468874E-01
2.301331E+00
1.845232E+00
4.577930E-01
1.589685E-01
1.981880E-01
3.927848E-02
2.789456E-01
5.304261E-02
3.897689E-01
9.413685E-02
9.140885E-01
2.822974E-01
1.887726E+00
7.592780E-01
1.841624E+00
1.680236E+00
4.059938E-01
1.562853E-01
2.897077E-01
8.393057E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.772980E-01
3.143459E-02
2.445051E-01
5.978276E-02
2.234657E-01
2.716625E-02
0.000000E+00
0.000000E+00
0.000000E+00
@ -353,24 +353,24 @@ tally 1:
0.000000E+00
4.783335E-02
2.288029E-03
8.744336E-01
2.611835E-01
1.754522E+00
1.459573E+00
3.549403E-01
6.365363E-02
1.843261E+00
9.517219E-01
1.009556E+00
4.611770E-01
1.540221E+00
7.611324E-01
8.030427E-01
2.255596E-01
5.606011E-01
1.607491E-01
1.908325E+00
1.505641E+00
1.255795E-01
1.541635E-02
7.395548E-01
2.274416E-01
5.986733E-01
9.576988E-02
1.095026E+00
4.872910E-01
1.043470E+00
3.153965E-01
1.004973E+00
6.046910E-01
1.664911E-01
2.771927E-02
1.724071E-01
2.775427E-02
0.000000E+00
0.000000E+00
0.000000E+00
@ -385,114 +385,114 @@ tally 1:
0.000000E+00
5.470039E-01
2.992133E-01
4.902613E-01
1.163360E-01
4.313918E-01
1.128703E-01
1.088037E+00
6.007745E-01
1.013469E+00
5.646900E-01
5.529904E-01
2.308161E-01
3.523818E-01
1.241729E-01
2.031343E-01
4.126355E-02
2.664927E+00
2.323850E+00
2.601166E+00
1.593528E+00
1.726764E+00
6.491827E-01
7.083106E-02
3.090181E-03
4.738741E-01
2.245567E-01
6.086518E-02
3.704570E-03
1.126662E+00
4.675322E-01
1.008459E+00
4.616352E-01
1.309592E+00
5.665211E-01
1.334050E+00
5.264819E-01
7.324996E-01
2.577124E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.227736E-01
1.041828E-01
9.218244E-01
2.000301E-01
2.119900E+00
1.123846E+00
4.366404E-02
1.015133E-03
0.000000E+00
0.000000E+00
2.213571E-01
2.455759E-02
1.471329E+00
7.001977E-01
6.348345E-01
1.603344E-01
0.000000E+00
0.000000E+00
7.240936E-01
1.701659E-01
1.280272E+00
3.624213E-01
9.794833E-01
4.411714E-01
2.309730E-01
2.570742E-02
1.270911E+00
4.617932E-01
1.107069E+00
4.574496E-01
1.269137E-01
1.610709E-02
1.120245E-01
1.254949E-02
2.207099E-01
4.871286E-02
9.075694E-02
8.236821E-03
1.046380E-01
7.875036E-03
2.836364E-01
3.508209E-02
4.509177E-01
7.393615E-02
1.077505E+00
3.209541E-01
1.204982E-02
1.451982E-04
0.000000E+00
0.000000E+00
8.072240E-01
2.272369E-01
3.032594E-01
5.076695E-02
7.265258E-01
1.893670E-01
2.902123E-01
3.514713E-02
0.000000E+00
0.000000E+00
2.502937E-01
5.035125E-02
1.367234E+00
5.128884E-01
5.563575E-01
2.510519E-01
3.174809E-01
1.007941E-01
9.152129E-01
2.609325E-01
9.040668E-01
2.263595E-01
7.868812E-01
2.436310E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.608891E-02
4.586038E-03
1.319691E+00
5.856361E-01
7.752262E-01
3.101466E-01
2.606443E-01
6.793546E-02
3.686503E-01
6.346661E-02
6.935008E-01
1.992259E-01
1.591064E+00
6.804860E-01
2.418186E-01
5.847626E-02
0.000000E+00
0.000000E+00
5.946272E-02
3.535814E-03
3.120645E-02
9.738427E-04
7.515014E-02
5.647544E-03
8.673868E-01
3.013632E-01
2.427216E-01
2.798481E-02
3.390904E-01
4.907887E-02
6.577001E-01
2.346964E-01
1.023735E-01
8.287220E-03
1.499600E-02
2.248801E-04
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.079687E-01
5.125039E-02
1.649743E+00
9.187008E-01
8.924207E-01
2.872456E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
8.137844E-02
5.984482E-03
1.395650E-01
1.947839E-02
1.040555E+00
3.043523E-01
1.426976E+00
6.218295E-01
8.342758E-01
2.539692E-01
3.101170E-01
9.617255E-02
6.319919E-02
3.629541E-03
1.292774E-01
8.674372E-03
0.000000E+00
0.000000E+00
0.000000E+00
@ -515,14 +515,14 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
8.400001E-01
3.766524E-01
2.709865E+00
1.701203E+00
2.332406E-01
4.591958E-02
9.579949E-02
9.177542E-03
9.587357E-01
4.992769E-01
1.756477E+00
7.884472E-01
2.541705E-01
4.325743E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
@ -551,10 +551,10 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
3.366480E-01
8.328374E-02
7.405305E-01
5.483854E-01
2.422913E-01
5.870510E-02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00

1634
tests/test_filter_mesh_3d/results_true.dat
View file

File diff suppressed because it is too large Load diff

18
tests/test_filter_universe/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
6.369043E+01
8.385422E+02
7.330762E+00
1.137874E+01
5.011385E+01
5.251276E+02
5.132453E+00
5.801019E+00
5.803437E+01
7.324822E+02
7.168085E+00
1.101286E+01
5.815302E+01
7.373019E+02
5.778029E+00
7.558661E+00

7
tests/test_fixed_source/results_true.dat
View file

@ -1,3 +1,6 @@
tally 1:
4.483337E+02
2.017057E+04
4.538791E+02
2.073271E+04
leakage:
9.830000E+00
9.663900E+00

4
tests/test_fixed_source/test_fixed_source.py
View file

@ -28,6 +28,10 @@ class FixedSourceTestHarness(TestHarness):
outstr += '\n'.join(results) + '\n'
tally_num += 1
outstr += 'leakage:\n'
outstr += '{0:12.6E}'.format(sp._global_tallies[3][0]) + '\n'
outstr += '{0:12.6E}'.format(sp._global_tallies[3][1]) + '\n'
return outstr

2
tests/test_infinite_cell/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
9.998895E-02 2.846817E-04
9.788797E-02 1.378250E-03

2
tests/test_lattice_mixed/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
9.605085E-01 8.180822E-03
9.922449E-01 1.281824E-02

2
tests/test_lattice_multiple/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02

214
tests/test_many_scores/results_true.dat
View file

@ -1,111 +1,111 @@
k-combined:
0.000000E+00 0.000000E+00
tally 1:
2.254760E+01
1.695439E+02
1.017400E+01
3.450623E+01
8.690000E+00
2.517771E+01
8.696000E+00
2.521248E+01
5.224881E-01
9.127522E-02
8.690000E+00
2.517771E+01
9.665236E-01
3.121401E-01
5.224881E-01
9.127522E-02
5.199964E-01
9.041516E-02
8.696000E+00
2.521248E+01
9.666447E-01
3.122176E-01
5.199964E-01
9.041516E-02
9.262308E+00
2.859762E+01
8.684000E+00
2.514296E+01
1.484000E+00
7.346160E-01
1.804688E+00
1.089100E+00
1.334833E+02
5.960573E+03
2.254760E+01
1.695439E+02
-6.230152E-02
1.023784E-02
-2.696818E-01
8.525477E-02
-1.670294E-01
4.366614E-02
1.331953E-02
8.064039E-05
-2.490295E-01
2.511425E-02
8.289978E-02
2.753446E-03
-1.390686E-01
1.592673E-02
1.600203E-01
1.344513E-02
1.017400E+01
3.450623E+01
-6.564301E-02
2.564587E-03
-1.245391E-01
1.143479E-02
-5.550309E-02
9.334547E-03
4.671277E-02
1.795787E-03
-9.901000E-02
4.609654E-03
2.619012E-02
3.011793E-04
-7.090362E-02
2.450750E-03
4.059073E-02
8.841300E-04
8.690000E+00
2.517771E+01
-2.685052E-02
2.769798E-04
4.047919E-03
3.384947E-03
-1.325869E-01
8.610169E-03
-1.594962E-02
5.528670E-04
-1.320550E-02
1.856248E-04
1.185337E-02
3.428536E-04
2.622196E-02
3.774884E-04
2.938463E-02
6.660370E-04
8.696000E+00
2.521248E+01
-2.634867E-02
2.718345E-04
4.042613E-03
3.367272E-03
-1.328903E-01
8.576852E-03
-1.632965E-02
5.577169E-04
-1.408565E-02
1.815109E-04
1.202519E-02
3.685250E-04
2.634054E-02
3.730853E-04
2.941341E-02
6.479457E-04
2.247257E+01
1.683779E+02
1.014000E+01
3.427460E+01
3.427342E+01
8.628000E+00
2.481430E+01
8.628000E+00
2.481430E+01
5.102293E-01
8.710841E-02
8.628000E+00
2.481430E+01
9.329009E-01
2.902534E-01
5.102293E-01
8.710841E-02
5.102293E-01
8.710841E-02
8.628000E+00
2.481430E+01
9.329009E-01
2.902534E-01
5.102293E-01
8.710841E-02
9.212024E+00
2.829472E+01
8.628000E+00
2.481430E+01
1.512000E+00
7.620560E-01
1.816851E+00
1.102658E+00
1.338067E+02
5.986137E+03
2.247257E+01
1.683779E+02
-1.512960E-01
2.623972E-02
-3.775020E-01
1.055377E-01
-1.916133E-01
4.680798E-02
2.754367E-02
3.320008E-04
-2.028374E-02
1.319357E-02
8.974271E-03
1.681081E-03
-1.658978E-01
1.520448E-02
2.878360E-01
5.645480E-02
1.014000E+01
3.427342E+01
-4.798897E-02
1.551226E-03
-1.818770E-01
1.492633E-02
-6.340651E-02
9.011305E-03
3.395308E-02
4.612818E-04
-2.640250E-02
6.434787E-04
-8.242639E-03
9.516540E-04
-8.378601E-02
2.645988E-03
9.567484E-02
7.262477E-03
8.628000E+00
2.481430E+01
-4.712248E-02
1.140942E-03
-6.431930E-02
4.290580E-03
-9.251642E-02
8.134201E-03
1.020119E-04
1.154184E-04
-2.994164E-02
3.079076E-04
2.128844E-02
2.046549E-04
1.637972E-02
1.459209E-04
4.629047E-02
7.823267E-04
8.628000E+00
2.481430E+01
-4.712248E-02
1.140942E-03
-6.431930E-02
4.290580E-03
-9.251642E-02
8.134201E-03
1.020119E-04
1.154184E-04
-2.994164E-02
3.079076E-04
2.128844E-02
2.046549E-04
1.637972E-02
1.459209E-04
4.629047E-02
7.823267E-04
1.014000E+01
3.427342E+01

2
tests/test_natural_element/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
1.053229E+00 9.241274E-02
1.013112E+00 2.551515E-02

2
tests/test_output/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.011726E-01 1.841644E-03
3.021779E-01 3.813358E-03

10
tests/test_particle_restart_eigval/results_true.dat
View file

@ -1,16 +1,16 @@
current batch:
1.200000E+01
9.000000E+00
current gen:
1.000000E+00
particle id:
6.160000E+02
5.550000E+02
run mode:
2.000000E+00
particle weight:
1.000000E+00
particle energy:
3.545295E-01
2.831611E-01
particle xyz:
3.516323E+01 -5.400148E+01 -1.588825E+01
4.973847E+01 6.971699E+00 -5.201827E+01
particle uvw:
4.129799E-01 7.649720E-01 4.942322E-01
6.945105E-01 6.295355E-01 -3.483393E-01

2
tests/test_particle_restart_eigval/test_particle_restart_eigval.py
View file

@ -6,5 +6,5 @@ from testing_harness import ParticleRestartTestHarness
if __name__ == '__main__':
harness = ParticleRestartTestHarness('particle_12_616.*')
harness = ParticleRestartTestHarness('particle_9_555.*')
harness.main()

10
tests/test_particle_restart_fixed/results_true.dat
View file

@ -1,16 +1,16 @@
current batch:
7.000000E+00
current gen:
0.000000E+00
1.000000E+00
particle id:
6.144000E+03
9.280000E+02
run mode:
1.000000E+00
particle weight:
1.000000E+00
particle energy:
5.749729E+00
4.412022E+00
particle xyz:
8.754675E+00 2.551620E+00 4.394350E-01
5.572639E+00 -9.139472E+00 -1.974824E+00
particle uvw:
-5.971721E-01 -4.845709E-01 6.391999E-01
1.410422E-01 5.330426E-01 -8.342498E-01

2
tests/test_particle_restart_fixed/test_particle_restart_fixed.py
View file

@ -6,5 +6,5 @@ from testing_harness import ParticleRestartTestHarness
if __name__ == '__main__':
harness = ParticleRestartTestHarness('particle_7_6144.*')
harness = ParticleRestartTestHarness('particle_7_928.*')
harness.main()

2
tests/test_ptables_off/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
2.984064E-01 3.464413E-03
2.998034E-01 5.227986E-03

2
tests/test_reflective_cone/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
2.281477E+00 5.206704E-03
2.269987E+00 4.469683E-03

2
tests/test_reflective_cylinder/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
2.276055E+00 3.186526E-03
2.272436E+00 7.831006E-04

2
tests/test_reflective_plane/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
2.281505E+00 3.254688E-03
2.276127E+00 4.678320E-03

2
tests/test_reflective_sphere/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
2.280812E+00 5.087968E-03
2.271012E+00 3.466350E-03

2
tests/test_resonance_scattering/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
6.452021E-02 1.738736E-03
6.842112E-02 8.480934E-04

2
tests/test_rotation/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
3.011726E-01 1.841644E-03
3.021779E-01 3.813358E-03

2
tests/test_salphabeta/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
8.418898E-01 5.633536E-03
8.538165E-01 6.355606E-03

2
tests/test_salphabeta_multiple/results_true.dat
View file

@ -1,2 +1,2 @@
k-combined:
9.686128E-01 5.056987E-02
9.141547E-01 2.728809E-02

42
tests/test_score_MT/results_true.dat
View file

@ -1,5 +1,5 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
@ -9,27 +9,27 @@ tally 1:
0.000000E+00
0.000000E+00
0.000000E+00
1.107640E-02
4.232516E-05
1.107640E-02
4.232516E-05
3.509526E-01
2.565440E-02
1.171523E+00
2.824070E-01
5.592851E-04
1.423656E-07
5.592851E-04
1.423656E-07
3.027322E-02
1.983924E-04
1.477313E-02
4.495633E-05
6.261679E-03
2.310296E-05
6.261679E-03
2.310296E-05
3.434524E-01
2.482946E-02
1.059004E+00
2.432555E-01
3.828838E-04
9.365499E-08
3.828838E-04
9.365499E-08
3.226036E-02
2.229037E-04
1.450446E-02
4.382852E-05
8.650696E-08
7.483455E-15
8.650696E-08
7.483455E-15
1.097403E-04
3.424171E-09
6.866942E-02
9.952916E-04
7.312640E-05
2.080857E-09
6.101318E-02
8.452067E-04

22
tests/test_score_absorption/results_true.dat
View file

@ -1,20 +1,20 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
2.267589E+00
1.054691E+00
1.476133E-02
4.485759E-05
3.306896E-01
2.302248E-02
2.006479E+00
8.808050E-01
1.444838E-02
4.349668E-05
2.926865E-01
1.944083E-02
tally 2:
0.000000E+00
0.000000E+00
2.290000E+00
1.070700E+00
2.030000E+00
8.879000E-01
0.000000E+00
0.000000E+00
3.800000E-01
3.820000E-02
4.000000E-01
4.240000E-02

2
tests/test_score_current/results_true.dat
View file

@ -1 +1 @@
8f7125c9686a94f0fdbd76b1667f84b9dab60c068e23d4f2b0105ce058fa533a8a4a221355614add10d54f3b18c14af41891e0fd0593154c73894e5dcb6fe5c2
1e6945632c55491d4584f4976cc6f5c7340874703cfaf739dd956b7124b4260955efb5b6ba041b32536f9a74572d071e0293dced55a41ea305223f698b734c2a

18
tests/test_score_events/results_true.dat
View file

@ -1,12 +1,12 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
5.627000E+01
6.565251E+02
4.892000E+01
5.009100E+02
5.328000E+01
6.082750E+02
5.594000E+01
6.836068E+02
tally 2:
1.458000E+01
4.391620E+01
1.269000E+01
3.363670E+01
1.372000E+01
4.018980E+01
1.474000E+01
4.802520E+01

18
tests/test_score_fission/results_true.dat
View file

@ -1,20 +1,20 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
1.089685E+00
2.438289E-01
9.432574E-01
1.970463E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.094323E-01
1.773090E-01
1.039280E+00
2.345809E-01
tally 2:
1.081920E+00
2.404883E-01
8.953531E-01
1.798609E-01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.945760E-01
2.073005E-01
9.923196E-01
2.067216E-01

26
tests/test_score_flux/results_true.dat
View file

@ -1,15 +1,15 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
3.403102E+01
2.393236E+02
1.143060E+01
2.704386E+01
5.691243E+01
6.685581E+02
2.924945E+01
1.782032E+02
9.965051E+00
2.067319E+01
4.931805E+01
5.050490E+02
3.224218E+01
2.217821E+02
1.079687E+01
2.494701E+01
5.252579E+01
5.951604E+02
3.325822E+01
2.414028E+02
1.145800E+01
2.880575E+01
5.605671E+01
6.804062E+02

890
tests/test_score_flux_yn/results_true.dat
View file

@ -1,448 +1,448 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
3.403102E+01
2.393236E+02
1.143060E+01
2.704386E+01
5.691243E+01
6.685581E+02
2.924945E+01
1.782032E+02
9.965051E+00
2.067319E+01
4.931805E+01
5.050490E+02
3.224218E+01
2.217821E+02
1.079687E+01
2.494701E+01
5.252579E+01
5.951604E+02
3.325822E+01
2.414028E+02
1.145800E+01
2.880575E+01
5.605671E+01
6.804062E+02
tally 2:
3.403102E+01
2.393236E+02
-3.055588E-01
3.594044E-02
3.171128E-01
3.805380E-01
-7.374877E-01
2.391349E-01
-1.621986E-01
5.578051E-02
4.710459E-01
1.850595E-01
-4.839599E-01
2.440735E-01
1.124798E-01
8.640649E-02
-4.976706E-01
8.943606E-02
1.323704E-01
6.713903E-02
-6.954103E-02
7.876989E-02
1.343675E-01
2.569171E-02
-5.205013E-01
7.985579E-02
2.037295E-01
1.501360E-01
-1.290735E-01
5.069422E-02
-2.045060E-01
2.197889E-01
1.307602E-01
5.938649E-02
-3.317004E-01
1.359233E-01
-4.347480E-01
7.528253E-02
7.685982E-02
4.559715E-02
-4.668784E-02
4.065614E-02
1.447872E-01
5.538956E-02
-2.027337E-01
7.095051E-02
-4.673353E-01
1.235371E-01
-3.321732E-01
6.905624E-02
-3.028859E-01
3.382353E-02
-2.320613E-01
2.845904E-02
8.038034E-02
3.466746E-02
1.766484E-01
9.996559E-02
3.385184E-02
5.510769E-02
6.952638E-02
2.741119E-02
3.630715E-01
5.336447E-02
2.876222E-01
9.938256E-02
-3.473240E-02
1.031412E-01
1.280013E-01
2.331384E-02
1.550774E-01
3.881493E-02
1.143060E+01
2.704386E+01
-1.955516E-01
9.789273E-03
4.070149E-02
3.530892E-02
-1.446334E-01
3.374966E-02
-4.072014E-02
1.407320E-02
1.056190E-01
2.144988E-02
-1.491728E-01
1.699068E-02
3.130540E-02
4.945783E-03
-1.706770E-01
6.943395E-03
6.278466E-02
4.522701E-03
-6.491915E-02
1.218265E-02
-6.226375E-02
8.920919E-03
-2.647477E-01
2.028918E-02
-9.179913E-03
1.186941E-02
-8.959288E-02
9.826050E-03
-4.800079E-02
3.500548E-02
9.919028E-02
1.278789E-02
-1.274511E-01
9.357943E-03
-9.037073E-02
4.686665E-03
-7.895199E-03
8.526809E-03
-2.071827E-02
5.844891E-03
-1.422877E-02
2.008170E-03
-7.787500E-02
9.642362E-03
-6.105504E-02
7.600775E-03
-1.249117E-01
1.314509E-02
-9.183605E-02
7.513215E-03
-9.171222E-02
4.220465E-03
8.798390E-02
5.037483E-03
4.366496E-02
8.411743E-03
1.605678E-02
6.994511E-03
-2.179875E-02
1.352098E-03
1.116755E-01
7.354424E-03
9.995739E-02
9.710026E-03
3.010324E-02
1.179938E-02
2.175815E-02
1.281441E-03
2.374886E-02
6.943432E-03
5.691243E+01
6.685581E+02
-7.130295E-01
1.971869E-01
1.138602E+00
1.071873E+00
-1.436194E+00
7.200870E-01
-2.384966E-01
1.129688E-01
8.003085E-01
5.809429E-01
-3.499868E-01
3.310435E-01
2.580763E-01
2.238385E-01
-3.576056E-01
1.437071E-01
2.560747E-01
4.492224E-02
-4.006295E-01
2.456809E-01
8.176054E-02
3.924768E-01
-8.105689E-01
2.397175E-01
3.777117E-01
4.811523E-01
-4.932440E-01
1.789106E-01
-3.705746E-02
5.517840E-01
7.049764E-01
4.243943E-01
-4.833188E-01
2.014189E-01
-5.420278E-01
1.201572E-01
7.606099E-02
1.451071E-01
2.244026E-01
7.628680E-02
2.495758E-01
1.882655E-01
-1.373888E-01
1.958724E-01
-3.799191E-01
2.258678E-01
-5.909624E-01
2.940322E-01
-6.336238E-01
1.272552E-01
-5.544697E-01
9.027652E-02
4.161023E-01
1.084744E-01
5.550548E-01
1.838427E-01
3.060033E-01
8.542356E-02
3.623817E-01
4.245612E-02
4.779023E-01
8.081477E-02
1.089872E-01
1.402042E-01
-1.107440E-01
3.140479E-01
2.333168E-01
1.035674E-01
-5.959476E-02
2.482914E-01
2.924945E+01
1.782032E+02
5.986908E-01
4.163713E-01
3.376452E-01
1.744277E-01
1.272993E-01
1.526084E-01
-2.551267E-01
7.284601E-02
4.755790E-02
6.832796E-02
-3.425133E-01
7.756266E-02
3.823240E-01
1.789133E-01
-3.767347E-01
1.676013E-01
1.655399E-01
1.153056E-01
1.432265E-01
2.956364E-02
-2.989439E-01
5.066565E-02
-3.226047E-02
5.960686E-02
-9.968550E-02
4.039270E-02
-1.235291E-02
2.805714E-01
-7.494860E-02
2.838999E-02
3.042503E-02
7.425538E-02
-1.849034E-01
8.599532E-02
-5.209242E-01
1.174277E-01
-2.314402E-01
4.065974E-02
-4.179560E-01
1.382716E-01
1.026654E-01
1.611177E-02
-2.270386E-01
5.157513E-02
7.884296E-02
2.677093E-02
7.284620E-01
1.871262E-01
1.360134E-01
9.437335E-03
-3.088199E-01
2.292832E-02
-2.543226E-01
1.264440E-01
6.256788E-02
7.604299E-02
-3.478518E-01
6.632508E-02
2.842847E-01
4.906462E-02
1.546241E-01
3.362908E-02
-2.247690E-01
1.912574E-02
-3.459210E-02
1.248892E-01
-2.204075E-01
4.281299E-02
2.145675E-01
5.593188E-02
9.965051E+00
2.067319E+01
2.381967E-01
4.786935E-02
-1.924308E-02
4.029344E-02
3.179175E-02
2.098033E-02
-1.652308E-02
9.391525E-03
-5.733039E-02
1.788395E-02
4.942155E-02
1.206623E-03
1.874861E-02
1.917059E-02
-1.179838E-01
2.185564E-02
-9.270508E-03
1.659002E-02
7.496511E-02
6.849931E-03
-1.004441E-01
4.570271E-03
-9.899601E-02
1.027435E-02
-9.963860E-02
4.686519E-03
-1.740638E-02
3.936514E-02
2.520349E-02
5.900520E-03
1.227791E-03
6.173783E-03
-3.068441E-02
6.867402E-03
-7.765654E-02
1.399383E-02
7.333551E-02
4.139444E-03
-3.508914E-02
1.421599E-02
4.334120E-02
4.918888E-03
-1.106116E-01
1.143785E-02
7.254477E-03
1.337965E-03
2.365948E-01
2.473485E-02
5.176646E-02
1.514750E-03
-8.723610E-02
2.395100E-03
-1.010207E-01
1.081070E-02
-5.113298E-02
6.053279E-03
-1.601115E-01
8.304230E-03
-5.105345E-03
5.079654E-03
-8.894549E-02
2.505026E-03
-1.023874E-01
4.311776E-03
7.972837E-02
1.178305E-02
-2.560942E-02
4.101793E-03
2.312203E-02
3.036209E-03
4.931805E+01
5.050490E+02
1.513494E+00
1.676567E+00
6.650656E-01
1.400732E+00
4.934620E-01
7.673322E-01
-2.964442E-02
2.262500E-01
-2.225956E-01
4.872938E-01
2.498269E-01
9.442256E-02
-1.454141E-01
2.683171E-01
-3.720260E-01
2.266316E-01
2.473591E-01
3.624408E-01
2.096941E-01
6.664835E-02
-4.700220E-01
9.789465E-02
-1.626165E-01
1.315362E-01
3.519966E-01
1.140516E-01
-1.185342E-01
5.708362E-01
-1.333083E-01
1.352294E-01
1.375803E-01
1.572243E-01
-6.099544E-02
2.042264E-01
-3.664336E-01
1.762506E-01
-1.717957E-01
5.864778E-02
-5.153204E-01
2.484316E-01
-2.258513E-01
8.839687E-02
-9.307565E-01
4.802462E-01
1.108577E-01
4.451363E-02
1.118965E+00
5.173979E-01
4.314600E-01
5.176751E-02
-2.864059E-01
9.159895E-02
-4.471918E-01
1.743046E-01
-7.850808E-02
1.317732E-01
-7.445394E-01
1.479808E-01
7.480873E-01
2.696289E-01
1.731449E-01
5.301367E-02
-3.663464E-01
3.430563E-02
3.553351E-01
2.051544E-01
-4.312304E-01
2.164291E-01
2.436515E-01
8.025128E-02
3.224218E+01
2.217821E+02
-4.524803E-01
1.326110E-01
-2.725152E-01
4.845463E-01
-4.030063E-01
1.825349E-01
-2.145346E-01
7.197305E-02
2.488985E-01
4.295894E-02
-6.797085E-01
2.125889E-01
1.698034E-01
3.750580E-02
-6.596628E-02
6.507545E-02
4.519597E-02
1.432412E-02
6.227959E-01
1.066624E-01
2.328135E-01
8.583522E-02
-8.899482E-01
1.982834E-01
3.269468E-01
5.392549E-02
1.684792E-01
1.461138E-01
-2.026161E-01
1.992812E-01
2.599678E-01
7.704699E-02
-3.329837E-01
8.761827E-02
-5.287177E-01
9.089410E-02
3.666982E-01
4.451200E-02
2.096594E-01
4.539545E-02
2.091762E-01
4.861771E-02
1.476292E-01
2.837569E-02
-4.091772E-01
1.141307E-01
-1.652694E-01
5.552090E-02
8.207679E-02
1.659512E-02
-1.203575E-01
6.408952E-02
-8.219791E-03
9.915873E-02
-3.768675E-02
8.535079E-02
1.305167E-01
1.917693E-02
1.823050E-01
1.945159E-02
2.999655E-01
4.587948E-02
3.687435E-01
1.382795E-01
1.220376E-01
6.297901E-02
3.899923E-01
5.894393E-02
-1.878957E-02
1.159448E-01
1.079687E+01
2.494701E+01
-1.067271E-01
9.075462E-03
-2.571244E-02
4.402172E-02
-3.646387E-02
2.948533E-02
-1.043828E-01
1.600528E-02
4.361032E-02
3.875958E-03
-3.096545E-01
3.184604E-02
8.723730E-02
4.599654E-03
6.541966E-03
9.605866E-03
5.837907E-02
3.026784E-03
1.523399E-01
7.846764E-03
2.778305E-02
8.106935E-03
-2.895648E-01
1.983532E-02
4.866247E-02
1.048055E-02
2.600097E-02
1.465465E-02
-6.958338E-02
3.287519E-02
1.690742E-01
1.204714E-02
-1.127534E-01
7.881362E-03
-1.662306E-01
9.412057E-03
1.674047E-01
6.690449E-03
-4.370664E-02
5.828314E-03
9.323095E-02
4.400433E-03
1.353164E-02
1.389416E-03
-1.062009E-01
1.297730E-02
-5.331064E-02
5.321055E-03
7.698466E-02
3.247975E-03
-2.095250E-02
4.788058E-03
-8.861817E-03
7.657174E-03
1.769365E-02
1.144058E-03
1.788678E-02
2.746804E-03
4.403056E-03
6.918351E-04
1.059091E-01
1.394217E-02
1.442934E-01
1.318500E-02
2.199047E-02
5.838245E-03
8.175546E-02
3.242012E-03
7.359786E-03
1.673277E-02
5.252579E+01
5.951604E+02
-8.293789E-01
2.900775E-01
6.285324E-01
1.487344E+00
-4.765860E-01
8.058451E-01
-3.101266E-01
1.597231E-01
3.333449E-01
1.162003E-01
-1.314742E+00
4.565068E-01
4.726655E-01
1.161194E-01
2.824987E-01
8.732802E-02
1.025468E-01
4.935830E-02
4.216767E-01
1.225700E-01
-4.022783E-02
1.178843E-01
-7.967853E-01
2.698001E-01
3.741686E-01
2.849501E-01
3.377745E-01
3.722806E-01
9.234355E-02
4.498496E-01
8.019475E-01
3.404218E-01
-3.086622E-01
1.857604E-01
-8.474654E-01
1.733363E-01
8.124627E-01
1.415201E-01
-3.075481E-01
3.345921E-02
7.126178E-01
2.478330E-01
2.983972E-01
2.913106E-02
-6.759256E-01
4.245338E-01
-3.491214E-01
1.017784E-01
1.637602E-01
7.414687E-02
-3.344942E-01
1.165491E-01
-6.038784E-02
1.791698E-01
2.645086E-01
8.129446E-02
-2.738265E-01
6.616303E-02
1.135151E-02
4.140154E-02
4.679883E-01
1.914750E-01
2.876984E-01
1.355086E-01
-2.741613E-01
1.060457E-01
5.972204E-01
2.144153E-01
-9.215131E-02
4.815135E-01
3.325822E+01
2.414028E+02
-3.072951E-02
3.385901E-01
9.729550E-01
3.745599E-01
-5.237870E-01
3.101477E-01
-5.994810E-01
1.696123E-01
-3.402784E-01
1.910685E-01
-3.779571E-01
2.598004E-01
-4.804007E-01
1.417261E-01
-5.444032E-01
1.756801E-01
-9.921809E-02
6.889191E-02
-3.768065E-01
8.931168E-02
5.814513E-02
1.660714E-02
-1.263664E-01
2.650287E-01
1.951860E-01
2.090892E-01
3.576168E-01
2.511782E-01
-9.547988E-02
7.900044E-03
5.460746E-01
1.319868E-01
-1.577655E-01
8.575762E-02
2.938792E-02
7.941653E-02
-7.505627E-01
2.530444E-01
-2.399786E-01
1.105633E-01
2.035740E-01
2.442135E-02
1.190368E-01
8.652216E-02
5.753449E-01
9.353813E-02
4.911025E-01
6.795258E-02
-2.967023E-03
7.063512E-03
-2.158307E-01
4.680921E-02
-1.811238E-01
7.008450E-02
1.809785E-01
3.414843E-02
-1.380825E-01
3.469381E-02
3.722241E-01
3.956589E-02
2.064516E-01
2.676312E-02
-3.889056E-01
6.737109E-02
1.625815E-01
1.356347E-02
-5.278455E-01
1.297993E-01
-1.269304E-02
1.253029E-02
1.145800E+01
2.880575E+01
6.201234E-02
4.410457E-02
2.694694E-01
6.262339E-02
-1.173739E-01
4.727885E-02
-1.785056E-01
1.745391E-02
-4.744392E-02
2.117667E-02
3.998416E-02
6.058339E-02
-1.927536E-01
1.710512E-02
-1.800560E-01
2.420071E-02
-1.136912E-01
1.049357E-02
-1.104415E-01
1.181394E-02
3.486245E-02
3.673104E-03
-8.343502E-02
2.270299E-02
-3.064437E-02
1.197968E-02
1.912081E-01
2.808487E-02
2.601014E-02
1.311244E-03
1.295023E-01
1.108983E-02
-1.158397E-01
1.033182E-02
8.316701E-02
2.061897E-02
-1.136796E-01
3.133021E-02
-1.722931E-02
2.168654E-02
1.464014E-01
8.570818E-03
5.859338E-02
6.076029E-03
1.776223E-01
1.345571E-02
1.582455E-01
9.719314E-03
9.088118E-03
1.681557E-03
-6.298727E-02
5.287397E-03
5.793395E-03
6.397495E-03
-1.835123E-02
4.506573E-03
-7.558061E-02
4.462082E-03
3.136371E-02
3.181332E-03
-2.787527E-02
1.765822E-03
-9.878307E-02
3.381183E-03
1.279690E-01
4.386067E-03
-1.043523E-01
5.381331E-03
1.000162E-02
6.707280E-04
5.605671E+01
6.804062E+02
1.431465E-01
1.700768E+00
2.010680E+00
2.307036E+00
-5.276592E-01
1.006050E+00
-9.123132E-01
4.824313E-01
-3.276490E-01
5.842232E-01
6.749445E-02
6.049661E-01
-1.317515E+00
4.738014E-01
-9.470910E-01
4.937085E-01
-3.338536E-01
2.196164E-01
-3.032729E-01
1.241029E-01
-3.285269E-01
5.276607E-02
-4.107012E-01
6.402858E-01
3.357626E-01
2.672801E-01
1.000507E+00
5.263322E-01
-1.834202E-01
9.495434E-02
8.133484E-01
3.011488E-01
-4.083516E-01
2.428630E-01
4.032934E-01
1.825835E-01
-8.291267E-01
6.773284E-01
-2.401557E-01
2.280331E-01
4.325658E-01
1.163540E-01
6.671759E-02
2.181769E-01
6.984462E-01
1.716985E-01
7.108347E-01
1.386394E-01
1.709795E-01
5.320122E-02
-2.354809E-01
7.105363E-02
1.600239E-01
9.179874E-02
2.802111E-02
6.149546E-02
-4.005049E-01
2.195356E-01
7.914623E-01
2.273026E-01
3.392453E-01
1.101813E-01
-7.821576E-01
1.637944E-01
4.736782E-01
1.531480E-01
-5.253391E-01
1.093913E-01
2.235961E-01
5.449263E-02

10
tests/test_score_kappafission/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.135627E+02
9.364189E+03
1.848110E+02
7.563211E+03
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.782510E+02
6.811324E+03
2.035912E+02
8.999693E+03

10
tests/test_score_nufission/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
2.879098E+00
1.700626E+00
2.486342E+00
1.368793E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.401881E+00
1.235627E+00
2.733038E+00
1.616903E+00

14
tests/test_score_nuscatter/results_true.dat
View file

@ -1,11 +1,11 @@
k-combined:
1.093844E+00 1.626801E-02
1.005983E+00 2.248579E-02
tally 1:
0.000000E+00
0.000000E+00
1.239000E+01
3.172630E+01
3.570000E+00
2.668900E+00
4.450000E+01
4.100412E+02
1.169000E+01
2.915330E+01
3.200000E+00
2.342600E+00
4.064000E+01
3.595168E+02

Some files were not shown because too many files have changed in this diff Show more