nojhan/paradiseo
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

GPU -> GPU-good

Browse source 
git-svn-id: svn://scm.gforge.inria.fr/svnroot/paradiseo@2177 331e1502-861f-0410-8da2-ba01fb791d7f
This commit is contained in:
liefooga 2011-03-11 10:48:32 +00:00
commit 842e7f83ad
158 changed files with 0 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

107
ParadisEO-GPU-good/src/cudaType/moCudaBitVector.h Normal file
View file

@ -0,0 +1,107 @@
/*
<moCudaBitVector.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaBitVector_H_
#define _moCudaBitVector_H_
#include <cudaType/moCudaVector.h>
/**
* Implementation of Bit vector representation on CUDA.
*/
template<class Fitness>
class moCudaBitVector: public moCudaVector<bool, Fitness> {
public:
using moCudaVector<bool, Fitness>::vect;
using moCudaVector<bool, Fitness>::N;
/**
* Default constructor.
*/
moCudaBitVector() :
moCudaVector<bool, Fitness> () {
}
/**
*Constructor.
*@param _size The neighborhood size.
*/
moCudaBitVector(unsigned _size) {
N = _size;
vect = new bool[_size];
create();
}
/**
*Constructor.
*@param _size The neighborhood size.
*@param _b Value to assign to vector.
*/
moCudaBitVector(unsigned _size, bool _b) {
N = _size;
vect = new bool[_size];
for (unsigned i = 0; i < _size; i++)
vect[i] = _b;
}
/**
*Initializer of random bit vector.
*/
void create() {
for (unsigned i = 0; i < N; i++) {
vect[i] = (bool) round((float) rng.rand() / RAND_MAX);
}
}
};
#endif

113
ParadisEO-GPU-good/src/cudaType/moCudaIntVector.h Normal file
View file

@ -0,0 +1,113 @@
/*
<moCudaIntVector.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaIntVector_H_
#define _moCudaIntVector_H_
#include <cudaType/moCudaVector.h>
/**
* Implementation of integer vector representation on CUDA.
*/
template<class Fitness>
class moCudaIntVector: public moCudaVector<int, Fitness> {
public:
using moCudaVector<int, Fitness>::vect;
using moCudaVector<int, Fitness>::N;
/**
* Default constructor.
*/
moCudaIntVector() :
moCudaVector<int, Fitness> () {
}
/**
*Constructor.
*@param _size The neighborhood size.
*/
moCudaIntVector(unsigned _size) {
N = _size;
vect = new int[_size];
create();
}
/**
*Assignment operator
*@param _vector The vector passed to the function determine the new content.
*@return a new vector.
*/
moCudaIntVector& operator=(const moCudaIntVector & _vector) {
N = _vector.N;
vect = new int[N];
for (unsigned i = 0; i < N; i++)
vect[i] = _vector.vect[i];
fitness(_vector.fitness());
return (*this);
}
/**
*Initializer of random integer vector.
*/
void create() {
unsigned random;
int temp;
for (unsigned i = 0; i < N; i++)
vect[i] = i;
// we want a random permutation so we shuffle
for (unsigned i = 0; i < N; i++) {
random = rng.rand() % (N - i) + i;
temp = vect[i];
vect[i] = vect[random];
vect[random] = temp;
}
}
};
#endif

159
ParadisEO-GPU-good/src/cudaType/moCudaVector.h Normal file
View file

@ -0,0 +1,159 @@
/*
<moCudaVector.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Thé Van LUONG, Karima Boufaras
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaVector_H_
#define _moCudaVector_H_
#include <eo>
/**
* Implementation of vector representation on CUDA.
*/
template<class ElemT, class Fitness>
class moCudaVector: public EO<Fitness> {
public:
typedef ElemT ElemType;
ElemType * vect;
unsigned N;
/**
* Default constructor.
*/
moCudaVector() :
N(0) {
}
/**
*Constructor.
*@param _size The neighborhood size.
*/
moCudaVector(unsigned _size) {
N = _size;
vect = new ElemType[_size];
create();
}
/**
* Destructor.
*/
~moCudaVector() {
if (N > 1)
delete[] vect;
}
/**
*How to fill the vector.
*/
virtual void create() {
}
/**
*Assignment operator
*@param _vector The vector passed to the function determine the new content.
*@return a new vector.
*/
virtual moCudaVector& operator=(const moCudaVector & _vector) {
N = _vector.N;
vect = new ElemType[N];
for (unsigned i = 0; i < N; i++)
vect[i] = _vector.vect[i];
fitness(_vector.fitness());
return (*this);
}
/**
*An accessor read only on the i'th element of the vector (function inline can be called from host or device).
*@param
*_i The i'th element of vector.
*@return
*The i'th element of the vector for read only
*/
inline __host__ __device__ const ElemType & operator[](unsigned _i) const {
return vect[_i];
}
/**
*An accessor read-write on the i'th element of the vector(function inline can be called from host or device).
*@param _i The i'th element of the vector.
*@return The i'th element of the vector for read-write
*/
inline __host__ __device__ ElemType & operator[](unsigned _i) {
return vect[_i];
}
/**
*Function inline to get the size of vector, called from host and device.
*@return The vector size's
*/
inline __host__ __device__ unsigned size() {
return N;
}
virtual void printOn(std::ostream& os) const {
EO<Fitness>::printOn(os);
os << ' ';
os << N << ' ';
unsigned int i;
for (i = 0; i < N; i++)
os << vect[i] << ' ';
}
};
#endif

136
ParadisEO-GPU-good/src/eval/moCudaEval.h Normal file
View file

@ -0,0 +1,136 @@
/*
<moCudaEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef moCudaEval_H
#define moCudaEval_H
#include <eval/moEval.h>
/**
* Abstract class for evaluation on GPU
*/
template<class Neighbor>
class moCudaEval: public moEval<Neighbor> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
typedef typename EOT::Fitness Fitness;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
*/
moCudaEval(unsigned int _neighborhoodSize) {
neighborhoodSize = _neighborhoodSize;
host_FitnessArray = new Fitness[neighborhoodSize];
cudaMalloc((void**) &device_FitnessArray, neighborhoodSize
* sizeof(Fitness));
kernel_Dim = neighborhoodSize / BLOCK_SIZE + ((neighborhoodSize
% BLOCK_SIZE == 0) ? 0 : 1);
}
/**
* Destructor
*/
~moCudaEval() {
delete[] host_FitnessArray;
cudaFree(device_FitnessArray);
cudaFree(&device_solution);
}
/**
* Set fitness of a solution neighbors
*@param _sol the solution which generate the neighborhood
*@param _neighbor the current neighbor
*/
void operator()(EOT & _sol, Neighbor & _neighbor) {
_neighbor.fitness(host_FitnessArray[_neighbor.index()]);
}
/**
* Compute fitness for all Kswap solution neighbors in device
* @param _sol the solution which generate the neighborhood
* @param _mapping the neighborhood mapping
* @param _Kswap the number of swap
*/
virtual void neighborhoodKswapEval(EOT & _sol, unsigned * _mapping,
unsigned _Kswap) {
}
/**
* Compute fitness for all Kflip solution neighbors in device
* @param _sol the solution which generate the neighborhood
* @param _mapping the neighborhood mapping
* @param _Kflip the number of bit to flip
*/
virtual void neighborhoodKflipEval(EOT & _sol, unsigned * _mapping,
unsigned _Kflip) {
}
/**
* Compute fitness for all solution neighbors in device
* @param _sol the solution which generate the neighborhood
*/
virtual void neighborhoodEval(EOT & _sol)=0;
protected:
//the host array to save all neighbors fitness
Fitness * host_FitnessArray;
//the device array to save neighbors fitness computed in device
Fitness * device_FitnessArray;
//the device solution
EOT device_solution;
//the size of neighborhood
unsigned int neighborhoodSize;
//Cuda kernel dimension
int kernel_Dim;
};
#endif

84
ParadisEO-GPU-good/src/eval/moCudaEvalFunc.h Normal file
View file

