10 #ifndef ROL_MINIMAX1_HPP
11 #define ROL_MINIMAX1_HPP
15 #include "ROL_Ptr.hpp"
42 Real f1 = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],4.0);
43 Real f2 = std::pow(2.0-(*xp)[0],2.0) + std::pow(2.0-(*xp)[1],2.0);
44 Real f3 = 2.0*std::exp(-(*xp)[0] + (*xp)[1]);
45 return std::max(f1,std::max(f2,f3));
51 Real f1 = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],4.0);
52 Real f2 = std::pow(2.0-(*xp)[0],2.0) + std::pow(2.0-(*xp)[1],2.0);
53 Real f3 = 2.0*std::exp(-(*xp)[0] + (*xp)[1]);
56 if( f1 >= std::max(f2,f3) ) {
57 (*gp)[0] = 2.0*(*xp)[0];
58 (*gp)[1] = 4.0*std::pow((*xp)[1],3.0);
60 else if ( f2 >= std::max(f1,f3) ) {
61 (*gp)[0] = 2.0*(*xp)[0]-4.0;
62 (*gp)[1] = 2.0*(*xp)[1]-4.0;
64 else if ( f3 >= std::max(f1,f2) ) {
65 (*gp)[0] = -2.0*std::exp(-(*xp)[0]+(*xp)[1]);
66 (*gp)[1] = 2.0*std::exp(-(*xp)[0]+(*xp)[1]);
77 return makePtr<Minimax1<Real>>();
81 Ptr<std::vector<Real> > x_ptr = makePtr<std::vector<Real>>(2, 0.0);
82 (*x_ptr)[0] = 1.0; (*x_ptr)[1] = -0.1;
83 return makePtr<StdVector<Real>>(x_ptr);
87 Ptr<std::vector<Real> > z_ptr = makePtr<std::vector<Real>>(2, 0.0);
88 (*z_ptr)[0] = 1.13904; (*z_ptr)[1] = 0.89956;
89 return makePtr<StdVector<Real>>(z_ptr);
Provides the interface to evaluate objective functions.
Ptr< vector > getVector(V &x)
Ptr< Vector< Real > > getSolution(const int i=0) const
Defines the linear algebra or vector space interface.
Contains definitions of test objective functions.
Ptr< Objective< Real > > getObjective(void) const
Ptr< Vector< Real > > getInitialGuess(void) const
Real value(const Vector< Real > &x, Real &tol)
Compute value.
Ptr< const vector > getVector(const V &x)
void gradient(Vector< Real > &g, const Vector< Real > &x, Real &tol)
Compute gradient.
std::vector< Real > vector