*** empty log message ***

minimize.cc

@@ -161,6 +161,18 @@ float diff2(const ColumnVector& x, const EvalFunction& func, int i,int j,float h
   return deriv;
 }
 
+ReturnMatrix gradient(const ColumnVector& x, const EvalFunction& func, float h,int errorord){
+  ColumnVector deriv(x.Nrows());
+  deriv = 0;
+
+  for(int i=1;i<=x.Nrows();i++){
+    deriv(i) = diff1(x,func,i,h,errorord);
+  }
+
+  deriv.Release();
+  return deriv;
+}
+
 ReturnMatrix hessian(const ColumnVector& x, const EvalFunction& func, float h,int errorord)
 { //evaluates the hessian of function "eval" at x in parameter space
   

minimize.h

@@ -76,6 +76,8 @@ float diff2(const ColumnVector& x, const EvalFunction& func, int i,float h,int e
 
 float diff2(const ColumnVector& x, const EvalFunction& func, int i,int j,float h,int errorord=4);// finite diff cross derivative
 
+ReturnMatrix gradient(const ColumnVector& x, const EvalFunction& func,float h,int errorord=4);// finite diff derivative vector 
+
 ReturnMatrix hessian(const ColumnVector& x, const EvalFunction& func,float h,int errorord=4);// finite diff hessian