68 typedef typename std::vector<Real>::size_type
uint;
77 Teuchos::RCP<const std::vector<Real> > ex
82 Real h = 1.0/((Real)(
dim_) + 1.0);
84 f = 2.0*(*ex)[i] + h*h*std::pow((*ex)[i] + (Real)(i+1)*h + 1.0,3.0)/2.0;
85 if ( i < (
dim_-1) ) { f -= (*ex)[i+1]; }
86 if ( i > 0 ) { f -= (*ex)[i-1]; }
93 Teuchos::RCP<std::vector<Real> > eg
95 Teuchos::RCP<const std::vector<Real> > ex
99 Real h = 1.0/((Real)(
dim_) + 1.0);
100 std::vector<Real> f(
dim_,0.0);
103 f[i] = 2.0*(*ex)[i] + h*h*std::pow((*ex)[i] + (Real)(i+1)*h + 1.0,3.0)/2.0;
104 if ( i < (
dim_-1) ) { f[i] -= (*ex)[i+1]; }
105 if ( i > 0) { f[i] -= (*ex)[i-1]; }
109 df = (2.0 + 3.0*h*h*std::pow((*ex)[i] + (Real)(i+1)*h + 1.0,2.0)/2.0)*f[i];
110 if ( i < (
dim_-1) ) { df -= f[i+1]; }
111 if ( i > 0 ) { df -= f[i-1]; }
117 Teuchos::RCP<std::vector<Real> > ehv
119 Teuchos::RCP<const std::vector<Real> > ev
121 Teuchos::RCP<const std::vector<Real> > ex
125 Real h = 1.0/((Real)(
dim_) + 1.0);
126 Real f = 0.0, df = 0.0, dfn = 0.0, hf = 0.0;
128 f = 2.0*(*ex)[i] + h*h*std::pow((*ex)[i] + (Real)(i+1)*h + 1.0,3.0)/2.0;
129 df = 2.0 + 3.0/2.0 * h*h * std::pow((*ex)[i] + (Real)(i+1)*h + 1.0,2.0);
130 hf = 3.0 * h*h * ((*ex)[i] + (Real)(i+1)*h + 1.0);
131 if ( i < (
dim_-2) ) {
132 (*ehv)[i] += 2.0*(*ev)[i+2];
134 if ( i < (
dim_-1) ) {
136 dfn = 2.0 + 3.0/2.0 * h*h * std::pow((*ex)[i+1] + (Real)(i+2)*h + 1.0,2.0);
137 (*ehv)[i] -= 2.0*(df + dfn)*(*ev)[i+1];
138 (*ehv)[i] += 2.0*(*ev)[i];
142 dfn = 2.0 + 3.0/2.0 * h*h * std::pow((*ex)[i-1] + (Real)(i)*h + 1.0,2.0);
143 (*ehv)[i] -= 2.0*(df + dfn)*(*ev)[i-1];
144 (*ehv)[i] += 2.0*(*ev)[i];
147 (*ehv)[i] += 2.0*(*ev)[i-2];
149 (*ehv)[i] += 2.0*(hf*f + df*df)*(*ev)[i];
164 Teuchos::RCP<std::vector<Real> > scale = Teuchos::rcp(
new std::vector<Real>(n,0.0));
187 Teuchos::RCP<std::vector<Real> > x0p = Teuchos::rcp(
new std::vector<Real>(n,0.0));
188 Real h = 1.0/((Real)n + 1.0);
189 for (
int i = 0; i < n; i++ ) {
190 (*x0p)[i] = (Real)(i+1)*h*((Real)(i+1)*h - 1.0);
195 Teuchos::RCP<std::vector<Real> > xp = Teuchos::rcp(
new std::vector<Real>(n,0.0));
196 (*xp)[0] = 1.2321000000000001e-01;
197 (*xp)[1] = 2.1743122909175336e-01;
198 (*xp)[2] = 2.8625218549543746e-01;
199 (*xp)[3] = 3.3309751851140840e-01;
200 (*xp)[4] = 3.6117201714254760e-01;
201 (*xp)[5] = 3.7342787212179440e-01;
202 (*xp)[6] = 3.7255212003706123e-01;
203 (*xp)[7] = 3.6096984201471016e-01;
204 (*xp)[8] = 3.4085861052124522e-01;
205 (*xp)[9] = 3.1417024791439530e-01;
206 (*xp)[10] = 2.8265678244892922e-01;
207 (*xp)[11] = 2.4789833165179542e-01;
208 (*xp)[12] = 2.1133139591375166e-01;
209 (*xp)[13] = 1.7427666644258599e-01;
210 (*xp)[14] = 1.3796594229036069e-01;
211 (*xp)[15] = 1.0356813245768780e-01;
212 (*xp)[16] = 7.2214621084083663e-02;
213 (*xp)[17] = 4.5024529114833199e-02;
214 (*xp)[18] = 2.3130648161534966e-02;
215 (*xp)[19] = 7.7070870882527927e-03;
222 Teuchos::RCP<std::vector<Real> > lp = Teuchos::rcp(
new std::vector<Real>);
223 Teuchos::RCP<std::vector<Real> > up = Teuchos::rcp(
new std::vector<Real>);
224 std::vector<Real> val(n,0.0);
226 val[1] = -0.1*0.4520;
227 val[2] = -0.1*0.6588;
228 val[3] = -0.1*0.8514;
229 val[4] = -0.1*1.0288;
230 val[5] = -0.1*1.1985;
231 val[6] = -0.1*1.3322;
232 val[7] = -0.1*1.4553;
233 val[8] = -0.1*1.5571;
234 val[9] = -0.1*1.6354;
235 val[10] = -0.1*1.6881;
236 val[11] = -0.1*1.7127;
237 val[12] = -0.1*1.7060;
238 val[13] = -0.1*1.6650;
239 val[14] = -0.1*1.5856;
240 val[15] = -0.1*1.4636;
241 val[16] = -0.1*1.2938;
242 val[17] = -0.1*1.0702;
243 val[18] = -0.1*0.7858;
244 val[19] = -0.1*0.4323;
245 for (
int i = 0; i < n; i++ ) {
247 lp->push_back(std::max(-0.2*(Real)(n),val[i]+0.1));
248 up->push_back(std::min( 0.2*(Real)(n),val[i]+1.1));
251 lp->push_back(-0.2*(Real)(n));
252 up->push_back( 0.2*(Real)(n));
255 Teuchos::RCP<Vector<Real> > l = Teuchos::rcp(
new StdVector<Real>(lp) );
256 Teuchos::RCP<Vector<Real> > u = Teuchos::rcp(
new StdVector<Real>(up) );
Provides the interface to evaluate objective functions.
void gradient(Vector< Real > &g, const Vector< Real > &x, Real &tol)
Compute gradient.
Real value(const Vector< Real > &x, Real &tol)
Compute value.
virtual void hessVec(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &x, Real &tol)
Apply Hessian approximation to vector.
virtual void zero()
Set to zero vector.
Defines the linear algebra or vector space interface.
std::vector< Real >::size_type uint
The discrete boundary value problem.
Provides the interface to apply upper and lower bound constraints.
Provides the std::vector implementation of the ROL::Vector interface that handles scalings in the inn...
void getBVP(Teuchos::RCP< Objective< Real > > &obj, Teuchos::RCP< BoundConstraint< Real > > &con, Teuchos::RCP< Vector< Real > > &x0, Teuchos::RCP< Vector< Real > > &x)