[Fortran 77] Fonction complexe

Version imprimable

08/08/2008, 19h16
feynman

[Fortran 77] Fonction complexe
Salut pour tous,
J'ai écrit une fonction qui calcule la fonction suivante : exp(i(kx+\bar(k)\bar(x))) (\bar c'est le conjugué complexe)
Mais quand j'utilise ma fonction, elle ne donne que la partie réelle de la réponse :

Code:

cmplx expr CEXP(-II*K*X-II*CONJG(K)..... truncated to real EXPP

Et voila ma fonction :
Code:

1 2 3 4 5 6 7 function expp(k,x) complex k,x,ii ii=cmplx(0.,1.) expp=cexp(-ii*k*x-ii*conjg(k)*conjg(x)) return end
Alors où est le problème ? J'ai essayé beaucoup, mais ....
Merci
09/08/2008, 21h51
François L.
Je rajouterai bien un:
Code:

1 2 COMPLEX expp
dans la déclaration des variables de la fonction.
Sinon, elle est déclarée implicitement comme un REAL.

(ou bien:
Code:

1 2 COMPLEX FUNCTION expp(...)
)

Entre parenthèses, mettre IMPLICIT NONE au début de toutes fonctions, sous-routines ou programme principal est une saine habitude qui t'aurait évité ce déboire. :P

on a utiliser comme tu as dit, mais toujours pas resolu , j ai eu ;

Code:

1
2
3
4
5
6
 
 C:\Documents and Settings\Toufic ARWADI\Desktop\eeee.f: In program `h':
C:\Documents and Settings\Toufic ARWADI\Desktop\eeee.f:5: 
         M=expp(II,ii)
           ^
Invalid declaration of or reference to symbol `expp' at (^) [initially seen at (^)]

et vvoila mon programme :

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
 
      program h
      implicit none
      complex ii,M
      ii=cmplx(0.,1.)
      M=expp(II,ii)
      print*,'M=',M
      end
 
 
 
 
 
 
 
 
 
        complex function expp(k,x)
        implicit none
        complex k,x,ii
        ii=cmplx(0.,1.)
        expp=cexp(-ii*k*x-ii*conjg(k)*conjg(x))
        return
        end

alors ou est tjs le pb?
merci

11/08/2008, 13h46
DaZumba
Citation:

Envoyé par feynman

alors ou est tjs le pb?

Comme d'habitude pour les fonctions, il faut declarer expp() comme etant de type complex la ou elle est utilisee. Donc:
Code:

1 2 3 4 5 6 7 8 9 program h implicit none complex ii,M complex expp ii=cmplx(0.,1.) M=expp(II,ii) print*,'M=',M end
11/08/2008, 13h53
feynman

Citation:

Envoyé par DaZumba

Comme d'habitude pour les fonctions, il faut declarer expp() comme etant de type complex la ou elle est utilisee. Donc:

Merde!! Merci bcp DaZumba, je ne savais pas ca!
11/08/2008, 14h00
feynman

Malgre cela, ca fait 2 ans que j utilise le F77, et j ai fait des centaines de fonctions et ca marchent toujours,!
11/08/2008, 17h53
François L.

Tu as sans doute toujours eu la chance de tomber sur des types implicites qui convenaient: si tes noms de fonctions commencent par a -> h ou o->z, elles sont implicitement définies comme des REAL.