From 53951d203e48fcc98a9515ec1862ff48db8b61fa Mon Sep 17 00:00:00 2001 From: Laura Blecha Date: Tue, 21 May 2024 14:05:23 -0400 Subject: [PATCH] fixed bug in GMR_Illustris (1+yy**mdelta instead of (1+yy)**mdelta) --- holodeck/sams/components.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/holodeck/sams/components.py b/holodeck/sams/components.py index 9572d0bb..81175556 100644 --- 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