44 #ifndef ROL_STDOBJECTIVE_DEF_H
45 #define ROL_STDOBJECTIVE_DEF_H
49 template<
typename Real>
55 template<
typename Real>
61 template<
typename Real>
67 template<
typename Real>
69 const unsigned size = x.size();
70 std::vector<Real> y; y.assign(x.begin(),x.end());
71 const Real cbrteps = std::cbrt(ROL::ROL_EPSILON<Real>()), one(1);
73 const Real val =
value(x,tol);
74 for (
unsigned i = 0; i < size; ++i) {
76 h = cbrteps * std::max(std::abs(xi),one) * sgn(xi);
80 g[i] = (
value(y,tol) - val)/h;
86 template<
typename Real>
93 template<
typename Real>
95 ROL_TEST_FOR_EXCEPTION(
true, std::invalid_argument,
96 ">>> ERROR (ROL::StdObjective): dirDeriv not implemented!");
99 template<
typename Real>
106 catch (std::exception &e) {
111 template<
typename Real>
113 ROL_TEST_FOR_EXCEPTION(
true, std::invalid_argument,
114 ">>> ERROR (ROL::StdObjective): hessVec not implemented!");
117 template<
typename Real>
125 catch (std::exception &e) {
130 template<
typename Real>
132 ROL_TEST_FOR_EXCEPTION(
true, std::invalid_argument,
133 ">>> ERROR (ROL::StdObjective): invHessVec not implemented!");
136 template<
typename Real>
144 template<
typename Real>
146 Pv.assign(v.begin(),v.end());
149 template<
typename Real>
157 template<
typename Real>
159 const Real
zero(0), one(1);
160 return (x <
zero ? -one : one);
virtual Real value(const std::vector< Real > &x, Real &tol)=0
virtual void update(const Vector< Real > &u, const Vector< Real > &z, bool flag=true, int iter=-1) override
virtual Real dirDeriv(const Vector< Real > &x, const Vector< Real > &d, Real &tol)
Compute directional derivative.
virtual void hessVec(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &x, Real &tol)
Apply Hessian approximation to vector.
Ptr< const std::vector< Element > > getVector() const
ROL::Objective_SimOpt value
Defines the linear algebra or vector space interface.
virtual void hessVec(std::vector< Real > &hv, const std::vector< Real > &v, const std::vector< Real > &x, Real &tol)
Objective_SerialSimOpt(const Ptr< Obj > &obj, const V &ui) z0_ zero()
virtual void gradient(std::vector< Real > &g, const std::vector< Real > &x, Real &tol)
virtual void invHessVec(std::vector< Real > &hv, const std::vector< Real > &v, const std::vector< Real > &x, Real &tol)
virtual void update(const std::vector< Real > &x, bool flag=true, int iter=-1)
virtual void precond(std::vector< Real > &Pv, const std::vector< Real > &v, const std::vector< Real > &x, Real &tol)
virtual Real dirDeriv(const std::vector< Real > &x, const std::vector< Real > &d, Real &tol)