F6.cpp

#include "F6.h"
#include <stdio.h>

/**
 * Single-group Shifted and m-rotated Ackley¡¯s Function
 *
 * as defined in "Benchmark Functions for the CEC'2010 Special Session
 * and Competition on Large-Scale Global Optimization" by Ke Tang,
 * Xiaodong Li, P. N. Suganthan, Zhenyu Yang, and Thomas Weise
 * published as technical report on January 8, 2010 at Nature Inspired
 * Computation and Applications Laboratory (NICAL), School of Computer
 * Science and Technology, University of Science and Technology of China,
 * Hefei, Anhui, China.
 */


F6::F6():Benchmarks(){
	m_havenextGaussian=0;
	Ovector = NULL;
	minX = -32;
	maxX = 32;
	ID = 6;
}

F6::~F6(){
	delete[] Ovector;
	delete[] Pvector;
	delete[] RotMatrix;
}

double F6::compute(double*x){
	int    m = nonSeparableGroupSize;
	int    i;
	double result;

	if(Ovector == NULL) {
		Ovector   = createShiftVector(dimension,minX,maxX);
		Pvector   = createPermVector(dimension);
		RotMatrix = createRotMatrix1D(nonSeparableGroupSize);
	}

	for(i = 0; i < dimension; i++) {
		anotherz[i] = x[i] - Ovector[i];
	}

	for(i = 0; i < m; i++) {
		anotherz1[i] = anotherz[Pvector[i]];
	}

	for(i = m; i < dimension; i++) {
		anotherz2[i - m] = anotherz[Pvector[i]];
	}

	result = rot_ackley(anotherz1,m) * 1e6 + ackley(anotherz2,dimension - m);
	return(result);
}

double F6::compute(vector<double> x){
	int    m = nonSeparableGroupSize;
	int    i;
	double result;

	if(Ovector == NULL) {
		Ovector   = createShiftVector(dimension,minX,maxX);
		Pvector   = createPermVector(dimension);
		RotMatrix = createRotMatrix1D(nonSeparableGroupSize);
	}

	for(i = 0; i < dimension; i++) {
		anotherz[i] = x[i] - Ovector[i];
	}

	for(i = 0; i < m; i++) {
		anotherz1[i] = anotherz[Pvector[i]];
	}

	for(i = m; i < dimension; i++) {
		anotherz2[i - m] = anotherz[Pvector[i]];
	}

	result = rot_ackley(anotherz1,m) * 1e6 + ackley(anotherz2,dimension - m);
	return(result);
}