rand() pour ifort

Bonjour,
j'ai récupéré un projet en fortran 77, que je compile avec gfortran (linux). Je voudrais (par curiosité intellectuelle et aussi car le débogueur est parait-il meilleur) le compiler avec ifort.
Mais voilà qu'avec ifort j'ai une erreur à cause d'une fonction rand() :

Code:

---

1 2	utils/UTR.F(62): error #6404: This name does not have a type, and must have an explicit type. [RAND] x1 = rand(s)

---

J'ai lu que le fortran 77 n'avait pas de standard en ce qui concerne la génération de pseudo-aléatoire, je pense donc que ça vient de là : confirmez-vous et que me conseilleriez-vous ?
Merci !

ps : voici le code de la fonction

Code:

---

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

DOUBLE PRECISION FUNCTION rn() IMPLICIT NONE DOUBLE PRECISION x1, x2, PI2, EPS PARAMETER ( PI2 = 6.283185307d0, EPS=1.d-15 ) INTEGER init_seed_linux, s EXTERNAL init_seed_linux INTRINSIC SQRT, LOG, COS   s = init_seed() x1 = rand(s) DO WHILE (x1 .LT. EPS) x1 = rand() ENDDO x2 = rand() rn = SQRT(-2.0*LOG(x1))*COS(PI2*x2) RETURN END

---

Bonjour

Citation:

---

J'ai lu que le fortran 77 n'avait pas de standard en ce qui concerne la génération de pseudo-aléatoire,

---

Exact. Ceci dit, certains éditeurs (Lahey Fortran par exemple, de mémoire) en proposent un.

Citation:

---

que me conseilleriez-vous ?

---

D'en écrire un. Voici par exemple un petit tout simple que j'ai utilisé il ya quelques années.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

FUNCTION RANDOM (N0)   C- On initialise N0 a une valeur quelconque ENTIERE C-RANDOM renvoie un REEL compris entre 0.0 et 1.0   REAL RANDOM INTEGER N0, K, CONST1, CONST2, CONST3, CONST4 PARAMETER (CONST1 = 127773) PARAMETER (CONST2 = 16807) PARAMETER (CONST3 = 2836) PARAMETER (CONST4 = 2147483647)   K = N0 / CONST1 N0 = CONST2 * (N0 - K * CONST1) - K * CONST3 IF (N0 .LT. 0) N0 = N0 + CONST4 RANDOM = REAL(N0) / REAL(CONST4)   END

---

Tu peux en trouver de plus complexes sur http://www.fortran.com/ dans la rubrique Tools, Libraries

Merci !

Bonjour

je viens de découvrir - par hasard - une autre solution.

Le compilateur g95 http://www.g95.org
qui tourne sous Windows, Mac OS X, Linux, Solaris, propose RAND dans ces fonctions intrinsèques supplémentaires.