OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Modified ComputeVelocityD3Q19 so zero is assigned inside solid

Browse Source
This commit is contained in:
James E McClure 2015-06-26 21:41:19 -04:00
parent 908b7fbb89
commit c8cbe211fe
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
gpu/D3Q19.cu
View File

@ -251,14 +251,19 @@ __global__ void dvc_ComputeVelocityD3Q19(char *ID, double *disteven, double *di
double f1,f2,f3,f4,f5,f6,f7,f8,f9;
double f10,f11,f12,f13,f14,f15,f16,f17,f18;
double vx,vy,vz;
char id;
N = Nx*Ny*Nz;
int S = N/NBLOCKS/NTHREADS + 1;
for (int s=0; s<S; s++){
//........Get 1-D index for this thread....................
n = S*blockIdx.x*blockDim.x + s*blockDim.x + threadIdx.x;
if (n<N && ID[n] > 0){
if (n<N){
id = ID[n];
if (id==0){
vel[n] = 0.0; vel[N+n] = 0.0; vel[2*N+n]=0.0;
}
else{
//........................................................................
// Registers to store the distributions
//........................................................................
@ -290,7 +295,7 @@ __global__ void dvc_ComputeVelocityD3Q19(char *ID, double *disteven, double *di
vel[N+n] = vy;
vel[2*N+n] = vz;
//........................................................................
}
}
}
}