C++

Template, OpenMP et clause reduction

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : Template, OpenMP et clause reduction

Vel-Ryphon

Salut,

J'ai une fonction template resolvant une équation de Poisson.
Dans cette méthode je fais appel des boucles que je parallelise avec OpenMP.
A un moment donné de l'algo, une boucle incrémente une variable extérieure à la boucle, donc je la passe en reduction.

Les problèmes arrivent lorsque mon template prendre des types non scalaire, par exemple un objet Coord_2D (pour gérer des coordonnées 2D).

Du coup la compile plante sur le pragma open mp reduction car la variable n'est plus de type scalaire.

J'ai donc besoin de commandes préprocesseur pour tester le type du template, et adapter les déclarations de variables / directives OpenMP en fonction.

Considérons la déclaration suivante :

Code :

template <class type> void Poisson_solver(..)

J'ai essayé ça :

Code :

typedef Coord_2D<float> coord_2d_float
...
#if type == coord_2d_
..
#endif

Mais cela ne semble pas fonctionner.

J'ai également essayé en passant un 2ème paramètre au template, comme ça :

Code :

template <class type, int dimension> void Poisson_solver(..)
..
#if dimension == 2
..
#endif

[/cpp]

Mais ça ne semble pas fonctionner non plus.

Est-il possible de tester à la précompile le type du template ?

Merci d'avance !

Publicité

Joel F

Real men use unique_ptr

Alors deja, on n'essaye aps des trucs au hasard en pensant que ca amrche.
Le preprocesseur est appelé en amont de la compilation, la résolution des types en plein mileiux, donc jamais ca amrche.

c'est pas du preprocesseur qu'il te faut mais meta-programmer deux variantes de ta boucle que tu spécialise selon le retour de la meta-fonction boost::is_fundamental qui renvoit un booleen statique si le type testé est un type atomique du langage.

en gros

Code :

template<class T, bool isFundamental> struct mon_solveur_impl;
template<class T> struct mon_solveur : mon_solveur_impl<T, typename boost::is_fundamental<T>::value> {};
template<class T> struct mon_solveur_impl<T,true>
{
void solve()
{
// T est fundamental, j'ecris du OpenMP sans probleme
}
};
template<class T> struct mon_solveur_impl<T,false>
{
void solve()
{
// T n'est pas fundamental, j'ecris autre chose
}
};

http://www.boost.org/doc/libs/1_38 [...] ental.html

Message édité par Joel F le 14-03-2009 à 09:21:06

Vel-Ryphon

Le truc c'est que je ne veux pas réécrire 2 fois l'agorithmes du solver, c'est pas beau à mon gout. Je veux une seule fonction template qui gère tout les cas.

---------------
http://www.ryphon-mechanics.com

Joel F

Real men use unique_ptr

C'est pas bien compliqué. Si la seule difference entre la version OpenMP et l'autre c'est les pragma openMP tu ecris :

Code :

template<class T, bool isFundamental> struct mon_solveur_impl;
template<class T> struct mon_solveur : mon_solveur_impl<T, typename boost::is_fundamental<T>::value> {};
template<class T> struct mon_solveur_impl<T,false>
{
void solve()
{
// T n'est pas fundamental, j'ecris autre chose
mon_solveur_qui_fait_le_vrai_calcul<T>();
}
};
template<class T> struct mon_solveur_impl<T,true>
{
void solve()
{
// T est fundamental, j'ecris du OpenMP sans probleme
#pragma omp parallel ????
mon_solveur_qui_fait_le_vrai_calcul<T>();
}
};

Aprés tu n'as pas le choix, le controle sur les types se fait là et pas avant.

Message édité par Joel F le 14-03-2009 à 10:24:46

FORUM HardWare.fr

Programmation

C++

Template, OpenMP et clause reduction

Sujets relatifs
Réalisation de la somme parallèle des éléments d’une matrice: openmp	foncteur et template
template instanciation explicite	template de déclaration d'un copy-constructeur en friend (combo \o/)
SVP, Comment inserer une template flash sur dreamweaver	XSLT template
template	template et functor
template et version des compilateurs	classe template avec pointeur membre sur la classe
Plus de sujets relatifs à : Template, OpenMP et clause reduction

Page générée en 0.050 secondes