From 197425408daa023f6383fae132c5233debe8d462 Mon Sep 17 00:00:00 2001
From: Henry <henry.korb@geo.uu.se>
Date: Thu, 5 Oct 2023 17:42:47 +0200
Subject: [PATCH] bug fix in collision kernel
---
.../LBM/GPUHelperFunctions/RunCollision.cuh | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh b/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
index ba1034ac5..6bb51f95c 100644
--- a/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
+++ b/src/gpu/VirtualFluids_GPU/LBM/GPUHelperFunctions/RunCollision.cuh
@@ -4,12 +4,15 @@
#include <cuda_runtime.h>
#include <basics/DataTypes.h>
+#include <basics/constants/NumericConstants.h>
#include <lbm/collision/CollisionParameter.h>
#include <lbm/collision/TurbulentViscosity.h>
#include "LBM/GPUHelperFunctions/KernelUtilities.h"
+using namespace vf::basics::constant;
+
namespace vf::gpu
{
@@ -58,9 +61,9 @@ __global__ void runCollision(CollisionFunctor collision, GPUCollisionParameter c
para.forceZ = (collisionParameter.forces[2] + collisionParameter.bodyForceZ[k_000]) * c1o2 * collisionParameter.forceFactor;
// Reset body force. To be used when not using round-off correction.
- collisionParameter.bodyForceX[k_000] = 0.0f;
- collisionParameter.bodyForceX[k_000] = 0.0f;
- collisionParameter.bodyForceX[k_000] = 0.0f;
+ collisionParameter.bodyForceX[k_000] = c0o1;
+ collisionParameter.bodyForceY[k_000] = c0o1;
+ collisionParameter.bodyForceZ[k_000] = c0o1;
////////////////////////////////////////////////////////////////////////////////////
//!> Round-off correction
--
GitLab