patch for vector valued problems

- fixes issues when pivot is not full rank

Author: maltezfaria

Project.toml

@@ -1,6 +1,6 @@
 name = "HMatrices"
 uuid = "8646bddf-ab1c-4fa7-9c51-ba187d647618"
-version = "0.2.11"
+version = "0.2.12"
 
 [deps]
 Distributed = "8ba89e20-285c-5b6f-9357-94700520ee1b"

src/compressor.jl

@@ -143,8 +143,8 @@ function _aca_partial(K, irange, jrange, atol, rmax, rtol, istart, buffer_ = not
                 break
             end
         else # δ != 0
-            iδ = inv(δ)
-            # rdiv!(b,δ) # b <-- b/δ
+            iδ = T <: AbstractArray ? pinv(δ; rtol = 1e-8) : inv(δ)
+            # b <-- b/δ
             for k in eachindex(b)
                 b[k] = b[k] * iδ
             end