tadd pressure set in jacobi loop - sphere - GPU-based 3D discrete element metho… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit ed357a684791e39fa86067c2e83825e288b7be6a | |
parent bebe373425e17094afa919530ad3c562435e06e0 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Fri, 19 Sep 2014 11:19:27 +0200 | |
add pressure set in jacobi loop | |
Diffstat: | |
M src/device.cu | 13 +++++++++++++ | |
1 file changed, 13 insertions(+), 0 deletions(-) | |
--- | |
diff --git a/src/device.cu b/src/device.cu | |
t@@ -1536,6 +1536,19 @@ __host__ void DEM::startTime() | |
&t_jacobiIterationNS); | |
checkForCudaErrorsIter("Post jacobiIterationNS", iter); | |
+ // set Dirichlet and Neumann BC at cells containing top wa… | |
+ if (walls.nw > 0 && walls.wmode[0] == 1) { | |
+ setNSepsilonAtTopWall<<<dimGridFluid, dimBlockFluid>>>( | |
+ dev_ns_epsilon, | |
+ dev_ns_epsilon_new, | |
+ wall0_iz, | |
+ epsilon_value, | |
+ dp_dz); | |
+ cudaThreadSynchronize(); | |
+ checkForCudaErrorsIter("Post setNSepsilonAtTopWall", | |
+ iter); | |
+ } | |
+ | |
// Flip flop: swap new and current array pointers | |
/*Float* tmp = dev_ns_epsilon; | |
dev_ns_epsilon = dev_ns_epsilon_new; |