From 9d349a7c709749635e6113d80b0011db45e9abd0 Mon Sep 17 00:00:00 2001
From: zy69guqi <richard.angersbach@fau.de>
Date: Thu, 24 Aug 2023 15:58:53 +0200
Subject: [PATCH] Fix cast to int in hydrodynamic 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 2b2843683..f76d200b5 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);
}
--
GitLab