Loic Hausammann
--- a/src/cooling/grackle/cooling.h
+++ b/src/cooling/grackle/cooling.h
    const struct phys_const* phys_const,
    struct cooling_function_data* cooling) {

+  fedisableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
--- a/src/cooling/grackle/cooling.h
+++ b/src/cooling/grackle/cooling.h
 */

 /* Some standard headers. */
+#include <fenv.h>
--- a/src/cooling/grackle/cooling.h
+++ b/src/cooling/grackle/cooling.h
+  fedisableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
+
  if (GRACKLE_NPART != 1)
    error("Grackle with multiple particles not implemented");

  /* read parameters */
-  cooling_read_parameters(parameter_file, cooling);
+  cooling_read_parameters(parameter_file, cooling, phys_const);

  /* Set up the units system. */
  cooling_init_units(us, cooling);

  /* Set up grackle */
  cooling_init_grackle(cooling);
+
+  feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
--- a/src/feedback/GEAR/feedback_iact.h 0 → 100644
+++ b/src/feedback/GEAR/feedback_iact.h 0 → 100644
+    return;
+  }
+
+  const float mj = hydro_get_mass(pj);
+  const float r = sqrtf(r2);
+
+  /* Get the kernel for hi. */
+  float hi_inv = 1.0f / hi;
+  float hi_inv_dim = pow_dimension(hi_inv);       /* 1/h^d */
+  float xi = r * hi_inv;
+  float wi, wi_dx;
+  kernel_deval(xi, &wi, &wi_dx);
+  wi *= hi_inv_dim;
+
+  /* Compute inverse enrichment weight */
+  const double si_inv_weight = si->feedback_data.enrichment_weight == 0?
--- a/src/feedback/GEAR/feedback_struct.h 0 → 100644
+++ b/src/feedback/GEAR/feedback_struct.h 0 → 100644
+  float delta_mass;
+
+  /*! specific energy received from supernovae */
+  float delta_u;
+
+  /*! Momemtum received from a supernovae */
+  float delta_p[3];
+
+};
+
+/**
+ * @brief Feedback fields carried by each star particles
+ */
+struct feedback_spart_data {
+
+  /*! Inverse of normalisation factor used for the enrichment. */