@ -0,0 +1,84 @@
/*
<moCudaEvalFunc.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Thé Van LUONG, Karima Boufaras
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __moCudaEvalFunc_H
#define __moCudaEvalFunc_H
/**
* Abstract class for CUDA evaluation of neighbor
*/
template<class Neighbor>
class moCudaEvalFunc {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
typedef typename EOT::ElemType T;
typedef typename EOT::Fitness Fitness;
/**
* Constructor
*/
moCudaEvalFunc() {
}
/**
* Destructor
*/
~moCudaEvalFunc() {
}
/**
*Virtual functor to compute fitness of a solution neighbor
*@param _solution the solution which generate the neighborhood
*@param _fitness the current solution fitness
*@param _index the index neighbor
*/
virtual inline __host__ __device__ Fitness operator() (EOT & _solution,Fitness _fitness, unsigned int * _index){
return _fitness;
};
};
#endif

212
ParadisEO-GPU-good/src/eval/moCudaKswapEval.h Normal file
View file

@ -0,0 +1,212 @@
/*
<moCudaKswapEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef moCudaKswapEval_H
#define moCudaKswapEval_H
#include <eval/moCudaEval.h>
#include <eval/moCudakernelEval.h>
/**
* class for the K-swap neighborhood evaluation
*/
template<class Neighbor, class IncrementEval>
class moCudaKswapEval: public moCudaEval<Neighbor> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
/**
* Define type of a vector corresponding to Solution
*/
typedef typename EOT::ElemType T;
/**
* Define type of a fitness corresponding to Solution
*/
typedef typename EOT::Fitness Fitness;
using moCudaEval<Neighbor>::neighborhoodSize;
using moCudaEval<Neighbor>::host_FitnessArray;
using moCudaEval<Neighbor>::device_FitnessArray;
using moCudaEval<Neighbor>::device_solution;
using moCudaEval<Neighbor>::kernel_Dim;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
* @param _incrEval the incremental evaluation
*/
moCudaKswapEval(unsigned int _neighborhoodSize, IncrementEval & _incrEval) :
moCudaEval<Neighbor> (_neighborhoodSize), incrEval(_incrEval) {
mutex = false;
mutex_kswap = false;
}
/**
* Destructor
*/
~moCudaKswapEval() {
if (mutex_kswap) {
cudaFree(&device_setSolution);
cudaFree(&device_tmp);
delete[] vect;
}
}
/**
* Compute fitness for all solution neighbors in device without specific mapping
* @param _sol the solution which generate the neighborhood
*/
virtual void neighborhoodEval(EOT & _sol) {
}
/**
* Compute fitness for all solution neighbors in device with K-swap mapping
* @param _sol the solution which generate the neighborhood
* @param _mapping the array of mapping indexes for K-swap neighborhood
* @param _Kswap the number of swap
*/
void neighborhoodKswapEval(EOT & _sol, unsigned * _mapping, unsigned _Kswap) {
// the solution size
unsigned _size = _sol.size();
// Get Current solution fitness
Fitness fitness = _sol.fitness();
//Case of Permutation
if (_Kswap == 1) {
if (!mutex) {
//Allocate the space for solution in the device global memory
cudaMalloc((void**) &device_solution.vect, _size * sizeof(T));
mutex = true;
}
//Copy the solution vector from the host to device
cudaMemcpy(device_solution.vect, _sol.vect, _size * sizeof(T),
cudaMemcpyHostToDevice);
//Launch the Kernel to compute all permutation neighbors fitness
kernelPermutation<EOT,Fitness,Neighbor,IncrementEval><<<kernel_Dim,BLOCK_SIZE >>>(incrEval,device_solution,device_FitnessArray,fitness,neighborhoodSize,_mapping,_size);
//Copy the result from device to host
cudaMemcpy(host_FitnessArray, device_FitnessArray, neighborhoodSize
* sizeof(Fitness), cudaMemcpyDeviceToHost);
}
//Case Kswap
else if (_Kswap > 1) {
if (!mutex_kswap) {
vect = new T[neighborhoodSize * _size];
//Allocate the space for set of solution in the device global memory
cudaMalloc((void**) &device_setSolution.vect, neighborhoodSize
* _size * sizeof(T));
//Allocate the space to save temporary EOT element to swap
cudaMalloc((void**) &device_tmp.vect, neighborhoodSize
* sizeof(T));
mutex_kswap = true;
}
for (int i = 0; i < neighborhoodSize; i++) {
for (int j = 0; j < _size; j++) {
vect[j + i * _size] = _sol.vect[j];
}
}
//Copy the set of solution from the host to device
cudaMemcpy(device_setSolution.vect, vect, neighborhoodSize * _size
* sizeof(T), cudaMemcpyHostToDevice);
//Launch the Kernel to compute all Kswap neighbors fitness
kernelKswap<EOT,Fitness,Neighbor,IncrementEval><<<kernel_Dim,BLOCK_SIZE >>>(incrEval,device_setSolution,device_tmp,device_FitnessArray,fitness,neighborhoodSize,_mapping,_Kswap,_size);
//Copy the result from device to host
cudaMemcpy(host_FitnessArray, device_FitnessArray, neighborhoodSize
* sizeof(Fitness), cudaMemcpyDeviceToHost);
}
}
/**
* Compute fitness for all solution neighbors(K-flip of binary solution) in device
* @param _sol the solution which generate the neighborhood
* @param _mapping the array of mapping indexes for k-flip neighborhood
* @param _Kflip the number of flip to do
*/
void neighborhoodKflipEval(EOT & _sol, unsigned * _mapping,
unsigned _Kflip) {
// the solution size
unsigned _size = _sol.size();
// Get Current solution fitness
Fitness fitness = _sol.fitness();
if (!mutex) {
//Allocate the space for solution in the device global memory
cudaMalloc((void**) &device_solution.vect, _size * sizeof(T));
mutex = true;
}
//Copy the solution vector from the host to device
cudaMemcpy(device_solution.vect, _sol.vect, _size * sizeof(T),
cudaMemcpyHostToDevice);
//Launch the Kernel to compute all flip neighbors fitness
kernelKflip<EOT,Fitness,Neighbor,IncrementEval><<<kernel_Dim,BLOCK_SIZE >>>(incrEval,device_solution,device_FitnessArray,fitness,neighborhoodSize,_mapping,_Kflip);
//Copy the result from device to host
cudaMemcpy(host_FitnessArray, device_FitnessArray, neighborhoodSize
* sizeof(Fitness), cudaMemcpyDeviceToHost);
}
protected:
IncrementEval & incrEval;
//NeighborhoodSize copy of solution
EOT device_setSolution;
//NeighborhoodSize element of EOT
EOT device_tmp;
//Vector of neighborhoodSize copy of solution
T * vect;
bool mutex_kswap;
bool mutex;
};
#endif

113
ParadisEO-GPU-good/src/eval/moCudaVectorEval.h Normal file
View file

