Fix #45 (#46)

* Fix #45

* Update Project.toml

* Update Project.toml

* Update Project.toml

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <solver@mac.com>"]
-version = "0.4.2"
+version = "0.4.3"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -29,8 +29,8 @@ ArrayLayouts = "0.7"
 BandedMatrices = "0.16.5"
 BlockArrays = "0.15, 0.16"
 BlockBandedMatrices = "0.10"
-ContinuumArrays = "0.8.1"
-DomainSets = "0.4, 0.5"
+ContinuumArrays = "0.8.5"
+DomainSets = "0.5"
 FFTW = "1.1"
 FastGaussQuadrature = "0.4.3"
 FastTransforms = "0.12"
@@ -39,9 +39,9 @@ HypergeometricFunctions = "0.3.4"
 InfiniteArrays = "0.11.1"
 InfiniteLinearAlgebra = "0.5.8"
 IntervalSets = "0.5"
-LazyArrays = "0.21"
+LazyArrays = "0.21.9"
 LazyBandedMatrices = "0.6"
-QuasiArrays = "0.6"
+QuasiArrays = "0.7"
 SpecialFunctions = "1.0"
 julia = "1.6"
 

test/test_chebyshev.jl

@@ -117,6 +117,13 @@ import ContinuumArrays: MappedWeightedBasisLayout, Map
         U = chebyshevu(0..1)
         @test (U*(U\x))[0.1] ≈ 0.1
         @test (U* (U \ exp.(x)))[0.1] ≈ exp(0.1)
+
+        @testset "Trivial map" begin
+            T = ChebyshevT()
+            x = Inclusion(-1..1)
+            @test T == T[x,:] == T[x,1:∞] == T[:,1:∞]
+            @test T[x,:] == T
+        end
     end
 
     @testset "weighted" begin

test/test_lanczos.jl

@@ -231,4 +231,13 @@ import ClassicalOrthogonalPolynomials: recurrencecoefficients, PaddedLayout, ort
         Q̃ = LanczosPolynomial( @.(inv(1+x̃^2)))
         @test all((x,w) .≈ golubwelsch(Q̃[:,Base.OneTo(10)]))
     end
+
+    @testset "ambiguity (#45)" begin
+        x = Inclusion(-1.0..1)
+        a = 1.5
+        ϕ = x.^4 - (a^2 + 1)*x.^2 .+ a^2
+        Pϕ = Normalized(LanczosPolynomial(ϕ))
+        P = Normalized(Legendre())
+        Cϕ = Pϕ\P
+    end
 end