diff --git a/src/black_holes/EAGLE/black_holes.h b/src/black_holes/EAGLE/black_holes.h index 9092e80a4761c62c8e70bc04d5aba8ba7ff40299..2dd7345e8cb50f2be713a94b9b0f8b60db447ad1 100644 --- a/src/black_holes/EAGLE/black_holes.h +++ b/src/black_holes/EAGLE/black_holes.h @@ -118,6 +118,8 @@ __attribute__((always_inline)) INLINE static void black_holes_predict_extra( } #endif + message("Reposition!"); + bp->x[0] = bp->reposition.x[0]; bp->x[1] = bp->reposition.x[1]; bp->x[2] = bp->reposition.x[2]; @@ -458,11 +460,9 @@ __attribute__((always_inline)) INLINE static void black_holes_prepare_feedback( } } -__attribute__((always_inline)) INLINE static void -black_holes_end_reposition(struct bpart* restrict bp, - const struct black_holes_props* props, - const struct phys_const* constants, - const struct cosmology* cosmo) { +__attribute__((always_inline)) INLINE static void black_holes_end_reposition( + struct bpart* restrict bp, const struct black_holes_props* props, + const struct phys_const* constants, const struct cosmology* cosmo) { const float potential = gravity_get_comoving_potential(bp->gpart); @@ -471,8 +471,13 @@ black_holes_end_reposition(struct bpart* restrict bp, if (potential < bp->reposition.min_potential || bp->subgrid_mass > props->max_reposition_mass) { + message("No repo!"); + /* No need to reposition */ bp->reposition.min_potential = FLT_MAX; + bp->reposition.x[0] = -FLT_MAX; + bp->reposition.x[1] = -FLT_MAX; + bp->reposition.x[2] = -FLT_MAX; } }