@ -0,0 +1,113 @@
/*
<moCudaVectorEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef moCudaVectorEval_H
#define moCudavectorEval_H
#include <eval/moCudaEval.h>
#include <eval/moCudakernelEval.h>
/**
* class for the parallel evaluation of neighborhood without mapping
*/
template<class Neighbor, class IncrementEval>
class moCudaVectorEval: public moCudaEval<Neighbor> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
/**
* Define type of a vector corresponding to Solution
*/
typedef typename EOT::ElemType T;
/**
* Define type of a fitness corresponding to Solution
*/
typedef typename EOT::Fitness Fitness;
using moCudaEval<Neighbor>::neighborhoodSize;
using moCudaEval<Neighbor>::host_FitnessArray;
using moCudaEval<Neighbor>::device_FitnessArray;
using moCudaEval<Neighbor>::device_solution;
using moCudaEval<Neighbor>::kernel_Dim;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
* @param _incrEval the incremental evaluation
*/
moCudaVectorEval(unsigned int _neighborhoodSize, IncrementEval & _incrEval) :
moCudaEval<Neighbor> (_neighborhoodSize), incrEval(_incrEval) {
}
/**
* Compute fitness for all solution neighbors in device
* @param _sol the solution which generate the neighborhood
*/
virtual void neighborhoodEval(EOT & _sol) {
// the solution vector size
unsigned _size = _sol.size();
// Get Current solution fitness
Fitness fitness = _sol.fitness();
//Allocate the space for solution in the global memory of device
cudaMalloc((void**) &device_solution.vect, _size * sizeof(T));
//Copy the solution vector from the host to device
cudaMemcpy(device_solution.vect, _sol.vect, _size * sizeof(T),
cudaMemcpyHostToDevice);
//Launch the Kernel to compute all neighbors fitness
kernelEval<EOT,Fitness,Neighbor,IncrementEval><<<kernel_Dim,BLOCK_SIZE >>>(incrEval,device_solution,device_FitnessArray,fitness,neighborhoodSize);
//Copy the result from device to host
cudaMemcpy(host_FitnessArray, device_FitnessArray, neighborhoodSize
* sizeof(Fitness), cudaMemcpyDeviceToHost);
}
protected:
IncrementEval & incrEval;
};
#endif

206
ParadisEO-GPU-good/src/eval/moCudakernelEval.h Normal file
View file

@ -0,0 +1,206 @@
/*
<moCudakernelEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __moCudakernelEval_H
#define __moCudakernelEval_H
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/**
* The kernel function called from the host and executed in device to compute all neighbors fitness at one time
* with linear mapping
* @param _eval how to evaluate each neighbor
* @param _solution the representation of solution( vector of int,float....)
* @param _allFitness the array of Fitness to save all neighbors fitness
* @param _fitness the current solution fitness
* @param _neighborhoodsize the size of the neighborhood
*/
template<class EOT, class Fitness, class Neighbor, class IncrementEval>
__global__ void kernelEval(IncrementEval _eval, EOT _solution, Fitness* _allFitness,
Fitness _fitness, unsigned _neighborhoodsize) {
// The thread identifier within a grid block's
int id = blockIdx.x * blockDim.x + threadIdx.x;
// array to save index to be changed
unsigned int index[1];
// In this representation each id identify one and only one neighbor in neighborhood
if (id < _neighborhoodsize) {
//Change the id'th element of solution
index[0]=id;
//Compute fitness for id'th neighbor
_allFitness[id] = _eval(_solution, _fitness,index);
}
}
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/**
* The kernel function called from the host and executed in device to compute all flip neighbors fitness at one time
* @param _eval how to evaluate each neighbor
* @param _solution representation of solution to flip
* @param _allFitness Array of Fitness type to save all neighbors fitness
* @param _fitness the current solution fitness
* @param _neighborhoodsize the size of the neighborhood
* @param _mapping the neighborhood mapping
* @param _Kflip the number of bit to flip
*/
template<class EOT, class Fitness, class Neighbor, class IncrementEval>
__global__ void kernelKflip(IncrementEval _eval, EOT _solution, Fitness* _allFitness,
Fitness _fitness, unsigned _neighborhoodsize, unsigned * _mapping,unsigned _Kflip) {
// The thread identifier within a grid block's
int id = blockIdx.x * blockDim.x + threadIdx.x;
//save temporary fitness
unsigned tmp_fitness;
//counter of number of flip to do
unsigned i;
// array to save index to be changed
unsigned index[1];
// In this representation each id identify one and only one neighbor in neighborhood
if (id < _neighborhoodsize) {
//Init fitness with fitness of solution
tmp_fitness=_fitness;
//Evaluate neighbor after Kflip
for(i=0;i<=_Kflip;i++){
//The designed index to flip
index[0]=_mapping[id + i * _neighborhoodsize];
//Evaluate the neighbor
tmp_fitness= _eval(_solution, tmp_fitness, index);
}
//The final fitness of the Id'th neighbor
_allFitness[id]=tmp_fitness;
}
}
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/**
* The kernel function called from the host and executed in device to compute all swap neighbors fitness at one time
* @param _eval how to evaluate each neighbor
* @param _solution representation ofsolution to swap
* @param _sol_tmp to save temporary a solution element to swap
* @param _allFitness Array of Fitness type to save all neighbors fitness
* @param _fitness the current solution fitness
* @param _neighborhoodsize the size of the neighborhood
* @param _mapping the neighborhood mapping
* @param _Kswap the number of swap to do
* @param _size the solution size
*/
template<class EOT,class Fitness, class Neighbor, class IncrementEval>
__global__ void kernelKswap(IncrementEval _eval,EOT _solution ,EOT _sol_tmp, Fitness* _allFitness,
Fitness _fitness, unsigned _neighborhoodsize, unsigned * _mapping,unsigned _Kswap,unsigned _size) {
// The thread identifier within a grid block's
int id = blockIdx.x * blockDim.x + threadIdx.x;
//save temporary fitness
int tmp_fitness;
//counter of number of swap to do
unsigned i;
// array to save index to be changed, solution size & thread id
unsigned index[4];
// In this representation each id identify one and only one neighbor in neighborhood
if (id < _neighborhoodsize) {
//the first index to swap
index[0]=_mapping[id];
//the second index to swap
index[1]=_mapping[id +_neighborhoodsize];
//the solution size
index[2]=_size;
//the thread id
index[3]=id;
//Init the temporary fitness with the initial solution fitness
tmp_fitness=_fitness;
//Evaluate neighbor after K-swap
for(i=2;i<=_Kswap+1;i++){
//Evaluate neighbor with index case
tmp_fitness=_eval(_solution, tmp_fitness, index);
//Permut the solution
_sol_tmp[id]=_solution[index[0]+id*index[2]];
_solution[index[0]+id*index[2]]=_solution[index[1]+id*index[2]];
_solution[index[1]+id*index[2]]=_sol_tmp[id];
//Init the next swap to do
index[0]=index[1];
index[1]=_mapping[id +i*_neighborhoodsize];
}
//save the final fitness of the id'th neighbor
_allFitness[id]=tmp_fitness;
}
}
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/**
* The kernel function called from the host and executed in device to compute all permutation neighbors fitness at one time
* @param _eval how to evaluate each neighbor
* @param _solution representation of solution
* @param _allFitness Array of Fitness type to save all neighbors fitness
* @param _fitness the current solution fitness
* @param _neighborhoodsize the size of the neighborhood
* @param _mapping the neighborhood mapping
* @param _size the solution size
*/
template<class EOT, class Fitness, class Neighbor, class IncrementEval>
__global__ void kernelPermutation(IncrementEval _eval, EOT _solution, Fitness* _allFitness,
Fitness _fitness, unsigned _neighborhoodsize, unsigned * _mapping,unsigned _size) {
// The thread identifier within a grid block's
int id = blockIdx.x * blockDim.x + threadIdx.x;
// array to save index to be changed, solution size
unsigned index[4];
// In this representation each id identify one and only one neighbor in neighborhood
if (id < _neighborhoodsize) {
//The first index of permutation
index[0]=_mapping[id];
//The second index of permutation
index[1]=_mapping[id +_neighborhoodsize];
//The solution size
index[2]=_size;
//Puch 0 in the 3 index
index[3]=0;
_allFitness[id]=_eval(_solution,_fitness,index);
}
}
#endif

81
ParadisEO-GPU-good/src/memory/moCudaAllocator.h Normal file
View file

@ -0,0 +1,81 @@
/*
<moCudaAllocator.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaAllocator_H_
#define _moCudaAllocator_H_
// CUDA includes
#include <cutil.h>
/**
* class for allocation data on GPU global memory
*/
class moCudaAllocator {
public:
/**
* Constructor
*/
moCudaAllocator() {
}
/**
*Allocate data on GPU global memory
*@param _data the data to allocate on GPU global memory
*@param _dataSize the size of data to allocate on GPU memory
*/
template<typename T>
void operator()(T* & _data, unsigned _dataSize) {
//Allocate data in GPU memory
CUDA_SAFE_CALL(cudaMalloc((void**) &_data, _dataSize * sizeof(T)));
// Check if data allocation is failed
CUT_CHECK_ERROR("Allocation of data on GPU global memory failed");
}
/**
* Destructor
*/
~moCudaAllocator() {
}
};
#endif

