From dcd594cb67d83ddc7dc120a1f9a3d6787df015d0 Mon Sep 17 00:00:00 2001
From: Filip Husko <filip.husko@durham.ac.uk>
Date: Tue, 19 Sep 2023 12:35:11 +0000
Subject: [PATCH] Fix new time-step for spin-jet scheme
---
src/black_holes/SPIN_JET/black_holes_spin.h | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/black_holes/SPIN_JET/black_holes_spin.h b/src/black_holes/SPIN_JET/black_holes_spin.h
index d7e611b0cb..30da69a95d 100644
--- a/src/black_holes/SPIN_JET/black_holes_spin.h
+++ b/src/black_holes/SPIN_JET/black_holes_spin.h
@@ -271,12 +271,12 @@ __attribute__((always_inline)) INLINE static float r_warp(
* @param constants Physical constants (in internal units).
* @param props Properties of the black hole scheme.
*/
-__attribute__((always_inline)) INLINE static float m_warp(
+__attribute__((always_inline)) INLINE static double m_warp(
struct bpart* bp, const struct phys_const* constants,
const struct black_holes_props* props) {
/* Define placeholder variable for the result */
- float Mw = -1.;
+ double Mw = -1.;
/* Gravitational radius */
const float R_G = R_gravitational(bp->subgrid_mass, constants);
@@ -363,12 +363,12 @@ __attribute__((always_inline)) INLINE static float m_warp(
* @param constants Physical constants (in internal units).
* @param props Properties of the black hole scheme.
*/
-__attribute__((always_inline)) INLINE static float j_warp(
+__attribute__((always_inline)) INLINE static double j_warp(
struct bpart* bp, const struct phys_const* constants,
const struct black_holes_props* props) {
/* Define placeholder variable for the result */
- float Jw = -1.;
+ double Jw = -1.;
/* Start branching depending on which accretion mode the BH is in */
if ((bp->accretion_mode == BH_thick_disc) ||
--
GitLab