25631610076

Committed 10 May 2026 02:44PM UTC coverage: 81.026% (-0.4%) from 81.388%

Build Type

github

Committed by

Commit Message

Merge debc5a921 into d56cda254

Coverage Stats

17769 of 25812 branches covered (68.84%)

Branch coverage included in aggregate %.

51 of 373 new or added lines in 25 files covered. (13.67%)

3 existing lines in 2 files now uncovered.

58805 of 68694 relevant lines covered (85.6%)

40257430.53 hits per line

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

0.0

/src/angle_energy.cpp

1	#include "openmc/angle_energy.h"
2
3	#include <cmath> // for sqrt
4
5	#include "openmc/random_lcg.h"
6
7	namespace openmc {
8
NEW 9	double get_jac_and_transform(	×
10	double E_in, double& mu, double& E_out, uint64_t* seed, double awr)
11	{
NEW 12	double E_com = E_in / ((awr + 1.0) * (awr + 1.0));	×
NEW 13	return get_jac_and_transform(E_in, mu, E_out, seed, awr, E_com);	×
14	}
15
NEW 16	double get_jac_and_transform(double E_in, double& mu, double& E_out,	×
17	uint64_t* seed, double awr, double E_com)
18	{
NEW 19	double E_cm = E_out;	×
NEW 20	double mu_lab = mu;	×
NEW 21	double D = mu_lab * mu_lab - 1.0 + E_cm / E_com;	×
NEW 22	if (D <= 0.0)	×
23	return 0.0;
NEW 24	D = std::sqrt(D);	×
25
NEW 26	if (mu_lab <= 0.0 && E_cm <= E_cm)	×
27	return 0.0;
NEW 28	double E_out1 = E_com * (mu_lab + D) * (mu_lab + D);	×
NEW 29	double mult;	×
NEW 30	if (E_cm > E_com) {	×
NEW 31	mult = 1.0;	×
NEW 32	E_out = E_out1;	×
33	} else {
NEW 34	mult = 2.0;	×
NEW 35	if (prn(seed) < 0.5) {	×
NEW 36	E_out = E_out1;	×
37	} else {
NEW 38	E_out = E_com * (mu_lab - D) * (mu_lab - D);	×
39	}
40	}
NEW 41	mu = mu_lab * std::sqrt(E_out / E_cm) - std::sqrt(E_com / E_cm);	×
NEW 42	return mult * E_out / (D * std::sqrt(E_cm * E_com));	×
43	}
44
45	} // namespace openmc