Commit 7c3bd5ed authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Removed the unnecessary additional viscosity limiter. A similar condition can...

Removed the unnecessary additional viscosity limiter. A similar condition can be aplied on dt when necessary, i.e. when cooling plays a role.


Former-commit-id: c07ff66fefc4874ae2c2ae2d1371ea1f49a6c2a7
parent 8776aab5
......@@ -406,11 +406,6 @@ __attribute__ ((always_inline)) INLINE static void runner_iact_force ( float r2
/* Apply balsara switch */
Pi_ij *= ( pi->force.balsara + pj->force.balsara );
/* Volker's modified viscosity */
/* dt_max = fmaxf(pi->dt, pj->dt);
if( dt_max > 0 && (wi_dr + wj_dr) < 0. )
Pi_ij = fminf( Pi_ij, 2.f * omega_ij / ( ( mi + mj ) * ( wi_dr + wj_dr ) * dt_max ) ); */
/* Get the common factor out. */
w = ri * ( ( POrho2i * wi_dr + POrho2j * wj_dr ) + 0.25f * Pi_ij * ( wi_dr + wj_dr ) );
......@@ -643,11 +638,6 @@ __attribute__ ((always_inline)) INLINE static void runner_iact_nonsym_force ( fl
/* Apply balsara switch */
Pi_ij *= ( pi->force.balsara + pj->force.balsara );
/* Volker's modified viscosity */
/* dt_max = fmaxf(pi->dt, pj->dt);
if(dt_max > 0 && (wi_dr + wj_dr) < 0.)
Pi_ij = fminf( Pi_ij, 2.f * omega_ij / ( ( mi + mj ) * ( wi_dr + wj_dr ) * dt_max ) ); */
/* Get the common factor out. */
w = ri * ( ( POrho2i * wi_dr + POrho2j * wj_dr ) + 0.25f * Pi_ij * ( wi_dr + wj_dr ) );
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment