From 603e40622b6913f2694e25c857cc4563ea1e333c Mon Sep 17 00:00:00 2001
From: Zedong Peng <peng_zedong@126.com>
Date: Thu, 6 Aug 2026 14:22:46 -0400
Subject: [PATCH] fix: add relative termination for power iteration

---
 src/utils.cu | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/utils.cu b/src/utils.cu
index 4d263b3..2f450e4 100644
--- a/src/utils.cu
+++ b/src/utils.cu
@@ -157,7 +157,7 @@ double estimate_maximum_singular_value(cusparseHandle_t sparse_handle,
         double residual_norm;
         CUBLAS_CHECK(cublasDnrm2_v2_64(blas_handle, m, next_eigenvector_d, 1, &residual_norm));
 
-        if (residual_norm < tolerance)
+        if (residual_norm < tolerance * fmin(1.0, sigma_max_sq))
             break;
     }