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

dfh force

Browse Source
This commit is contained in:
James E McClure
2018-07-18 13:05:58 -04:00
parent 1e3f3cd466
commit 6662698836
2 changed files with 12 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
cpu/dfh.cpp
View File

@@ -519,9 +519,9 @@ extern "C" void ScaLBL_D3Q19_AAeven_DFH(int *neighborList, double *dist, double
//........................................................................
// write the velocity
ux = jx / rho0;
uy = jy / rho0;
uz = jz / rho0;
ux = (jx + force_x) / rho0;
uy = (jy + force_y) / rho0;
uz = (jz + force_z) / rho0;
//Velocity[n] = ux;
//Velocity[Np+n] = uy;
//Velocity[2*Np+n] = uz;
@@ -1109,9 +1109,9 @@ extern "C" void ScaLBL_D3Q19_AAodd_DFH(int *neighborList, double *dist, double *
dist[nread] = fq;
// write the velocity
ux = jx / rho0;
uy = jy / rho0;
uz = jz / rho0;
ux = (jx + force_x) / rho0;
uy = (jy + force_y) / rho0;
uz = (jz + force_z) / rho0;
//Velocity[n] = ux;
//Velocity[Np+n] = uy;
//Velocity[2*Np+n] = uz;

12
gpu/dfh.cu
View File

@@ -536,9 +536,9 @@ __global__ void dvc_ScaLBL_D3Q19_AAeven_DFH(int *neighborList, double *dist, dou
//........................................................................
// write the velocity
ux = jx / rho0;
uy = jy / rho0;
uz = jz / rho0;
ux = (jx + force_x) / rho0;
uy = (jy + force_y) / rho0;
uz = (jz + force_z) / rho0;
//Velocity[n] = ux;
//Velocity[Np+n] = uy;
//Velocity[2*Np+n] = uz;
@@ -1129,9 +1129,9 @@ __global__ void dvc_ScaLBL_D3Q19_AAodd_DFH(int *neighborList, double *dist, doub
dist[nread] = fq;
// write the velocity
ux = jx / rho0;
uy = jy / rho0;
uz = jz / rho0;
ux = (jx + force_x) / rho0;
uy = (jy + force_y) / rho0;
uz = (jz + force_z) / rho0;
//Velocity[n] = ux;
//Velocity[Np+n] = uy;
//Velocity[2*Np+n] = uz;