diff --git a/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H b/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H index 99c2b0b89df..d36eaceaebd 100644 --- a/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H +++ b/Source/Particles/Collision/BinaryCollision/Bremsstrahlung/BremsstrahlungFunc.H @@ -298,9 +298,9 @@ public: // k will be between k_im1 and k_i amrex::ParticleReal const f_im1 = kdsigdk_im1/koT1_grid_im1; - amrex::ParticleReal const fi = kdsigdk_i/m_koT1_grid[i]; - amrex::ParticleReal const x = (std::sqrt(f_im1*f_im1 + 2._prt*(fi - f_im1)*(random_number*sigma_total - sigma)/dk) - - f_im1)/(fi - f_im1); + amrex::ParticleReal const f_i = kdsigdk_i/m_koT1_grid[i]; + amrex::ParticleReal const x = (std::sqrt(f_im1*f_im1 + 2._prt*(f_i - f_im1)*(random_number*sigma_total - sigma)/dk) + - f_im1)/(f_i - f_im1); amrex::ParticleReal const result = x*dk + koT1_grid_im1; return result;