Discussion :

[Ladgalen]

Bonjour

D'après ce post il semble qu'en fortran 90 il soit possible d'écrire une function qui renvoi un vecteur. J'ai écris cette fonction qui me renvoit une erreur que je n'arrive pas à comprendre :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
program test
  implicit none
  integer :: k
  double precision, dimension(3) :: vect_product
  double precision, dimension(3) :: u,v,w
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
end program test
 
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3)             :: w
    double precision                           :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product

Voici l'erreur :

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
 In file test.f90:11
 
  w(:) = vect_product(u(:),v(:))
                     1
Warning: Extension: REAL array index at (1)
 In file test.f90:11
 
  w(:) = vect_product(u(:),v(:))
                          1
Warning: Extension: REAL array index at (1)
 In file test.f90:11
 
  w(:) = vect_product(u(:),v(:))
                    1
Error: Rank mismatch in array reference at (1) (2/1)

Merci pour votre aide

[FR119492]

Salut!
Je ne possède pas de copie de la norme Fortran 90 ou 95. En revanche, le manuel de référence de mon Compaq Fortran affirme d'une part qu'il est compatible Fortran 90 ou 95 et d'autre part qu'une fonction ne peut retourner qu'une seule valeur.
Jean-Marc Blanc

[Sylvain Bergeron]

Il n'y a pas de problème à retourner un vecteur en F90. Le problème de ton programme est simplement qu'il ne voie pas la fonction. Tu as besoin d'une interface.

Je n'ai pas le temps de tester, mais le programme suivant devrait fonctionner :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
 
program test
  implicit none
  integer :: k
!  double precision, dimension(3) :: vect_product
  double precision, dimension(3) :: u,v,w
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
!end program test
contains
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3)             :: w
    double precision                           :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product
end program test

Évidemment que tu peux alternativement utiliser un bloc interface ou un module plutôt que de placer la fonction en « contains » de ton programme.

[Ladgalen]

Il n'y a pas de problème à retourner un vecteur en F90. Le problème de ton programme est simplement qu'il ne voie pas la fonction. Tu as besoin d'une interface.

Je n'ai pas le temps de tester, mais le programme suivant devrait fonctionner :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
 
program test
  implicit none
  integer :: k
!  double precision, dimension(3) :: vect_product
  double precision, dimension(3) :: u,v,w
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
!end program test
contains
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3)             :: w
    double precision                           :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product
end program test

Évidemment que tu peux alternativement utiliser un bloc interface ou un module plutôt que de placer la fonction en « contains » de ton programme.

Bonjour

Ton code compile et fonctionne correctement. J'ai rajouté une interface (je n'y connais pas grand chose en interface). Le code suivant compile mais plante à l'éxécution :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
program test
  implicit none
  integer :: k
  double precision, dimension(3) :: u,v,w
 
  interface
    function vect_product(u,v)
      double precision,dimension(3),intent(in) :: u,v
    end function vect_product
  end interface
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
 
end program test
 
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
  function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
    implicit none
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3)             :: w
    double precision                           :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product

Le mesage d'erreur est le suivant :
Fortran runtime error: Array bound mismatch for dimension 1 of array 'vect_product' (in file 'test.f90', around line 26)

J'ai compilé avec gfortran.

En dehors de ce problème je ne comprend pas pourquoi une interface est nécessaire. J'ai l'habitude, pour des petits codes, de mettre toutes les fonctions et subroutines dans un seul fichier les unes après les autres, sans utiliser de module ou d'interface et ça ne pose pas de problème.

Pourquoi ce problème apparait quand j'ai une fonction qui renvoie un vecteur !

Merci

[Sylvain Bergeron]

Bonjour

Ton code compile est fonctionne correctement. J'ai rajouté une interface (je n'y connais pas grand chose en interface). Le code suivant compile mais plante à l'éxécution :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
program test
  implicit none
  integer :: k
  double precision, dimension(3) :: u,v,w
 
  interface
    function vect_product(u,v)
      double precision,dimension(3),intent(in) :: u,v
    end function vect_product
  end interface
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
 
end program test
 
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
  function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
    implicit none
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3)             :: w
    double precision                           :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product

Le mesage d'erreur est le suivant :
Fortran runtime error: Array bound mismatch for dimension 1 of array 'vect_product' (in file 'test.f90', around line 26)

J'ai compilé avec gfortran.

Il manque la ligne :

double precision, dimension(3) :: vect_product

dans l'interface. C'est tout.

