From 411574a39c53b16d7056883dd3b8d207176d70e6 Mon Sep 17 00:00:00 2001 From: Germain Clavier Date: Fri, 26 Apr 2024 14:46:30 +0200 Subject: [PATCH] Changed alpha_init initialization to avoid infinite loop with 0 starting value. --- src/min_linesearch.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/min_linesearch.cpp b/src/min_linesearch.cpp index 24ba4c5c23..4dfeee6823 100644 --- a/src/min_linesearch.cpp +++ b/src/min_linesearch.cpp @@ -690,7 +690,7 @@ int MinLineSearch::linemin_forcezero(double eoriginal, double &alpha) // choosing the initial alpha that we'll use // rough estimate that'll decrease energy to 1/10 - alpha_init = 0.1*fabs(eoriginal)/fdothall; + alpha_init = MAX(EPS_QUAD, 0.1*fabs(eoriginal)/fdothall); // initialize aplha to 0.0