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

Refactoring gpu/Color.cu: MassColorCollideD3Q7

Browse Source
This commit is contained in:
James E McClure
2015-06-24 18:21:17 -04:00
parent 2e4822385a
commit fe07dfc2ad
1 changed files with 7 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
gpu/Color.cu
View File

@@ -1085,24 +1085,23 @@ __global__ void dvc_ColorCollideOpt( char *ID, double *disteven, double *distod
__global__ void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_odd, double *B_even, double *B_odd, __global__ void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_odd, double *B_even, double *B_odd,
double *Den, double *Phi, double *ColorGrad, double *Velocity, double beta, int N, bool pBC) double *Den, double *Phi, double *ColorGrad, double *Velocity, double beta, int N, bool pBC)
{ {
int idx,n,q,Cqx,Cqy,Cqz; char id;
// int sendLoc; int n;
double f0,f1,f2,f3,f4,f5,f6; double f0,f1,f2,f3,f4,f5,f6;
double na,nb; // density values double na,nb; // density values
double ux,uy,uz; // flow velocity double ux,uy,uz; // flow velocity
double nx,ny,nz,C; // color gradient components double nx,ny,nz,C; // color gradient components
double a1,a2,b1,b2; double a1,a2,b1,b2;
double sp,delta; double delta;
//double feq[6]; // equilibrium distributions
// Set of Discrete velocities for the D3Q19 Model
//int D3Q7[3][3]={{1,0,0},{0,1,0},{0,0,1}};
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){
//.....Load the Color gradient......... //.....Load the Color gradient.........
nx = ColorGrad[n]; nx = ColorGrad[n];