19795099684

Committed 30 Nov 2025 06:38AM UTC coverage: 81.545% (-0.5%) from 82.055%

Build # 19795099684

Build Type

Pull #3339

github

Committed by

web-flow

Commit Message

Merge 47e5a263d into ef22558f4

Pull Request Pull Request #3339: Enable Scikit-Build-Core Support

Run Details

16899 of 23552 branches covered (71.75%)

Branch coverage included in aggregate %.

34 of 34 new or added lines in 3 files covered. (100.0%)

383 existing lines in 8 files now uncovered.

54550 of 64067 relevant lines covered (85.15%)

42133974.15 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

56.0

/src/random_dist.cpp

#include "openmc/random_dist.h"

#include <cmath>

#include "openmc/constants.h"
#include "openmc/random_lcg.h"

namespace openmc {

double uniform_distribution(double a, double b, uint64_t* seed)
{
  return a + (b - a) * prn(seed);
}

int64_t uniform_int_distribution(int64_t a, int64_t b, uint64_t* seed)
{
  return a + static_cast<int64_t>(prn(seed) * (b - a + 1));
}

double maxwell_spectrum(double T, uint64_t* seed)
{
  // Set the random numbers
  double r1 = prn(seed);
  double r2 = prn(seed);
  double r3 = prn(seed);

  // determine cosine of pi/2*r
  double c = std::cos(PI / 2. * r3);

  // Determine outgoing energy
  return -T * (std::log(r1) + std::log(r2) * c * c);
}

double watt_spectrum(double a, double b, uint64_t* seed)
{
  double w = maxwell_spectrum(a, seed);
  return w + 0.25 * a * a * b +
         uniform_distribution(-1., 1., seed) * std::sqrt(a * a * b * w);
}

double normal_variate(double mean, double standard_deviation, uint64_t* seed)
{
  // Sample a normal variate using Marsaglia's polar method
  double x, y, r2;
  do {
    x = uniform_distribution(-1., 1., seed);
    y = uniform_distribution(-1., 1., seed);
    r2 = x * x + y * y;
  } while (r2 > 1 || r2 == 0);
  double z = std::sqrt(-2.0 * std::log(r2) / r2);
  return mean + standard_deviation * z * x;
}

} // namespace openmc

1	#include "openmc/random_dist.h"
2
3	#include <cmath>
4
5	#include "openmc/constants.h"
6	#include "openmc/random_lcg.h"
7
8	namespace openmc {
9
10	double uniform_distribution(double a, double b, uint64_t* seed)	681,334,984✔
11	{
12	return a + (b - a) * prn(seed);	681,334,984✔
13	}
14
15	int64_t uniform_int_distribution(int64_t a, int64_t b, uint64_t* seed)	5,247,000✔
16	{
17	return a + static_cast<int64_t>(prn(seed) * (b - a + 1));	5,247,000✔
18	}
19
20	double maxwell_spectrum(double T, uint64_t* seed)	13,542,519✔
21	{
22	// Set the random numbers
23	double r1 = prn(seed);	13,542,519✔
24	double r2 = prn(seed);	13,542,519✔
25	double r3 = prn(seed);	13,542,519✔
26
27	// determine cosine of pi/2*r
28	double c = std::cos(PI / 2. * r3);	13,542,519✔
29
30	// Determine outgoing energy
31	return -T * (std::log(r1) + std::log(r2) * c * c);	13,542,519✔
32	}
33
34	double watt_spectrum(double a, double b, uint64_t* seed)	13,394,283✔
35	{
36	double w = maxwell_spectrum(a, seed);	13,394,283✔
37	return w + 0.25 * a * a * b +	13,394,283✔
38	uniform_distribution(-1., 1., seed) * std::sqrt(a * a * b * w);	13,394,283✔
39	}
40
UNCOV 41	double normal_variate(double mean, double standard_deviation, uint64_t* seed)	×
42	{
43	// Sample a normal variate using Marsaglia's polar method
44	double x, y, r2;
45	do {
UNCOV 46	x = uniform_distribution(-1., 1., seed);	×
UNCOV 47	y = uniform_distribution(-1., 1., seed);	×
UNCOV 48	r2 = x * x + y * y;	×
UNCOV 49	} while (r2 > 1 \|\| r2 == 0);	×
UNCOV 50	double z = std::sqrt(-2.0 * std::log(r2) / r2);	×
UNCOV 51	return mean + standard_deviation * z * x;	×
52	}
53
54	} // namespace openmc

openmc-dev / openmc / 19795099684

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous