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

Try using 10**x instead of logspace in source test

Browse source
This commit is contained in:
Paul Romano 2023-02-13 09:00:33 -06:00
parent 4783b34fe2
commit 893c4d0648
1 changed files with 8 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

13
tests/regression_tests/source/test.py
View file

@ -39,15 +39,15 @@ class SourceTestHarness(PyAPITestHarness):
spatial2 = openmc.stats.Box([-4., -4., -4.], [4., 4., 4.])
spatial3 = openmc.stats.Point([1.2, -2.3, 0.781])
spatial4 = openmc.stats.SphericalIndependent(r_dist, cos_theta_dist,
phi_dist,
phi_dist,
origin=(1., 1., 0.))
spatial5 = openmc.stats.CylindricalIndependent(r_dist, phi_dist,
spatial5 = openmc.stats.CylindricalIndependent(r_dist, phi_dist,
z_dist,
origin=(1., 1., 0.))
spatial6 = openmc.stats.SphericalIndependent(r_dist2, cos_theta_dist,
phi_dist,
phi_dist,
origin=(1., 1., 0.))
spatial7 = openmc.stats.CylindricalIndependent(r_dist1, phi_dist,
spatial7 = openmc.stats.CylindricalIndependent(r_dist1, phi_dist,
z_dist,
origin=(1., 1., 0.))
@ -57,7 +57,10 @@ class SourceTestHarness(PyAPITestHarness):
angle2 = openmc.stats.Monodirectional(reference_uvw=[0., 1., 0.])
angle3 = openmc.stats.Isotropic()
E = np.logspace(0, 7)
# Note that the definition for E is equivalent to logspace(0, 7) but we
# manually take powers because of last-digit differences that may cause
# test failures with different versions of numpy
E = np.array([10**x for x in np.linspace(0, 7)])
p = np.sin(np.linspace(0., pi))
p /= sum(np.diff(E)*p[:-1])
energy1 = openmc.stats.Maxwell(1.2895e6)