Remove Rmath sampling methods

Author: devmotion

Project.toml

@@ -1,6 +1,6 @@
 name = "StatsFuns"
 uuid = "4c63d2b9-4356-54db-8cca-17b64c39e42c"
-version = "1.5.0"
+version = "2.0.0"
 
 [deps]
 HypergeometricFunctions = "34004b35-14d8-5ef3-9330-4cdb6864b03a"

src/rmath.jl

@@ -34,17 +34,9 @@ function _import_rmath(rname::Symbol, jname::Symbol, pargs)
     invlogcdf = Symbol(jname, "invlogcdf")
     invlogccdf = Symbol(jname, "invlogccdf")
 
-    is_tukey = rname == :tukey
-
-    rand = Symbol(jname, "rand")
-    has_rand = true
-    if rname == :nbeta || rname == :nf || rname == :nt || is_tukey
-        has_rand = false
-    end
-
     # arguments & argument types
     _pts = fill(:Cdouble, length(pargs))
-
+    is_tukey = rname == :tukey
     if is_tukey
         dtypes = Expr(:tuple, :Cdouble, :Cdouble, _pts..., :Cint)
         ptypes = Expr(:tuple, :Cdouble, :Cdouble, _pts..., :Cint, :Cint)
@@ -117,11 +109,6 @@ function _import_rmath(rname::Symbol, jname::Symbol, pargs)
             T = float(Base.promote_typeof($(pargs...), lq))
             return convert(T, ccall(($qfun, libRmath), Float64, $qtypes, lq, $(pargs...), 0, 1))
         end
-
-        if $has_rand
-            $rand($(pdecls...)) =
-                ccall(($rfun, libRmath), Float64, $rtypes, $(pargs...))
-        end
     end
 end
 

test/rmath.jl

@@ -85,13 +85,6 @@ function rmathcomp(basename, params, X::AbstractArray, rtol = _default_rtol(para
         rmath_invlogccdf = get_rmathfun(invlogccdf)
     end
 
-    # tackle rand specially
-    has_rand = true
-    if basename == "nbeta" || basename == "nfdist" || basename == "ntdist" || basename == "srdist"
-        has_rand = false
-    end
-    rmath_rand = has_rand ? get_rmathfun(rand) : nothing
-
     if has_pdf
         @testset "pdf with x=$x" for x in X
             check_rmath(
@@ -163,10 +156,7 @@ function rmathcomp(basename, params, X::AbstractArray, rtol = _default_rtol(para
         end
     end
 
-    # make sure that rand works
-    return if has_rand
-        rmath_rand(params...)
-    end
+    return nothing
 end
 
 function rmathcomp_tests(basename, configs)