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

don't compute density inside solid (just reverting a debugging change)

Browse Source
This commit is contained in:
James E McClure 2015-07-04 11:29:15 -04:00
parent d6cf96f24c
commit 0eb4221b57
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
cpu/D3Q7.cpp
View File

@ -150,8 +150,8 @@ extern "C" void ComputeDensityD3Q7(char *ID, double *disteven, double *distodd,
int N = Nx*Ny*Nz;
for (n=0; n<N; n++){
//id = ID[n];
//if (id > 0 ){
id = ID[n];
if (id > 0 ){
// Read the distributions
f0 = disteven[n];
f2 = disteven[N+n];
@ -162,6 +162,6 @@ extern "C" void ComputeDensityD3Q7(char *ID, double *disteven, double *distodd,
f5 = distodd[2*N+n];
// Compute the density
Den[n] = f0+f1+f2+f3+f4+f5+f6;
//}
}
}
}