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;
}
