Discussion :

[hammah]

Bonjour,

Je suis débutant dans en Fortran, et j'aimerais savoir si une fonction en Fortran 77 peut rendre un tableau?

Mon cas est le suivant: j'ai une fonction qui fait le produit Matrice (M) Vecteur (X), de dimension (dim), et me rend un vecteur MX. Lors de la compilation j'ai l'erreur suivante:

FUNCTION MX(M,X, dim)
...............1................

Si j'ai bien compris le compilateur n'aime pas trop le nom de ma fonction... il y a quelquechose qui cloche. J'ai tout essayé (changé de nom, vérification des types ...) mais je ne trouve pas. Mon code:

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
      FUNCTION MX(M,X, dim)
      REAL*8 MX
      INTEGER dim
      REAL*8 M(dim, dim), X(dim)
      INTEGER i,j
 
      DO i =1,dim
      MX(i)=0.d0
        DO j=1, dim
          MX(i)= MX(i) + M(i,j)*X(j)
        ENDDO
      ENDDO
 
      RETURN
      END

J'ai donc conclu qu'une fonction ne peut pas rendre un tableau. J'ai déjà fais du C, et la notion de pointeur m'est venu à l'esprit...

[Ehouarn]

Bonjour,

Si tu comptes en rester au Fortran 77 (pourquoi?), tu devrais utiliser une subroutine (en Fortran, contrairement au C, les routines travaillent sur les variables et non pas des copies):

Code :

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

1
2
3
4
5
6
7
8
9
 
      subroutine produitMV(M,X, MX,dim)
      implicit none
      INTEGER dim
      REAL*8 M(dim, dim), X(dim),MX(dim)
c
... tes calculs ici
c
      end

Bonne continuation.

[Ladgalen]

Salut

Petite remarque sur le 77 et le 90 : si tu passes au 90 il existe une fonction toute faite qui fait un produit de matrice :

Code :

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

VEC(:) = matmul( A(:,:) , U(:) )

Plus besoin de la programmer et en plus dans un même programme de fortran 90 tu peux très bien avoir des bout en 77 ça ne gène pas la compilation.

Ceci dis en fortran tu as 2 sous programmes :
SUBROUTINE
FUNCTION

La fonction sert à renvoyer une valeur (comme cosinus par exemple) stockée dans le nom de la fonction. La subroutine fait tout le reste. Elle renvoit rien ou plusieurs valeurs (dans ses arguments). C'est pour ça que ça marche pas, même si ton vecteur est une seule variable tu peux pas déclarer une fonction avec une dimension !

[FR119492]

Salut !

rester au Fortran 77 (pourquoi?)

Pourquoi pas ?

Jean-Marc Blanc

[feynman]

Mon code:

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
      FUNCTION MX(M,X, dim)
      REAL*8 MX
      INTEGER dim
      REAL*8 M(dim, dim), X(dim)
      INTEGER i,j
 
      DO i =1,dim
      MX(i)=0.d0
        DO j=1, dim
          MX(i)= MX(i) + M(i,j)*X(j)
        ENDDO
      ENDDO
 
      RETURN
      END

Il ne faut pas mettre le nom de ton tableau analogue au nom de ta fonction, tu appelle ta fonction autrement par exemple Produit(M,X,dim) et ensuite tu ecrit :

Code :

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

    MX=  Produit(M,X,dim)

A mon avis les subroutines sont mieux

[Sylvain Bergeron]

Bonjour,

Je suis débutant dans en Fortran, et j'aimerais savoir si une fonction en Fortran 77 peut rendre un tableau?

Non. Comme d'autre l'on dit, si tu dois te limiter à du F77, tu dois utiliser une subroutine.

Mon cas est le suivant: j'ai une fonction qui fait le produit Matrice (M) Vecteur (X), de dimension (dim), et me rend un vecteur MX. Lors de la compilation j'ai l'erreur suivante:

FUNCTION MX(M,X, dim)
...............1................

Si j'ai bien compris le compilateur n'aime pas trop le nom de ma fonction... il y a quelquechose qui cloche. J'ai tout essayé (changé de nom, vérification des types ...) mais je ne trouve pas. Mon code:

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
      FUNCTION MX(M,X, dim)
      REAL*8 MX
      INTEGER dim
      REAL*8 M(dim, dim), X(dim)
      INTEGER i,j
 
      DO i =1,dim
      MX(i)=0.d0
        DO j=1, dim
          MX(i)= MX(i) + M(i,j)*X(j)
        ENDDO
      ENDDO
 
      RETURN
      END

J'ai donc conclu qu'une fonction ne peut pas rendre un tableau. J'ai déjà fais du C, et la notion de pointeur m'est venu à l'esprit...

Comme dit par d'autre, cette fonction est native en F90. Toujours en F90, tu peux cependant l'écrire toi-même. Elle aura la forme :

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
 
function MX(M,X)
   real M(:,:)
   real X(:)
   real MX(size(X))
 
      INTEGER dim
      INTEGER i,j
 
      dim = size(X)
      DO i =1,dim
      MX(i)=0.d0
        DO j=1, dim
          MX(i)= MX(i) + M(i,j)*X(j)
        ENDDO
      ENDDO
 
      RETURN
      END

[hammah]

Merci pour les réponses.

En fait, je commence avec le Fortran 77, et une fois que j'ai compris comment tout ça marche, je passe au Fortran 90.

