# \#specialfunctions

**URL:** https://discourse.julialang.org/tag/specialfunctions/671.md

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

---

## [Using SpecialFunctions spherical bessel functions with Reactant](https://discourse.julialang.org/t/using-specialfunctions-spherical-bessel-functions-with-reactant/131993)

<div class="topic-metadata">

**Author:** [@dameka](https://discourse.julialang.org/u/dameka)\
**Replies:** 3\
**Last updated:** [August 31, 2025, 6:06pm UTC](https://discourse.julialang.org/t/using-specialfunctions-spherical-bessel-functions-with-reactant/131993 "2025-08-31T18:06:56Z")

</div>

I’m trying to apply Reactant to my code which uses spherical Bessel functions from the SpecialFunctions package. I’m running into this error; has anyone been able to use spherical Bessel functions with Reactant? @avikpal …

---

## [How to notify JuliaMath of a PR to one of their (almost) unwatched packages](https://discourse.julialang.org/t/how-to-notify-juliamath-of-a-pr-to-one-of-their-almost-unwatched-packages/131655)

<div class="topic-metadata">

**Author:** [@PeterSimon](https://discourse.julialang.org/u/PeterSimon)\
**Replies:** 4\
**Last updated:** [August 18, 2025, 7:41am UTC](https://discourse.julialang.org/t/how-to-notify-juliamath-of-a-pr-to-one-of-their-almost-unwatched-packages/131655 "2025-08-18T07:41:23Z")

</div>

I recently submitted a PR to FunctionZeros.jl to implement new functionality I need for a package I’m developing. FunctionZeros.jl is in the JuliaMath organization, and I’m not sure of the proper protocol, if any, to not…

---

## [Why functions in SpecialFunctions package work on CUDA arrays?](https://discourse.julialang.org/t/why-functions-in-specialfunctions-package-work-on-cuda-arrays/128194)

<div class="topic-metadata">

**Author:** [@singularity](https://discourse.julialang.org/u/singularity)\
**Replies:** 7\
**Last updated:** [April 20, 2025, 7:04am UTC](https://discourse.julialang.org/t/why-functions-in-specialfunctions-package-work-on-cuda-arrays/128194 "2025-04-20T07:04:45Z")

</div>

HI. I am confused how the behavior of functions inside the SpecialFunctions package on CUDA arrays. I am confused why the code below works and doesn’t throw an error at the last line. using CUDA using SpecialFunctions…

---

## [SpecialFunctions + Turing, AD support for quantile function of Gamma distribution gamma\_inv\_cdf](https://discourse.julialang.org/t/specialfunctions-turing-ad-support-for-quantile-function-of-gamma-distribution-gamma-inv-cdf/113663)

<div class="topic-metadata">

**Author:** [@mleprovost](https://discourse.julialang.org/u/mleprovost)\
**Replies:** 0\
**Last updated:** [April 30, 2024, 5:14pm UTC](https://discourse.julialang.org/t/specialfunctions-turing-ad-support-for-quantile-function-of-gamma-distribution-gamma-inv-cdf/113663 "2024-04-30T17:14:39Z")

</div>

Hello, I would like to use AD in Turing.jl to differentiate a custom density that involves the inverse cdf (i.e., quantile function) of a Gamma distribution. Unfortunately, the quantile function of the Gamma distributi…

---

## [Fast implementation of incomplete extended gamma function, QuadGK?](https://discourse.julialang.org/t/fast-implementation-of-incomplete-extended-gamma-function-quadgk/112239)

<div class="topic-metadata">

**Author:** [@mleprovost](https://discourse.julialang.org/u/mleprovost)\
**Replies:** 4\
**Last updated:** [March 29, 2024, 3:10am UTC](https://discourse.julialang.org/t/fast-implementation-of-incomplete-extended-gamma-function-quadgk/112239 "2024-03-29T03:10:42Z")

</div>

Hello, I am interested in developing a fast routine to compute the following extended incomplete gamma function for various pairs (x, q). \\gamma(\\alpha, \\kappa, p, q)(x) = \\int\_{0}^x t^{\\alpha-1} \\exp( -\\frac{t^\\kappa}…

---

## [Problem with elliptic integrals](https://discourse.julialang.org/t/problem-with-elliptic-integrals/99922)

<div class="topic-metadata">

**Author:** [@Eben60](https://discourse.julialang.org/u/Eben60)\
**Replies:** 1\
**Last updated:** [June 6, 2023, 6:09am UTC](https://discourse.julialang.org/t/problem-with-elliptic-integrals/99922 "2023-06-06T06:09:30Z")

</div>

using Plots using SpecialFunctions # using Elliptic K = ellipk # Computes Complete Elliptic Integral of 1st kind K(m) E = ellipe # Computes Complete Elliptic Integral of 2nd kind E(m) psi0(k) = (K(k) - E(k)) / k @sho…

---

## [Numerically stable derivative of logerf without catastrophic cancellation](https://discourse.julialang.org/t/numerically-stable-derivative-of-logerf-without-catastrophic-cancellation/92992)

<div class="topic-metadata">

**Author:** [@user22](https://discourse.julialang.org/u/user22)\
**Replies:** 4\
**Last updated:** [January 15, 2023, 4:56pm UTC](https://discourse.julialang.org/t/numerically-stable-derivative-of-logerf-without-catastrophic-cancellation/92992 "2023-01-15T16:56:08Z")

</div>

I need to evaluate the derivative of a logerf function fun with respect to x within the widest possible range of x. At the moment, the derivative is calculated through the function dfun. using LogExpFunctions using Fini…

---

## [InvalidIRError with CuArray broadcast when complex values are involved in Bessel functions](https://discourse.julialang.org/t/invalidirerror-with-cuarray-broadcast-when-complex-values-are-involved-in-bessel-functions/91312)

<div class="topic-metadata">

**Author:** [@jinml](https://discourse.julialang.org/u/jinml)\
**Replies:** 2\
**Last updated:** [December 8, 2022, 12:47am UTC](https://discourse.julialang.org/t/invalidirerror-with-cuarray-broadcast-when-complex-values-are-involved-in-bessel-functions/91312 "2022-12-08T00:47:51Z")

</div>

Hi, I’m receiving runtime errors from julia when playing with the following simple code \> \> julia\> using CUDA, SpecialFunctions \> \> julia\> function cutest(x::Float64,f::Float64) \> c=besselj0(f); \> return…

---

## [\[ANN\] Bessels.jl - Bessel, Hankel, and Airy functions](https://discourse.julialang.org/t/ann-bessels-jl-bessel-hankel-and-airy-functions/87941)

<div class="topic-metadata">

**Author:** [@heltonmc](https://discourse.julialang.org/u/heltonmc)\
**Replies:** 12\
**Last updated:** [October 12, 2022, 11:26pm UTC](https://discourse.julialang.org/t/ann-bessels-jl-bessel-hankel-and-airy-functions/87941 "2022-10-12T23:26:36Z")

</div>

We are excited to announce Bessels.jl which is a special function library written in Julia to compute Bessel and Modified Bessel functions, Spherical Bessel functions, Airy functions, and Hankel functions. Currently the …

---

## [Elliptic theta functions](https://discourse.julialang.org/t/elliptic-theta-functions/65742)

<div class="topic-metadata">

**Author:** [@Gabrielmlando](https://discourse.julialang.org/u/Gabrielmlando)\
**Replies:** 15\
**Last updated:** [July 13, 2022, 1:44pm UTC](https://discourse.julialang.org/t/elliptic-theta-functions/65742 "2022-07-13T13:44:26Z")

</div>

Are elliptic theta functions implemented in Julia? They exist in Mathematica as EllipticTheta, but I need them for heavy numerical calculations.

---

## [How to work around this error in expint?](https://discourse.julialang.org/t/how-to-work-around-this-error-in-expint/76133)

<div class="topic-metadata">

**Author:** [@amrods](https://discourse.julialang.org/u/amrods)\
**Replies:** 1\
**Last updated:** [February 11, 2022, 8:28pm UTC](https://discourse.julialang.org/t/how-to-work-around-this-error-in-expint/76133 "2022-02-11T20:28:16Z")

</div>

I ran into this error while evaluating the Gamma function from SpecialFunctions: julia\> using SpecialFunctions julia\> gamma(exp(19.39)+1, 0.0125\*0.6) ERROR: ArgumentError: Unsupported order |ν| \> 50 off the positive re…

---

## [The endless packages to add as as a dependency or to resolve](https://discourse.julialang.org/t/the-endless-packages-to-add-as-as-a-dependency-or-to-resolve/70450)

<div class="topic-metadata">

**Author:** [@Redmusk888](https://discourse.julialang.org/u/Redmusk888)\
**Replies:** 2\
**Last updated:** [October 27, 2021, 12:08pm UTC](https://discourse.julialang.org/t/the-endless-packages-to-add-as-as-a-dependency-or-to-resolve/70450 "2021-10-27T12:08:00Z")

</div>

1.I’m using julia 1.6.3 as my jupyternote kernel. 2.I was learning Linear Algebra on juliaacademy.com .https://juliaacademy.com/courses/937702/lectures/17339511 3.i got the problem below 4.i have been add and reso…

---

## [Hypergeometric functions](https://discourse.julialang.org/t/hypergeometric-functions/63621)

<div class="topic-metadata">

**Author:** [@marcobonici](https://discourse.julialang.org/u/marcobonici)\
**Replies:** 6\
**Last updated:** [June 29, 2021, 9:48am UTC](https://discourse.julialang.org/t/hypergeometric-functions/63621 "2021-06-29T09:48:29Z")

</div>

Hello, I need to evaluate the Gauss Hypergeometric functions \_2F\_1(a,b;c;z). Here are some specifications for the arguments of this function a has compex values. The real part is \\sim2 (and is strictly positive), whi…

---

## [Differences in MATLAB's versus Julia's beta(x, y) for large x or y?](https://discourse.julialang.org/t/differences-in-matlabs-versus-julias-beta-x-y-for-large-x-or-y/60138)

<div class="topic-metadata">

**Author:** [@bashonubuntu](https://discourse.julialang.org/u/bashonubuntu)\
**Replies:** 16\
**Last updated:** [May 1, 2021, 7:20pm UTC](https://discourse.julialang.org/t/differences-in-matlabs-versus-julias-beta-x-y-for-large-x-or-y/60138 "2021-05-01T19:20:06Z")

</div>

Hi everyone, I was comparing a colleague’s MATLAB code to my version (we are coding in parallel, and I am using Julia. I am also trying to convince him to move to Julia. Making progress… :slight\_smile: ), and I noticed …

---

## [gamma(::BigFloat) overflows on Windows](https://discourse.julialang.org/t/gamma-bigfloat-overflows-on-windows/60169)

<div class="topic-metadata">

**Author:** [@DNF](https://discourse.julialang.org/u/DNF)\
**Replies:** 7\
**Last updated:** [April 29, 2021, 2:36pm UTC](https://discourse.julialang.org/t/gamma-bigfloat-overflows-on-windows/60169 "2021-04-29T14:36:24Z")

</div>

jl\> p = big(3.6e13) 3.6e+13 jl\> q = 1/big(0.563483398)+1 2.774675178628776489804929213602311938188089959295862682297487944939957397047059 jl\> beta(p, q) NaN ?? :confused: