diff --git a/src/runner_iact.h b/src/runner_iact.h index 05c73c047b7b627ad582ca22a34e296c9e7a072c..c9338c9e2adf3d2cc576629a02819d0a393d9d69 100644 --- a/src/runner_iact.h +++ b/src/runner_iact.h @@ -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 ) );