Discussion :

[shadock0]

Bonjour,

Alors voici mon problème. Je définis une interface pour une fonction à l'intérieur d'un module ainsi:

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
 
 
module utils
 
  use blablabla1
  use blablabla2
 
  implicit none
 
  contains
 
  [plusieurs fonctions et sous-routines]
 
interface
 
   function convert2matrX(tab)
     real,dimension(:,:,:,:,16):: tab
     real,dimension(:,:,:,:,4,4):: conv_scat
 
     [... instruction 1 ...]
 
   end function convert2matrX
 
   function convert2matrX(tab)
     real,dimension(:,:,7):: tab
     real,dimension(:,:,4,4):: conv_ext
 
     [...instructions 2...]
 
   end function convert2matrX
 
end interface

Cette fonction transforme une liste de 7 ou ben 16 nombre en une matrice 4x4. Bien sur les instructions ne sont pas les même selon la taille de la liste, j'ai donc opté pour cette solution d'une interface avec 2 fonctions à l'intérieur.

Cependant, il semblerait que g95 n'aime pas beaucoup car j'ai le droit à:

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 utils.f90:573
 
interface
         1
Error: Unexpected INTERFACE statement in CONTAINS section at (1)
In file utils.f90:579
 
     real,dimension(:,:,:,:,16):: tab
                              1
Error: Bad specification for deferred shape array at (1)
In file utils.f90:580
 
     real,dimension(:,:,:,:,4,4):: conv_scat
                               1

Pourquoi ce problème? L,interface est-elle mal déclarée? Me conseillez-vous une autre solution (par exemple un case dans la fonction) qui serait plus rapide à l'exécution (car cette fonction risque d'être utilisée de nombreuses fois) ?

merci beaucoup pour votre précieuse aide

[shadock0]

Alors pour le problème de "Error: Bad specification for deferred shape array", j'ai trouvé la solution.

je ne spécifie pas le tableau et je le met en allocatable, puis je l'alloue dans la fonction ainsi:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
   function convert2matrX(tab)
     implicit none
     real,dimension(:,:,:,:,16),intent(in):: tab
     real,allocatable,dimension(:,:,:,:,:,:):: convert2matrX
     integer::sz
 
     allocate(convert2matrX(size(tab(:,1,1,1,1)),size(tab(1,:,1,1,1)),&
          size(tab(1,1,:,1,1)),size(tab(1,1,1,:,1)),4,4))
 
     ... instructions ...
 
   end function convert2matrX

Cependant le problème avec la structure "interface" demeure... Je ne sais pas comment faire pour que le compilateur en veuille...

[shadock0]

J'ai également essayé de bouger tout le module interface plus haut dans le module, avant CONTAINS.

Désormais j'ai l'erreur suivant lorsque j'essaie de récupérées les données des différents tableaux:

Code :

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

1
2
3
4
 
       conv16(:,:,:,:,1,:)=tab(:,:,:,:,1:4)
                                           1
Error: Unexpected assignment statement at (1) in INTERFACE body

où conv16 est la fonction chargée les listes de 16 chiffres et donc qui renvoie un tableau (:,:,:,:,4,4) et tab, est la liste de valeurs en entrée...

faut-il que j'utilise plutôt une subroutine ? d'où vient le pb ?

[bobbyboy]

Bonjour,

Si tu veux utiliser une interface regroupant tes deux fonctions, ces deux fonctions doivent avoir un nom différent dans l'interface. La syntaxe suivante devrai mieux marcher :

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
 
module utils
 
   use blabla1
   use blabla 2
 
   implicit none
 
!  Déclaration de l'interface avec un nom générique
!  ------------------------------------------------
   interface convert2matriX
 
      function convert2matriX_num1 (...)
         ! Déclarations et instructions
      end function convert2matriX
 
      function convert2matriX_num2 (...)
         ! Déclarations et instructions
      end function convert2matriX
 
   end interface convert2matriX
 
end module utils

Puis dans ton programme principal, tu utilise le nom de l'interface comme une fonction :

Code :

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

1
2
3
4
5
6
7
 
program main
   use utils
   implicit none
 
   convert2matriX(...)
end

[shadock0]

@bobbyboy

Non ce n'est pas ça. j'y aie bien sûr pensé et j,ai modifié le code en conséquence mais l'erreur est toujours la même...

Peut-être cela vient-il du fait que des fonctions ne peuvent pas renvoyer de tabeaux? Je ne vois pas sinon...

J'ai bougé hors de du bloc "CONTAINS", j'ai changé les noms des fonctions, j'ai fait en sorte qu'elles soient parfaitement définies et qu'il n'y aie pas de "ALLOCATE" à l'interrieur (ce qui apparemment n'est pas possible). Bref... je ne sais plus...

[Sylvain Bergeron]

la ligne :

! déclarations et instructions

