using std::cos;

using boost::math::constants::two_pi;

Real A = 10;

Real y = 10 * v.size();

for (auto x : v) {

y += x * x - A * cos(two_pi<Real>() * x);

}

return y;

std::cout << "Running random search on Rastrigin function (global minimum = (0,0,...,0))\n";

signal(SIGINT, ctrl_c_handler);

using ArgType = std::vector<double>;

auto params = cma_es_parameters<ArgType>();

params.lower_bounds.resize(50, -5.12);

params.upper_bounds.resize(50, 5.12);

std::random_device rd;

std::mt19937_64 gen(rd());

// By definition, the value of the function which a target value is provided must be <= target_value.

double target_value = 1e-3;

std::atomic<double> current_minimum_cost;

std::cout << "Hit ctrl-C to gracefully terminate the optimization." << std::endl;

auto f = [&]() {

return cma_es(rastrigin<double>, params, gen, target_value, &cancel, &current_minimum_cost);

};

auto future = std::async(std::launch::async, f);

std::future_status status = future.wait_for(3ms);

while (!cancel && (status != std::future_status::ready)) {

status = future.wait_for(3ms);

std::cout << "Current cost is " << current_minimum_cost << "\r";

}

auto local_minima = future.get();

std::cout << "Local minimum is {";

for (size_t i = 0; i < local_minima.size() - 1; ++i) {

std::cout << local_minima[i] << ", ";

}

std::cout << local_minima.back() << "}.\n";

std::cout << "Final cost: " << current_minimum_cost << "\n";