tadd Neumann BCs to iterative solution - sphere - GPU-based 3D discrete element… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 833df6a19e3c16b0bf1e1cbbe760841b03ec3c54 | |
parent 04e358abff4061084ea4a73a0652eae4dfa42af1 | |
Author: Anders Damsgaard Christensen <[email protected]> | |
Date: Thu, 11 Aug 2016 11:43:05 -0700 | |
add Neumann BCs to iterative solution | |
Diffstat: | |
M src/darcy.cuh | 9 ++++++++- | |
1 file changed, 8 insertions(+), 1 deletion(-) | |
--- | |
diff --git a/src/darcy.cuh b/src/darcy.cuh | |
t@@ -1785,9 +1785,16 @@ __global__ void updateDarcySolution( | |
Float p_zp = dev_darcy_p[d_idx(x,y,z+1)]; | |
// Neumann BCs | |
+ if (x == 0 && bc_xn == 1) | |
+ p_xn = p; | |
+ if (x == nx-1 && bc_xp == 1) | |
+ p_xp = p; | |
+ if (y == 0 && bc_yn == 1) | |
+ p_yn = p; | |
+ if (y == ny-1 && bc_yp == 1) | |
+ p_yp = p; | |
if (z == 0 && bc_bot == 1) | |
p_zn = p; | |
- //if ((z == nz-1 && bc_top == 1) || z >= wall0_iz) | |
if (z == nz-1 && bc_top == 1) | |
p_zp = p; | |