fonction et matrice?

Version imprimable

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:

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

08/04/2011, 07h54
Ehouarn

Bonjour,

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

Bonne continuation.

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:

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 ...