En dehors de ce problème je ne comprend pas pourquoi une interaface est nécessaire. J'ai l'habitude, pour des petits codes, de mettre toutes les fonctions et subroutines dans un seul fichier les unes après les autres, sans utiliser de module ou d'interface et ça ne pose pas de problème.

Pourquoi ce problème apparait quand j'ai une fonction qui renvoit un vecteur !

Merci

Première règle : le concept de fichier source n'existe pas en Fortran. Que tu mettes toutes tes routines dans le même fichier, dans n'importe quel ordre, ou dans des fichiers distincts, n'a aucun effet sur la visibilité que le compilateur a des différentes composantes : le compilateur compile un sous-programme à la fois de façon indépendante.

Deuxième règle : en absence d'interface, le compilateur ne peut que faire des déductions sur le protocole à utiliser pour invoquer la fonction :

• Le nom de la fonction (si c'est une fonction) est t'il déclaré ? Sinon, y a t'il une règle implicite ? implicite none ?
• Du point précédent, le compilateur « connait » maintenant le type de résultat, qui ne peut être que scalaire, puisqu'il n'y a pas moyen de déduire que le résultat est vectoriel.
• Reste maintenant à déduire le type des arguments. Seule source d'info disponible : les arguments utilisés.
• Si tu passes un vecteur, l'argument dummy correspondant ne peut être que « assumed-size », soit (*), puisqu'il n'y a pas moyen de savoir si l'argument est déclaré « assumed-shape », soit (:).

En gros, lorsqu'il n'y a pas d'interface, tu es limité à du protocole F77.

De ces déductions, rien ne garantie qu'elles seront justes. Il est donc possible qu'il y ait plantage à l'exécution...

[Ladgalen]

Merci, effectivement ça compile et ça marche comme ça :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
program test
  implicit none
  integer :: k
  double precision, dimension(3) :: u,v,w
 
  interface
    function vect_product(u,v)
      double precision, dimension(3)           :: vect_product
      double precision,dimension(3),intent(in) :: u,v
    end function vect_product
  end interface
 
  u(:) = (/ 1.d0, 0.d0, 0.d0 /)
  v(:) = (/ 0.d0, 1.d0, 0.d0 /)
  w(:) = (/ 0.d0, 0.d0, 0.d0 /)
 
  w(:) = vect_product(u(:),v(:))
 
  write(*,"(3F5.2)")(w(k),k=1,3)
 
 
end program test
 
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
  function vect_product( u, v)
    ! calcule le produit vectoriel ( u x v )
    ! retourne un vecteur norme
    implicit none
 
    ! la fonction :
    double precision, dimension(3)             :: vect_product
    ! arguments :
    double precision, dimension(3), intent(in) :: u, v
    ! locales :
    double precision, dimension(3) :: w
    double precision               :: norme
 
    w(1) = u(2)*v(3) - u(3)*v(2)
    w(2) = u(3)*v(1) - u(1)*v(3)
    w(3) = u(1)*v(2) - u(2)*v(1)
 
    norme = sqrt( sum( (w(:))**2 ) )
 
    if( norme > 1.d-15 ) then
      vect_product(:) = w(:) / norme
    else
      vect_product(:) = 0.d0
      write(*,"('ATTENTION : norme du produit vectoriel nulle')")
    end if
 
    return
  end function vect_product

D'après ce que tu me dis l'interface est un équivalent du prototype en C qui permet au compilateur de vérifier les types et dimensions des arguments lors de l'appel ?

En conclusion un module devrait donc ressembler à (dis moi si je me trompe) :

Code :

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

1
2
3
4
5
6
7
8
 
module mon_module
  ! variables communes au fonctions et routine du modules
 
  contains
  ! les fonctions et subroutine du modules
 
end module mon_module

Edit : En fait dans les modules l'interface ne sert à rien (et donne une erreur de compilation). Quand on place un sous programme dans un module, et que ce module est inclu dans les programmes appelant ces sous programme, l'interface est crée automatiquement. Une interface n'est nécessaire (conseillé en fait, obligatoire uniquement dans certains cas) que s'il s'agit d'un sous programme externe.

[Ladgalen]

Je pensais naivement qu'un vecteur ou une matrice pouvait entrer dans la définition "ne renvoit qu'une valeur". En fait il faudrait préciser "ne renvoi qu'un scalaire".

Il semble qu'en fortran 90 il y ait un moyen de bidouillé en définissant une structure qui contient trois réel et déclarer une fonction qui renvois cette strucutre ... je n'ai pas testé.

De toute façon c'est pas bien grave je ferai une subroutine à la place.