From 15d00d4f8fc6ea0c22259e1094bdf4ed07538e26 Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Mon, 19 Oct 2020 23:07:10 +0200
Subject: [PATCH] initialise the rays in black_holes_init_bpart()

---
 src/black_holes/EAGLE/black_holes.h | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/black_holes/EAGLE/black_holes.h b/src/black_holes/EAGLE/black_holes.h
index 0cb18d6ed1..5311b1bb85 100644
--- a/src/black_holes/EAGLE/black_holes.h
+++ b/src/black_holes/EAGLE/black_holes.h
@@ -30,6 +30,7 @@
 #include "minmax.h"
 #include "physical_constants.h"
 #include "random.h"
+#include "rays.h"
 
 /* Standard includes */
 #include <float.h>
@@ -180,6 +181,9 @@ __attribute__((always_inline)) INLINE static void black_holes_init_bpart(
   bp->f_visc = FLT_MAX;
   bp->accretion_boost_factor = -FLT_MAX;
   bp->mass_at_start_of_step = bp->mass; /* bp->mass may grow in nibbling mode */
+
+  /* Reset the rays carried by this BH */
+  ray_init(bp->rays, eagle_blackhole_number_of_rays);
 }
 
 /**
-- 
GitLab