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

Modified ComputePressureD3Q19 so zero is assigned inside solid

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

11
gpu/D3Q19.cu
View File

@ -302,15 +302,18 @@ __global__ void dvc_ComputePressureD3Q19(char *ID, double *disteven, double *di
// distributions // distributions
double f0,f1,f2,f3,f4,f5,f6,f7,f8,f9; double f0,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;
chard id;
N = Nx*Ny*Nz; N = Nx*Ny*Nz;
int S = N/NBLOCKS/NTHREADS + 1; int S = N/NBLOCKS/NTHREADS + 1;
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) Pressure[n] = 0.0;
else{
//.......................................................................
// Registers to store the distributions // Registers to store the distributions
//........................................................................ //........................................................................
f0 = disteven[n]; f0 = disteven[n];
@ -336,7 +339,7 @@ __global__ void dvc_ComputePressureD3Q19(char *ID, double *disteven, double *di
//.................Compute the velocity................................... //.................Compute the velocity...................................
Pressure[n] = 0.3333333333333333*(f0+f2+f1+f4+f3+f6+f5+f8+f7+f10+ Pressure[n] = 0.3333333333333333*(f0+f2+f1+f4+f3+f6+f5+f8+f7+f10+
f9+f12+f11+f14+f13+f16+f15+f18+f17); f9+f12+f11+f14+f13+f16+f15+f18+f17);
}
} }
} }
} }