Use expm1

rubenseyer · devmotion · web-flow · commit 8f68d065590a · 2024-01-08T15:47:27.000+01:00
Co-authored-by: David Widmann &lt;devmotion@users.noreply.github.com&gt;
diff --git a/src/univariate/continuous/gumbel.jl b/src/univariate/continuous/gumbel.jl
@@ -102,4 +102,4 @@ logcdf(d::Gumbel, x::Real) = -exp(-zval(d, x))
 
 quantile(d::Gumbel, p::Real) = d.μ - d.θ * log(-log(p))
 
-gradlogpdf(d::Gumbel, x::Real) = (-1 + exp(-zval(d, x))) / d.θ
+gradlogpdf(d::Gumbel, x::Real) = expm1(-zval(d, x)) / d.θ

Original file line number	Diff line number	Diff line change
`@@ -102,4 +102,4 @@ logcdf(d::Gumbel, x::Real) = -exp(-zval(d, x))`
`102`	`102`
`103`	`103`	`quantile(d::Gumbel, p::Real) = d.μ - d.θ * log(-log(p))`
`104`	`104`
`105`		`-gradlogpdf(d::Gumbel, x::Real) = (-1 + exp(-zval(d, x))) / d.θ`
	`105`	`+gradlogpdf(d::Gumbel, x::Real) = expm1(-zval(d, x)) / d.θ`