Handle limit case in LM

Inria-Empenn · Nov 12, 2018 · 8157338 · 8157338
1 parent 2a21515
commit 8157338
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/Anima/math-tools/optimizers/animaBoundedLevenbergMarquardtOptimizer.cxx b/Anima/math-tools/optimizers/animaBoundedLevenbergMarquardtOptimizer.cxx
@@ -1,6 +1,7 @@
 #include <animaBoundedLevenbergMarquardtOptimizer.h>
 #include <animaBaseTensorTools.h>
 #include <vnl_qr.h>
+#include <limits>
 
 namespace anima
 {
@@ -243,7 +244,7 @@ bool BoundedLevenbergMarquardtOptimizer::CheckConditions(unsigned int numIterati
     if (numIterations == m_NumberOfIterations)
         return true;
 
-    if (m_LambdaParameter == 0.0)
+    if ((m_LambdaParameter == 0.0)||(m_LambdaParameter >= std::numeric_limits<double>::max() / 2.0))
         return true;
 
     double normDiffParams = 0.0;