RastriginFunction::RastriginFunction(
const arma::uword numberOfDimensions)
: BlackBoxOptimisationBenchmark(numberOfDimensions),
parameterConditioning_(getParameterConditioning(std::sqrt(10.0))) {
if (numberOfDimensions_ < 2) {
throw std::domain_error("RastriginFunction: The number of dimensions must be greater than 1.");
} else if (!isRepresentableAsFloatingPoint(numberOfDimensions_)) {
throw std::overflow_error("RastriginFunction: The number of elements must be representable as a floating point.");
}
setParameterTranslation(getRandomParameterTranslation());
setObjectiveFunctions({{[this](
const arma::vec& parameter_) {
assert(parameter_.n_elem == numberOfDimensions_);
const arma::vec& z = parameterConditioning_ % getAsymmetricParameter(0.2, getOscillatedParameter(parameter_));
return 10.0 * (static_cast<double>(numberOfDimensions_) - arma::accu(arma::cos(2.0 * arma::datum::pi * z))) + std::pow(arma::norm(z), 2.0);
},
"BBOB Rastrigin Function (f3)"}});
}
SchaffersF7FunctionIllConditioned::SchaffersF7FunctionIllConditioned(
const arma::uword numberOfDimensions)
: BlackBoxOptimisationBenchmark(numberOfDimensions),
parameterConditioning_(getParameterConditioning(std::sqrt(1000.0))),
rotationQ_(randomRotationMatrix(numberOfDimensions_)) {
setParameterTranslation(getRandomParameterTranslation());
setParameterRotation(randomRotationMatrix(numberOfDimensions_));
#if defined(SUPPORT_MPI)
MPI_Bcast(rotationQ_.memptr(), static_cast<int>(rotationQ_.n_elem), MPI_DOUBLE, 0, MPI_COMM_WORLD);
#endif
setObjectiveFunction(
[this](
const arma::Col<double>& parameter_) {
assert(parameter_.n_elem == numberOfDimensions_);
const arma::Col<double>& s = arma::square(parameterConditioning_ % (rotationQ_ * getAsymmetricParameter(0.5, parameter_)));
const arma::Col<double>& z = arma::pow(s.head(s.n_elem - 1) + s.tail(s.n_elem - 1), 0.25);
return std::pow(arma::mean(z % (1.0 + arma::square(arma::sin(50.0 * arma::pow(z, 0.4))))), 2.0);
},
"BBOB Schaffers F7 Function, ill-conditioned (f18)");
}
SchaffersF7Function::SchaffersF7Function(
const arma::uword numberOfDimensions)
: BlackBoxOptimisationBenchmark(numberOfDimensions),
parameterConditioning_(getParameterConditioning(std::sqrt(10.0))),
rotationQ_(synchronise(randomRotationMatrix(numberOfDimensions_))) {
if (numberOfDimensions_ < 2) {
throw std::domain_error("SchaffersF7Function: The number of dimensions must be greater than 1.");
}
setParameterTranslation(getRandomParameterTranslation());
setParameterRotation(randomRotationMatrix(numberOfDimensions_));
setObjectiveFunctions({{[this](
const arma::vec& parameter_) {
assert(parameter_.n_elem == numberOfDimensions_);
const arma::vec& s = arma::square(parameterConditioning_ % (rotationQ_ * getAsymmetricParameter(0.5, parameter_)));
const arma::vec& z = arma::pow(s.head(s.n_elem - 1) + s.tail(s.n_elem - 1), 0.25);
return std::pow(arma::mean(z % (1.0 + arma::square(arma::sin(50.0 * arma::pow(z, 0.4))))), 2.0);
},
"BBOB Schaffers F7 Function (f17)"}});
}
