xdgv1.cpp

#include"xdgv1.hpp"

using std::abs;

XDGv1::XDGv1() {}

XDGv1::XDGv1(const XDGv1& o) : XDG(o) {}

XDGv1::~XDGv1() {}

string XDGv1::info() {
	return "Extended Differential Grouping version 1 (" + sinfo() + ")";
}

string XDGv1::sinfo() {
	return "XDGv1";
}

tuple<vector<size_t>, vector<vector<size_t>>> XDGv1::run(BoundedObjectiveFunction<double>& ifunc) {
	_epsilon = calc_epsilon(ifunc);
	return XDG::run(ifunc);
}

double XDGv1::calc_epsilon(BoundedObjectiveFunction<double>& fitf) {
	auto x = vector<vector<double>>(np, vector<double>(fitf.dim()));
	auto xf = vector<double>(np);
	std::uniform_real_distribution<double> dist(0.0,1.0);
	for (int i = 0; i < np; i++) for (int j = 0; j < fitf.dim(); j++) x[i][j] = fitf.x_bound_max(j) - fitf.x_bound_min(j) * dist(rand_gen) + fitf.x_bound_min(j);
	for (int i = 0; i < np; i++) xf[i] = fitf(x[i].data());
	auto minf = abs(xf[0]);
	for (int i = 1; i < np; i++) if (minf > abs(xf[i])) minf = abs(xf[i]);
	auto epsilon = alpha * minf;
	return alpha * minf;
}

double XDGv1::epsilon(double y1, double y2, double y3, double y4, size_t dim) {
	return _epsilon;
}