SpecialFunctions.jl

SpecialFunctions.jl copied to clipboard

WIP: Partial derivatives of betainc

3

Chain rules for partial derivatives of incomplete beta. Code is by @arzwa and taken from https://github.com/arzwa/IncBetaDer/blob/main/src/beta_inc_grad.jl Related to effort to provide julia implementations of pdf/cdf like functions (https://github.com/JuliaStats/StatsFuns.jl/pull/113) and associated...

bdeonovic

It would be good (to go with #364 ) to have a 1.x release that include the deprecation for `factorial(::Number)`, which was removed in https://github.com/JuliaMath/SpecialFunctions.jl/pull/297

oxinabox


cossio

Use `RealDot.realdot` for more efficient `real(dot(x, y))`

2

This PR uses `RealDot.realdot` to perform computations of `real(dot(x, y))` more efficiently. [RealDot.jl](https://github.com/JuliaMath/RealDot.jl) is a very lightweight package that only defines `realdot` and only depends on LinearAlgebra: https://github.com/JuliaMath/RealDot.jl/blob/main/src/RealDot.jl

devmotion

Handle arguments of 1 in `beta` and `logabsbeta`

4

This PR adds special handling for arguments of 1 in `beta` and `logabsbeta`. This fixes one of the issues in https://github.com/JuliaStats/StatsFuns.jl/pull/126 (I checked locally that the problem is fixed) but...

devmotion

Ensure all functions correctly handle Float16

2

It was noted that `gamma` and `lgamma` did not handle `Float16` in #131, which was subsequently fixed in #132. I'm opening this issue to track correct handling in the remaining...

ararslan

What the title says. The problem seems to be with the Newton's method implementation in lines 432-437 of `erf.jl`. ```julia while true # Newton iterations Δx = sqrtπhalf * (erf(x)...

aquohn

Derivatives of the incomplete gamma functions with respect to parameter `a` are missing

13

Copied from https://github.com/JuliaMath/SpecialFunctions.jl/pull/305: > I added ChainRules definitions for gamma(a, x), loggamma(a, x), and gamma_inc according to https://functions.wolfram.com/GammaBetaErf/GammaRegularized/introductions/Gammas/ShowAll.html . Similar to the Bessel functions, derivatives with respect to the first...

devmotion

Bessel functions only support BigFloat arguments for integer orders

1

Hi! I'm using the besselj and bessely ( )functions on a function to solve their zeros, as is shown in the source code of boostlibrary (`https://github.com/Ceboc/BesselZeros`). As i want to...

Ceboc

Factorial with big float throws domain error

3

As stated [here](https://github.com/JuliaMath/SpecialFunctions.jl/blob/master/src/gamma.jl#L888), if `x` is a float, instead of `factorial(x)` it should use the `gamma(x+1)` function instead. This is not the what is happening, as seen by this output:...

jakewilliami