Simplify feasibility loop

docs/src/index.md

@@ -27,7 +27,7 @@ A concrete example: a 3×3 matrix whose rows and columns correspond to physical
 variables with different units (position in meters, velocity in m/s, force
 in N):
 
-```jldoctest coverones
+```jldoctest coverones; filter = r"(\d+\.\d{6})\d+" => s"\1"
 julia> using ScaleInvariantAnalysis
 
 julia> A = [1e6 1e3 1.0; 1e3 1.0 1e-3; 1.0 1e-3 1e-6];

src/covers.jl

@@ -151,21 +151,11 @@ function cover(A::AbstractMatrix; kwargs...)
     for j in axes(A, 2)
         for i in axes(A, 1)
             Aij, ai, bj = abs(A[i, j]), a[i], b[j]
-            if iszero(bj)
-                if !iszero(ai)
-                    b[j] = Aij / ai
-                else
-                    a[i] = b[j] = sqrt(Aij)
-                end
-            elseif iszero(ai)
-                a[i] = Aij / bj
-            else
-                aprod = ai * bj
-                aprod >= Aij && continue
-                s = sqrt(Aij / aprod)
-                a[i] = s * ai
-                b[j] = s * bj
-            end
+            aprod = ai * bj
+            aprod >= Aij && continue
+            s = sqrt(Aij / aprod)
+            a[i] = s * ai
+            b[j] = s * bj
         end
     end
     return tighten_cover!(a, b, A; kwargs...)