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

gpu interfaces not moving :/

Browse Source
This commit is contained in:
James E McClure 2015-07-04 14:23:47 -04:00
parent 465853410b
commit 9f2840ccdc
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
gpu/Color.cu
View File

@ -1147,7 +1147,7 @@ __global__ void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_od
nab = 1.0/(na+nb); nab = 1.0/(na+nb);
//........................................................................ //........................................................................
delta = beta*na*nb*nab*0.1111111111111111*nx; delta = beta*na*nb*nab*0.1111111111111111*nx;
if (!(na*nb*nab>0)) delta=0; if (na*nb*nab<0.0)) delta=0.0;
a1 = na*(0.1111111111111111*(1+4.5*ux))+delta; a1 = na*(0.1111111111111111*(1+4.5*ux))+delta;
b1 = nb*(0.1111111111111111*(1+4.5*ux))-delta; b1 = nb*(0.1111111111111111*(1+4.5*ux))-delta;
@ -1163,7 +1163,8 @@ __global__ void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_od
// q = 2 // q = 2
// Cq = {0,1,0} // Cq = {0,1,0}
delta = beta*na*nb*nab*0.1111111111111111*ny; delta = beta*na*nb*nab*0.1111111111111111*ny;
if (!(na*nb*nab>0)) delta=0; if (na*nb*nab<0.0)) delta=0.0;
a1 = na*(0.1111111111111111*(1+4.5*uy))+delta; a1 = na*(0.1111111111111111*(1+4.5*uy))+delta;
b1 = nb*(0.1111111111111111*(1+4.5*uy))-delta; b1 = nb*(0.1111111111111111*(1+4.5*uy))-delta;
a2 = na*(0.1111111111111111*(1-4.5*uy))-delta; a2 = na*(0.1111111111111111*(1-4.5*uy))-delta;
@ -1177,7 +1178,8 @@ __global__ void dvc_MassColorCollideD3Q7(char *ID, double *A_even, double *A_od
// q = 4 // q = 4
// Cq = {0,0,1} // Cq = {0,0,1}
delta = beta*na*nb*nab*0.1111111111111111*nz; delta = beta*na*nb*nab*0.1111111111111111*nz;
if (!(na*nb*nab>0)) delta=0; if (na*nb*nab<0.0)) delta=0.0;
a1 = na*(0.1111111111111111*(1+4.5*uz))+delta; a1 = na*(0.1111111111111111*(1+4.5*uz))+delta;
b1 = nb*(0.1111111111111111*(1+4.5*uz))-delta; b1 = nb*(0.1111111111111111*(1+4.5*uz))-delta;
a2 = na*(0.1111111111111111*(1-4.5*uz))-delta; a2 = na*(0.1111111111111111*(1-4.5*uz))-delta;