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

Compute velocity instead of momentum

Browse Source
This commit is contained in:
James E McClure
2017-11-26 07:47:00 -05:00
parent 53fbe986d6
commit ff31f21337
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
cpu/Color.cpp
View File

@@ -1299,9 +1299,9 @@ extern "C" void ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven, doubl
distodd[7*N+n] = f15;
distodd[8*N+n] = f17;
//...Store the Velocity..........................
Velocity[n] = jx;
Velocity[N+n] = jy;
Velocity[2*N+n] = jz;
Velocity[n] = jx/rho_0;
Velocity[N+n] = jy/rho_0;
Velocity[2*N+n] = jz/rho_0;
//***************************************************************
} // check if n is in the solid
} // loop over n

6
gpu/Color.cu
View File

@@ -1300,9 +1300,9 @@ __global__ void dvc_ScaLBL_D3Q19_ColorCollide_gen( char *ID, double *disteven,
distodd[7*N+n] = f15;
distodd[8*N+n] = f17;
//...Store the Velocity..........................
Velocity[n] = jx;
Velocity[N+n] = jy;
Velocity[2*N+n] = jz;
Velocity[n] = jx/rho0;
Velocity[N+n] = jy/rho0;
Velocity[2*N+n] = jz/rho0;
//***************************************************************
}// check if n is in the solid