Test KronTrav bug

Author: dlfivefifty

Project.toml

@@ -36,7 +36,7 @@ FillArrays = "1.0"
 HarmonicOrthogonalPolynomials = "0.6"
 InfiniteArrays = "0.15"
 InfiniteLinearAlgebra = "0.9"
-LazyArrays = "2.0"
+LazyArrays = "2.3.1"
 LazyBandedMatrices = "0.11"
 QuasiArrays = "0.11"
 SpecialFunctions = "1, 2"

test/test_rect.jl

@@ -1,4 +1,4 @@
-using MultivariateOrthogonalPolynomials, ClassicalOrthogonalPolynomials, StaticArrays, LinearAlgebra, BlockArrays, FillArrays, Base64, Test
+using MultivariateOrthogonalPolynomials, ClassicalOrthogonalPolynomials, StaticArrays, LinearAlgebra, BlockArrays, FillArrays, Base64, LazyBandedMatrices, Test
 using ClassicalOrthogonalPolynomials: expand
 using MultivariateOrthogonalPolynomials: weaklaplacian
 using ContinuumArrays: plotgridvalues
@@ -146,4 +146,12 @@ using ContinuumArrays: plotgridvalues
         f = expand(P, splat((x,y) -> exp(cos(x^2*y))))
         @test sum(f) ≈ 10.546408460894801 # empirical
     end
+
+    @testset "KronTrav bug" begin
+        W = Weighted(Ultraspherical(3/2))
+        D² = diff(W)'diff(W)
+        M = W'W
+        Dₓ = KronTrav(D²,M)
+        @test Dₓ[Block.(1:1),Block.(1:1)] == Dₓ[Block(1,1)]
+    end
 end