void rSQPOptimizer::find_optimum() { //------------------------------------------------------------------ // Solve the problem. //------------------------------------------------------------------ namespace mmp = MemMngPack; namespace mp = MoochoPack; using mp::MoochoSolver; Cout << "\nCalling rSQP++ to solve the NLP! ...\n"; if(model_->dv()) { Cout << "\nWarning!!! The dakota model has " << model_->dv() << " discrete variables " << "associated with it that will be ignored by this optimizer!"; } try { MoochoSolver solver; /* nlp_.initialize(true); // Terminate early! Cout << "\nrSQPOptimizer::find_optimum(): Terminate early!\n"; bestVariablesArray.front().continuous_variables(model_->cv()); RealVector dummy_func_values(model_->num_functions()); bestResponseArray.front().function_values(dummy_func_values); return; // Don't do anything! */ solver.set_nlp( mmp::rcp(&nlp_,false) ); std::ofstream Moocho_out("MoochoConsole.out"); solver.set_console_out(mmp::rcp(&Moocho_out,false)); solver.set_error_handling(true,mmp::rcp(&Moocho_out,false)); const MoochoSolver::ESolutionStatus solution_status = solver.solve_nlp(); } catch(const std::exception& excpt) { std::cerr << "\nrSQPOptimizer::find_optimum(): Caught a std::exception: " << excpt.what() << std::endl; } catch(...) { std::cerr << "\nrSQPOptimizer::find_optimum(): Caught an unknown exception: \n"; } //------------------------------------------------------------------ // Save the final results. //------------------------------------------------------------------ bestVariablesArray.front().continuous_variables(nlp_.dakota_x()); // Just use the last point computed bestResponseArray.front().function_values(nlp_.dakota_functions()); // "" }
MoochoPack::MoochoSolver::ESolutionStatus NLPInterfacePack::ExampleNLPFirstOrderRun( const VectorSpace& vec_space ,value_type xo ,bool has_bounds ,bool dep_bounded ,std::ostream* console_out ,std::ostream* error_out ,bool throw_solve_exception ,std::ostream* algo_out ,std::ostream* summary_out ,std::ostream* journal_out ) { using std::endl; using std::setw; namespace rcp = MemMngPack; using Teuchos::RCP; namespace ofsp = OptionsFromStreamPack; using ofsp::OptionsFromStream; namespace rsqp = MoochoPack; using rsqp::MoochoSolver; using rsqp::NLPAlgoConfigMamaJama; MoochoSolver::ESolutionStatus solve_return = MoochoSolver::SOLVE_RETURN_EXCEPTION; int prec = 8; if(console_out) *console_out << std::setprecision(prec) << std::scientific << "*************************************************\n" << "*** Running Tests on ExampleNLPFirstOrder ***\n" << "*************************************************\n" << "\nUsing a vector space of type \'" << typeName(vec_space) << "\'" << "\nwith a dimension of vec_space.dim() = " << vec_space.dim() << std::endl; // Create the nlp ExampleNLPFirstOrder nlp(VectorSpace::space_ptr_t(&vec_space,false),xo,has_bounds,dep_bounded); // Create the solver object and set it up MoochoSolver solver; solver.set_nlp(Teuchos::rcp(&nlp,false)); // Set nlp // set up outputting solver.set_error_handling( throw_solve_exception ,Teuchos::rcp(error_out,false) ); solver.set_console_out(Teuchos::rcp(console_out,false)); solver.set_summary_out(Teuchos::rcp(summary_out,false)); solver.set_journal_out(Teuchos::rcp(journal_out,false)); solver.set_algo_out( Teuchos::rcp(algo_out,false) ); // Run MOOCHO using the MamaJama configuration solve_return = solver.solve_nlp(); return solve_return; }