| 89 |
|
|
| 90 |
|
// conjugate gradient search direction |
| 91 |
|
direction = getUpdatedDirection(P, gold2, prevGradient); |
| 92 |
– |
|
| 92 |
|
sddiff = direction - lineSearch_->searchDirection(); |
| 93 |
|
lineSearch_->searchDirection() = direction; |
| 94 |
|
// Now compute accuracy and check end criteria |
| 96 |
|
fnew = P.functionValue(); |
| 97 |
|
fdiff = 2.0*std::fabs(fnew-fold) / |
| 98 |
|
(std::fabs(fnew) + std::fabs(fold) + NumericConstant::epsilon); |
| 100 |
– |
std::cerr << "fdiff = " << fdiff << "ftol = " << ftol << "\n"; |
| 99 |
|
if (fdiff < ftol || |
| 100 |
|
endCriteria.checkMaxIterations(iterationNumber_, ecType)) { |
| 101 |
|
endCriteria.checkStationaryFunctionValue(0.0, 0.0, |
| 105 |
|
} |
| 106 |
|
P.setCurrentValue(x_); // update problem current value |
| 107 |
|
++iterationNumber_; // Increase iteration number |
| 110 |
– |
std::cerr << "in = " << iterationNumber_ << "\n"; |
| 108 |
|
first_time = false; |
| 109 |
|
} else { |
| 110 |
|
done = true; |