tdecrease factor to be added - 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 ddf90fec75f449401838ada68084b5d2c4aea6a5
(DIR) parent 8170eb9480b254ce94b2d0d2fefa461f386cf335
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Fri, 20 Mar 2015 11:40:03 +0100
decrease factor to be added
Diffstat:
M src/darcy.cuh | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
---
(DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
t@@ -444,32 +444,32 @@ __global__ void findDarcyPorositiesLinear(
s = weight(x3,
X + MAKE_FLOAT3(-0.5*dx, 0.0, 0.0),
dx, dy, dz);
- v_p_xn += s*vol_p*v3.x / (s*vol_p + 1.0e-8);
+ v_p_xn += s*vol_p*v3.x/(s*vol_p + 1.0e-16);
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-8);
+ 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),
dx, dy, dz);
- v_p_yn += s*vol_p*v3.y / (s*vol_p + 1.0e-8);
+ v_p_yn += s*vol_p*v3.y/(s*vol_p + 1.0e-16);
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-8);
+ 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),
dx, dy, dz);
- v_p_zn += s*vol_p*v3.z / (s*vol_p + 1.0e-8);
+ v_p_zn += s*vol_p*v3.z/(s*vol_p + 1.0e-16);
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-8);
+ v_p_zp += s*vol_p*v3.z/(s*vol_p + 1.0e-16);
}
}
}