117
ParadisEO-GPU-good/src/memory/moCudaCopy.h Normal file
View file

@ -0,0 +1,117 @@
/*
<moCudaCopy.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaCopy_H_
#define _moCudaCopy_H_
// CUDA includes
#include <cutil.h>
/**
* class to copy data from CPU memory to GPU global memory and vice versa
*/
class moCudaCopy {
public:
/**
* Constructor
*/
moCudaCopy() {
}
/**
*Copy data from CPU memory to GPU global memory (default copy)
*@param _data the data representation where the data will be copied
*@param _dataTocpy the data to copy from CPU memory to GPU memory
*@param _dataSize the size of data to copy
*/
template<typename T>
void operator()(T* & _data, T * & _dataTocpy, unsigned _dataSize) {
//copy data from CPU memory to GPU memory
CUDA_SAFE_CALL(cudaMemcpy(_data, _dataTocpy, _dataSize * sizeof(T),
cudaMemcpyHostToDevice));
// Check if the copy of data is failed
CUT_CHECK_ERROR("Copy of data from CPU to GPU global memory failed");
}
/**
*Copy device data from GPU global memory to global variable declared in device
*@param _dev_data the device global variable
*@param _dataTocpy the data to copy GPU global memory to GPU global variable
*/
template<typename T>
void operator()(T* & _dev_data, T * & _dataTocpy) {
//Copy n bytes from the memory area pointed to by _dataTocpy to the memory area pointed to by offset bytes from the start of symbol _dev_data
cudaMemcpyToSymbol(_dev_data, &_dataTocpy, sizeof(_dataTocpy));
}
/**
*Copy data from CPU memory to GPU global memory and vice versa
*@param _data the data representation where the data will be copied
*@param _dataTocpy the data to copy from CPU memory to GPU memory and vice versa
*@param _dataSize the size of data to copy
*@param _HostToDevice the direction of copy(true if copy will be done from CPU memory to GPU memory)
*/
template<typename T>
void operator()(T* & _data, T * & _dataTocpy, unsigned _dataSize,
bool _HostToDevice) {
if (_HostToDevice) {
//copy data from CPU memory to GPU global memory
CUDA_SAFE_CALL(cudaMemcpy(_data, _dataTocpy, _dataSize * sizeof(T),
cudaMemcpyHostToDevice));
}
else {
//copy data from GPU global memory to GPU memory
CUDA_SAFE_CALL(cudaMemcpy(_data, _dataTocpy, _dataSize * sizeof(T),
cudaMemcpyDeviceToHost));
}
}
~moCudaCopy() {
}
};
#endif

77
ParadisEO-GPU-good/src/memory/moCudaDesallocator.h Normal file
View file

@ -0,0 +1,77 @@
/*
<moCudaDesallocator.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaDesallocator_H_
#define _moCudaDesallocator_H_
// CUDA includes
#include <cutil.h>
/**
* class for Desallocation of data from GPU global memory
*/
class moCudaDesallocator {
public:
/**
* Constructor
*/
moCudaDesallocator() {
}
/**
*Desallocate data on GPU global memory
*@param _data the data to desallocate from GPU global memory
*/
template<typename T>
void operator()(T* & _data) {
//Desallocate data from GPU global memory
cudaFree(_data);
}
/**
* Destructor
*/
~moCudaDesallocator() {
}
};
#endif

91
ParadisEO-GPU-good/src/memory/moCudaObject.h Normal file
View file

@ -0,0 +1,91 @@
/*
<moCudaObject.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaObject_H_
#define _moCudaObject_H_
// CUDA includes
#include <cutil.h>
#include <memory/moCudaAllocator.h>
#include <memory/moCudaDesallocator.h>
#include <memory/moCudaCopy.h>
/**
* class of managment of data on GPU global memory (allocation,desallocation & copy)
*/
class moCudaObject {
public:
/**
*Copy data from CPU memory to GPU global memory (default copy)
*@param _data the data to allocate on GPU
*@param _dataTocpy the data to copy from CPU memory to _data on GPU memory
*@param _dataSize the size of data to copy
*/
template<typename T>
void memCopy(T* & _data, T * & _dataTocpy, unsigned _dataSize) {
malloc(_data, _dataSize);
copy(_data, _dataTocpy, _dataSize);
}
/**
*Copy device data from GPU global memory to global variable declared in device
*@param _dev_data the device global variable
*@param _dataTocpy the data to copy GPU global memory to GPU global variable
*/
template<typename T>
void memCopyGlobalVariable(T* & _dev_data, T * & _dataTocpy) {
copy(_dev_data, _dataTocpy);
}
/**
*Desallocate data on GPU global memory
*@param _data the data to desallocate from GPU global memory
*/
template<typename T>
void memFree(T* & _data) {
free(_data);
}
public:
moCudaAllocator malloc;
moCudaCopy copy;
moCudaDesallocator free;
};
#endif

85
ParadisEO-GPU-good/src/memory/moCudaSpecificData.h Normal file
View file

@ -0,0 +1,85 @@
/*
<moCudaSpecificData.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaSpecificData_H_
#define _moCudaSpecificData_H_
#include <memory/moCudaObject.h>
/**
* class of managment of specific data problem
*/
class moCudaSpecificData {
public:
/*
* Default constructor
*/
moCudaSpecificData() {
sizeData = 0;
}
/*
* Load data from file given in argument
* @param _file the name of data file
*/
virtual void load(char* _file)=0;
/*
* Return the size
* @return the size of data to load
*/
virtual unsigned int getSize() {
return sizeData;
}
/*
* Set the size of data to load
*/
virtual void setSize(unsigned int _size) {
sizeData = _size;
}
public:
unsigned int sizeData;
moCudaObject cudaObject;
};
#endif

121
ParadisEO-GPU-good/src/neighborhood/moCudaBitNeighbor.h Normal file
View file

@ -0,0 +1,121 @@
/*
<moCudaBitNeighbor.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Jerémie Humeau, Boufaras Karima, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaBitNeighbor_h
#define _moCudaBitNeighbor_h
#include <neighborhood/moBackableNeighbor.h>
#include <neighborhood/moIndexNeighbor.h>
/**
* Neighbor related to a solution vector of EOT (type)
*/
template<class EOT, class Fitness = typename EOT::Fitness>
class moCudaBitNeighbor: public moBackableNeighbor<EOT, Fitness> ,
public moIndexNeighbor<EOT, Fitness> {
public:
using moBackableNeighbor<EOT, Fitness>::fitness;
using moIndexNeighbor<EOT, Fitness>::key;
/**
* move the solution
* @param
* _solution the solution to move
*/
virtual void move(EOT & _solution) {
_solution[key] = !_solution[key];
}
/**
* move back the solution (useful for the evaluation by modif)
* @param
* _solution the solution to move back
*/
virtual void moveBack(EOT & _solution) {
_solution[key] = !_solution[key];
}
/**
* Return the class name.
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaBitNeighbor";
}
/**
* Read object.\
* Calls base class, just in case that one had something to do.
* The read and print methods should be compatible and have the same format.
* In principle, format is "plain": they just print a number
* @param _is a std::istream.
* @throw runtime_std::exception If a valid object can't be read.
*/
virtual void readFrom(std::istream& _is) {
std::string fitness_str;
int pos = _is.tellg();
_is >> fitness_str;
if (fitness_str == "INVALID") {
throw std::runtime_error("invalid fitness");
} else {
Fitness repFit;
_is.seekg(pos);
_is >> repFit;
_is >> key;
fitness(repFit);
}
}
/**
* Write object. Called printOn since it prints the object _on_ a stream.
* @param _os A std::ostream.
*/
virtual void printOn(std::ostream& _os) const {
_os << fitness() << ' ' << key << std::endl;
}
};
#endif

