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(
"General").set(
"Inexact Hessian-Times-A-Vector",
true);
84 parlist->sublist(
"General").set(
"Inexact Hessian-Times-A-Vector",
false);
86 parlist->sublist(
"Step").set(
"Type",
"Line Search");
91 ROL::Ptr<ROL::Vector<RealT> > x0;
92 std::vector<ROL::Ptr<ROL::Vector<RealT> > > z;
93 ROL::Ptr<ROL::OptimizationProblem<RealT> > problem;
94 ROL::GetTestProblem<RealT>(problem,x0,z,objFunc);
102 int dim = x0->dimension();
103 parlist->sublist(
"General").sublist(
"Krylov").set(
"Iteration Limit", 2*
dim);
106 ROL::Ptr<ROL::Vector<RealT> > e = x0->clone();
109 parlist->sublist(
"Step").sublist(
"Line Search").sublist(
"Descent Method").set(
"Type",
ROL::EDescentToString(desc));
122 solver.
solve(*outStream);
126 for (
int i = 0; i < static_cast<int>(z.size()); ++i) {
133 err = std::min(err,e->norm());
136 *outStream << std::endl <<
"Norm of Error: " << err << std::endl;
142 catch (std::logic_error& err) {
143 *outStream << err.what() << std::endl;
148 std::cout <<
"End Result: TEST FAILED" << std::endl;
150 std::cout <<
"End Result: TEST PASSED" << std::endl;