Merge pull request #148 from JuliaSymbolics/s/specfun

Special Functions

Author: shashi

.travis.yml

@@ -8,8 +8,7 @@ branches:
 os:
   - linux
 julia:
-  - 1.3.1
-  - 1.4.1
+  - 1
   - nightly
 matrix:
   allow_failures:

src/methods.jl

@@ -1,6 +1,22 @@
-const monadic = [deg2rad, rad2deg, transpose, -, conj, asind, log1p, acsch, acos, asec, acosh, acsc, cscd, log, tand, log10, csch, asinh, abs2, cosh, sin, cos, atan, cospi, cbrt, acosd, acoth, acotd, asecd, exp, acot, sqrt, sind, sinpi, asech, log2, tan, exp10, sech, coth, asin, cotd, cosd, sinh, abs, csc, tanh, secd, atand, sec, acscd, cot, exp2, expm1, atanh, real]
-
-const diadic = [+, -, max, min, *, /, \, hypot, atan, mod, rem, ^, copysign]
+import SpecialFunctions: gamma, loggamma, erf, erfc, erfcinv, erfi, erfcx,
+                         dawson, digamma, trigamma, invdigamma, polygamma,
+                         airyai, airyaiprime, airybi, airybiprime, besselj0,
+                         besselj1, bessely0, bessely1, besselj, bessely, besseli,
+                         besselk, hankelh1, hankelh2, polygamma, beta, logbeta
+
+const monadic = [deg2rad, rad2deg, transpose, -, conj, asind, log1p, acsch,
+                 acos, asec, acosh, acsc, cscd, log, tand, log10, csch, asinh,
+                 abs2, cosh, sin, cos, atan, cospi, cbrt, acosd, acoth, acotd,
+                 asecd, exp, acot, sqrt, sind, sinpi, asech, log2, tan, exp10,
+                 sech, coth, asin, cotd, cosd, sinh, abs, csc, tanh, secd,
+                 atand, sec, acscd, cot, exp2, expm1, atanh, real, gamma,
+                 loggamma, erf, erfc, erfcinv, erfi, erfcx, dawson, digamma,
+                 trigamma, invdigamma, polygamma, airyai, airyaiprime, airybi,
+                 airybiprime, besselj0, besselj1, bessely0, bessely1]
+
+const diadic = [+, -, max, min, *, /, \, hypot, atan, mod, rem, ^, copysign,
+                besselj, bessely, besseli, besselk, hankelh1, hankelh2,
+                polygamma, beta, logbeta]
 
 const previously_declared_for = Set([])
 

test/fuzzlib.jl

@@ -92,7 +92,8 @@ function gen_rand_expr(inputs;
     try
         return f(args...)
     catch err
-        if err isa DomainError || err isa DivideError || err isa MethodError
+        if err isa DomainError || err isa DivideError || err isa MethodError ||
+            err isa SymbolicUtils.SpecialFunctions.AmosException
             return gen_rand_expr(inputs,
                                  spec=spec,
                                  leaf_prob=leaf_prob,