Update beta_inc.jl

Author: devmotion

test/beta_inc.jl

@@ -218,10 +218,13 @@
         @test beta_inc(a, b, x, 1 - x) isa Tuple{T,T}
         @test beta_inc(a, b, x) == beta_inc(a, b, x, 1 - x)
     end
-    a, b = randexp(2)
-    x = rand()
-    @test beta_inc(a, Float32(b), x) == beta_inc(a, b, x)
-    @test beta_inc(a, Float32(b), Float32(x), 1-x) == beta_inc(a, b, x, 1-x)
+    a = randexp()
+    b = randexp(Float32)
+    x = rand(Float32)
+    @test beta_inc(a, b, x) isa Tuple{Float64,Float64}
+    @test beta_inc(a, b, x) == beta_inc(a, Float64(b), Float64(x))
+    @test beta_inc(a, b, x, 1 - x) isa Tuple{Float64,Float64}
+    @test beta_inc(a, b, x, 1 - x) == beta_inc(a, Float64(b), Float64(x), 1 - x)
 
     @test SpecialFunctions.loggammadiv(13.89, 21.0001) ≈ log(gamma(big(21.0001))/gamma(big(21.0001)+big(13.89)))
     @test SpecialFunctions.stirling_corr(11.99, 100.1) ≈ SpecialFunctions.stirling_error(11.99) + SpecialFunctions.stirling_error(100.1) - SpecialFunctions.stirling_error(11.99 + 100.1)