From 7c3bd5edcf52e644efd139b3f09ac32d571f46b7 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Sun, 19 May 2013 22:56:00 +0000
Subject: [PATCH] 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
---
src/runner_iact.h | 10 ----------
1 file changed, 10 deletions(-)
diff --git a/src/runner_iact.h b/src/runner_iact.h
index 05c73c047b..c9338c9e2a 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 ) );
--
GitLab