From b4550e32de76f66b6d26765d34b89d81500af596 Mon Sep 17 00:00:00 2001
From: zy69guqi <richard.angersbach@fau.de>
Date: Thu, 24 Aug 2023 15:51:19 +0200
Subject: [PATCH] Fix cast to int in collision stress calculation
---
apps/showcases/ChargedParticles/ChargedParticles.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/apps/showcases/ChargedParticles/ChargedParticles.cpp b/apps/showcases/ChargedParticles/ChargedParticles.cpp
index 251e38553..2b2843683 100644
--- a/apps/showcases/ChargedParticles/ChargedParticles.cpp
+++ b/apps/showcases/ChargedParticles/ChargedParticles.cpp
@@ -1211,9 +1211,9 @@ int main(int argc, char** argv)
ps->forEachParticle(
useOpenMP, mesa_pd::kernel::SelectLocal(), *accessor,
[&collisionStress, &diameter](const size_t idx1, auto& ac) {
- collisionStress[0] += (abs(ac.getForce(idx1)[0])) / ((math::pi) *diameter * diameter);
- collisionStress[1] += (abs(ac.getForce(idx1)[1])) / ((math::pi) *diameter * diameter);
- collisionStress[2] += (abs(ac.getForce(idx1)[2])) / ((math::pi) *diameter * diameter);
+ collisionStress[0] += real_c(fabs(ac.getForce(idx1)[0])) / ((math::pi) * diameter * diameter);
+ collisionStress[1] += real_c(fabs(ac.getForce(idx1)[1])) / ((math::pi) * diameter * diameter);
+ collisionStress[2] += real_c(fabs(ac.getForce(idx1)[2])) / ((math::pi) * diameter * diameter);
},
*accessor);
}
--
GitLab