Discussion :

[feynman]

Bonjour,
Je veux bien tracer une courbe 1D à une varible donner comme suit :
gamma(x)=2 si x entre -0.5 et 0.5
gamma(x)=1 ailleurs.
J ai taper un petit code sur mon PC avec windows, ca marche tres bien, mais quand je suis passer à une machine lunix avec gfortran tout est boulverser, alors ou est le problème?

Voila ce que j ai recu sur lunix :

Code :

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

1
2
3
 
 ht=  5.00000000000000028E-002
 gamma=   9.5135076986687324

et voilà mon petit 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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
 
C     fonction test
 
        real function gamma(x)
 
        double precision x
 
        if(x.GT.-0.5.AND.0.5.GT.x) then
 
             gamma =2.d0
 
        else
 
             gamma =1.d0
 
        end if
 
        return
 
        end
 
 
        program ghg
       integer i,j,f,j1,j2,ss,M,p
 
       parameter (M=2**5)
 
       parameter(p=22)
        double precision ht
        ht=0.05d0
        print*,'ht=',ht
        print*,'gamma=',gamma(2.d0*ht)
       open(2,file='gammanexx.dat',status='unknown')
 
        do i=-p,p
 
            write(2,*) i*ht, gamma(i*ht)
 
        end do
 
        close(2)
        end

merci

[Grame]

J'obtiens le même résultat.

En ajoutant des print *, pour déboguer, je m'aperçois
que le programme n'entre pas dans la fonction gamma.

Conclusion : le fortran unix que j'utilise (gfortran) contient une fonction intrinsèque gamma.

En renommant ta fonction zzzgamma, j'obtiens bien ce que tu veux :

ht= 5.00000000000000028E-002
zzzgamma= 2.0000000

Voili, voila.
Sois prudent en utilisant des réels simple et double précision dans la même
routine, il y a toujours un moment où on risque de s'emmêler les pinceaux.