| 25 |
|
|
| 26 |
|
gplus = calcGrad(tempX); |
| 27 |
|
|
| 28 |
< |
hessian.Cloume(i) = (gplus - currentGrad) / hi; |
| 28 |
> |
//hessian.Colume(i) = (gplus - currentGrad) / hi; |
| 29 |
|
|
| 30 |
|
//restore tempX to its original value |
| 31 |
|
tempX[i] -= hi; |
| 66 |
|
double ConcreteNLMode1::calcF(){ |
| 67 |
|
|
| 68 |
|
currentF = (*objfunc)(currentX, currentGrad); |
| 69 |
< |
numOfProc++; |
| 69 |
> |
numOfFunEval ++; |
| 70 |
|
|
| 71 |
< |
return curretF; |
| 71 |
> |
return currentF; |
| 72 |
|
} |
| 73 |
|
|
| 74 |
< |
double ConcreteNLMode1::calcF(const vector<double>& x){ |
| 74 |
> |
double ConcreteNLMode1::calcF(vector<double>& x){ |
| 75 |
|
|
| 76 |
|
vector<double> tempGrad(x.size()); |
| 77 |
|
|
| 84 |
|
} |
| 85 |
|
|
| 86 |
|
vector<double> ConcreteNLMode1::calcGrad(){ |
| 87 |
< |
|
| 87 |
> |
|
| 88 |
|
currentF = (*objfunc)(currentX, currentGrad); |
| 89 |
|
|
| 90 |
|
return currentGrad; |
| 92 |
|
} |
| 93 |
|
|
| 94 |
|
vector<double> ConcreteNLMode1::calcGrad(vector<double>& x){ |
| 95 |
+ |
vector<double> tempGrad(x.size()); |
| 96 |
|
|
| 96 |
– |
vector<double> tempGrad(x.szie()); |
| 97 |
– |
|
| 97 |
|
double tempF; |
| 98 |
|
|
| 99 |
|
tempF = (*objfunc)(x, tempGrad); |
| 101 |
|
return tempGrad; |
| 102 |
|
} |
| 103 |
|
|
| 104 |
+ |
/* |
| 105 |
|
SymMatrix ConcreteNLMode1::calcHessian(){ |
| 106 |
|
calcGrad(currentX); |
| 107 |
|
return FiniteHessian(currentX); |
| 110 |
|
SymMatrix ConcreteNLMode1::calcHessian(vector<double>& x){ |
| 111 |
|
return FiniteHessian(x); |
| 112 |
|
} |
| 113 |
+ |
*/ |