PERF: rvcoef subtlety

ANTsX · Nov 3, 2024 · 57531cc · 57531cc
1 parent c60cbc8
commit 57531cc
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/R/multiscaleSVDxpts.R b/R/multiscaleSVDxpts.R
@@ -4415,7 +4415,8 @@ rvcoef <- function(X, Y) {
   C <- t(X_norm) %*% Y_norm
 
   # Compute SVD
-  svd_C <- ba_svd(C, dividebymax=TRUE )
+  # svd_C <- corpcor::fast.svd( C )
+  svd_C <- svd( C + diag(1e-8 , nrow(C))  )
 
   # Compute RV coefficient
   sigma_sq <- sum(svd_C$d^2)