44 #ifndef ROL_MINIMAX3_HPP 45 #define ROL_MINIMAX3_HPP 51 #include "Teuchos_RCP.hpp" 64 Teuchos::RCP<const vector>
getVector(
const V& x ) {
65 using Teuchos::dyn_cast;
70 using Teuchos::dyn_cast;
82 Real F = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],2.0) + 2.0*std::pow((*xp)[2],2.0)
83 + std::pow((*xp)[3],2.0) - 5.0*(*xp)[0] - 5.0*(*xp)[1] - 21.0*(*xp)[2] + 7.0*(*xp)[3];
84 Real g2 = -std::pow((*xp)[0],2.0)-std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)-std::pow((*xp)[3],2.0)
85 -(*xp)[0]+(*xp)[1]-(*xp)[2]+(*xp)[3]+8.0;
86 Real g3 = -std::pow((*xp)[0],2.0)-2.0*std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)
87 -2.0*std::pow((*xp)[3],2.0)+(*xp)[0]+(*xp)[3]+10.0;
88 Real g4 = -std::pow((*xp)[0],2.0)-std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)
89 -2.0*(*xp)[0]+(*xp)[1]+(*xp)[3]+5.0;
90 Real a2 = 10.0, a3 = 10.0, a4 = 10.0;
91 return std::max(F,std::max(F-a2*g2,std::max(F-a3*g3,F-a4*g4)));
100 Real F = std::pow((*xp)[0],2.0) + std::pow((*xp)[1],2.0) + 2.0*std::pow((*xp)[2],2.0)
101 + std::pow((*xp)[3],2.0) - 5.0*(*xp)[0] - 5.0*(*xp)[1] - 21.0*(*xp)[2] + 7.0*(*xp)[3];
102 Real g2 = -std::pow((*xp)[0],2.0)-std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)-std::pow((*xp)[3],2.0)
103 -(*xp)[0]+(*xp)[1]-(*xp)[2]+(*xp)[3]+8.0;
104 Real g3 = -std::pow((*xp)[0],2.0)-2.0*std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)
105 -2.0*std::pow((*xp)[3],2.0)+(*xp)[0]+(*xp)[3]+10.0;
106 Real g4 = -std::pow((*xp)[0],2.0)-std::pow((*xp)[1],2.0)-std::pow((*xp)[2],2.0)
107 -2.0*(*xp)[0]+(*xp)[1]+(*xp)[3]+5.0;
108 Real a2 = 10.0, a3 = 10.0, a4 = 10.0;
110 (*gp)[0] = 2.0*(*xp)[0] - 5.0;
111 (*gp)[1] = 2.0*(*xp)[1] - 5.0;
112 (*gp)[2] = 4.0*(*xp)[2] - 21.0;
113 (*gp)[3] = 2.0*(*xp)[3] + 7.0;
114 if ( F-a2*g2 >= std::max(F,std::max(F-a3*g3,F-a4*g4)) ) {
115 (*gp)[0] += a2*(2.0*(*xp)[0] + 1.0);
116 (*gp)[1] += a2*(2.0*(*xp)[1] - 1.0);
117 (*gp)[2] += a2*(2.0*(*xp)[2] + 1.0);
118 (*gp)[3] += a2*(2.0*(*xp)[3] - 1.0);
120 else if ( F-a3*g3 >= std::max(F,std::max(F-a2*g2,F-a4*g4)) ) {
121 (*gp)[0] += a2*(2.0*(*xp)[0] - 1.0);
122 (*gp)[1] += a2*(4.0*(*xp)[1]);
123 (*gp)[2] += a2*(2.0*(*xp)[2]);
124 (*gp)[3] += a2*(4.0*(*xp)[3] - 1.0);
126 else if ( F-a4*g4 >= std::max(F,std::max(F-a2*g2,F-a3*g3)) ) {
127 (*gp)[0] += a2*(2.0*(*xp)[0] + 2.0);
128 (*gp)[1] += a2*(2.0*(*xp)[1] - 1.0);
129 (*gp)[2] += a2*(2.0*(*xp)[2]);
130 (*gp)[3] += a2*(-1.0);
Provides the interface to evaluate objective functions.
std::vector< Real > vector
Defines the linear algebra or vector space interface.
Teuchos::RCP< const vector > getVector(const V &x)
void gradient(Vector< Real > &g, const Vector< Real > &x, Real &tol)
Compute gradient.
Teuchos::RCP< vector > getVector(V &x)
Real value(const Vector< Real > &x, Real &tol)
Compute value.