119
ParadisEO-GPU-good/src/neighborhood/moCudaKflipNeighborhood.h Normal file
View file

@ -0,0 +1,119 @@
/*
<moKswapNeighborhood.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can ue,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaKswapNeighborhood_h
#define _moCudaKswapNeighborhood_h
#include <neighborhood/moKswapNeighborhood.h>
#include <eval/moCudaEval.h>
/**
* K-flip Neighborhood
*/
template<class N>
class moCudaKflipNeighborhood: public moKswapNeighborhood<N> {
public:
/**
* Define a Neighbor and type of a solution corresponding
*/
typedef N Neighbor;
typedef typename Neighbor::EOT EOT;
using moKswapNeighborhood<Neighbor>::neighborhoodSize;
using moKswapNeighborhood<Neighbor>::currentIndex;
using moKswapNeighborhood<Neighbor>::indices;
using moKswapNeighborhood<Neighbor>::mapping;
using moKswapNeighborhood<Neighbor>::Kswap;
using moKswapNeighborhood<Neighbor>::size;
using moKswapNeighborhood<Neighbor>::mutex;
/**
* Constructor
* @param _size the size of solution
* @param _size the solution size
* @param _Kflip the number of bit to flip
* @param _eval show how to evaluate neighborhood of a solution at one time
*/
moCudaKflipNeighborhood(unsigned int _size, unsigned int _Kflip,
moCudaEval<Neighbor>& _eval) :
moKswapNeighborhood<Neighbor> (_size, _Kflip), eval(_eval) {
sendMapping = true;
cudaMalloc((void**) &device_Mapping, sizeof(unsigned int)
* neighborhoodSize * (_Kflip + 1));
}
;
/**
*Destructor
*/
~moCudaKflipNeighborhood() {
cudaFree(device_Mapping);
}
/**
* Initialization of the neighborhood
* @param _solution the solution to explore
* @param _current the first neighbor
*/
virtual void init(EOT& _solution, Neighbor& _current) {
moKswapNeighborhood<Neighbor>::init(_solution, _current);
if (sendMapping) {
cudaMemcpy(device_Mapping, mapping, (Kswap + 1) * neighborhoodSize
* sizeof(unsigned int), cudaMemcpyHostToDevice);
sendMapping = false;
}
//Compute all neighbors fitness at one time
eval.neighborhoodKflipEval(_solution, device_Mapping, Kswap);
}
/**
* Return the class Name
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaKSwapNeighborhood";
}
protected:
moCudaEval<Neighbor>& eval;
bool sendMapping;
unsigned int * device_Mapping;
};
#endif

119
ParadisEO-GPU-good/src/neighborhood/moCudaKswapNeighborhood.h Normal file
View file

@ -0,0 +1,119 @@
/*
<moKswapNeighborhood.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can ue,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaKswapNeighborhood_h
#define _moCudaKswapNeighborhood_h
#include <neighborhood/moKswapNeighborhood.h>
#include <eval/moCudaEval.h>
/**
* K-Swap Neighborhood
*/
template<class N>
class moCudaKswapNeighborhood: public moKswapNeighborhood<N> {
public:
/**
* Define a Neighbor and type of a solution corresponding
*/
typedef N Neighbor;
typedef typename Neighbor::EOT EOT;
using moKswapNeighborhood<Neighbor>::neighborhoodSize;
using moKswapNeighborhood<Neighbor>::currentIndex;
using moKswapNeighborhood<Neighbor>::indices;
using moKswapNeighborhood<Neighbor>::mapping;
using moKswapNeighborhood<Neighbor>::Kswap;
using moKswapNeighborhood<Neighbor>::size;
using moKswapNeighborhood<Neighbor>::mutex;
/**
* Constructor
* @param _size the size of solution
* @param _Kswap the number of swap
* @param _eval show how to evaluate neighborhood of a solution at one time
*/
moCudaKswapNeighborhood(unsigned int _size, unsigned int _Kswap,
moCudaEval<Neighbor>& _eval) :
moKswapNeighborhood<Neighbor> (_size, _Kswap), eval(_eval) {
sendMapping = true;
cudaMalloc((void**) &device_Mapping, sizeof(unsigned int)
* neighborhoodSize * (Kswap + 1));
}
;
/**
*Destructor
*/
~moCudaKswapNeighborhood() {
cudaFree(device_Mapping);
}
/**
* Initialization of the neighborhood
* @param _solution the solution to explore
* @param _current the first neighbor
*/
virtual void init(EOT& _solution, Neighbor& _current) {
moKswapNeighborhood<Neighbor>::init(_solution, _current);
//Send Mapping for only the first time
if (sendMapping) {
cudaMemcpy(device_Mapping, mapping, (Kswap + 1) * neighborhoodSize
* sizeof(unsigned int), cudaMemcpyHostToDevice);
sendMapping = false;
}
//Compute all neighbors fitness at one time
eval.neighborhoodKswapEval(_solution, device_Mapping, Kswap);
}
/**
* Return the class Name
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaKSwapNeighborhood";
}
protected:
moCudaEval<Neighbor>& eval;
bool sendMapping;
unsigned int * device_Mapping;
};
#endif

96
ParadisEO-GPU-good/src/neighborhood/moCudaOrderNeighborhood.h Normal file
View file

@ -0,0 +1,96 @@
/*
<moCudaOrderNeighborhood.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Jerémie Humeau, Thé Van LUONG, Karima Boufaras
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaOrderNeighborhood_h
#define _moCudaOrderNeighborhood_h
#include <neighborhood/moOrderNeighborhood.h>
#include <eval/moCudaEval.h>
/**
* An ordered neighborhood with parallel evaluation
*/
template<class N>
class moCudaOrderNeighborhood: public moOrderNeighborhood<N> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef N Neighbor;
typedef typename Neighbor::EOT EOT;
using moOrderNeighborhood<Neighbor>::neighborhoodSize;
using moOrderNeighborhood<Neighbor>::currentIndex;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
* @param _eval show how to evaluate neighborhood of a solution at one time
*/
moCudaOrderNeighborhood(unsigned int _neighborhoodSize,
moCudaEval<Neighbor>& _eval) :
moOrderNeighborhood<Neighbor> (_neighborhoodSize), eval(_eval) {
}
/**
* Initialization of the neighborhood
*@param _solution the solution to explore
*@param _neighbor the first neighbor
*/
virtual void init(EOT & _solution, Neighbor & _neighbor) {
moOrderNeighborhood<Neighbor>::init(_solution, _neighbor);
//Compute all neighbors fitness at one time
eval.neighborhoodEval(_solution);
}
/**
* Return the class name.
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaOrderNeighborhood";
}
protected:
moCudaEval<Neighbor>& eval;
};
#endif

88
ParadisEO-GPU-good/src/neighborhood/moCudaRndWithReplNeighborhood.h Normal file
View file

@ -0,0 +1,88 @@
/*
<moCudaRndWithReplNeighborhood.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Jerémie Humeau, Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaRndWithReplNeighborhood_h
#define _moCudaRndWithReplNeighborhood_h
#include <neighborhood/moRndWithReplNeighborhood.h>
#include <eval/moCudaEval.h>
/**
* A Random With replacement Neighborhood with parallel evaluation
*/
template<class Neighbor>
class moCudaRndWithReplNeighborhood: public moRndWithReplNeighborhood<Neighbor> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
using moRndWithReplNeighborhood<Neighbor>::neighborhoodSize;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
* @param _eval show how to evaluate neighborhood of a solution at one time
*/
moCudaRndWithReplNeighborhood(unsigned int _neighborhoodSize, moCudaEval<
Neighbor>& _eval) :
moRndWithReplNeighborhood<Neighbor> (_neighborhoodSize), eval(_eval) {
}
/**
* Initialization of the neighborhood
* @param _solution the solution to explore
* @param _neighbor the first neighbor
*/
virtual void init(EOT & _solution, Neighbor & _neighbor) {
moRndWithReplNeighborhood<Neighbor>::init(_solution, _neighbor);
//Compute all neighbors fitness at one time
eval.neighborhoodEval(_solution);
}
/**
* Return the class Name
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaRndWithReplNeighborhood";
}
protected:
moCudaEval<Neighbor>& eval;
};
#endif

90
ParadisEO-GPU-good/src/neighborhood/moCudaRndWithoutReplNeighborhood.h Normal file
View file

@ -0,0 +1,90 @@
/*
<moCudaRndWithoutReplNeighborhood.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Jerémie Humeau, Boufaras Karima, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaRndWithoutReplNeighborhood_h
#define _moCudaRndWithoutReplNeighborhood_h
#include <neighborhood/moRndWithoutReplNeighborhood.h>
#include <eval/moCudaEval.h>
/**
* A Random without replacement Neighborhood with parallel evaluation
*/
template<class Neighbor>
class moCudaRndWithoutReplNeighborhood: public moRndWithoutReplNeighborhood<Neighbor> {
public:
/**
* Define type of a solution corresponding to Neighbor
*/
typedef typename Neighbor::EOT EOT;
using moRndWithoutReplNeighborhood<Neighbor>::neighborhoodSize;
using moRndWithoutReplNeighborhood<Neighbor>::maxIndex;
using moRndWithoutReplNeighborhood<Neighbor>::indexVector;
/**
* Constructor
* @param _neighborhoodSize the size of the neighborhood
* @param _eval show how to evaluate neighborhood of a solution at one time
*/
moCudaRndWithoutReplNeighborhood(unsigned int _neighborhoodSize,moCudaEval<
Neighbor>& _eval) :
moRndWithoutReplNeighborhood<Neighbor> (_neighborhoodSize),eval(_eval) {
for (unsigned int i = 0; i < neighborhoodSize; i++)
indexVector.push_back(i);
}
/**
* Initialization of the neighborhood
* @param _solution the solution to explore
* @param _neighbor the first neighbor
*/
virtual void init(EOT & _solution, Neighbor & _neighbor) {
moRndWithoutReplNeighborhood<Neighbor>::init(_solution, _neighbor);
//Compute all neighbors fitness at one time
eval.neighborhoodEval(_solution);
}
/**
* Return the class Name
* @return the class name as a std::string
*/
virtual std::string className() const {
return "moCudaRndWithoutReplNeighborhood";
}
protected:
moCudaEval<Neighbor>& eval;
};
#endif

116
ParadisEO-GPU-good/src/performance/moCudaTimer.h Normal file
View file

@ -0,0 +1,116 @@
/*
<moCudaTimer.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _moCudaTimer_H_
#define _moCudaTimer_H_
// CUDA Utility Tools
#include <cutil.h>
/**
* To compute execution time
*/
class moCudaTimer {
public:
unsigned int timer;
/**
* Constructor
*/
moCudaTimer() {
timer = 0;
//Create timer
cutCreateTimer(&timer);
}
/**
* Destructor
*/
~moCudaTimer() {
}
/**
* to start compute execution time
*/
void start() {
// Start timer
cutStartTimer(timer);
}
/**
* to stop compute execution time
*/
void stop() {
// Stop timer
cutStopTimer(timer);
}
/**
* to get timer value in ms
* @return execution time in ms
*/
double getTime() {
// get time
return cutGetTimerValue(timer);
}
/**
* to delete timer
*/
void deleteTimer() {
//Delete timer
cutDeleteTimer(timer);
}
};
#endif

178
ParadisEO-GPU-good/src/problems/data/PPPData.h Normal file
View file

@ -0,0 +1,178 @@
/*
<PPPData.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _PPPData_H_
#define _PPPData_H_
#include <memory/moCudaSpecificData.h>
template<class ElemType>
class PPPData: public moCudaSpecificData {
public:
using moCudaSpecificData::cudaObject;
/**
* Default Constructor
*/
PPPData() :
moCudaSpecificData() {
(*this).load();
}
/**
* Constructor by copy
* @param _pppData the specific data of PPP
*/
PPPData(const PPPData & _pppData) {
a_h = new int[Md * Nd];
H_h = new int[Nd + 1];
for (int i = 0; i < Md; i++)
for (int j = 0; j < Nd; j++) {
a_h[i * Nd + j] = _pppData.a_h[i * Nd + j];
}
for (int k = 0; k <= Nd; k++) {
H_h[k] = _pppData.H_h[k];
}
cudaObject.memCopy(a_d, a_h, Nd * Md);
cudaObject.memCopy(H_d, H_h, Nd + 1);
}
/**
* Assignement operator
* @param _pppData the specific data of PPP
* @return a PPP Data
*/
PPPData & operator=(const PPPData & _pppData) {
a_h = new int[Md * Nd];
H_h = new int[Nd + 1];
for (int i = 0; i < Md; i++)
for (int j = 0; j < Nd; j++) {
a_h[i * Nd + j] = _pppData.a_h[i * Nd + j];
}
for (int k = 0; k <= Nd; k++) {
H_h[k] = _pppData.H_h[k];
}
cudaObject.memCopy(a_d, a_h, Nd * Md);
cudaObject.memCopy(H_d, H_h, Nd + 1);
return (*this);
}
/*
* Destructor
*/
~PPPData() {
cudaObject.memFree(a_d);
cudaObject.memFree(H_d);
delete[] a_h;
delete[] H_h;
}
/*
*Load PPP data
*/
void load(char * _fileName) {
}
void load() {
int *v = new int[Nd];
int *s = new int[Md];
a_h = new int[Md * Nd];
H_h = new int[Nd + 1];
for (int i = 0; i < Md; i++) {
for (int j = 0; j < Nd; j++) {
if (rng.rand() % 2)
a_h[i * Nd + j] = 1;
else
a_h[i * Nd + j] = -1;
}
}
for (int i = 0; i < Nd; i++) {
if (rng.rand() % 2)
v[i] = 1;
else
v[i] = -1;
}
for (int i = 0; i < Md; i++) {
s[i] = 0;
for (int j = 0; j < Nd; j++)
s[i] += a_h[i * Nd + j] * v[j];
if (s[i] < 0) {
for (int k = 0; k < Nd; k++)
a_h[i * Nd + k] = -a_h[i * Nd + k];
s[i] = -s[i];
}
}
for (int i = 0; i <= Nd; i++)
H_h[i] = 0;
for (int i = 0; i < Md; i++)
H_h[s[i]]++;
//Allocate and copy QAP data from CPU memory to GPU global memory
cudaObject.memCopy(a_d, a_h, Nd * Md);
cudaObject.memCopy(H_d, H_h, Nd + 1);
delete[] v;
delete[] s;
}
public:
ElemType* a_h;
ElemType* H_h;
ElemType* a_d;
ElemType* H_d;
};
#endif

167
ParadisEO-GPU-good/src/problems/data/QAPData.h Normal file
View file

@ -0,0 +1,167 @@
/*
<QAPData.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _QAPData_H_
#define _QAPData_H_
#include <memory/moCudaSpecificData.h>
template<class ElemType>
class QAPData: public moCudaSpecificData {
public:
using moCudaSpecificData::sizeData;
using moCudaSpecificData::cudaObject;
/**
* Default Constructor
*/
QAPData() :
moCudaSpecificData() {
}
/**
* Constructor
* @param _fileName the data file name
*/
QAPData(char* _fileName) {
(*this).load(_fileName);
}
/**
* Constructor by copy
* @param _qapData the specific data of QAP
*/
QAPData(const QAPData & _qapData) {
sizeData = _qapData.sizeData;
a_h = new int[sizeData * sizeData];
b_h = new int[sizeData * sizeData];
for (int i = 0; i < sizeData; i++)
for (int j = 0; j < sizeData; j++) {
a_h[i * sizeData + j] = _qapData.a_h[i * sizeData + j];
b_h[i * sizeData + j] = _qapData.b_h[i * sizeData + j];
}
cudaObject.memCopy(a_d, a_h, sizeData * sizeData);
cudaObject.memCopy(b_d, b_h, sizeData * sizeData);
}
/**
* Assignement operator
* @param _qapData the specific data of QAP
* @return a QAP Data
*/
QAPData & operator=(const QAPData & _qapData) {
sizeData = _qapData.sizeData;
a_h = new int[sizeData * sizeData];
b_h = new int[sizeData * sizeData];
for (int i = 0; i < sizeData; i++)
for (int j = 0; j < sizeData; j++) {
a_h[i * sizeData + j] = _qapData.a_h[i * sizeData + j];
b_h[i * sizeData + j] = _qapData.b_h[i * sizeData + j];
}
cudaObject.memCopy(a_d, a_h, sizeData * sizeData);
cudaObject.memCopy(b_d, b_h, sizeData * sizeData);
return (*this);
}
/*
* Destructor
*/
~QAPData() {
cudaObject.memFree(a_d);
cudaObject.memFree(b_d);
delete[] a_h;
delete[] b_h;
}
/*
*Load QAP data from file name
*@param _fileName the data file name to load
*/
void load(char* _fileName) {
fstream file(_fileName, ios::in);
if (!file) {
string str = "QAPData: Could not open file [" + (string) _fileName
+ "].";
throw runtime_error(str);
}
unsigned i, j;
file >> sizeData;
a_h = new ElemType[sizeData * sizeData];
b_h = new ElemType[sizeData * sizeData];
for (i = 0; i < sizeData; i++)
for (j = 0; j < sizeData; j++)
file >> a_h[i * sizeData + j];
for (i = 0; i < sizeData; i++)
for (j = 0; j < sizeData; j++)
file >> b_h[i * sizeData + j];
//Allocate and copy QAP data from CPU memory to GPU global memory
cudaObject.memCopy(a_d, a_h, sizeData * sizeData);
cudaObject.memCopy(b_d, b_h, sizeData * sizeData);
}
public:
ElemType* a_h;
ElemType* b_h;
ElemType* a_d;
ElemType* b_d;
};
#endif

