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

Correction to maxwell dist sampling

Browse source
This commit is contained in:
Patrick Shriwise 2022-04-05 13:24:18 -05:00
parent 7484618c7c
commit e856c6a270
1 changed files with 2 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

4
openmc/stats/univariate.py
View file

@ -452,13 +452,13 @@ class Maxwell(Univariate):
def sample(self, n_samples=1, seed=None):
np.random.seed(seed)
self.sample_maxwell(self.theta, n_samples)
return self.sample_maxwell(self.theta, n_samples)
@staticmethod
def sample_maxwell(t, n_samples):
r1, r2, r3 = np.random.rand(3, n_samples)
c = np.power(np.cos(0.5 * np.pi * r3), 2)
return t * (np.log(r1) + np.log(r2) * c)
return -t * (np.log(r1) + np.log(r2) * c)
def to_xml_element(self, element_name):
"""Return XML representation of the Maxwellian distribution