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

mix gradient compiles

Browse Source
This commit is contained in:
JamesEMcclure 2020-10-06 14:05:26 -04:00
parent d64e9f9bb0
commit db5579530d
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
cpu/MixedGradient.cpp
View File

@ -29,7 +29,7 @@ extern "C" void ScaLBL_D3Q19_MixedGradient(int *Map, double *Phi, double *Gradie
vp = Phi[np]; vp = Phi[np];
vp2 = Phi[np2]; vp2 = Phi[np2];
vm = Phi[nm]; vm = Phi[nm];
grad += 0.25*(5.0*vp1-vp2-3.0*v-vm); grad += 0.25*(5.0*vp-vp2-3.0*v-vm);
} }
for (int q=6; q<18; q++){ for (int q=6; q<18; q++){
int iqx = D3Q19[q][0]; int iqx = D3Q19[q][0];
@ -41,7 +41,7 @@ extern "C" void ScaLBL_D3Q19_MixedGradient(int *Map, double *Phi, double *Gradie
vp = Phi[np]; vp = Phi[np];
vp2 = Phi[np2]; vp2 = Phi[np2];
vm = Phi[nm]; vm = Phi[nm];
grad += 0.125*(5.0*vp1-vp2-3.0*v-vm); grad += 0.125*(5.0*vp-vp2-3.0*v-vm);
} }
Gradient[n] = grad; Gradient[n] = grad;
} }