From 19e43ed1022e24c482867c579709d15a2a5d7ff7 Mon Sep 17 00:00:00 2001 From: GuySten Date: Sun, 10 May 2026 18:23:54 +0300 Subject: [PATCH] wip --- src/angle_energy.cpp | 6 +++--- src/tallies/next_event_scoring.cpp | 5 ++--- 2 files changed, 5 insertions(+), 6 deletions(-) diff --git a/src/angle_energy.cpp b/src/angle_energy.cpp index abbe1af896..1157b4c33c 100644 --- a/src/angle_energy.cpp +++ b/src/angle_energy.cpp @@ -10,11 +10,11 @@ double get_jac_and_transform( double E_in, double& mu, double& E_out, uint64_t* seed, double awr) { double E_com = E_in / ((awr + 1.0) * (awr + 1.0)); - return get_jac_and_transform(E_in, mu, E_out, seed, awr, E_com); + return get_jac_and_transform_impl(E_in, mu, E_out, seed, awr, E_com); } -double get_jac_and_transform(double E_in, double& mu, double& E_out, - uint64_t* seed, double awr, double E_com) +double get_jac_and_transform_impl( + double E_com, double& mu, double& E_out, uint64_t* seed, double awr) { double E_cm = E_out; double mu_lab = mu; diff --git a/src/tallies/next_event_scoring.cpp b/src/tallies/next_event_scoring.cpp index 185c7c537c..0a235f729c 100644 --- a/src/tallies/next_event_scoring.cpp +++ b/src/tallies/next_event_scoring.cpp @@ -17,11 +17,10 @@ void score_point_tally_elastic( // Velocity of center-of-mass Direction v_cm = (v_n + awr * v_t) / (awr + 1.0); + auto u_cm = v_cm / v_cm.norm(); double E_com = v_cm.dot(v_cm); double E_out = (v_n - v_cm).dot(v_n - v_cm); - double E_in = p.E(); - auto u_cm = v_cm / v_cm.norm(); auto& d = rx.products_[i_product].distribution_[0]; auto d_ = dynamic_cast(d.get()); @@ -30,7 +29,7 @@ void score_point_tally_elastic( double mu = u.dot(u_cm); E = E_out; double jac = - get_jac_and_transform(E_in, mu, E, p.current_seed(), awr, E_com); + get_jac_and_transform_impl(E_com, mu, E, p.current_seed(), awr); if (!d_->angle().empty()) { return jac * d_->angle().evaluate(p.E(), mu) / (2.0 * PI); } else {