Discussion :

[feynman]

Bonjour je veux declarer une matrice (2x2) qui depend d'une variable xi,
j'ai ecrit le code suivant, apres compilation le compilateur me dit qu'on ne peut pas redimentionner la fonction, pouvez vous m'aider svp?
merci
Voila 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
16
17
18
 
      function A(xi)
      double precision xi, A(2,2)
      double precision u,u0,l,ky,lb,delta,del,e
      del=5.0d0
	  U0=50.0d0
      L=200.0d0
      ky=0.03d0
      lb=4.5d-4
	  e=6.8d-2
	  delta=0.0d0
      U=(U0/2.0d0)*(tanh(((-xi*L)-(L/2.0d0)/del)+tanh(((xi*L)-(L/2.0d0))/del)+2.0d0)
      A(1,1)=-L*(ky-(1.0d0/(lb**2))*L*xi)
	  A(1,2)=L*(e-U+delta)
	  A(2,1)=-L*(e-U-delta)
	  A(2,2)=L*(ky-(1.0d0/(lb**2))*L*xi)
      return
      end

[Ehouarn]

Bonjour,

Une fonction ne renvoie qu'une seule valeur.
Pour en récupérer plus, il faut passer par une subroutine.

Bonne continuation.

[Sylvain Bergeron]

Sous réserve d'une erreur de parenthèse à la ligne commençant par "U=(U0/2.0d0)*", ton code fonctionne de mon côté.

La seule chose que j'ai fait est d'insérer la fonction dans une section "contains" dans le programme appelant. Ton problème est donc un problème d'interface. Tes pistes de solution sont donc :

• Section contains comme j'ai fait
• Bloc interface
• Module

Voici mon code fonctionnel :

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
program test
   implicit none
   real (8)  Xi
 
   Xi = 1.0
   write (*,*) A(Xi)
contains
      function A(xi)
      double precision xi, A(2,2)
      double precision u,u0,l,ky,lb,delta,del,e
      del=5.0d0
	  U0=50.0d0
      L=200.0d0
      ky=0.03d0
      lb=4.5d-4
	  e=6.8d-2
	  delta=0.0d0
      U=(U0/2.0d0)*tanh(((-xi*L)-(L/2.0d0)/del)+tanh(((xi*L)-(L/2.0d0))/del)+2.0d0)
      A(1,1)=-L*(ky-(1.0d0/(lb**2))*L*xi)
	  A(1,2)=L*(e-U+delta)
	  A(2,1)=-L*(e-U-delta)
	  A(2,2)=L*(ky-(1.0d0/(lb**2))*L*xi)
      return
      end function
end program

Commentaire éditorial pour terminer : Ta fonction ne contient pas de « implicit none », ce qui est une mauvaise habitude de travail ...