fixed bug in GMR_Illustris (1+yy**mdelta instead of (1+yy)**mdelta)

nanograv · May 21, 2024 · 53951d2 · 53951d2
1 parent 91d4ca0
commit 53951d2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/holodeck/sams/components.py b/holodeck/sams/components.py
@@ -475,7 +475,7 @@ def __call__(self, mtot, mrat, redz):
         xx = (mtot/self._mref)
         mt = np.power(xx, malpha)
         yy = mtot/self._mref_delta
-        mp1t = np.power(1.0 + yy, mdelta)
+        mp1t = 1.0 + np.power(yy, mdelta)
         qt = np.power(mrat, qgamma)
 
         rate = norm * mt * mp1t * qt