🤖 Format .jl files

Author: dpo

src/utils.jl

@@ -18,28 +18,29 @@ function opnorm_eig(B; max_attempts::Int = 3)
   ncv = max(20, 2 * nev + 1)
 
   while !(have_eig || attempt >= max_attempts)
-      attempt += 1
-      try
-          # Estimate largest eigenvalue in absolute value
-          d, nconv, niter, nmult, resid = eigs(B; nev = nev, ncv = ncv, which = :LM, ritzvec = false, check = 1)
-          
-          # Check if eigenvalue has converged
-          have_eig = nconv == 1
-          if have_eig
-              λ = abs(d[1])  # Take absolute value of the largest eigenvalue
-              break  # Exit loop if successful
-          else
-              # Increase NCV for the next attempt if convergence wasn't achieved
-              ncv = min(2 * ncv, n)
-          end
-      catch e
-          if occursin("XYAUPD_Exception", string(e)) && ncv < n
-              @warn "Arpack error: $e. Increasing NCV to $ncv and retrying."
-              ncv = min(2 * ncv, n)  # Increase NCV but don't exceed matrix size
-          else
-              rethrow(e)  # Re-raise if it's a different error
-          end
+    attempt += 1
+    try
+      # Estimate largest eigenvalue in absolute value
+      d, nconv, niter, nmult, resid =
+        eigs(B; nev = nev, ncv = ncv, which = :LM, ritzvec = false, check = 1)
+
+      # Check if eigenvalue has converged
+      have_eig = nconv == 1
+      if have_eig
+        λ = abs(d[1])  # Take absolute value of the largest eigenvalue
+        break  # Exit loop if successful
+      else
+        # Increase NCV for the next attempt if convergence wasn't achieved
+        ncv = min(2 * ncv, n)
+      end
+    catch e
+      if occursin("XYAUPD_Exception", string(e)) && ncv < n
+        @warn "Arpack error: $e. Increasing NCV to $ncv and retrying."
+        ncv = min(2 * ncv, n)  # Increase NCV but don't exceed matrix size
+      else
+        rethrow(e)  # Re-raise if it's a different error
       end
+    end
   end
 
   return λ, have_eig
@@ -54,28 +55,28 @@ function opnorm_svd(J; max_attempts::Int = 3)
   ncv = 10
 
   while !(have_svd || attempt >= max_attempts)
-      attempt += 1
-      try
-          # Estimate largest singular value
-          s, nconv, niter, nmult, resid = svds(J; nsv = nsv, ncv = ncv, ritzvec = false, check = 1)
-          
-          # Check if singular value has converged
-          have_svd = nconv >= 1
-          if have_svd
-              σ = maximum(s.S)  # Take the largest singular value
-              break  # Exit loop if successful
-          else
-              # Increase NCV for the next attempt if convergence wasn't achieved
-              ncv = min(2 * ncv, n)
-          end
-      catch e
-          if occursin("XYAUPD_Exception", string(e)) && ncv < n
-              @warn "Arpack error: $e. Increasing NCV to $ncv and retrying."
-              ncv = min(2 * ncv, n)  # Increase NCV but don't exceed matrix size
-          else
-              rethrow(e)  # Re-raise if it's a different error
-          end
+    attempt += 1
+    try
+      # Estimate largest singular value
+      s, nconv, niter, nmult, resid = svds(J; nsv = nsv, ncv = ncv, ritzvec = false, check = 1)
+
+      # Check if singular value has converged
+      have_svd = nconv >= 1
+      if have_svd
+        σ = maximum(s.S)  # Take the largest singular value
+        break  # Exit loop if successful
+      else
+        # Increase NCV for the next attempt if convergence wasn't achieved
+        ncv = min(2 * ncv, n)
+      end
+    catch e
+      if occursin("XYAUPD_Exception", string(e)) && ncv < n
+        @warn "Arpack error: $e. Increasing NCV to $ncv and retrying."
+        ncv = min(2 * ncv, n)  # Increase NCV but don't exceed matrix size
+      else
+        rethrow(e)  # Re-raise if it's a different error
       end
+    end
   end
 
   return σ, have_svd