60int main(
int argc,
char *argv[]) {
62 Teuchos::GlobalMPISession mpiSession(&argc, &argv);
65 int iprint = argc - 1;
66 ROL::Ptr<std::ostream> outStream;
69 outStream = ROL::makePtrFromRef(std::cout);
71 outStream = ROL::makePtrFromRef(bhs);
79 std::string filename =
"input.xml";
81 auto parlist = ROL::getParametersFromXmlFile( filename );
82 parlist->sublist(
"Step").set(
"Type",
"Trust Region");
89 ROL::Ptr<ROL::Vector<RealT> > x0;
90 std::vector<ROL::Ptr<ROL::Vector<RealT> > > z;
91 ROL::Ptr<ROL::OptimizationProblem<RealT> > problem;
92 ROL::GetTestProblem<RealT>(problem,x0,z,objFunc);
101 int dim = x0->dimension();
102 parlist->sublist(
"General").sublist(
"Krylov").set(
"Iteration Limit", 5*
dim);
105 ROL::Ptr<ROL::Vector<RealT> > e = x0->clone();
109 parlist->sublist(
"Step").sublist(
"Trust Region").set(
"Subproblem Solver", ETrustRegionToString(tr));
115 solver.
solve(*outStream);
119 for (
int i = 0; i < static_cast<int>(z.size()); ++i) {
126 err = std::min(err,e->norm());
129 *outStream << std::endl <<
"Norm of Error: " << err << std::endl;
135 catch (std::logic_error& err) {
136 *outStream << err.what() << std::endl;
141 std::cout <<
"End Result: TEST FAILED" << std::endl;
143 std::cout <<
"End Result: TEST PASSED" << std::endl;