Move RecurrenceArray to RecurrenceRelationships.jl (#48)

* Move RecurrenceArray to RecurrenceRelationships.jl

* Update Project.toml

* Update SingularIntegrals.jl

* tests pass

* Support C^(-1/2)

* Update Project.toml

* Update ci.yml

* v0.3.3

* Update ci.yml

Author: dlfivefifty

.github/workflows/ci.yml

@@ -20,15 +20,16 @@ jobs:
       fail-fast: false
       matrix:
         version:
-          - '1.10'
+          - 'lts'
+          - '1'
         os:
           - ubuntu-latest
           - macOS-latest
           - windows-latest
         arch:
           - x64
     steps:
-      - uses: actions/checkout@v2
+      - uses: actions/checkout@v4
       - uses: julia-actions/setup-julia@v2
         with:
           version: ${{ matrix.version }}
@@ -48,4 +49,5 @@ jobs:
       - uses: julia-actions/julia-processcoverage@v1
       - uses: codecov/codecov-action@v4
         with:
+          token: ${{ secrets.CODECOV_TOKEN }}
           file: lcov.info

Project.toml

@@ -1,7 +1,7 @@
 name = "SingularIntegrals"
 uuid = "d7440221-8b5e-42fc-909c-0567823f424a"
 authors = ["Sheehan Olver <solver@mac.com>"]
-version = "0.3.2"
+version = "0.3.3"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -16,20 +16,24 @@ LazyArrays = "5078a376-72f3-5289-bfd5-ec5146d43c02"
 LazyBandedMatrices = "d7e5e226-e90b-4449-9968-0f923699bf6f"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 QuasiArrays = "c4ea9172-b204-11e9-377d-29865faadc5c"
+RecurrenceRelationships = "807425ed-42ea-44d6-a357-6771516d7b2c"
+RecurrenceRelationshipArrays = "b889d2dc-af3c-4820-88a8-238fa91d3518"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 
 [compat]
 ArrayLayouts = "1.4"
 BandedMatrices = "1"
-ClassicalOrthogonalPolynomials = "0.12, 0.13"
-ContinuumArrays = "0.17.2, 0.18"
-FastTransforms = "0.15, 0.16"
+ClassicalOrthogonalPolynomials = "0.13.7"
+ContinuumArrays = "0.18"
+FastTransforms = "0.16"
 FillArrays = "1"
 HypergeometricFunctions = "0.3.4"
-InfiniteArrays = "0.13, 0.14"
-LazyArrays = "1.10, 2"
-LazyBandedMatrices = "0.9, 0.10"
+InfiniteArrays = "0.14"
+LazyArrays = "2"
+LazyBandedMatrices = "0.10"
 QuasiArrays = "0.11"
+RecurrenceRelationships = "0.1"
+RecurrenceRelationshipArrays = "0.1"
 SpecialFunctions = "1, 2"
 julia = "1.10"
 

src/SingularIntegrals.jl

@@ -2,18 +2,18 @@ module SingularIntegrals
 using ClassicalOrthogonalPolynomials, ContinuumArrays, QuasiArrays, LazyArrays, LazyBandedMatrices, FillArrays, BandedMatrices, LinearAlgebra, SpecialFunctions, HypergeometricFunctions, InfiniteArrays
 using ContinuumArrays: @simplify, Weight, AbstractAffineQuasiVector, inbounds_getindex, broadcastbasis, MappedBasisLayouts, MemoryLayout, MappedWeightLayout, AbstractWeightLayout, ExpansionLayout, demap, basismap, AbstractBasisLayout, SubBasisLayout
 using QuasiArrays: AbstractQuasiMatrix, BroadcastQuasiMatrix, LazyQuasiArrayStyle, AbstractQuasiVecOrMat
-import ClassicalOrthogonalPolynomials: AbstractJacobiWeight, WeightedBasis, jacobimatrix, orthogonalityweight, recurrencecoefficients, _p0, Clenshaw, chop, initiateforwardrecurrence, MappedOPLayouts, unweighted, WeightedOPLayout, MappedOPLayout
+import ClassicalOrthogonalPolynomials: AbstractJacobiWeight, AbstractJacobi, WeightedBasis, jacobimatrix, orthogonalityweight, recurrencecoefficients, _p0, chop, initiateforwardrecurrence, MappedOPLayouts, unweighted, WeightedOPLayout, MappedOPLayout
 using LazyBandedMatrices: Tridiagonal, SymTridiagonal, subdiagonaldata, supdiagonaldata, diagonaldata, ApplyLayout
 import LazyArrays: AbstractCachedMatrix, AbstractCachedArray, paddeddata, arguments, resizedata!, cache_filldata!, zero!, cacheddata, LazyArrayStyle
 import Base: *, +, -, /, \, Slice, axes, getindex, sum, ==, oneto, size, broadcasted, copy, tail, view
 import LinearAlgebra: dot
 using BandedMatrices: _BandedMatrix
-using FastTransforms: _forwardrecurrence!, _forwardrecurrence_next
+using RecurrenceRelationshipArrays
+using RecurrenceRelationshipArrays: Clenshaw
 
 export associated, stieltjes, logkernel, powerkernel, complexlogkernel
 
 
-include("recurrence.jl")
 include("stieltjes.jl")
 include("logkernel.jl")
 include("power.jl")

src/logkernel.jl

@@ -166,6 +166,8 @@ function complexlogkernel_recurrence(wP::Weighted{<:Any,<:Legendre})
 end
 complexlogkernel(P::Legendre, z...) = complexlogkernel(Weighted(P), z...)
 logkernel(P::Legendre, x...) = logkernel(Weighted(P), x...)
+complexlogkernel(J::AbstractJacobi{T}, z...) where T = complexlogkernel(Legendre{T}(), z...) * (Legendre{T}() \ J)
+logkernel(J::AbstractJacobi{T}, z...) where T = logkernel(Legendre{T}(), z...) * (Legendre{T}() \ J)
 
 
 

src/recurrence.jl

@@ -114,8 +114,6 @@ function stieltjes(wP::Weighted{<:Any,<:OrthogonalPolynomial})
     (-A[1]*sum(w))*[zero(axes(P,1)) Q] + stieltjes(w) .* P
 end
 
-stieltjes(P::Legendre) = stieltjes(Weighted(P))
-
 
 ##
 # OffStieltjes
@@ -199,6 +197,7 @@ sqrtx2(x::Real) = sign(x)*sqrt(x^2-1)
 
 
 stieltjes(P::Legendre, z...) = stieltjes(Weighted(P), z...)
+stieltjes(J::AbstractJacobi{T}, z...) where T = stieltjes(Legendre{T}(), z...) * (Legendre{T}() \ J)
 
 @simplify function *(S::StieltjesPoints, wP::Weighted)
     z = S.args[1].args[1] # vector of points to eval at

src/stieltjes.jl

@@ -4,7 +4,6 @@ using SingularIntegrals: Stieltjes, StieltjesPoint, ChebyshevInterval, associate
 using LazyArrays: MemoryLayout, PaddedLayout, colsupport, rowsupport, paddeddata
 using LazyBandedMatrices: blockcolsupport, Block, BlockHcat, blockbandwidths
 
-include("test_recurrence.jl")
 
 @testset "Associated" begin
     T = ChebyshevT()

test/runtests.jl

@@ -24,7 +24,7 @@ using ClassicalOrthogonalPolynomials: affine
 
     @testset "expand" begin
         @test complexlogkernel(exp.(x), 2 + im) ≈ sum(log.((2+im) .- x) .* exp.(x))
-        @test_throws ErrorException complexlogkernel(Jacobi(0.1,0.2), 2+im)
+        @test_throws InexactError complexlogkernel(Jacobi(0.1,0.2), 2+im)
     end
 end