dans les interfaces, ça contient quoi ?

[shadock0]

Le code se présente comme ça en entier:

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
 
module utils
 
  use machin
  use truc
  use bidule
 
  implicit none
 
!************************************************************
 
  interface convert2matrX
     ! this function transform the list read in 4x4 matrix
 
     function conv16(tab)
       !convert diffusion matrix size: 16
       implicit none
       real,dimension(nphio,ntheo,nphii,nthei,16),intent(in):: tab
       real,dimension(nphio,ntheo,nphii,nthei,4,4):: conv16
       integer::sz
 
       conv16(:,:,:,:,1,:)=tab(:,:,:,:,1:4)
       conv16(:,:,:,:,2,:)=tab(:,:,:,:,5:8)
       conv16(:,:,:,:,3,:)=tab(:,:,:,:,9:12)
       conv16(:,:,:,:,4,:)=tab(:,:,:,:,13:16)
     end function conv16
 
     function conv7(tab)
       !convert extinction matrix size: 7
       implicit none
       real,dimension(nphii,nthei,7),intent(in):: tab
       real,dimension(nphii,nthei,4,4):: conv7
       integer::sz
 
       conv7(:,:,1,:)=tab(:,:,1:4)
       conv7(:,:,2,3:4)=tab(:,:,5:6)
       conv7(:,:,3,4)=tab(:,:,6)
       conv7(:,:,2,2)=tab(:,:,1)
       conv7(:,:,3,3)=tab(:,:,1)
       conv7(:,:,4,4)=tab(:,:,1)
       conv7(:,:,2:4,:)=conv7(:,:,:,2:4)
       conv7(:,:,3,2)=-conv7(:,:,2,3)
       conv7(:,:,4,2:3)=-conv7(:,:,2:3,4)
 
     end function conv7
  end interface
 
  CONTAINS
 
  ... autres fonctions ...
 
end module utils

[Sylvain Bergeron]

Il y a confusion ici.

Il y a plusieurs utilisations à l'instruction interface. Tu peux :

• définir une fonction générique
• définir la fonction à appeler pour résoudre un opérateur ou une assignation
• définir comment appeler une fonction (liste et types des arguments).

La 3ième utilisation est automatiquement générée par le compilateur à la compilation d'un module (l'interface est stockée dans le fichier .mod). Une interface, ça sert à dire au compilateur, lors de la compilation du code client (code utilisant une fonction) comment appeler la fonction. Ce n'est pas le code de la fonction, c'est uniquement comment l'appeler.

Dans ton cas, tu semble vouloir définir une fonction générique (convert2matrX) implantée en conv16 ou conv7 selon que l'argument est de rang 5 ou 3.

Tu dois donc avoir :

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
54
55
56
 
module utils
 
  use machin
  use truc
  use bidule
 
  implicit none
 
  interface convert2matrX
    module procedure conv16
    module procedure conv7
  end interface
 
  CONTAINS
 
     ! this function transform the list read in 4x4 matrix
 
     function conv16(tab)
       !convert diffusion matrix size: 16
       implicit none
       real,dimension(nphio,ntheo,nphii,nthei,16),intent(in):: tab
 
!? d'où viennent nphio,ntheo,nphii,nthei ? D'un des use ?
 
       real,dimension(nphio,ntheo,nphii,nthei,4,4):: conv16
       integer::sz
 
       conv16(:,:,:,:,1,:)=tab(:,:,:,:,1:4)
       conv16(:,:,:,:,2,:)=tab(:,:,:,:,5:8)
       conv16(:,:,:,:,3,:)=tab(:,:,:,:,9:12)
       conv16(:,:,:,:,4,:)=tab(:,:,:,:,13:16)
     end function conv16
 
     function conv7(tab)
       !convert extinction matrix size: 7
       implicit none
       real,dimension(nphii,nthei,7),intent(in):: tab
       real,dimension(nphii,nthei,4,4):: conv7
       integer::sz
 
       conv7(:,:,1,:)=tab(:,:,1:4)
       conv7(:,:,2,3:4)=tab(:,:,5:6)
       conv7(:,:,3,4)=tab(:,:,6)
       conv7(:,:,2,2)=tab(:,:,1)
       conv7(:,:,3,3)=tab(:,:,1)
       conv7(:,:,4,4)=tab(:,:,1)
       conv7(:,:,2:4,:)=conv7(:,:,:,2:4)
       conv7(:,:,3,2)=-conv7(:,:,2,3)
       conv7(:,:,4,2:3)=-conv7(:,:,2:3,4)
 
     end function conv7
 
  ... autres fonctions ...
 
end module utils

[bobbyboy]

Bonjour,

la ligne :

! déclarations et instructions

dans les interfaces, ça contient quoi ?

Désolé c'est une erreur, je pensais qu'on pouvais définir les fonction directement dans le bloc de l'interface, mais non.