fix tests

Author: dlfivefifty

Project.toml

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

src/choleskyQR.jl

@@ -26,7 +26,7 @@ orthogonalityweight(Q::ConvertedOrthogonalPolynomial) = Q.weight
 # transform to P * inv(U) if needed for differentiation, etc.
 arguments(::ApplyLayout{typeof(*)}, Q::ConvertedOrthogonalPolynomial) = Q.P, ApplyArray(inv, Q.U)
 
-OrthogonalPolynomial(w::AbstractQuasiVector) = OrthogonalPolynomial(w, orthogonalpolynomial(axes(w,1), singularities(w)))
+OrthogonalPolynomial(w::AbstractQuasiVector) = OrthogonalPolynomial(w, orthogonalpolynomial_axis(axes(w,1), singularities(w)))
 function OrthogonalPolynomial(w::AbstractQuasiVector, P::AbstractQuasiMatrix)
     Q = normalized(P)
     X,U = cholesky_jacobimatrix(w, Q)
@@ -39,8 +39,8 @@ orthogonalpolynomial(w::SubQuasiArray) = orthogonalpolynomial(parent(w))[parenti
 OrthogonalPolynomial(w::Function, P::AbstractQuasiMatrix) = OrthogonalPolynomial(w.(axes(P,1)), P)
 orthogonalpolynomial(w::Function, P::AbstractQuasiMatrix) = orthogonalpolynomial(w.(axes(P,1)), P)
 
-orthogonalpolynomial(ax, ::NoSingularities) = legendre(ax)
-orthogonalpolynomial(ax, w) = orthogonalpolynomial(w)
+orthogonalpolynomial_axis(ax, ::NoSingularities) = legendre(ax)
+orthogonalpolynomial_axis(ax, w) = orthogonalpolynomial(w)
 resizedata!(P::ConvertedOrthogonalPolynomial, ::Colon, n::Int) = resizedata!(P.X.dv, n)
 
 

src/classical/jacobi.jl

@@ -92,7 +92,7 @@ hasboundedendpoints(w::AbstractJacobiWeight) = w.a ≥ 0 && w.b ≥ 0
 singularities(a::AbstractAffineQuasiVector) = singularities(a.x)
 
 
-singularities(w::JacobiWeight) = w
+singularities(w::AbstractJacobiWeight) = w
 
 ## default is to just assume no singularities
 singularitiesbroadcast(_...) = NoSingularities()

src/lanczos.jl

@@ -165,7 +165,7 @@ function LanczosPolynomial(w_in::AbstractQuasiVector, P::AbstractQuasiMatrix)
     LanczosPolynomial(w, Q, LanczosData(w, Q))
 end
 
-LanczosPolynomial(w::AbstractQuasiVector) = LanczosPolynomial(w, normalized(orthogonalpolynomial(axes(w,1), singularities(w))))
+LanczosPolynomial(w::AbstractQuasiVector) = LanczosPolynomial(w, normalized(orthogonalpolynomial_axis(axes(w,1), singularities(w))))
 
 ==(A::LanczosPolynomial, B::LanczosPolynomial) = A.w == B.w
 ==(::LanczosPolynomial, ::OrthogonalPolynomial) = false # TODO: fix

test/test_monic.jl

@@ -48,7 +48,6 @@ using ClassicalOrthogonalPolynomials: Monic, _p0, orthogonalityweight, recurrenc
             Hermite() _Hermite
         ]
         for (P, _P) in eachrow(Ps)
-            @show P
             Q = Monic(P)
             @test Q[0.2, 1] == 1.0
             @test Q[0.25, 2] ≈ _P(0.25, 1)