Support Jacobi(-1,0)\Ultraspherical(-1/2)

Author: DanielVandH

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.15.3"
+version = "0.15.4"
 
 
 [deps]

src/classical/ultraspherical.jl

@@ -199,8 +199,16 @@ function \(A::Ultraspherical, B::Jacobi)
     Diagonal(Ã[1,:]./A[1,:]) * (Ã\B)
 end
 function \(A::Jacobi, B::Ultraspherical)
-    B̃ = Jacobi(B)
-    (A\B̃)*Diagonal(B[1,:]./B̃[1,:])
+    if isone(-2B.λ) && (A.a, A.b) == (-1, 0)
+        T = promote_type(eltype(A), eltype(B))
+        n = -2one(T) ./ (2 .* (2:∞) .- one(T))
+        dv = Vcat(one(T), -2one(T), n)
+        ev = Vcat(-one(T), n)
+        LazyBandedMatrices.Bidiagonal(dv, ev, :U)
+    else
+        B̃ = Jacobi(B)
+        (A\B̃)*Diagonal(B[1,:]./B̃[1,:])
+    end 
 end
 
 function \(wA::Weighted{<:Any,<:Ultraspherical}, wB::Weighted{<:Any,<:Jacobi})

test/test_ultraspherical.jl

@@ -207,4 +207,14 @@ using ClassicalOrthogonalPolynomials: grammatrix
     @testset "ChebyshevInterval constructior" begin
         @test ultraspherical(2,ChebyshevInterval()) ≡ Ultraspherical(2)
     end
+end
+
+@testset "Jacobi(-1, 0) \ Ultraspherical(-1/2)" begin 
+    A = Jacobi(-1, 0)
+    B = Ultraspherical(-1/2)
+    R = A \ B 
+    AR = A * R 
+    lhs = AR[0.3918, 1:100]
+    rhs = B[0.3918, 1:100]
+    @test lhs ≈ rhs
 end