Quel algorithme pour rand ?

Bonjour,

Je m'intéresse au fonctionnement interne de rand, et j'aimerais savoir quel algorithme est utilisé dans le standard C (et si la norme ne dit rien, qu'en est-il en pratique avec gcc par exemple ?).

J'ai déjà fait quelques recherches, et suis tombé sur les congruences linéaires et Mersenne Twister, mais pas moyen de savoir ce qui est utilisé dans rand.

Quelqu'un peut éclairer ma lanterne ?

Merci,
Sébastien

Parce que ce qui est utilisé dans rand est spécifique à chaque implémentation. La norme propose (pas impose) cependant un algorithme portable :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

/* #define RAND_MAX 32767 */   static unsigned long int next = 1;   void srand(unsigned int seed) { next = seed; }   int rand(void) { next = next * 1103515245 + 12345; return (unsigned int)(next/65536) % 32768; }

---

Et de façon générale, il ne faut JAMAIS présupposer de l'algorithme de rand, notamment pour une pseudo-cryptographie. Si tu as besoin d'un générateur aléatoire particulier, tu DOIS l'implémenter toi-même et ne pas compter sur rand pour ça, même si l'implémentation que tu utilises est "la bonne". En effet, cela peut changer à tout moment sans préavis suivant les patchs du compilateur, et/ou ne pas être la même fonction d'une plate-forme à l'autre.

Tu peux également aller lire cette discussion, ainsi que celle-ci, à titre d'information.

Citation:

---

Envoyé par Melem

Parce que ce qui est utilisé dans rand est spécifique à chaque implémentation. La norme propose (pas impose) cependant un algorithme portable :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

/* #define RAND_MAX 32767 */   static unsigned long int next = 1;   void srand(unsigned int seed) { next = seed; }   int rand(void) { next = next * 1103515245 + 12345; return (unsigned int)(next/65536) % 32768; }

---

---

Si c'est vraiment comme ça dans la pratique c'est très moche...

Très moche tu veux dire ...

Tu trouveras un certain nombre d'information sur des algorithmes de génération de nombres pseudo-aléatoire au chapitre 7 de Numerical Recipes. Cette question aurait plus sa place dans le forum algo que C.

Meilleures salutations

Thierry