Discussion :

[feynman]

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 :

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

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

Et voila ma fonction :

Code :

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

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

[François L.]

Je rajouterai bien un:

Code :

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

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 :

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

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.

[feynman]

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

Code :

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

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 :

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

[DaZumba]

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 :

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

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

[feynman]

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!

[feynman]

Malgre cela, ca fait 2 ans que j utilise le F77, et j ai fait des centaines de fonctions et ca marchent toujours,!