diff --git a/openmc/data/thermal.py b/openmc/data/thermal.py index bc2e1d9819..0366af1a66 100644 --- a/openmc/data/thermal.py +++ b/openmc/data/thermal.py @@ -19,12 +19,12 @@ from . import HDF5_VERSION, HDF5_VERSION_MAJOR, endf from .data import K_BOLTZMANN, ATOMIC_SYMBOL, EV_PER_MEV, NATURAL_ABUNDANCE from .ace import Table, get_table, Library from .angle_energy import AngleEnergy -from .correlated import CorrelatedAngleEnergy from .function import Tabulated1D, Function1D from .njoy import make_ace_thermal from .thermal_angle_energy import (CoherentElasticAE, IncoherentElasticAE, IncoherentElasticAEDiscrete, - IncoherentInelasticAEDiscrete) + IncoherentInelasticAEDiscrete, + IncoherentInelasticAEContinuous) _THERMAL_NAMES = { @@ -673,7 +673,7 @@ class ThermalScattering(EqualityMixin): breakpoints = [n_energy] interpolation = [2] energy = inelastic_xs.x - distribution = CorrelatedAngleEnergy( + distribution = IncoherentInelasticAEContinuous( breakpoints, interpolation, energy, energy_out, mu_out) table.inelastic = ThermalScatteringReaction( diff --git a/openmc/data/thermal_angle_energy.py b/openmc/data/thermal_angle_energy.py index cf9587cbe8..e5b77b597a 100644 --- a/openmc/data/thermal_angle_energy.py +++ b/openmc/data/thermal_angle_energy.py @@ -1,6 +1,7 @@ import numpy as np from .angle_energy import AngleEnergy +from .correlated import CorrelatedAngleEnergy class CoherentElasticAE(AngleEnergy): @@ -134,3 +135,6 @@ class IncoherentInelasticAEDiscrete(AngleEnergy): mu_out = group['mu_out'][()] skewed = bool(group['skewed']) return cls(energy_out, mu_out, skewed) + +class IncoherentInelasticAEContinuous(CorrelatedAngleEnergy): + pass