Tests for type-inference in kronecker operator indexing (#288)

* Tests for type-inference in kronecker operator indexing

* Import BlockArrays

* Import BlockBandedMatrices

Author: jishnub

Project.toml

@@ -27,7 +27,7 @@ Static = "aedffcd0-7271-4cad-89d0-dc628f76c6d3"
 ApproxFunOrthogonalPolynomialsStaticExt = "Static"
 
 [compat]
-ApproxFunBase = "0.8.39, 0.9"
+ApproxFunBase = "0.8.47, 0.9.2"
 ApproxFunBaseTest = "0.1"
 Aqua = "0.6"
 BandedMatrices = "0.16, 0.17"

test/JacobiTest.jl

@@ -11,6 +11,8 @@ using ApproxFunBaseTest: testbandedbelowoperator, testbandedoperator, testspace,
                     testfunctional
 using ApproxFunOrthogonalPolynomials: jacobip
 using BandedMatrices
+using BlockArrays
+using BlockBandedMatrices
 using StaticArrays: SVector
 using Static
 using HalfIntegers
@@ -841,12 +843,17 @@ include("testutils.jl")
             testspaces(d1, d2, r1, r2)
         end
 
-        d1, r1 = Legendre(), Jacobi(2,2)
-        d2, r2 = Chebyshev(), Chebyshev()
-        K = (Operator(I, d1) ⊗ Operator(I, d2)) → (r1 ⊗ r2)
-        M = K[1:20, 1:20]
-        for ind in CartesianIndices(M)
-            @test K[ind] ≈ M[ind]
+        @testset "Indexing" begin
+            d1, r1 = Legendre(), Jacobi(2,2)
+            d2, r2 = Chebyshev(), Chebyshev()
+            K = (Operator(I, d1) ⊗ Operator(I, d2)) → (r1 ⊗ r2)
+            M = @inferred K[1:20, 1:20]
+            for ind in CartesianIndices(M)
+                @test K[ind] ≈ M[ind]
+            end
+            Kv = view(K, Block(1):Block(1), Block(1):Block(1))
+            Y = @inferred BandedBlockBandedMatrix(Kv)
+            @test Y == reshape([K[1,1]], 1, 1)
         end
     end
 end