diff --git a/apps/showcases/ChargedParticles/ChargedParticles.cpp b/apps/showcases/ChargedParticles/ChargedParticles.cpp index 2b284368336f71795c23f31c9b955846fd30a152..f76d200b5de4b83e9a86befef5fbba580f972e36 100644 --- a/apps/showcases/ChargedParticles/ChargedParticles.cpp +++ b/apps/showcases/ChargedParticles/ChargedParticles.cpp @@ -1277,11 +1277,11 @@ int main(int argc, char** argv) useOpenMP, mesa_pd::kernel::SelectLocal(), *accessor, [&hydrodynamicStress, diameter, gravitationalForce](const size_t idx1, auto& ac) { hydrodynamicStress[0] += - (abs(ac.getHydrodynamicForce(idx1)[0] + gravitationalForce[0])) / ((math::pi) *diameter * diameter); + real_c(fabs(ac.getHydrodynamicForce(idx1)[0] + gravitationalForce[0])) / ((math::pi) * diameter * diameter); hydrodynamicStress[1] += - (abs(ac.getHydrodynamicForce(idx1)[1] + gravitationalForce[1])) / ((math::pi) *diameter * diameter); + real_c(fabs(ac.getHydrodynamicForce(idx1)[1] + gravitationalForce[1])) / ((math::pi) * diameter * diameter); hydrodynamicStress[2] += - (abs(ac.getHydrodynamicForce(idx1)[2] + gravitationalForce[2])) / ((math::pi) *diameter * diameter); + real_c(fabs(ac.getHydrodynamicForce(idx1)[2] + gravitationalForce[2])) / ((math::pi) * diameter * diameter); }, *accessor); }