Le fichier source

5.1.1.2.2 : Le fichier source

Écrivons le fichier function_vectorize.cpp :

Nous commençons par inclure le header de la bibliothèque standard de math :

1 2	#include <cmath> #include "phoenix_intrinsics.h"

Ensuite, nous incluons notre header :

1	#include "function_vectorize.h"

Ensuite nous implémentons l'appel d'une fonction std::COMPUTE_FUNCTION_DEF définie par une macro en prennant en compte le fait que nous devons indiquer au compilateur que les tableaux que nous utiliserons seront bien alignés en mémoire sur PLIB_VECTOR_SIZE_BYTE_FLOAT octets avec la fonction __builtin_assume_aligned :

///Do the function call
/**	@param[out] ptabResult : table of results of function(tabX)
 * 	@param ptabX : input table
 * 	@param nbElement : number of elements in the tables
*/
void function_vectorize(float* __restrict__ ptabResult, const float* __restrict__ ptabX, size_t nbElement){
	float* tabResult = (float*)__builtin_assume_aligned(ptabResult, PLIB_VECTOR_SIZE_BYTE_FLOAT);
	const float* tabX = (const float*)__builtin_assume_aligned(ptabX, PLIB_VECTOR_SIZE_BYTE_FLOAT);
	
	for(size_t i(0lu); i < nbElement; ++i){
		tabResult[i] = std::COMPUTE_FUNCTION_DEF(tabX[i]);
	}
}

Le fichier function_vectorize.cpp complet :

/***************************************
	Auteur : Pierre Aubert
	Mail : aubertp7@gmail.com
	Licence : CeCILL-C
****************************************/

#include <cmath>
#include "phoenix_intrinsics.h"

#include "function_vectorize.h"

///Do the function call
/**	@param[out] ptabResult : table of results of function(tabX)
 * 	@param ptabX : input table
 * 	@param nbElement : number of elements in the tables
*/
void function_vectorize(float* __restrict__ ptabResult, const float* __restrict__ ptabX, size_t nbElement){
	float* tabResult = (float*)__builtin_assume_aligned(ptabResult, PLIB_VECTOR_SIZE_BYTE_FLOAT);
	const float* tabX = (const float*)__builtin_assume_aligned(ptabX, PLIB_VECTOR_SIZE_BYTE_FLOAT);
	
	for(size_t i(0lu); i < nbElement; ++i){
		tabResult[i] = std::COMPUTE_FUNCTION_DEF(tabX[i]);
	}
}

Vous pouvez le télécharger ici.