--- a/OptimLib3.mod
+++ b/OptimLib3.mod
@@ -65,16 +65,12 @@
FROM LMathLib IMPORT MachEps,MachEpsR4,GoldSec,CGold,Funktion1,FunktionN,sqr;
IMPORT TIO;
-FROM MatLib IMPORT ENorm;
-FROM LibDBlasM2 IMPORT dcopy,ddot;
-
-
CONST reltest = 10.0; (* A relative size used to check equality of *)
(* numbers. Numbers x and y are considered equal if *)
(* the floating-point representation of reltest+x *)
(* equals that of reltest+y. *)
- debug = FALSE;
+CONST debug = FALSE;
TermOut = TRUE;
@@ -124,7 +120,7 @@
END;
x:=x + d;
fu := f(x); dfu := df(x); INC(cnt);
- IF (dfu >= 0) OR ((fu >= fa) AND (dfa <= 0)) THEN
+ IF (dfu >= 0.0) OR ((fu >= fa) AND (dfa <= 0.0)) THEN
b := x; fb := fu; dfb := dfu;
ELSE
a := x; fa := fu; dfa := dfu;
@@ -444,7 +440,7 @@
IF (count < n) THEN
newstep := 2.0*((f - Fmin) - gradproj*steplength);
- IF (newstep > 0) THEN
+ IF (newstep > 0.0) THEN
newstep := -gradproj*sqr(steplength)/newstep;
FOR i:=0 TO n-1 DO
Bvec[i] := X[i] + newstep*t^[i];