54 int main(
int argc,
char *argv[]) {
56 Teuchos::GlobalMPISession mpiSession(&argc, &argv);
59 int iprint = argc - 1;
60 Teuchos::RCP<std::ostream> outStream;
61 Teuchos::oblackholestream bhs;
63 outStream = Teuchos::rcp(&std::cout,
false);
65 outStream = Teuchos::rcp(&bhs,
false);
79 Teuchos::ParameterList list;
80 list.sublist(
"SimOpt").sublist(
"Solve").set(
"Residual Tolerance",1.e2*
ROL::ROL_EPSILON);
83 Teuchos::RCP<std::vector<RealT> > z_rcp = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
84 Teuchos::RCP<std::vector<RealT> > gz_rcp = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
85 Teuchos::RCP<std::vector<RealT> > yz_rcp = Teuchos::rcp(
new std::vector<RealT> (nx+2, 1.0) );
86 for (
int i=0; i<nx+2; i++) {
93 Teuchos::RCP<ROL::Vector<RealT> > zp = Teuchos::rcp(&z,
false);
94 Teuchos::RCP<ROL::Vector<RealT> > gzp = Teuchos::rcp(&z,
false);
95 Teuchos::RCP<ROL::Vector<RealT> > yzp = Teuchos::rcp(&yz,
false);
97 Teuchos::RCP<std::vector<RealT> > u_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
98 Teuchos::RCP<std::vector<RealT> > gu_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
99 Teuchos::RCP<std::vector<RealT> > yu_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
100 for (
int i=0; i<nx; i++) {
107 Teuchos::RCP<ROL::Vector<RealT> > up = Teuchos::rcp(&u,
false);
108 Teuchos::RCP<ROL::Vector<RealT> > gup = Teuchos::rcp(&u,
false);
109 Teuchos::RCP<ROL::Vector<RealT> > yup = Teuchos::rcp(&yu,
false);
111 Teuchos::RCP<std::vector<RealT> > c_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
112 Teuchos::RCP<std::vector<RealT> > l_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
128 Teuchos::RCP<std::vector<RealT> > p_rcp = Teuchos::rcp(
new std::vector<RealT> (nx, 1.0) );
130 Teuchos::RCP<ROL::Vector<RealT> > pp = Teuchos::rcp(&p,
false);
131 Teuchos::RCP<ROL::Objective_SimOpt<RealT> > pobj = Teuchos::rcp(&obj,
false);
132 Teuchos::RCP<ROL::EqualityConstraint_SimOpt<RealT> > pcon = Teuchos::rcp(&con,
false);
139 std::string filename =
"input.xml";
140 Teuchos::RCP<Teuchos::ParameterList> parlist = Teuchos::rcp(
new Teuchos::ParameterList() );
141 Teuchos::updateParametersFromXmlFile( filename, parlist.ptr() );
142 parlist->sublist(
"Status Test").set(
"Gradient Tolerance",1.e-14);
143 parlist->sublist(
"Status Test").set(
"Constraint Tolerance",1.e-14);
144 parlist->sublist(
"Status Test").set(
"Step Tolerance",1.e-16);
145 parlist->sublist(
"Status Test").set(
"Iteration Limit",1000);
147 Teuchos::RCP<ROL::Algorithm<RealT> > algo;
151 RealT zerotol = std::sqrt(ROL::ROL_EPSILON);
153 con.
solve(c,u,z,zerotol);
155 algo->run(x, g, l, c, obj, con,
true, *outStream);
156 Teuchos::RCP<ROL::Vector<RealT> > zCS = z.
clone();
162 algo->run(z,robj,
true,*outStream);
165 Teuchos::RCP<ROL::Vector<RealT> > err = z.
clone();
166 err->set(*zCS); err->axpy(-1.,z);
167 errorFlag += ((err->norm()) > 1.e-8) ? 1 : 0;
169 catch (std::logic_error err) {
170 *outStream << err.what() <<
"\n";
175 std::cout <<
"End Result: TEST FAILED\n";
177 std::cout <<
"End Result: TEST PASSED\n";
virtual void setSolveParameters(Teuchos::ParameterList &parlist)
Set solve parameters.
Defines the linear algebra or vector space interface for simulation-based optimization.
virtual std::vector< std::vector< Real > > checkApplyAdjointJacobian(const Vector< Real > &x, const Vector< Real > &v, const Vector< Real > &c, const Vector< Real > &ajv, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS)
Finite-difference check for the application of the adjoint of constraint Jacobian.
void solve(ROL::Vector< Real > &c, ROL::Vector< Real > &u, const ROL::Vector< Real > &z, Real &tol)
Given , solve for .
virtual void zero()
Set to zero vector.
virtual std::vector< std::vector< Real > > checkGradient(const Vector< Real > &x, const Vector< Real > &d, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
Finite-difference gradient check.
virtual std::vector< std::vector< Real > > checkApplyJacobian(const Vector< Real > &x, const Vector< Real > &v, const Vector< Real > &jv, const std::vector< Real > &steps, const bool printToStream=true, std::ostream &outStream=std::cout, const int order=1)
Finite-difference check for the constraint Jacobian application.
Provides the interface to evaluate simulation-based reduced objective functions.
Provides the std::vector implementation of the ROL::Vector interface.
Provides an interface to run optimization algorithms.
virtual Teuchos::RCP< Vector< Real > > clone() const
Clone to make a new (uninitialized) vector.
int main(int argc, char *argv[])
virtual std::vector< std::vector< Real > > checkHessVec(const Vector< Real > &x, const Vector< Real > &v, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
Finite-difference Hessian-applied-to-vector check.
virtual std::vector< std::vector< Real > > checkApplyAdjointHessian(const Vector< Real > &x, const Vector< Real > &u, const Vector< Real > &v, const Vector< Real > &hv, const std::vector< Real > &step, const bool printToScreen=true, std::ostream &outStream=std::cout, const int order=1)
Finite-difference check for the application of the adjoint of constraint Hessian. ...
static const double ROL_EPSILON
Platform-dependent machine epsilon.