Changeset 12476 for issm/trunk-jpl/src/c/shared/Numerics/BrentSearch.cpp

Timestamp:

06/20/12 12:15:29 (13 years ago)

Author:

utke

Message:

type renames

File:

• issm/trunk-jpl/src/c/shared/Numerics/BrentSearch.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/shared/Numerics/BrentSearch.cpp

@@ -16 +16 @@
 #include <float.h>
 
-void BrentSearch(double* psearch_scalar,double* pJ,OptPars* optpars,double (*f)(double,OptArgs*), OptArgs* optargs){
+void BrentSearch(IssmDouble* psearch_scalar,IssmDouble* pJ,OptPars* optpars,IssmDouble (*f)(IssmDouble,OptArgs*), OptArgs* optargs){
 
         /* This routine is optimizing a given function using Brent's method
@@ -22 +22 @@
 
         /*Intermediary*/
-        double si,gold,intervalgold,oldintervalgold;
-        double parab_num,parab_den;
-        double distance,cm_jump;
-        double fxmax,fxmin,fxbest;
-        double fx,fx1,fx2;
-        double xmax,xmin,xbest;
-        double x,x1,x2,xm;
-        double tol1,tol2,seps;
-        double tolerance=1.e-4;
+        IssmDouble si,gold,intervalgold,oldintervalgold;
+        IssmDouble parab_num,parab_den;
+        IssmDouble distance,cm_jump;
+        IssmDouble fxmax,fxmin,fxbest;
+        IssmDouble fx,fx1,fx2;
+        IssmDouble xmax,xmin,xbest;
+        IssmDouble x,x1,x2,xm;
+        IssmDouble tol1,tol2,seps;
+        IssmDouble tolerance=1.e-4;
         int    maxiter,iter;
         bool   loop=true,goldenflag;
@@ -58 +58 @@
 
         /*initialize optimization variables*/
-        seps=sqrt(DBL_EPSILON);    //precision of a double
+        seps=sqrt(DBL_EPSILON);    //precision of a IssmDouble
         distance=0.0;              //new_x=old_x + distance
         gold=0.5*(3.0-sqrt(5.0));  //gold = 1 - golden ratio