crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Make sure to use std::fabs

Browse source
This commit is contained in:
Paul Romano 2018-09-04 07:26:57 -05:00
parent 8af84064ca
commit 1b42fa20c6
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/mesh.cpp
View file

@ -432,9 +432,9 @@ void RegularMesh::bins_crossed(const Particle* p, std::vector<int>& bins,
}
for (j = 0; j < n; ++j) {
if (std::abs(u[j]) < FP_PRECISION) {
if (std::fabs(u[j]) < FP_PRECISION) {
d(j) = INFTY;
} else if (u[j] > 0) {
} else if (u[j] > 0.0) {
double xyz_cross = lower_left_[j] + ijk0(j) * width_[j];
d(j) = (xyz_cross - r0[j]) / u[j];
} else {
@ -473,7 +473,7 @@ void RegularMesh::bins_crossed(const Particle* p, std::vector<int>& bins,
// surface.
xt::xtensor<double, 1> d = xt::zeros<double>({n});
for (int j = 0; j < n; ++j) {
if (abs(u[j]) < FP_PRECISION) {
if (std::fabs(u[j]) < FP_PRECISION) {
d(j) = INFTY;
} else if (u[j] > 0) {
double xyz_cross = lower_left_[j] + ijk0(j) * width_[j];