From 805cdcecd7e958388429fdfdac73586d795af7ab Mon Sep 17 00:00:00 2001 From: Cristian Lalescu <Cristian.Lalescu@ds.mpg.de> Date: Wed, 1 Jul 2015 10:21:22 +0200 Subject: [PATCH] make Kolmogorov force *actually* A*sin(y) --- src/fluid_solver.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/fluid_solver.cpp b/src/fluid_solver.cpp index 0ebae211..d476691b 100644 --- a/src/fluid_solver.cpp +++ b/src/fluid_solver.cpp @@ -218,12 +218,12 @@ void fluid_solver<R>::add_forcing(\ if (this->cd->myrank == this->cd->rank[this->fmode]) \ { \ cindex = ((this->fmode - this->cd->starts[0]) * this->cd->sizes[1])*this->cd->sizes[2]; \ - field[cindex*3+2][1] -= this->famplitude*factor/2; \ + field[cindex*3+2][0] -= this->famplitude*factor/2; \ } \ if (this->cd->myrank == this->cd->rank[this->cd->sizes[0] - this->fmode]) \ { \ cindex = ((this->cd->sizes[0] - this->fmode - this->cd->starts[0]) * this->cd->sizes[1])*this->cd->sizes[2]; \ - field[cindex*3+2][1] += this->famplitude*factor/2; \ + field[cindex*3+2][0] -= this->famplitude*factor/2; \ } \ } \ \ -- GitLab