From c639aa34a8729b48295d8c03f44c9d8bf011024a Mon Sep 17 00:00:00 2001
From: Lizzzka007 <gashchuk2011@mail.ru>
Date: Thu, 17 Oct 2024 16:51:15 +0300
Subject: [PATCH] Arithmetic fabs fix
---
includeCU/sfx-model-compute-subfunc.cuh | 10 +++++-----
srcCU/sfx-sheba.cu | 2 +-
srcCXX/sfx-sheba.cpp | 2 +-
3 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/includeCU/sfx-model-compute-subfunc.cuh b/includeCU/sfx-model-compute-subfunc.cuh
index 0b45bf9..9bf794e 100644
--- a/includeCU/sfx-model-compute-subfunc.cuh
+++ b/includeCU/sfx-model-compute-subfunc.cuh
@@ -18,7 +18,7 @@ FUCNTION_DECLARATION_SPECIFIER void get_convection_lim(T &zeta_lim, T &Rib_lim,
f_m = zeta_lim / h0_m;
f_h = zeta_lim / h0_t;
- if (fabs(B) < 1.0e-10) f_h = f_m;
+ if (fabsf(B) < 1.0e-10) f_h = f_m;
f_m = powf(1.0 - param.alpha_m * f_m, 0.25);
f_h = sqrtf(1.0 - param.alpha_h_fix * f_h);
@@ -67,7 +67,7 @@ FUCNTION_DECLARATION_SPECIFIER void get_psi_convection(T &psi_m, T &psi_h, T &ze
{
zeta0_m = zeta / h0_m;
zeta0_h = zeta / h0_t;
- if (fabs(B) < 1.0e-10)
+ if (fabsf(B) < 1.0e-10)
zeta0_h = zeta0_m;
f0_m = powf(1.0 - param.alpha_m * zeta0_m, 0.25);
@@ -97,7 +97,7 @@ FUCNTION_DECLARATION_SPECIFIER void get_psi_neutral(T &psi_m, T &psi_h, T &zeta,
zeta = 0.0;
psi_m = logf(h0_m);
psi_h = logf(h0_t) / param.Pr_t_0_inv;
- if (fabs(B) < 1.0e-10)
+ if (fabsf(B) < 1.0e-10)
psi_h = psi_m / param.Pr_t_0_inv;
}
@@ -112,7 +112,7 @@ FUCNTION_DECLARATION_SPECIFIER void get_psi_semi_convection(T &psi_m, T &psi_h,
psi_m = logf(h0_m);
psi_h = logf(h0_t);
- if (fabs(B) < 1.0e-10)
+ if (fabsf(B) < 1.0e-10)
psi_h = psi_m;
zeta = Rib * param.Pr_t_0_inv * psi_m * psi_m / psi_h;
@@ -121,7 +121,7 @@ FUCNTION_DECLARATION_SPECIFIER void get_psi_semi_convection(T &psi_m, T &psi_h,
{
zeta0_m = zeta / h0_m;
zeta0_h = zeta / h0_t;
- if (fabs(B) < 1.0e-10)
+ if (fabsf(B) < 1.0e-10)
zeta0_h = zeta0_m;
f_m = powf(1.0 - param.alpha_m * zeta, 0.25e0);
diff --git a/srcCU/sfx-sheba.cu b/srcCU/sfx-sheba.cu
index cb00722..3fff68a 100644
--- a/srcCU/sfx-sheba.cu
+++ b/srcCU/sfx-sheba.cu
@@ -79,7 +79,7 @@ __global__ void sfx_kernel::compute_flux(sfxDataVecTypeC sfx,
if (U > 0.0)
Cm = Udyn / U;
Ct = 0.0;
- if (fabs(dT) > 0.0)
+ if (fabsf(dT) > 0.0)
Ct = Tdyn / dT;
// --- define eddy viscosity & inverse Prandtl number
diff --git a/srcCXX/sfx-sheba.cpp b/srcCXX/sfx-sheba.cpp
index 7901bfd..bf389dc 100644
--- a/srcCXX/sfx-sheba.cpp
+++ b/srcCXX/sfx-sheba.cpp
@@ -82,7 +82,7 @@ void FluxSheba<T, memIn, memOut, MemType::CPU>::compute_flux()
if (U > 0.0)
Cm = Udyn / U;
Ct = 0.0;
- if (fabs(dT) > 0.0)
+ if (fabsf(dT) > 0.0)
Ct = Tdyn / dT;
// --- define eddy viscosity & inverse Prandtl number
--
GitLab