diff --git a/test/navier-stokes/howto/taylor-green.cc b/test/navier-stokes/howto/taylor-green.cc index d88bf544520ef179c4e90bfb2aa5ad51be895ed4..218947d3ad879a2c0f5bd1799d9f03ecd7715e83 100644 --- a/test/navier-stokes/howto/taylor-green.cc +++ b/test/navier-stokes/howto/taylor-green.cc @@ -254,9 +254,9 @@ void taylor_green (const GV& gv, const Dune::ParameterTree& configuration, std:: while (time < time_end - dt_min*0.5){ osm.apply(time,dt,xold,x); - // Correct pressure after each step. Implementation not suited for - // parallel computations so it's turned off for periodic boundary - // conditions. + // Correct pressure after each step. Without this pressure + // correction the velocity will be ok but the pressure will be + // shifted by a constant. Dune::PDELab::integrateGridFunction(pdgf,pressure_integral,2); pressure_integral = gv.comm().sum(pressure_integral); std::cout << gv.comm().rank() << " palpo pressure_integral before: " << pressure_integral << std::endl;