Product of QuasiDiagonals fails

[jagot]

I discovered the following while reviving some untested functionality in CompactBases (see also JuliaApproximation/QuasiArrays.jl#32):

julia> using IntervalSets, QuasiArrays, ContinuumArrays

julia> r = Inclusion(0.0..1.0)
Inclusion(0.0..1.0)

julia> Q = QuasiDiagonal(r)
QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}

julia> Q*Q
ERROR: MethodError: no method matching Int64(::Infinities.InfiniteCardinal{1})
Closest candidates are:
  (::Type{T})(::T) where T<:Number at boot.jl:763
  (::Type{T})(::AbstractChar) where T<:Union{Int32, Int64} at char.jl:51
  (::Type{T})(::AbstractChar) where T<:Union{AbstractChar, Number} at char.jl:50
  ...
Stacktrace:
  [1] _similar_for(c::UnitRange{Int64}, #unused#::Type{Float64}, itr::Inclusion{Float64, ClosedInterval{Float64}}, #unused#::Base.HasLength)
    @ Base ./array.jl:576
  [2] _collect(cont::UnitRange{Int64}, itr::Inclusion{Float64, ClosedInterval{Float64}}, #unused#::Base.HasEltype, isz::Base.HasLength)
    @ Base ./array.jl:609
  [3] collect(itr::Inclusion{Float64, ClosedInterval{Float64}})
    @ Base ./array.jl:603
  [4] map
    @ ./tuple.jl:214 [inlined]
  [5] QuasiArray(a::ApplyQuasiMatrix{Float64, typeof(*), Tuple{QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}}})
    @ QuasiArrays ~/.julia/dev/QuasiArrays/src/quasiarray.jl:84
  [6] QuasiArray(M::ArrayLayouts.Mul{QuasiArrays.QuasiArrayLayout, QuasiArrays.QuasiArrayLayout, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}})
    @ QuasiArrays ~/.julia/dev/QuasiArrays/src/lazyquasiarrays.jl:214
  [7] _quasi_mul(M::ArrayLayouts.Mul{QuasiArrays.QuasiArrayLayout, QuasiArrays.QuasiArrayLayout, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}}, #unused#::Tuple{Inclusion{Float64, ClosedInterval{Float64}}, Inclusion{Float64, ClosedInterval{Float64}}})
    @ QuasiArrays ~/.julia/dev/QuasiArrays/src/matmul.jl:94
  [8] copy(M::ArrayLayouts.Mul{QuasiArrays.QuasiArrayLayout, QuasiArrays.QuasiArrayLayout, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}})
    @ QuasiArrays ~/.julia/dev/QuasiArrays/src/matmul.jl:96
  [9] materialize(M::ArrayLayouts.Mul{QuasiArrays.QuasiArrayLayout, QuasiArrays.QuasiArrayLayout, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}})
    @ ArrayLayouts ~/.julia/packages/ArrayLayouts/wBWqq/src/mul.jl:105
 [10] mul
    @ ~/.julia/packages/ArrayLayouts/wBWqq/src/mul.jl:106 [inlined]
 [11] *(A::QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}}, B::QuasiDiagonal{Float64, Inclusion{Float64, ClosedInterval{Float64}}})
    @ QuasiArrays ~/.julia/dev/QuasiArrays/src/matmul.jl:23
 [12] top-level scope
    @ REPL[10]:1

This is on QuasiArrays 0.5.1 and dl/QuasiArrays-v0.5 of ContinuumArrays.