79
ParadisEO-GPU-good/src/problems/eval/EvalOneMax.h Normal file
View file

@ -0,0 +1,79 @@
/*
<EvalOneMax.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __EvalOneMax
#define __EvalOneMax
/**
* Full Evaluation of the solution
*/
template<class EOT>
class EvalOneMax: public eoEvalFunc<EOT> {
public:
/**
* Constructor
*/
EvalOneMax() {
}
/**
* Destructor
*/
~EvalOneMax(void) {
}
/**
* Full evaluation of the solution
* @param _bitVector the solution to evaluate
*/
void operator()(EOT & _bitVector) {
unsigned sum = 0;
for (unsigned i = 0; i < _bitVector.size(); i++)
sum += _bitVector[i];
//set the solution fitness
_bitVector.fitness(sum);
}
};
#endif

90
ParadisEO-GPU-good/src/problems/eval/OneMaxIncrEval.h Normal file
View file

@ -0,0 +1,90 @@
/*
<OneMaxIncrEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __OneMaxIncrEval_H
#define __OneMaxIncrEval_H
#include <eval/moCudaEvalFunc.h>
/**
* Incremental Evaluation of OneMax
*/
template<class Neighbor>
class OneMaxIncrEval: public moCudaEvalFunc<Neighbor> {
public:
typedef typename Neighbor::EOT EOT;
typedef typename EOT::Fitness Fitness;
/**
* Constructor
*/
OneMaxIncrEval() {
}
/**
* Destructor
*/
~OneMaxIncrEval() {
}
/**
* Incremental evaluation of the solution(function inline can be called from host or device)
* @param _bitVector the solution to evaluate
* @param _fitness the fitness of the current solution
* @param _index the index of solution neighbor
*/
inline __host__ __device__ Fitness operator() (EOT & _bitVector,Fitness _fitness, unsigned int * _index) {
Fitness tmp;
if (_bitVector[_index[0]] == 0)
tmp= _fitness+1;
else
tmp= _fitness-1;
return tmp;
}
};
#endif