J'ai finalement utilisé une subroutine, et ça marche. je croyais avoir tout compris sur les subroutines et les fonctions, mais ce n'était évidemment pas le cas... ça l'est maintenant.

Encore merci!

le code qui marche:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
SUBROUTINE MX(M,X,Y, dim)
      INTEGER dim
      REAL*8 M(dim, dim), X(dim), Y(dim)
      INTEGER i,j
 
      DO i =1,dim
      Y(i)=0.d0
        DO j=1, dim
          Y(i)= Y(i) + M(i,j)*X(j)
        ENDDO
      ENDDO
 
      RETURN
      END

[Sylvain Bergeron]

...je commence avec le Fortran 77, et une fois que j'ai compris comment tout ça marche, je passe au Fortran 90...

Il est à mon avis beaucoup plus facile de commencer en F90 avec des fonctions et routines « contains » à cause des interfaces automatiques qui simplifient beaucoup beaucoup le débogage.

Ensuite, on peut passer à F77.

[Ehouarn]

Je pense de même qu'il vaut mieux commencer par du Fortran 90 puis s'intéresser au "sous-ensemble" qu'est le Fortran 77; ou éventuellement attaquer les deux simultanément.
Si le but final est de faire du Fortran 90, je te conseille de ne pas trop t'attarder sur le 77.

rester au Fortran 77 (pourquoi?)

Pourquoi pas ?

Mouais, et pourquoi pas du Fortran 66 (voir du Fortran IV) ?

[FR119492]

Salut !

pourquoi pas du Fortran 66

T'en fais pas, Ehouarn, par ordre chronologique, le Fortran II n'a été que le troisième langage avec lequel j'ai travaillé (READ INPUT FILE, etc.).

Pour répondre plus directement à la question, parce qu'il existe une masse de bibliothèques de calcul écrites en Fortran 77 et qui ont l'avantage d'être anciennes, donc nettoyées de tous les bogues qu'elles pouvaient contenir initialement.

Une idée serait que Developpez.com initie une campagne de tests comparatifs pour la vitesse d'exécution de divers algorithmes de calcul numérique (systèmes linéaires, équations différentielles ordinaires, etc.) programmés en différents langages. Personnellement, je serais partant, et on verra ce qu'on verra.

Jean-Marc Blanc

[Ehouarn]

Ta réponse ne tient pas Jean-Marc; un programme Fortran 90 peut faire appel à des librairies Fortan 77. En déduire qu'il faut donc travailler en Fortran 77 relève du mauvais syllogisme.

Pour ce qui est d'une campagne comparative; c'est à mon avis simplement impossible et donc bien futile. Comment définir la 'vitesse' d'un programme d'un langage donnée alors que tant dépend du compilateur et de la plateforme sur laquelle on tourne ?

Ehouarn

[FR119492]

Salut Ehouarn !

'vitesse' d'un programme d'un langage donnée alors que tant dépend du compilateur et de la plateforme

Il est bien clair que des tests effectués sur ces machines différentes ne sont pas comparables. Il serait indispensable de transférer soit tous les sources, soit tous les .exe sur la même machine.

un programme Fortran 90 peut faire appel à des librairies Fortan 77

Pratiquement, c'est ce que je fais, en particulier pour utiliser le format EN. L'allocation dynamique constitue une autre raison. Mais, pour des raisons d'habitude, je préfère le format fixe.

Jean-Marc Blanc

[MABROUKI]

Voici un exemple .Tourne sur compilteur Silverfrost Fnt95 et Compaq Visual Fortran 6.6.
Ceci pour montrer que Fortran permet autant que le C ou VB.
Bon courage pour les "Fortraners".
De la à dire qu' il faut utiliser un "subroutine"


PROGRAM MAIN
REAL, DIMENSION(2) :: V, W
REAL, DIMENSION(2,2) :: A

REAL, DIMENSION(2) :: T, U
REAL, DIMENSION(2,2) :: B

!exemple 1-remplit V & A
W(1)=0.0
W(2)=0.0
A(1,1)=1.
A(1,2)=2.
A(2,1)=3.
A(2,2)=4.
V(1)=1.
V(2)=2.
!exemple 2-remplit T & B
B(1,1)=2.
B(1,2)=4.
B(2,1)=6.
B(2,2)=8.
T(1)=0.0
T(2)=0.0
U(1)=1.
U(2)=2.

! print *, W ! imprime colonnes
W= MatrixVector(A, V)
DO J=1,2
PRINT *, "W = ", W(J)
END DO
PAUSE
! print *, T ! imprime colonnes
T = MatrixVector(B, U)
DO J=1,2
PRINT *, "T = ", T(J)
END DO
PAUSE
CONTAINS

FUNCTION MatrixVector(M, X)
REAL, DIMENSION(:, :: M
REAL, DIMENSION( :: X
INTEGER ::N,I,J
REAL :: MatrixVector(SIZE(X))

N=SIZE(X)
DO I=1, N
MatrixVector(I)=0.
DO J=1, N
MatrixVector(I) = MatrixVector(I) + M(I,J)*X(J)
ENDDO
ENDDO
! print *,MatrixVector ! imprime colonnes

DO J=1,2
PRINT *, "Function MatrixVector = ", MatrixVector(J)
END DO
PAUSE
RETURN
END FUNCTION MatrixVector

END PROGRAM MAIN