@@ -23,7 +23,7 @@ import Base: conj, sin, cos, tan, sinh, cosh, tanh,
2323 exp2, exp10, expm1, cbrt, sqrt,
2424 significand, frexp, modf
2525import LinearAlgebra: norm, ⋅ , dot, promote_leaf_eltypes # norm1, norm2, normInf
26- import SpecialFunctions: gamma, lgamma , lfact
26+ import SpecialFunctions: gamma, logabsgamma , lfact
2727using Statistics
2828import Statistics: middle, _mean_promote
2929
@@ -213,17 +213,21 @@ const unaryOps = (:~, :conj, :abs,
213213 :asind , :atand , :rad2deg , :deg2rad ,
214214 :log , :log2 , :log10 , :log1p , :exponent , :exp ,
215215 :exp2 , :exp10 , :expm1 , :cbrt , :sqrt ,
216- :significand , :lgamma ,
216+ :significand ,
217217 :gamma , :lfact , :frexp , :modf ,
218218 :(SpecialFunctions. erf), :(SpecialFunctions. erfc),
219219 :(SpecialFunctions. erfcx), :(SpecialFunctions. erfi), :(SpecialFunctions. dawson),
220- :(SpecialFunctions. airy), :(SpecialFunctions . airyai),
221- :(SpecialFunctions. airyprime), :(SpecialFunctions . airyaiprime), :(SpecialFunctions. airybi), :(SpecialFunctions. airybiprime),
220+ :(SpecialFunctions. airyai),
221+ :(SpecialFunctions. airyaiprime), :(SpecialFunctions. airybi), :(SpecialFunctions. airybiprime),
222222 :(SpecialFunctions. besselj0), :(SpecialFunctions. besselj1), :(SpecialFunctions. bessely0), :(SpecialFunctions. bessely1),
223223 :(SpecialFunctions. eta), :(SpecialFunctions. zeta), :(SpecialFunctions. digamma))
224224for op in unaryOps
225225 @eval ($ op)(c:: AbstractGray ) = Gray ($ op (gray (c)))
226226end
227+ function logabsgamma (c:: AbstractGray )
228+ lagc, s = logabsgamma (gray (c))
229+ return Gray (lagc), s
230+ end
227231
228232"""
229233 y = complement(x)
0 commit comments