104
ParadisEO-GPU-good/src/problems/eval/PPPEval.h Normal file
View file

@ -0,0 +1,104 @@
/*
<PPPEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __PPPEval
#define __PPPEval
#include <problems/data/PPPData.h>
template<class EOT, class ElemType = typename EOT::ElemType>
class PPPEval: public eoEvalFunc<EOT> {
public:
/**
* Constructor
* @param _qapData the specific data problem useful to evalute solution(flow & distance matrices of QAP problem)
*/
PPPEval(PPPData<ElemType> & _pppData) {
pppData = _pppData;
}
/**
* Destructor
*/
~PPPEval() {
}
/**
* Full evaluation of the solution
* @param _sol the solution to evaluate
*/
void operator()(EOT & _sol) {
int *H;
int *S;
int tmp_1=0;
int tmp_2=0;
H=new int[Nd+1];
S=new int[Md];
for (int i=0; i<Md; i++){
S[i]=0;
for (int j=0; j<Nd; j++){
S[i]+=pppData.a_h[i*Nd+j]*_sol[j];
}
tmp_1+=abs(S[i])-S[i];
if(S[i]>=0)
H[S[i]]++;
}
for (int j=0; j<=Nd; j++){
tmp_2+=abs(pppData.H_h[j]-H[j]);
}
_sol.fitness(ca*tmp_1+cb*tmp_2);
delete[] H;
delete[] S;
}
protected:
PPPData<ElemType> pppData;
};
#endif

105
ParadisEO-GPU-good/src/problems/eval/PPPIncrEval.h Normal file
View file

@ -0,0 +1,105 @@
/*
<PPPIncrEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __PPPIncrEval_H
#define __PPPIncrEval_H
#include <eval/moCudaEvalFunc.h>
/**
* Incremental Evaluation of PPP
*/
template<class Neighbor>
class PPPIncrEval: public moCudaEvalFunc<Neighbor> {
public:
typedef typename Neighbor::EOT EOT;
typedef typename EOT::Fitness Fitness;
/**
* Constructor
*/
PPPIncrEval() {
}
/**
* Destructor
*/
~PPPIncrEval() {
}
/**
* functor of incremental evaluation of the solution(function inline can be called from host or device)
* @param _sol the solution to evaluate
* @param _fitness the fitness of the current solution
* @param _index the set of information to compute fitness neighbor
*/
inline __host__ __device__ Fitness operator() (EOT & _sol,Fitness _fitness, unsigned int *_index) {
int H[Nd+1];
int S[Md];
int tmp_1=0;
int tmp_2=0;
for (int i=0; i<Md; i++){
S[i]=0;
for (int j=0; j<Nd; j++){
if(_sol[j])
S[i]=S[i]+dev_a[i*Nd+j];
else
S[i]=S[i]-dev_a[i*Nd+j];
}
tmp_1=tmp_1+abs(S[i])-S[i];
if(S[i]>=0)
H[S[i]]=H[S[i]]+1;
}
for (int j=0; j<=Nd; j++){
tmp_2=tmp_2+abs(dev_h[j]-H[j]);
}
return ca*tmp_1+cb*tmp_2;
}
};
#endif

84
ParadisEO-GPU-good/src/problems/eval/QAPEval.h Normal file
View file

@ -0,0 +1,84 @@
/*
<QAPEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __QAPEval
#define __QAPEval
#include <problems/data/QAPData.h>
template<class EOT, class ElemType = typename EOT::ElemType>
class QAPEval: public eoEvalFunc<EOT> {
public:
/**
* Constructor
* @param _qapData the specific data problem useful to evalute solution(flow & distance matrices of QAP problem)
*/
QAPEval(QAPData<ElemType> & _qapData) {
qapData = _qapData;
}
/**
* Destructor
*/
~QAPEval() {
}
/**
* Full evaluation of the solution
* @param _sol the solution to evaluate
*/
void operator()(EOT & _sol) {
int cost = 0;
unsigned int size = qapData.getSize();
for (unsigned int i = 0; i < size; i++)
for (unsigned int j = 0; j < size; j++) {
cost += qapData.a_h[i * size + j] * qapData.b_h[_sol[i] * size
+ _sol[j]];
}
_sol.fitness(cost);
}
protected:
QAPData<ElemType> qapData;
};
#endif

