Skip to content

Commit 6b84662

Browse files
committed
Use Speve's version
1 parent 4d64e3c commit 6b84662

File tree

1 file changed

+5
-10
lines changed

1 file changed

+5
-10
lines changed

src/gamma_inc.jl

Lines changed: 5 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -1015,18 +1015,13 @@ function _gamma_inc_inv(a::Float64, p::Float64, q::Float64)
10151015
end
10161016

10171017
function _gamma_inc_inv(a::T, p::T, q::T) where {T <: Union{Float16, Float32}}
1018+
if p + q != one(T)
1019+
throw(ArgumentError("p + q must equal one but was $(p + q)"))
1020+
end
10181021
p64, q64 = if p < q
1019-
_q = 1 - Float64(p)
1020-
if q != T(_q)
1021-
throw(ArgumentError("p + q must equal one but was $(p + q)"))
1022-
end
1023-
Float64(p), _q
1022+
(Float64(p), 1 - Float64(p))
10241023
else
1025-
_p = 1 - Float64(q)
1026-
if p != T(_p)
1027-
throw(ArgumentError("p + q must equal one but was $(p + q)"))
1028-
end
1029-
_p, Float64(q)
1024+
(1 - Float64(q), Float64(q))
10301025
end
10311026
return T(_gamma_inc_inv(Float64(a), p64, q64))
10321027
end

0 commit comments

Comments
 (0)