Version imprimable
Bonjour,
quelqu'un pourrait-il m'expliquer clairement comment on utilisa la fonction RAN et SECNDS en fortran?
En gros je dois l'utiliser dans un programme pour lancer deux des.. Donc la valeur qu'elel doit me retourner doit etre comprise entre 1 et 6..
Merci d'avance et merci aux admin qui vont tourner ce site, il est vraiment très bien fait!
Salut
En fortran 90 tu as
random_number(plop)
plop est un real (je crois que tu peux mettre un tableau 1D)
Je sais pas si c'est cette fonction que tu cherches ?
Pour avoir un nombre entre 1 et 6, sachant que plop est entre 0 et 1, cette instruction te donne un entier entre 1 et 6 (nint c'est l'entier le plus proche) :
Ceci dit dans random_number tu ne maitrise rien et notamment tu ne connais pas la graine (c'est un entier qui permet d'initialiser le tirage des nombres aléatoires). Je te donne une fonction issue de numerical recipes qui donne des nombres aléatoires. La graine c'est idum ! il faut que idum soit un entier négatif assez grand, du genre idum = - 35463364 (d'après mon clavier).Code:
2
3
En premiere ligne de ton programme tu initialises idum et TU NE CHANGES PLUS SA VALEUR !!!!
idum = - 35463364
ensuite dans ton programmes tu mets
val_des = nint( 5.d0 * ran2(idum) + 1.d0)
et tu as une valeur entre 1 et 6 !
La fonction est en F77 mais normalement elle devrait marcher si tu compiles tu F90. Tu peux la convertir en F90 il n'y a pas grand chose à changer.
Problème pour toi : à une graine correspond une série de nombre aléatoire donc si tu ne changes pas la graine les dés donneront toujours les mêmes valeurs dans le même ordre chaque fois que tu relanceras ton programme.Code:
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35FUNCTION ran2(idum) INTEGER idum,IM1,IM2,IMM1,IA1,IA2,IQ1,IQ2,IR1,IR2,NTAB,NDIV REAL ran2,AM,EPS,RNMX PARAMETER (IM1=2147483563,IM2=2147483399,AM=1./IM1,IMM1=IM1-1, *IA1=40014,IA2=40692,IQ1=53668,IQ2=52774,IR1=12211,IR2=3791, *NTAB=32,NDIV=1+IMM1/NTAB,EPS=1.2e-7,RNMX=1.-EPS) INTEGER idum2,j,k,iv(NTAB),iy SAVE iv,iy,idum2 DATA idum2/123456789/, iv/NTAB*0/, iy/0/ if (idum.le.0) then idum=max(-idum,1) idum2=idum do 11 j=NTAB+8,1,-1 k=idum/IQ1 idum=IA1*(idum-k*IQ1)-k*IR1 if (idum.lt.0) idum=idum+IM1 if (j.le.NTAB) iv(j)=idum 11 continue iy=iv(1) endif k=idum/IQ1 idum=IA1*(idum-k*IQ1)-k*IR1 if (idum.lt.0) idum=idum+IM1 k=idum2/IQ2 idum2=IA2*(idum2-k*IQ2)-k*IR2 if (idum2.lt.0) idum2=idum2+IM2 j=1+iy/NDIV iy=iv(j)-idum2 iv(j)=idum if(iy.lt.1)iy=iy+IMM1 ran2=min(AM*iy,RNMX) return END C (C) Copr. 1986-92 Numerical Recipes Software 0?#<
La fonction ci dessous donne en premier argument un entier généralement grand qui peut être utilisé comme graine. Il donne un entier lié au bip de l'horloge de l'ordinateur.
system_clock(count,count_rate,count_max)
les trois variables sont des sorties de la subroutine qui sont des integer. Le premier est la valeur de l'horloge quand tu appelles la subroutine(que tu peux prendre comme graine), le deuxième c'est la fréquence de l'horloge et le troisème est la valeur max de l'horloge.
Conseil :
Pour débeuguer le programme il est souvent préférable d'avoir une graine fixe. De cette manière la série de nombre aléatoire est toujours la même et tu peux donc comparer les résultats issus de deux exécutions successives.
Je ne connais pas SECNDS !