tchange sign on every 2nd value - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
(HTM) git clone git://src.adamsgaard.dk/sphere
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) LICENSE
---
(DIR) commit 4ac1331a468e9543d1b7ec944b75aebe8d229065
(DIR) parent ddf90fec75f449401838ada68084b5d2c4aea6a5
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Fri, 20 Mar 2015 11:40:42 +0100
change sign on every 2nd value
Diffstat:
M src/darcy.cuh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
t@@ -449,7 +449,7 @@ __global__ void findDarcyPorositiesLinear(
s = weight(x3,
X + MAKE_FLOAT3( 0.5*dx, 0.0, 0.0),
dx, dy, dz);
- v_p_xp += s*vol_p*v3.x/(s*vol_p + 1.0e-16);
+ v_p_xp += s*vol_p*v3.x/(s*vol_p - 1.0e-16);
s = weight(x3,
X + MAKE_FLOAT3( 0.0, -0.5*dy, 0.0),
t@@ -459,7 +459,7 @@ __global__ void findDarcyPorositiesLinear(
s = weight(x3,
X + MAKE_FLOAT3( 0.0, 0.5*dy, 0.0),
dx, dy, dz);
- v_p_yp += s*vol_p*v3.y/(s*vol_p + 1.0e-16);
+ v_p_yp += s*vol_p*v3.y/(s*vol_p - 1.0e-16);
s = weight(x3,
X + MAKE_FLOAT3( 0.0, 0.0, -0.5*dz),
t@@ -469,7 +469,7 @@ __global__ void findDarcyPorositiesLinear(
s = weight(x3,
X + MAKE_FLOAT3( 0.0, 0.0, 0.5*dz),
dx, dy, dz);
- v_p_zp += s*vol_p*v3.z/(s*vol_p + 1.0e-16);
+ v_p_zp += s*vol_p*v3.z/(s*vol_p - 1.0e-16);
}
}
}