C++

RNG et multithreading

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : RNG et multithreading

Lan Wezel

Bonjour,

J'utilise boost::mt19937 comme RNG.
Avant de multithreader mon code je l'avait mis dans un singleton pour le seeder une seule fois pour toutes.

Code :

class RandomGenerator : public boost::mt19937{
static ml::RandomGenerator m_instance;
RandomGenerator(){
seed( time(NULL) * getpid() );
}
~RandomGenerator(){}
public:
static RandomGenerator& instance(){
return m_instance;
}
};

Maintenant je suis en train de multithreader mon code avec OPENMP, et je voudrais faire ça un peu propremment. Comme je voudrais avoir le même RNG dans tous les threads, je suis confronté au problème suivant :
L'appel est donc de la forme :

Code :

boost::bernoulli_distribution<> distrib(threshold);
event = distrib(RandomGenerator::instance());

sauf que j'ai cru comprendre que boost::mt19937 n'est pas thread-safe. Donc je peux avoir deux threads qui appellent le RNG en même temps et qui vont donc avoir le même résultat. Et je ne vois pas trop comment m'en sortir. En effet, ce qui devrait être locké, ce n'est pas le fait de demander l'instance du RNG, mais l'utilisation qu'en fait la distribution.
Est-ce que je vais devoir me rabattre sur un RNG par thread (comme j'utilise OpenMP je sais pas trop comment ça va être possible) ou y a-t-il un solution à ce problème ?

Message édité par Lan Wezel le 28-10-2010 à 18:19:57

Publicité

Lan Wezel

J'ai trouvé une solution, mais je ne sais pas si c'est très propre.
Comme l'operateur () des random::distribution est template de l'Engine (RNG), j'ai réimplémenté l'operateur () du RNG en mettant une section critique.

Code :

class RandomGenerator : public boost::mt19937{
static ml::RandomGenerator m_instance;
RandomGenerator(){
seed(time(NULL) * getpid() );
}
~RandomGenerator(){}
public:
static RandomGenerator& instance(){
return m_instance;
}
boost::mt19937::result_type operator()(){
boost::mt19937::result_type r;
#pragma omp critical
{
r = boost::mt19937::operator()();
}
return r;
}
};

Je suis ouvert à toute critique (acerbe ou non) sur la propreté de cette solution car l'operateur boost::mt19937::operator()() n'est pas "virtual". Il aurait peut-être été plus propre de le wrapper complètement plutôt que d'en hériter.

Message édité par Lan Wezel le 28-10-2010 à 18:15:34

el muchacho

Comfortably Numb

Ton problème est le sujet de cette discussion: http://www.wilmott.com/messageview [...] TARTPAGE=1
En page 3, sont mentionnées des impléms particulières comme http://trng.berlios.de/

Message édité par el muchacho le 30-10-2010 à 00:11:11

---------------
Les aéroports où il fait bon attendre, voila un topic qu'il est bien

FORUM HardWare.fr

Programmation

C++

RNG et multithreading

Sujets relatifs
librairie multithreading	[J2EE][Multithreading] séparer les traitements et resynchroniser
UDP/TCP multithreading/sockets asynchrones	Programme Multithreading
Multithreading pour les nuls	select() + multithreading
[Java] Multithreading et gestion des exceptions	exceptions et multithreading ?
MultiThreading C/C++ + VBS	[VB] Multithreading ?
Plus de sujets relatifs à : RNG et multithreading

Page générée en 0.043 secondes