Split chop! and resize! (#101)

Author: dlfivefifty

Project.toml

@@ -1,6 +1,6 @@
 name = "InfiniteLinearAlgebra"
 uuid = "cde9dba0-b1de-11e9-2c62-0bab9446c55c"
-version = "0.6.6"
+version = "0.6.7"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/InfiniteLinearAlgebra.jl

@@ -74,10 +74,11 @@ function choplength(c::AbstractVector, tol)
     return 0
 end
 
-_chop!(_, c::AbstractVector, tol::Real) = resize!(c, choplength(c, tol))
-_chop!(ax::BlockedUnitRange, c::AbstractVector, tol::Real) = resize!(c, findblock(ax, choplength(c, tol)))
-
-chop!(c::AbstractVector{T}, tol::Real=zero(real(T))) where T = _chop!(axes(c,1), c, tol)
+# resize! to nearest block
+compatible_resize!(_, c::AbstractVector, n) = resize!(c, n)
+compatible_resize!(ax::BlockedUnitRange, c::AbstractVector, n) = resize!(c, iszero(n) ? Block(0) : findblock(ax, n))
+compatible_resize!(c, n) = compatible_resize!(axes(c,1), c, n)
+chop!(c::AbstractVector{T}, tol::Real=zero(real(T))) where T = compatible_resize!(c, choplength(c, tol))
 
 function chop(A::AbstractMatrix{T}, tol::Real=zero(real(T))) where T
     for k = size(A,1):-1:1