Skip to content
Snippets Groups Projects
Commit dcd594cb authored by Filip Husko's avatar Filip Husko Committed by Matthieu Schaller
Browse files

Fix new time-step for spin-jet scheme

parent adf9aee1
No related branches found
No related tags found
2 merge requests!1766added read Vz factor from the yml files.,!1762Fix new time-step for spin-jet scheme
...@@ -271,12 +271,12 @@ __attribute__((always_inline)) INLINE static float r_warp( ...@@ -271,12 +271,12 @@ __attribute__((always_inline)) INLINE static float r_warp(
* @param constants Physical constants (in internal units). * @param constants Physical constants (in internal units).
* @param props Properties of the black hole scheme. * @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, struct bpart* bp, const struct phys_const* constants,
const struct black_holes_props* props) { const struct black_holes_props* props) {
/* Define placeholder variable for the result */ /* Define placeholder variable for the result */
float Mw = -1.; double Mw = -1.;
/* Gravitational radius */ /* Gravitational radius */
const float R_G = R_gravitational(bp->subgrid_mass, constants); const float R_G = R_gravitational(bp->subgrid_mass, constants);
...@@ -363,12 +363,12 @@ __attribute__((always_inline)) INLINE static float m_warp( ...@@ -363,12 +363,12 @@ __attribute__((always_inline)) INLINE static float m_warp(
* @param constants Physical constants (in internal units). * @param constants Physical constants (in internal units).
* @param props Properties of the black hole scheme. * @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, struct bpart* bp, const struct phys_const* constants,
const struct black_holes_props* props) { const struct black_holes_props* props) {
/* Define placeholder variable for the result */ /* Define placeholder variable for the result */
float Jw = -1.; double Jw = -1.;
/* Start branching depending on which accretion mode the BH is in */ /* Start branching depending on which accretion mode the BH is in */
if ((bp->accretion_mode == BH_thick_disc) || if ((bp->accretion_mode == BH_thick_disc) ||
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment