Permalink
Browse files

potential bugfix for bounce back in CPU lb

  • Loading branch information...
1 parent 8d50f0d commit 6ad3e767e8763eab4be6da8a9673ff0bb458ec2f @lahnerml lahnerml committed Jan 30, 2017
Showing with 3 additions and 3 deletions.
  1. +3 −3 src/core/lb-boundaries.cpp
@@ -558,12 +558,12 @@ void lb_bounce_back() {
population_shift = 0;
for (l = 0; l < 3; l++) {
population_shift -=
- lbpar.agrid * lbpar.agrid * lbpar.agrid * lbpar.agrid *
- lbpar.agrid * lbpar.rho[0] * 2 * lbmodel.c[i][l] *
- lbmodel.w[i] *
+ lbpar.agrid * lbpar.agrid * lbpar.agrid * lbpar.rho[0] * 2 *
+ lbmodel.c[i][l] * lbmodel.w[i] *
lb_boundaries[lbfields[k].boundary - 1].velocity[l] /
lbmodel.c_sound_sq;
}
+
if (x - lbmodel.c[i][0] > 0 &&
x - lbmodel.c[i][0] < lblattice.grid[0] + 1 &&
y - lbmodel.c[i][1] > 0 &&

0 comments on commit 6ad3e76

Please sign in to comment.