BUG: in autoscaling

ANTsX · Aug 9, 2024 · e78d73a · e78d73a
1 parent 4291af0
commit e78d73a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/multiscaleSVDxpts.R b/R/multiscaleSVDxpts.R
@@ -3504,7 +3504,7 @@ simlr <- function(
       )
     }
     myorthEnergy = invariant_orthogonality_defect( myenergysearchv )
-    if ( last_energy > 0 )
+    if ( last_energy > .Machine$double.eps & myorthEnergy > .Machine$double.eps)
       myorthEnergy = 0.5 * myorthEnergy*(0.5 * last_energy/myorthEnergy)
     if (ccaEnergy) {
       # ( v'*X'*Y )/( norm2(X*v ) * norm2( u ) )