fixed linesearch equation

jacobwilliams · Dec 14, 2021 · e7c6f5c · e7c6f5c
1 parent a9b0a86
commit e7c6f5c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/nlesolver_module.f90 b/src/nlesolver_module.f90
@@ -783,7 +783,7 @@ subroutine backtracking_linesearch(me,xold,p,fjac,x,f,fvec)
             end if
         end if
 
-        if ((f - 0.5_wp * ftmp >= alpha*t) .or. min_alpha_reached) then
+        if (((f - ftmp) / 2.0_wp >= alpha*t) .or. min_alpha_reached) then
             if (min_alpha_reached) then
                 write(me%iunit,'(A)') '        Minimum alpha reached'
             end if