/**
* @brief Allocates a problem constructed by stacking two COCO problems.
*
* This is particularly useful for generating multi-objective problems, e.g. a bi-objective problem from two
* single-objective problems. The stacked problem must behave like a normal COCO problem accepting the same
* input. The region of interest in the decision space is defined by parameters smallest_values_of_interest
* and largest_values_of_interest, which are two arrays of size equal to the dimensionality of both problems.
*
* @note Regions of interest in the decision space must either agree or at least one of them must be NULL.
* @note Best parameter becomes somewhat meaningless, but the nadir value make sense now.
*/
static coco_problem_t *coco_problem_stacked_allocate(coco_problem_t *problem1,
coco_problem_t *problem2,
const double *smallest_values_of_interest,
const double *largest_values_of_interest) {
const size_t number_of_variables = coco_problem_get_dimension(problem1);
const size_t number_of_objectives = coco_problem_get_number_of_objectives(problem1)
+ coco_problem_get_number_of_objectives(problem2);
const size_t number_of_constraints = coco_problem_get_number_of_constraints(problem1)
+ coco_problem_get_number_of_constraints(problem2);
size_t i;
char *s;
coco_problem_stacked_data_t *data;
coco_problem_t *problem; /* the new coco problem */
assert(coco_problem_get_dimension(problem1) == coco_problem_get_dimension(problem2));
problem = coco_problem_allocate(number_of_variables, number_of_objectives, number_of_constraints);
s = coco_strconcat(coco_problem_get_id(problem1), "__");
problem->problem_id = coco_strconcat(s, coco_problem_get_id(problem2));
coco_free_memory(s);
s = coco_strconcat(coco_problem_get_name(problem1), " + ");
problem->problem_name = coco_strconcat(s, coco_problem_get_name(problem2));
coco_free_memory(s);
problem->evaluate_function = coco_problem_stacked_evaluate_function;
if (number_of_constraints > 0)
problem->evaluate_constraint = coco_problem_stacked_evaluate_constraint;
assert(smallest_values_of_interest);
assert(largest_values_of_interest);
for (i = 0; i < number_of_variables; ++i) {
problem->smallest_values_of_interest[i] = smallest_values_of_interest[i];
problem->largest_values_of_interest[i] = largest_values_of_interest[i];
}
if (problem->best_parameter) /* logger_bbob doesn't work then anymore */
coco_free_memory(problem->best_parameter);
problem->best_parameter = NULL;
/* Compute the ideal and nadir values */
assert(problem->best_value);
assert(problem->nadir_value);
problem->best_value[0] = problem1->best_value[0];
problem->best_value[1] = problem2->best_value[0];
coco_evaluate_function(problem1, problem2->best_parameter, &problem->nadir_value[0]);
coco_evaluate_function(problem2, problem1->best_parameter, &problem->nadir_value[1]);
/* setup data holder */
data = (coco_problem_stacked_data_t *) coco_allocate_memory(sizeof(*data));
data->problem1 = problem1;
data->problem2 = problem2;
problem->data = data;
problem->problem_free_function = coco_problem_stacked_free;
return problem;
}
/* The gateway function */
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
size_t *ref;
coco_problem_t *pb = NULL;
const char *res;
/* check for proper number of arguments */
if(nrhs!=1) {
mexErrMsgIdAndTxt("cocoProblemGetName:nrhs","One input required.");
}
/* get the problem */
ref = (size_t *)mxGetData(prhs[0]);
pb = (coco_problem_t *)(*ref);
/* call coco_problem_get_name(...) */
res = coco_problem_get_name(pb);
/* prepare the return value */
plhs[0] = mxCreateString(res);
}
/*
* Class: CocoJNI
* Method: cocoProblemGetName
* Signature: (J)Ljava/lang/String;
*/
JNIEXPORT jstring JNICALL Java_org_moeaframework_problem_BBOB2016_CocoJNI_cocoProblemGetName
(JNIEnv *jenv, jclass interface_cls, jlong jproblem_pointer) {
coco_problem_t *problem = NULL;
const char *result;
jstring jresult;
/* This test is both to prevent warning because interface_cls was not used and to check for exceptions */
if (interface_cls == NULL) {
jclass Exception = (*jenv)->FindClass(jenv, "java/lang/Exception");
(*jenv)->ThrowNew(jenv, Exception, "Exception in cocoProblemGetName\n");
}
problem = (coco_problem_t *) jproblem_pointer;
result = coco_problem_get_name(problem);
jresult = (*jenv)->NewStringUTF(jenv, result);
return jresult;
}
/**
* Initializes the toy logger.
*/
static coco_problem_t *logger_toy(coco_observer_t *observer, coco_problem_t *problem) {
logger_toy_t *logger;
coco_problem_t *self;
FILE *output_file;
if (problem->number_of_objectives != 1) {
coco_warning("logger_toy(): The toy logger shouldn't be used to log a problem with %d objectives", problem->number_of_objectives);
}
logger = coco_allocate_memory(sizeof(*logger));
logger->observer = observer;
logger->next_target = 0;
logger->number_of_evaluations = 0;
output_file = ((observer_toy_t *) logger->observer->data)->log_file;
fprintf(output_file, "\n%s, %s\n", coco_problem_get_id(problem), coco_problem_get_name(problem));
self = coco_transformed_allocate(problem, logger, NULL);
self->evaluate_function = logger_toy_evaluate;
return self;
}
/**
* Return a problem that stacks the output of two problems, namely
* of coco_evaluate_function and coco_evaluate_constraint. The accepted
* input remains the same and must be identical for the stacked
* problems.
*
* This is particularly useful to generate multiobjective problems,
* e.g. a biobjective problem from two single objective problems.
*
* Details: regions of interest must either agree or at least one
* of them must be NULL. Best parameter becomes somewhat meaningless.
*/
coco_problem_t *coco_stacked_problem_allocate(coco_problem_t *problem1,
coco_problem_t *problem2,
void *userdata,
coco_stacked_problem_free_data_t free_data) {
const size_t number_of_variables = coco_problem_get_dimension(problem1);
const size_t number_of_objectives = coco_problem_get_number_of_objectives(problem1)
+ coco_problem_get_number_of_objectives(problem2);
const size_t number_of_constraints = coco_problem_get_number_of_constraints(problem1)
+ coco_problem_get_number_of_constraints(problem2);
size_t i;
char *s;
const double *smallest, *largest;
coco_stacked_problem_data_t *data;
coco_problem_t *problem; /* the new coco problem */
assert(coco_problem_get_dimension(problem1) == coco_problem_get_dimension(problem2));
problem = coco_problem_allocate(number_of_variables, number_of_objectives, number_of_constraints);
s = coco_strconcat(coco_problem_get_id(problem1), "__");
problem->problem_id = coco_strconcat(s, coco_problem_get_id(problem2));
coco_free_memory(s);
s = coco_strconcat(coco_problem_get_name(problem1), " + ");
problem->problem_name = coco_strconcat(s, coco_problem_get_name(problem2));
coco_free_memory(s);
problem->evaluate_function = coco_stacked_problem_evaluate;
if (number_of_constraints > 0)
problem->evaluate_constraint = coco_stacked_problem_evaluate_constraint;
/* set/copy "boundaries" and best_parameter */
smallest = problem1->smallest_values_of_interest;
if (smallest == NULL)
smallest = problem2->smallest_values_of_interest;
largest = problem1->largest_values_of_interest;
if (largest == NULL)
largest = problem2->largest_values_of_interest;
for (i = 0; i < number_of_variables; ++i) {
if (problem2->smallest_values_of_interest != NULL)
assert(smallest[i] == problem2->smallest_values_of_interest[i]);
if (problem2->largest_values_of_interest != NULL)
assert(largest[i] == problem2->largest_values_of_interest[i]);
if (smallest != NULL)
problem->smallest_values_of_interest[i] = smallest[i];
if (largest != NULL)
problem->largest_values_of_interest[i] = largest[i];
if (problem->best_parameter) /* bbob2009 logger doesn't work then anymore */
coco_free_memory(problem->best_parameter);
problem->best_parameter = NULL;
if (problem->best_value)
coco_free_memory(problem->best_value);
problem->best_value = NULL; /* bbob2009 logger doesn't work */
}
/* setup data holder */
data = coco_allocate_memory(sizeof(*data));
data->problem1 = problem1;
data->problem2 = problem2;
data->data = userdata;
data->free_data = free_data;
problem->data = data;
problem->free_problem = coco_stacked_problem_free; /* free self->data and coco_problem_free(self) */
return problem;
}
