Fix extremely large yields from Bremsstrahlung (#3386)

Co-authored-by: Paul Romano <paul.k.romano@gmail.com>
This commit is contained in:
Hunter Belanger 2025-05-03 05:10:18 +00:00 committed by GitHub
parent d0354dbd71
commit a921280fa9
No known key found for this signature in database
GPG key ID: B5690EEEBB952194

View file

@ -640,7 +640,7 @@ void Material::init_bremsstrahlung()
// Allocate arrays for TTB data
ttb->pdf = xt::zeros<double>({n_e, n_e});
ttb->cdf = xt::zeros<double>({n_e, n_e});
ttb->yield = xt::empty<double>({n_e});
ttb->yield = xt::zeros<double>({n_e});
// Allocate temporary arrays
xt::xtensor<double, 1> stopping_power_collision({n_e}, 0.0);