Merge pull request #213 from JuliaArrays/ChrisRackauckas-patch-1

Improve banded matrix sparsity calculations

Author: ChrisRackauckas

src/ArrayInterface.jl

@@ -887,7 +887,7 @@ function __init__()
                     (startinds[blockcolors[i]]:endinds[blockcolors[i]])[1:cols[i]]
                     for i = 1:nblock
                 ]
-                return vcat(colors...)
+                return reduce(vcat,colors)
             end
 
             function matrix_colors(A::BlockBandedMatrices.BandedBlockBandedMatrix)
@@ -909,7 +909,7 @@ function __init__()
                     _cycle(startinds[blockcolors[i]]:endinds[blockcolors[i]], cols[i])
                     for i = 1:nblock
                 ]
-                return vcat(colors...)
+                return reduce(vcat,colors)
             end
         end
     end