From 69855169b9c007953963b88602d4c45eec05c133 Mon Sep 17 00:00:00 2001
From: Henry <henry.korb@geo.uu.se>
Date: Thu, 12 Oct 2023 15:48:14 +0200
Subject: [PATCH] fix wrong type

---
 src/lbm/collision/TurbulentViscosity.h | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/src/lbm/collision/TurbulentViscosity.h b/src/lbm/collision/TurbulentViscosity.h
index 0122dda77..8ec8130b7 100644
--- a/src/lbm/collision/TurbulentViscosity.h
+++ b/src/lbm/collision/TurbulentViscosity.h
@@ -71,7 +71,7 @@ inline __host__ __device__ real calcTurbulentViscositySmagorinsky(real Cs, real
 }
 
 template <typename T>
-__host__ __device__ int max( T a, T b )
+__host__ __device__ T max( T a, T b )
 {
     return ( a > b ) ? a : b;
 }
@@ -87,6 +87,11 @@ inline __host__ __device__ real calcTurbulentViscosityQR(real C, real dxux, real
     return C * max(R, c0o1) / Q;
 }
 
+inline __host__ __device__ real calculateOmegaWithturbulentViscosity(const real omega_old, real turbulenceViscosity)
+{
+    return omega_old / (c1o1 + c3o1 * omega_old * turbulenceViscosity);
+}
+
 } // namespace vf::lbm
 
 #endif //TURBULENT_VISCOSITY_H
-- 
GitLab