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

Debugging TestCommD3Q19

Browse Source
This commit is contained in:
James E McClure 2015-06-25 15:45:25 -04:00
parent 4dc71449a7
commit 79384a42bb
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
gpu/D3Q19.cu
View File

@ -102,6 +102,7 @@ __global__ void dvc_SwapD3Q19(char *ID, double *disteven, double *distodd, int
{ {
int i,j,k,n,nn,N; int i,j,k,n,nn,N;
// distributions // distributions
char id;
double f1,f2,f3,f4,f5,f6,f7,f8,f9; double f1,f2,f3,f4,f5,f6,f7,f8,f9;
double f10,f11,f12,f13,f14,f15,f16,f17,f18; double f10,f11,f12,f13,f14,f15,f16,f17,f18;
@ -111,7 +112,9 @@ __global__ void dvc_SwapD3Q19(char *ID, double *disteven, double *distodd, int
for (int s=0; s<S; s++){ for (int s=0; s<S; s++){
//........Get 1-D index for this thread.................... //........Get 1-D index for this thread....................
n = S*blockIdx.x*blockDim.x + s*blockDim.x + threadIdx.x; 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){
//.......Back out the 3-D indices for node n.............. //.......Back out the 3-D indices for node n..............
k = n/(Nx*Ny); k = n/(Nx*Ny);
j = (n-Nx*Ny*k)/Nx; j = (n-Nx*Ny*k)/Nx;
@ -231,6 +234,7 @@ __global__ void dvc_SwapD3Q19(char *ID, double *disteven, double *distodd, int
//........................................................................ //........................................................................
} }
}
} }
} }