BandedBlockBandedMatrix function barrier (#376)

* Disable empty vector in findmindomain/findmaxrange

* Function barrier in BandedBlockBandedMatrix

Author: jishnub

src/Operators/SubOperator.jl

@@ -263,23 +263,25 @@ function BandedBlockBandedMatrix(::Type{Zeros}, S::SubOperator)
                                 rows,cols, (l,u), (λ-jsh,μ+ksh))
 end
 
-function BandedBlockBandedMatrix(::Type{Zeros}, S::SubOperator{<:Any,<:Any,Tuple{BlockRange1,BlockRange1}})
-    KR,JR = parentindices(S)
-    KO = parent(S)
-    l,u = blockbandwidths(KO)::Tuple{Int,Int}
-    λ,μ = subblockbandwidths(KO)::Tuple{Int,Int}
-
-    rt = rangespace(KO)
-    dt = domainspace(KO)
+function _BandedBlockBandedMatrixZeros(::Type{T}, KR, JR, (l,u), (λ,μ), rt, dt) where {T}
     J = first(JR)
     K = first(KR)
     bl_sh = Int(J) - Int(K)
 
     KBR = blocklengthrange(rt,KR)
     KJR = blocklengthrange(dt,JR)
 
-    BandedBlockBandedMatrix(Zeros{eltype(KO)}(sum(KBR),sum(KJR)),
-                                AbstractVector{Int}(KBR),AbstractVector{Int}(KJR), (l+bl_sh,u-bl_sh), (λ,μ))
+    BandedBlockBandedMatrix(Zeros{T}(sum(KBR),sum(KJR)),
+                                convert(AbstractVector{Int}, KBR),
+                                convert(AbstractVector{Int}, KJR),
+                                (l+bl_sh,u-bl_sh), (λ,μ))
+end
+function BandedBlockBandedMatrix(::Type{Zeros}, S::SubOperator{<:Any,<:Any,Tuple{BlockRange1,BlockRange1}})
+    KR,JR = parentindices(S)
+    KO = parent(S)
+    l,u = blockbandwidths(KO)::Tuple{Int,Int}
+    λ,μ = subblockbandwidths(KO)::Tuple{Int,Int}
+    _BandedBlockBandedMatrixZeros(eltype(KO), KR, JR, (l,u), (λ,μ), rangespace(KO), domainspace(KO))
 end
 
 

src/Operators/general/algebra.jl

@@ -439,10 +439,10 @@ for TYP in (:Matrix, :BandedMatrix, :RaggedMatrix)
                     return $TYP(BandedMatrix(V))
                 end
             elseif isbandedblockbanded(P)
-                N = block(rangespace(P), last(parentindices(V)[1]))
-                M = block(domainspace(P), last(parentindices(V)[2]))
-                B = P[Block(1):N, Block(1):M]
-                return $TYP(view(B, parentindices(V)...), _colstops(V))
+                N = block(rangespace(P), last(parentindices(V)[1]))::Block{1,Int}
+                M = block(domainspace(P), last(parentindices(V)[2]))::Block{1,Int}
+                B = P[Block(1):N, Block(1):M]::AbstractMatrix{T}
+                return $TYP(view(B, parentindices(V)...), _colstops(V))::$TYP{T}
             end
 
             kr, jr = parentindices(V)