Tests #182 (#188)

* Tests #182

* v0.13.1

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "ClassicalOrthogonalPolynomials"
 uuid = "b30e2e7b-c4ee-47da-9d5f-2c5c27239acd"
 authors = ["Sheehan Olver <[email protected]>"]
-version = "0.13"
+version = "0.13.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"
@@ -39,7 +39,7 @@ HypergeometricFunctions = "0.3.4"
 InfiniteArrays = " 0.14"
 InfiniteLinearAlgebra = "0.8"
 IntervalSets = "0.7"
-LazyArrays = "2"
+LazyArrays = "2.0.3"
 LazyBandedMatrices = "0.10"
 QuasiArrays = "0.11"
 SpecialFunctions = "1.0, 2"

test/test_normalized.jl

@@ -286,4 +286,10 @@ import ContinuumArrays: MappedWeightedBasisLayout
         g = ApplyQuasiArray(*, W, [1:3; zeros(∞)])
         @test P \ g ≈ transform(P, x -> g[x])
     end
+
+    @testset "inv bug (#182)" begin
+        P = Jacobi(2.0, 0.5)
+        Q = Jacobi(3.0, 0.5)
+        @test (P \ Normalized(Q))[1:10,1:10] ≈ inv((Normalized(Q) \ P)[1:10,1:10])
+    end
 end