/*------------------------------------------*/
int main ( int argc , char ** argv )
{
// display:
NOMAD::Display out ( NOMAD::rout ); //zhenghua
out.precision ( NOMAD::DISPLAY_PRECISION_STD );
std::string error;
{
// NOMAD initializations:
NOMAD::begin ( argc , argv );
// usage:
if ( argc < 2 ) {
NOMAD::display_usage ( out); //zhenghua
NOMAD::end();
return EXIT_FAILURE;
}
// parameters file:
std::string param_file_name = argv[1];
std::string opt = param_file_name;
NOMAD::toupper ( opt );
// display version if option '-v' has been specified:
if ( opt == "-U" ) {
NOMAD::display_usage ( out );
NOMAD::end();
return EXIT_SUCCESS;
}
// display version if option '-v' has been specified:
if ( opt == "-V" || opt =="-VERSION") {
NOMAD::display_version ( out );
NOMAD::end();
return EXIT_SUCCESS;
}
// display info if option '-i' has been specified:
if ( opt == "-I" || opt == "-INFO" ) {
NOMAD::display_info ( out );
NOMAD::display_usage ( out ); //zhenghua
NOMAD::end();
return EXIT_SUCCESS;
}
// parameters creation:
NOMAD::Parameters p ( out );
// display help on parameters if option '-h' has been specified:
if ( opt == "-H" || opt == "-HELP" ) {
p.help ( argc , argv );
NOMAD::end();
return EXIT_SUCCESS;
}
// display developer help on parameters if option '-d' has been specified:
if ( opt == "-D" ) {
p.help ( argc , argv,true );
NOMAD::end();
return EXIT_SUCCESS;
}
// check the number of processess:
#ifdef USE_MPI
if ( NOMAD::Slave::get_nb_processes() < 2 ) {
NOMAD::rout << "ERROR: Incorrect command to run with MPI." << std::endl; //zhenghua
NOMAD::display_usage ( out ); //zhenghua
NOMAD::end();
return EXIT_FAILURE;
}
#endif
try {
// read parameters file:
p.read ( param_file_name );
// parameters check:
p.check();
// display NOMAD info:
if ( p.get_display_degree() > NOMAD::MINIMAL_DISPLAY)
NOMAD::display_info ( out );
// parameters display:
if ( NOMAD::Slave::is_master() &&
p.get_display_degree() == NOMAD::FULL_DISPLAY )
out << std::endl
<< NOMAD::open_block ( "parameters" ) << std::endl
<< p
<< NOMAD::close_block();
// algorithm creation and execution:
NOMAD::Mads mads ( p , NULL );
if ( p.get_nb_obj() == 1 )
mads.run();
else
mads.multi_run();
#ifdef MODEL_STATS
mads.display_model_stats ( out );
#endif
}
catch ( std::exception & e ) {
if ( NOMAD::Slave::is_master() ) {
error = std::string ( "NOMAD has been interrupted: " ) + e.what();
NOMAD::rout << std::endl << error << std::endl << std::endl; //zhenghua
}
}
NOMAD::Slave::stop_slaves ( out );
NOMAD::end();
}
#ifdef MEMORY_DEBUG
NOMAD::display_cardinalities ( out );
#endif
return ( error.empty() ) ? EXIT_SUCCESS : EXIT_FAILURE;
}
/**
- Basic version.
\param p Parameters -- \b IN.
\param ev A pointer to the evaluator -- \b IN
-- \b optional (default = \c NULL).
*/
Mads ( NOMAD::Parameters & p , NOMAD::Evaluator * ev = NULL )
: _p ( p ) ,
_stats ( p.get_sgte_cost() ) ,
_ev_control ( p , _stats , ev , NULL , NULL ) ,
_ev_control_for_sorting( p , _stats , _ev_control.get_evaluator() , &(_ev_control.get_cache()) , &(_ev_control.get_sgte_cache()) ) ,
_true_barrier ( p , NOMAD::TRUTH ) ,
_sgte_barrier ( p , NOMAD::SGTE ) ,
_mesh ( p.get_signature()->get_mesh() ) ,
_pareto_front ( NULL ) ,
_user_search ( NULL ) ,
_model_search1 ( NULL ) ,
_model_search2 ( NULL ) ,
_VNS_search ( NULL ) ,
_cache_search ( NULL ) ,
_L_curve ( NULL ) ,
_extended_poll ( NULL ) ,
_user_ext_poll ( false ) { init(); }
//Constructor
matlabEval(const NOMAD::Parameters &p, usrFcn *_obj, int _nobj, usrFcn *_con, int _ncon, iter_fun_data *_iterF) : NOMAD::Evaluator(p)
{
fun = _obj;
hasSur = p.has_sgte();
nobj = _nobj;
con = _con;
ncon = _ncon;
iterF = _iterF;
citer = 0;
}
/**
- Advanced version.
\param p Parameters -- \b IN.
\param ev A pointer to the evaluator -- \b IN (may be \c NULL).
\param extended_poll A pointer to a NOMAD::Extended_Poll object
-- \b IN (may be \c NULL).
\param cache A pointer to a cache -- \b IN (may be \c NULL).
\param sgte_cache A pointer to a cache for surrogates
-- \b IN (may be \c NULL).
*/
Mads ( NOMAD::Parameters & p ,
NOMAD::Evaluator * ev , // may be NULL
NOMAD::Extended_Poll * extended_poll , // may be NULL
NOMAD::Cache * cache , // may be NULL
NOMAD::Cache * sgte_cache ) // may be NULL
: _p ( p ) ,
_stats ( p.get_sgte_cost() ) ,
_ev_control ( p , _stats , ev , cache , sgte_cache ) ,
_ev_control_for_sorting( p , _stats , _ev_control.get_evaluator() , cache , sgte_cache ) ,
_true_barrier ( p , NOMAD::TRUTH ) ,
_sgte_barrier ( p , NOMAD::SGTE ) ,
_mesh ( p.get_signature()->get_mesh() ) ,
_pareto_front ( NULL ) ,
_user_search ( NULL ) ,
_model_search1 ( NULL ) ,
_model_search2 ( NULL ) ,
_VNS_search ( NULL ) ,
_cache_search ( NULL ) ,
_L_curve ( NULL ) ,
_extended_poll ( extended_poll ) ,
_user_ext_poll ( (extended_poll!=NULL) ) { init(); }
/**
\param p Parameters -- \b IN.
\param et Truth or surrogate -- \b IN.
*/
Barrier(const NOMAD::Parameters &p , NOMAD::eval_type et)
: _p(p) ,
_eval_type(et) ,
_h_max(p.get_h_max_0()) ,
_best_feasible(NULL) ,
_ref(NULL) ,
_rho_leaps(0) ,
_poll_center(NULL) ,
_sec_poll_center(NULL) ,
_peb_changes(0) ,
_peb_filter_reset(0) ,
_one_eval_succ(NOMAD::UNSUCCESSFUL) ,
_success(NOMAD::UNSUCCESSFUL) {}
//
// Set options specified in R object opts in IpoptApplcation app.
// If we set the option to approximate the Hessian, then IpoptRNLP
// needs to know that, so it can return false when calling eval_h,
// instead of trying to find an R function that evaluates the Hessian.
//
// SEXP opts is an R list containing three sub-lists, with names
// integer, string and numeric. These sub-lists contain the actual
// options and there values that were specified by the user.
// Passing the options in this way makes it easier to call different
// SetValue function in IpoptApplication of the different types.
//
void setApplicationOptions(NOMAD::Parameters & p, SEXP opts ) {
std::stringbuf nomadbuf;
std::iostream fp(&nomadbuf);
fp.precision(15);
fp.seekg(0, std::ios::beg);
// extract the sub-lists with options of the different types into separate lists
SEXP opts_integer = getListElement(opts, "integer");
SEXP opts_numeric = getListElement(opts, "numeric");
SEXP opts_string = getListElement(opts, "string");
// showArgs1(opts_integer);
// showArgs1(opts_numeric);
// showArgs1(opts_string);
// loop over the integer options and set them
SEXP opts_integer_names;
opts_integer_names = getAttrib(opts_integer, R_NamesSymbol);
for (int list_cnt=0;list_cnt<length( opts_integer );list_cnt++) {
SEXP opt_value;
PROTECT(opt_value = AS_INTEGER(VECTOR_ELT(opts_integer, list_cnt)));
fp << CHAR(STRING_ELT(opts_integer_names, list_cnt)) << "\t" << INTEGER(opt_value)[0] << std::endl;
UNPROTECT(1);
}
// loop over the numeric options and set them
SEXP opts_numeric_names;
opts_numeric_names = getAttrib(opts_numeric, R_NamesSymbol);
for (int list_cnt=0;list_cnt<length( opts_numeric );list_cnt++) {
SEXP opt_value;
PROTECT(opt_value = VECTOR_ELT(opts_numeric, list_cnt));
fp << CHAR(STRING_ELT(opts_numeric_names, list_cnt)) << "\t" << REAL(opt_value)[0] << std::endl;
UNPROTECT(1);
}
// loop over the string options and set them
SEXP opts_string_names;
opts_string_names = getAttrib(opts_string, R_NamesSymbol);
for (int list_cnt=0;list_cnt<length( opts_string );list_cnt++) {
// opt_value will contain the first (should be the only one) element of the list
SEXP opt_value;
PROTECT(opt_value = STRING_ELT(VECTOR_ELT(opts_string, list_cnt),0));
fp << CHAR(STRING_ELT(opts_string_names, list_cnt)) << "\t" << CHAR(opt_value) << std::endl;
UNPROTECT(1);
}
fp.seekg (0, std::ios::beg);
p.read(fp);
return;
}
/*-----------------------------------*/
int main ( int argc , char ** argv ) {
// MPI initialization:
MPI_Init ( &argc , &argv );
int rank , np;
MPI_Comm_rank ( MPI_COMM_WORLD, &rank );
MPI_Comm_size ( MPI_COMM_WORLD, &np );
// check the arguments and the number of processes:
if ( np <= 2 || argc != 4 ) {
if ( rank == 0 )
cerr << "usage: mpirun -np p " << argv[0]
<< " param_file bbe ns, with p>2,"
<< " bbe>1, and 1<=ns<=n."
<< endl;
MPI_Finalize();
return 1;
}
// display:
Display out ( cout );
out.precision ( 16 );
// parameters:
NOMAD::Parameters p ( out );
int bbe = atoi ( argv[2] );
int ns = atoi ( argv[3] );
try {
// read the parameters file:
p.read ( argv[1] );
// check the parameters:
p.check();
if ( ns < 1 || ns > p.get_dimension() )
throw Exception ( __FILE__ , __LINE__ ,
"Bad value for ns the number of free variables for each process" );
if ( p.get_nb_obj() > 1 )
throw Exception ( __FILE__ , __LINE__ ,
"PSD-MADS is not designed for multi-objective optimization" );
}
catch ( exception & e ) {
if ( rank == 0 )
cerr << "error with parameters" << endl;
MPI_Finalize();
return 1;
}
// start the master:
Master_Slaves master_slaves ( rank , np , bbe , ns , p , DEBUG );
master_slaves.start();
// cache server:
Cache_Server cache ( out ,
rank ,
np ,
p.get_h_min() ,
p.get_max_bb_eval() ,
false , // ALLOW_MULTIPLE_EVALS
DEBUG );
// start the cache server:
if ( rank == np-1 ) {
if ( !DEBUG )
out << endl << "TIME\tBBE\tOBJ" << endl << endl;
cache.start();
}
// slaves: algorithm creation and execution:
if ( rank != 0 && rank != np-1 ) {
// MADS run:
master_slaves.mads_run ( cache );
// stop the master:
master_slaves.stop();
}
// stop the cache server:
cache.stop();
// display the final solution:
if ( !DEBUG && rank == np-1 )
cache.display_best_points ( out );
// MPI finalization:
MPI_Finalize();
return 0;
}