121
ParadisEO-GPU-good/src/problems/eval/QAPIncrEval.h Normal file
View file

@ -0,0 +1,121 @@
/*
<QAPIncrEval.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef __QAPIncrEval_H
#define __QAPIncrEval_H
#include <eval/moCudaEvalFunc.h>
/**
* Incremental Evaluation of QAP
*/
template<class Neighbor>
class QAPIncrEval: public moCudaEvalFunc<Neighbor> {
public:
typedef typename Neighbor::EOT EOT;
typedef typename EOT::Fitness Fitness;
/**
* Constructor
*/
QAPIncrEval() {
}
/**
* Destructor
*/
~QAPIncrEval() {
}
/**
* functor of incremental evaluation of the solution(function inline can be called from host or device)
* @param _sol the solution to evaluate
* @param _fitness the fitness of the current solution
* @param _index the set of information to compute fitness neighbor
*/
inline __host__ __device__ Fitness operator() (EOT & _sol,Fitness _fitness, unsigned int *_index) {
Fitness tmp;
/**
* dev_a & dev_b are global device variable, data specific to QAP problem (flow & distance matices)
* _index[0] the first position of swap
* _index[1] the second position of swap
* _index[2] the solution size
* _index[3] the id of neighbor
*/
tmp=_fitness+compute_delta(dev_a,dev_b,_sol,_index[0],_index[1],_index[2],_index[3]);
return tmp;
}
/**
* compute the new fitness of the solution after permutation of pair(i,j)(function inline called from host device)
* @param _a the flow matrix of size*size (specific data of QAP problem must be declared as global device variable)
* @param _b the distance matrix of size*size (specific data of QAP problem must be declared as global device variable)
* @param _sol the solution to evaluate
* @param _i the first position of swap
* @param _j the second position of swap
* @param _size the id of neighbor
* @param _id the neighbor identifier
*/
inline __host__ __device__ int compute_delta(int * _a,int * _b,EOT & _sol, int _i, int _j,int _size, int _id){
int d;
int k;
d = (_a[_i*_size+_i]-_a[_j*_size+_j])*(_b[_sol[_j+_id*_size]*_size+_sol[_j+_id*_size]]-_b[_sol[_i+_id*_size]*_size+_sol[_i+_id*_size]]) +
(_a[_i*_size+_j]-_a[_j*_size+_i])*(_b[_sol[_j+_id*_size]*_size+_sol[_i+_id*_size]]-_b[_sol[_i+_id*_size]*_size+_sol[_j+_id*_size]]);
for (k = 0; k < _size; k=k+1)
if (k!=_i && k!=_j)
d = d + (_a[k*_size+_i]-_a[k*_size+_j])*(_b[_sol[k+_id*_size]*_size+_sol[_j+_id*_size]]-_b[_sol[k+_id*_size]*_size+_sol[_i+_id*_size]]) +
(_a[_i*_size+k]-_a[_j*_size+k])*(_b[_sol[_j+_id*_size]*_size+_sol[k+_id*_size]]-_b[_sol[_i+_id*_size]*_size+_sol[k+_id*_size]]);
return(d);
}
};
#endif

103
ParadisEO-GPU-good/src/problems/neighborhood/PPPNeighbor.h Normal file
View file

@ -0,0 +1,103 @@
/*
<PPPNeighbor.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Boufaras Karima, Thé Van Luong
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can ue,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _PPPNeighbor_h
#define _PPPNeighbor_h
#include <neighborhood/moBackableNeighbor.h>
#include <neighborhood/moIndexSwapNeighbor.h>
/**
* PPP Neighbor
*/
template<class EOT, class Fitness = typename EOT::Fitness>
class PPPNeighbor: public moBackableNeighbor<EOT> ,
public moIndexSwapNeighbor<EOT> {
public:
using moIndexSwapNeighbor<EOT>::indices;
using moIndexSwapNeighbor<EOT>::Kswap;
using moIndexSwapNeighbor<EOT>::size;
/**
*Default Constructor
*/
PPPNeighbor() :
moIndexSwapNeighbor<EOT> () {
}
/**
* Constructor
* @param _Kflip the number of bit to flip
*/
PPPNeighbor(unsigned int _Kflip) :
moIndexSwapNeighbor<EOT> (_Kflip) {
}
/**
* Apply the K-Flip in solution
* @param _solution the solution to move
*/
virtual void move(EOT& _solution) {
size = _solution.size();
if (!Kswap) {
_solution[indices[Kswap]] = -_solution[indices[Kswap]];
} else {
for (unsigned int i = 0; i <= Kswap; i++) {
_solution[indices[i]] = -_solution[indices[i]];
}
}
_solution.invalidate();
}
/**
* apply the K-Flip to restore the solution (use by moFullEvalByModif)
* @param _solution the solution to move back
*/
virtual void moveBack(EOT& _solution) {
move(_solution);
}
/**
* Return the class name.
* @return the class name as a std::string
*/
virtual std::string className() const {
return "PPPNeighbor";
}
};
#endif

109
ParadisEO-GPU-good/src/problems/types/PPPSolution.h Normal file
View file

@ -0,0 +1,109 @@
/*
<PPPSolution.h>
Copyright (C) DOLPHIN Project-Team, INRIA Lille - Nord Europe, 2006-2010
Karima Boufaras, Thé Van LUONG
This software is governed by the CeCILL license under French law and
abiding by the rules of distribution of free software. You can use,
modify and/ or redistribute the software under the terms of the CeCILL
license as circulated by CEA, CNRS and INRIA at the following URL
"http://www.cecill.info".
As a counterpart to the access to the source code and rights to copy,
modify and redistribute granted by the license, users are provided only
with a limited warranty and the software's author, the holder of the
economic rights, and the successive licensors have only limited liability.
In this respect, the user's attention is drawn to the risks associated
with loading, using, modifying and/or developing or reproducing the
software by the user in light of its specific status of free software,
that may mean that it is complicated to manipulate, and that also
therefore means that it is reserved for developers and experienced
professionals having in-depth computer knowledge. Users are therefore
encouraged to load and test the software's suitability as regards their
requirements in conditions enabling the security of their systems and/or
data to be ensured and, more generally, to use and operate it in the
same conditions as regards security.
The fact that you are presently reading this means that you have had
knowledge of the CeCILL license and that you accept its terms.
ParadisEO WebSite : http://paradiseo.gforge.inria.fr
Contact: paradiseo-help@lists.gforge.inria.fr
*/
#ifndef _PPPSolution_H_
#define _PPPSolution_H_
#include <cudaType/moCudaVector.h>
/**
* Implementation of PPP vector representation on CUDA.
*/
template<class Fitness>
class PPPSolution: public moCudaVector<int, Fitness> {
public:
using moCudaVector<int, Fitness>::vect;
using moCudaVector<int, Fitness>::N;
/**
* Default constructor.
*/
PPPSolution() :
moCudaVector<int, Fitness> () {
}
/**
*Constructor.
*@param _size The neighborhood size.
*/
PPPSolution(unsigned _size) {
N = _size;
vect = new int[_size];
create();
}
/**
*Assignment operator
*@param _vector The vector passed to the function determine the new content.
*@return a new vector.
*/
PPPSolution& operator=(const PPPSolution & _vector) {
N = _vector.N;
vect = new int[N];
for (unsigned i = 0; i < N; i++)
vect[i] = _vector.vect[i];
fitness(_vector.fitness());
return (*this);
}
/**
*Initializer of random PPP vector.
*/
void create() {
for (int i = 0; i <N; i++){
if((rng.rand()%2)==0)
vect[i] = -1;
else
vect[i] = 1;
}
}
};
#endif