Discussion :

[salseropom]

Bonjour, dans une fonction que j'ai intégrée à une librairie dynamique, je tire des nombres aléatoires suivant une loi de Weibull (fonction utilisant des nombres aléatoires suivant une loi uniforme sur [0, 1]). Quand je fais

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
for(i = 0 ; i < 10 ; i++)
{
  srand(10);
  d = wblrnd(...);
}

j'obtiens 10 valeurs différentes ! En revanche, si je mets le srand(10) DANS la fonction qui est dans la librairie, j'ai bien les mêmes nombres.

Ma question est : pourquoi est-ce qeu srand() n'est pas connu dans la librairie ?

Merci d'avance

[Emmanuel Delahaye]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
for(i = 0 ; i < 10 ; i++)
{
  srand(10);
  d = wblrnd(...);
}

Non.

http://emmanuel-delahaye.developpez....=Page7#LXXXVII

[salseropom]

Salut Emmanuel,

dans le code que j'ai écrit, je désirais avoir 10 fois le même nombre aléatoire, d'où la présence du srand(10) DANS la boucle for... ceci n'est qu'un test. Mais malgré ce srand(10) dans la boucle for, j'obtiens 10 nombres différents. Je rappelle que wblrnd() est une fonction que j'ai écrite et qui est dans une librairie.

Ma question est donc : pourquoi malgré ce srand(10) DANS la boucle for j'ai 10 valeurs différentes ?

[Emmanuel Delahaye]

Salut Emmanuel,

dans le code que j'ai écrit, je désirais avoir 10 fois le même nombre aléatoire,

Que peut bien signifier 'avoir 10 fois le même nombre aléatoire. Si c'est le même, tu fais un tirage et tu l'utilise 10 fois. Il ne faut pas faire 10 tirages...

d'où la présence du srand(10) DANS la boucle for... ceci n'est qu'un test. Mais malgré ce srand(10) dans la boucle for, j'obtiens 10 nombres différents. Je rappelle que wblrnd() est une fonction que j'ai écrite et qui est dans une librairie.

Ma question est donc : pourquoi malgré ce srand(10) DANS la boucle for j'ai 10 valeurs différentes ?

Alors ça vient de la fonction.

Je commence à comprendre ce que tu veux faire. Tu es sûr que ta fonction est basée sur rand() ?

[diogene]

n1256 :

7.20.2 Pseudo-random sequence generation functions
7.20.2.1 The rand function
Synopsis
1 #include <stdlib.h>
int rand(void);
....
3 The implementation shall behave as if no library function calls the rand function.
....
7.20.2.2 The srand function
Synopsis
1 #include <stdlib.h>
void srand(unsigned int seed);
....
3 The implementation shall behave as if no library function calls the srand function.

[salseropom]

Bonjour à vous deux et merci de vos réponses.

Emmanuel, voici mon code :

mon fichie aleatoire.c (qui fait partie de ma librairie)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
 
/* 
   RAND is a macro which returns a pseudo-random numbers from a uniform
   distribution on the interval [0 1]
 */
#define RAND ((rand())/((double) RAND_MAX))
 
/*******************************************************************************
 * wblrnd : returns a pseudo-random number chosen from the Weibull distribution*
 * with scale parameter lambda > 0, shape parameter k > 0 and localisation     *
 * parameter theta                                                             *
 *******************************************************************************/
EXTERN_C RISOTO_API double __stdcall wblrnd(double k,double lambda,double theta)
{
  double d;
  do
  {
    d = RAND;
  }
  while(d < DBL_EPSILON);
 
  return lambda * pow(-log(d),1./k) + theta;
}

Dans une autre fonction en gros je calcule

wblrnd(k,lambda,theta) ou bien
theta + wblrnd(k,lambda,theta) ou bien
theta - wblrnd(k,lambda,theta)

donc si theta == 0 je suis sensé obtenir les mêmes réponses, ce qui n'est pas le cas... Pourtant j'ai bien forcé à avoir les mêmes nombres avec mon srand(10) ==> mais je suis tout à fait d'accrod avec toi que mettre le srand(10) dans la boucle "for" n'a aucun intérêt, sauf pour mon test...

Si je mets mon srand(10) dans ma fonction wblrnd() qui est dans la librairie, j'ai bien les même valeurs... ce qui me rassure.

Diogène :

que dois-je faire alors ?

Merci encore de vos réponses