26 ROL::Ptr<CauchyPoint<Real> >
cpt_;
28 ROL::Ptr<Vector<Real> >
s_;
29 ROL::Ptr<Vector<Real> >
Hp_;
37 cpt_ = ROL::makePtr<CauchyPoint<Real>>(parlist);
42 cpt_->initialize(x,s,g);
53 Real tol = std::sqrt(ROL_EPSILON<Real>());
54 const Real
zero(0), half(0.5), one(1), two(2);
60 Real sNnorm =
s_->norm();
61 Real gsN = -
s_->dot(s);
62 bool negCurv = (gsN >
zero ?
true :
false);
66 cpt_->run(s,snorm,iflag,iter,del,model);
83 Real gnorm = s.
norm();
84 Real gnorm2 = gnorm*gnorm;
86 Real gamma = gnorm2/gBg;
87 if ( gamma*gnorm >= del || gBg <=
zero ) {
97 Real a = sNnorm*sNnorm + two*gamma*gsN + gamma*gamma*gnorm2;
98 Real b = -gamma*gsN - gamma*gamma*gnorm2;
99 Real c = gamma*gamma*gnorm2 - del*del;
100 alpha = (-b + sqrt(b*b - a*c))/a;
101 beta = gamma*(one-alpha);
107 pRed_ = (alpha*(half*alpha-one)*gsN - half*beta*beta*gBg + beta*(one-alpha)*gnorm2);
virtual const Vector & dual() const
Return dual representation of , for example, the result of applying a Riesz map, or change of basis...
virtual void scale(const Real alpha)=0
Compute where .
virtual ROL::Ptr< Vector > clone() const =0
Clone to make a new (uninitialized) vector.
virtual void axpy(const Real alpha, const Vector &x)
Compute where .
ROL::Ptr< Vector< Real > > Hp_
virtual void initialize(const Vector< Real > &x, const Vector< Real > &s, const Vector< Real > &g)
Contains definitions of custom data types in ROL.
Provides interface for and implements trust-region subproblem solvers.
Provides the interface to evaluate trust-region model functions.
Defines the linear algebra or vector space interface.
virtual const Ptr< const Vector< Real > > getGradient(void) const
Objective_SerialSimOpt(const Ptr< Obj > &obj, const V &ui) z0_ zero()
virtual void dualTransform(Vector< Real > &tv, const Vector< Real > &v)
void run(Vector< Real > &s, Real &snorm, int &iflag, int &iter, const Real del, TrustRegionModel< Real > &model)
void setPredictedReduction(const Real pRed)
virtual void hessVec(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &s, Real &tol)
Apply Hessian approximation to vector.
virtual void invHessVec(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &s, Real &tol)
Apply inverse Hessian approximation to vector.
ROL::Ptr< CauchyPoint< Real > > cpt_
ROL::Ptr< Vector< Real > > s_
Provides interface for dog leg trust-region subproblem solver.
virtual void set(const Vector &x)
Set where .
virtual Real norm() const =0
Returns where .
DogLeg(ROL::ParameterList &parlist)
void initialize(const Vector< Real > &x, const Vector< Real > &s, const Vector< Real > &g)
virtual void primalTransform(Vector< Real > &tv, const Vector< Real > &v)