diff --git a/Project.toml b/Project.toml index cc18e35..9f18b9b 100644 --- a/Project.toml +++ b/Project.toml @@ -1,6 +1,6 @@ name = "BlockBandedMatrices" uuid = "ffab5731-97b5-5995-9138-79e8c1846df0" -version = "0.13.5" +version = "0.13.6" [deps] ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a" diff --git a/src/BlockSkylineMatrix.jl b/src/BlockSkylineMatrix.jl index 633ffa6..833c6fd 100644 --- a/src/BlockSkylineMatrix.jl +++ b/src/BlockSkylineMatrix.jl @@ -5,7 +5,7 @@ checkbandwidths(N, M, l::AbstractVector{Int}, u::AbstractVector{Int}) = #### Routines for BandedSizes function bb_blockstarts(ax, l::AbstractVector{Int}, u::AbstractVector{Int}) N,M = blocksize.(ax,1) - L,U = maximum(l), maximum(u) + L,U = maximum(l; init=-720), maximum(u; init=-720) b_start = BandedMatrix{Int}(undef, (N, M), (L, U)) -L > U && return b_start @@ -28,7 +28,7 @@ end function bb_blockstrides(b_axes, l::AbstractVector{Int}, u::AbstractVector{Int}) N, M = blocksize.(b_axes,1) - L,U = maximum(l), maximum(u) + L,U = maximum(l; init=-720), maximum(u; init=-720) checkbandwidths(N, M, l, u) b_strides = Vector{Int}(undef, M) for J=1:M @@ -368,7 +368,7 @@ axes(A::BlockSkylineMatrix) = A.block_sizes.axes MemoryLayout(::Type{<:BlockSkylineMatrix}) = BlockBandedColumnMajor() colblockbandwidths(A::BlockSkylineMatrix) = (A.block_sizes.l, A.block_sizes.u) -blockbandwidths(A::BlockSkylineMatrix) = maximum.(colblockbandwidths(A)) +blockbandwidths(A::BlockSkylineMatrix) = maximum.(colblockbandwidths(A); init=-720) BroadcastStyle(::Type{<:BlockSkylineMatrix}) = BlockSkylineStyle() BroadcastStyle(::Type{<:BlockBandedMatrix}) = BlockBandedStyle()