From b41060a4b588b76e77c918324d42535ec53d0282 Mon Sep 17 00:00:00 2001
From: lhausamm <loic_hausammann@hotmail.com>
Date: Fri, 12 Jan 2018 13:57:04 +0100
Subject: [PATCH] Remove diffusion/particle mixing functions

---
 src/cooling.h                      | 8 ++++----
 src/cooling/grackle/cooling_iact.h | 6 +++---
 src/hydro/Gadget2/hydro_iact.h     | 2 --
 3 files changed, 7 insertions(+), 9 deletions(-)

diff --git a/src/cooling.h b/src/cooling.h
index 06085e6794..72caff33ea 100644
--- a/src/cooling.h
+++ b/src/cooling.h
@@ -29,16 +29,16 @@
 
 /* Import the right cooling definition */
 #if defined(COOLING_NONE)
-#include "./cooling/none/cooling_iact.h"
+//#include "./cooling/none/cooling_iact.h"
 #include "./cooling/none/cooling.h"
 #elif defined(COOLING_CONST_DU)
-#include "./cooling/const_du/cooling_iact.h"
+//#include "./cooling/const_du/cooling_iact.h"
 #include "./cooling/const_du/cooling.h"
 #elif defined(COOLING_CONST_LAMBDA)
-#include "./cooling/const_lambda/cooling_iact.h"
+//#include "./cooling/const_lambda/cooling_iact.h"
 #include "./cooling/const_lambda/cooling.h"
 #elif defined(COOLING_GRACKLE)
-#include "./cooling/grackle/cooling_iact.h"
+//#include "./cooling/grackle/cooling_iact.h"
 #include "./cooling/grackle/cooling.h"
 #else
 #error "Invalid choice of cooling function."
diff --git a/src/cooling/grackle/cooling_iact.h b/src/cooling/grackle/cooling_iact.h
index 502858800d..5b465bf5ca 100644
--- a/src/cooling/grackle/cooling_iact.h
+++ b/src/cooling/grackle/cooling_iact.h
@@ -25,8 +25,8 @@
 __attribute__((always_inline)) INLINE static void cooling_density_iact(
     float wi, float wj, struct part *pi, struct part *pj) {
 
-  float diff = pi->cooling_data.He_density - pj->cooling_data.He_density;
+  /* float diff = pi->cooling_data.He_density - pj->cooling_data.He_density; */
 
-  pi->cooling_data.He_density += wj * diff;
-  pj->cooling_data.He_density += wi * diff;
+  /* pi->cooling_data.He_density += wj * diff; */
+  /* pj->cooling_data.He_density += wi * diff; */
 };
diff --git a/src/hydro/Gadget2/hydro_iact.h b/src/hydro/Gadget2/hydro_iact.h
index 29a86da3f9..6f839a4710 100644
--- a/src/hydro/Gadget2/hydro_iact.h
+++ b/src/hydro/Gadget2/hydro_iact.h
@@ -32,7 +32,6 @@
  * Gadget-2 tree-code neighbours search.
  */
 
-#include "cooling.h"
 #include "cache.h"
 #include "minmax.h"
 
@@ -105,7 +104,6 @@ __attribute__((always_inline)) INLINE static void runner_iact_density(
   pj->density.rot_v[1] += facj * curlvr[1];
   pj->density.rot_v[2] += facj * curlvr[2];
 
-  cooling_density_iact(wi, wj, pi, pj);
 
 #ifdef DEBUG_INTERACTIONS_SPH
   /* Update ngb counters */
-- 
GitLab