FastTransforms v0.17 (#223)

* FastTransforms v0.17

* Update jacobi.jl

* v0.14.1

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.14.0"
+version = "0.14.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -42,7 +42,7 @@ DomainSets = "0.6, 0.7"
 DynamicPolynomials = "0.6"
 FFTW = "1.1"
 FastGaussQuadrature = "1"
-FastTransforms = "0.16.6"
+FastTransforms = "0.16.6, 0.17"
 FillArrays = "1"
 HypergeometricFunctions = "0.3.4"
 InfiniteArrays = " 0.14, 0.15"
@@ -53,7 +53,7 @@ LazyBandedMatrices = "0.10, 0.11"
 MutableArithmetics = "1"
 QuasiArrays = "0.11"
 RecurrenceRelationshipArrays = "0.1.2"
-RecurrenceRelationships = "0.1.1"
+RecurrenceRelationships = "0.1.1, 0.2"
 SpecialFunctions = "1.0, 2"
 julia = "1.10"
 

src/classical/jacobi.jl

@@ -217,8 +217,8 @@ basis_singularities(w::JacobiWeight) = Weighted(Jacobi(w.a, w.b))
 computes the `n`-th Jacobi polynomial, orthogonal with
 respec to `(1-x)^a*(1+x)^b`, at `z`.
 """
-jacobip(n::Integer, a, b, z) = Base.unsafe_getindex(Jacobi{polynomialtype(typeof(a), typeof(b), typeof(z))}(a,b), z, n+1)
-normalizedjacobip(n::Integer, a, b, z) = Base.unsafe_getindex(Normalized(Jacobi{polynomialtype(typeof(a), typeof(b), typeof(z))}(a,b)), z, n+1)
+jacobip(n::Integer, a, b, z) = Base.unsafe_getindex(Jacobi{polynomialtype(promote_type(typeof(a), typeof(b)), typeof(z))}(a,b), z, n+1)
+normalizedjacobip(n::Integer, a, b, z) = Base.unsafe_getindex(Normalized(Jacobi{polynomialtype(promote_type(typeof(a), typeof(b)), typeof(z))}(a,b)), z, n+1)
 
 OrthogonalPolynomial(w::JacobiWeight) = Jacobi(w.a, w.b)
 orthogonalityweight(P::Jacobi) = JacobiWeight(P.a, P.b)