Discussion :

[feynman]

Bonjour,
j'ai un truc très bizarre qui passe avec moi.
J'ai écrit un programme qui calcule une fonction, ensuite met ses valeurs dans un fichier data, mais malgré que tout mon programme est écrit en double précision, quand j'ai affiché mon fichier data, il est affiché en real ! ce qui m'étonne, alors où est le problème ?
Merci

Voilà 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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
 
      program rien
      integer i
      double precision h
      h=1.d-2
       open(1,file='gamma8.dat',status='unknown')
          do i=-100,100
              write(1,*) i*h, gammasuppcompact(i*h,100,500)
          end do
        close(1)
       end
 
 
      function gammasuppcompact(x,n,m)
      integer n,m,i
      double precision f1,f2,h,pi,x,gammasuppcompact
      pi=4.d0*datan(1.d0)
      f1=0.d0
      f2=0.d0
      h=2.d0/(dble(n)*dble(m))
      do i=1,m-1
         f1=f1+gamma(x+1.d0/dble(n)-i*h)*dble(n)*
&              (dcos((n*pi/2.d0)*(-1.d0/dble(n)+i*h)))**2
      enddo
      do i=0,m-1
         f2=f2+gamma(x+1.d0/dble(n)-i*h-h/2.d0)*n*
&                (dcos((dble(n)*pi/2.d0)*(-1.d0/dble(n)+i*h+h/2.d0)))**2
      enddo
      gammasuppcompact=(h/6.d0)*(2.d0*f1+4.d0*f2)
      return
      end
 
      function gamma(x)
      double precision x,gamma
      gamma=1.d0+(1.d0+signe(-x**2/4.d-2+1.d0))/2.d0
      return
      end
 
      function signe(x)
      double precision x,signe
      if(x.ge.0.d0) then
         signe=1.d0
      else
         signe=-1.d0
      end if
      return
      end

[genteur slayer]

utilisez les format les gars!!!!

dans ton write tu lui demande de mettre le format par défaut ( * ) donc, comme de normalement il t'affiche les double comme cela l'arrange....

Code :

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

1
2
 
write(1,'(2d25.15)')i*h, gammasuppcompact(i*h,100,500)

là, il t'affichera un truc du genre:

Code :

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

1
2
 
0.582325236512548D+00 0.882322876512548D+08

tu peux lui dire de mettre un E à la place du D en utilisant le format:
'(2e25.15)'

[feynman]

Merci genteur slayer , t a totalement raison, j ai oublié de mettre le format,
par contre peut me detailler un peut cette ecriture (2d25.15)?
merci

[genteur slayer]

pour plus de précision, va voir le cours
http://www.idris.fr/data/cours/lang/...tran_cours.pdf

et le chapitre format (car il y en a plein d'autre)

mais pour le 2d25.15 cela signifi, dans l'ordre: 2 objets / d :double précision / 25 caractère / 15 chiffres après la virgule

dans les caractère, il faut compter : l'espace devant, le signe, le "0.", les chiffres significatif (ici 15), la lettre "D" et l'exposant (+00 par exemple)

[feynman]

bon mais quand je met write(1,'(2d25.20)') ou plus que 20 ( do i=1,74
write(*,'(2d25.20)'),'rho(',i,')=',PP(i))
enddo, des choses bizares apparaisent :
*************************

[Sylvain Bergeron]

Bonjour M. l'irréductible du implicit none.

L'ensemble de tes problèmes de précision, d'affichage (même avec * comme format), ... découlent du fait que tu ne déclare pas en double precision les fonctions dans les appelants. Les fonctions retournent donc des doubles, mais ces résultats sont interprétés comme des simples par les appelants.

Avec des implicit none, tu aurais trouvé l'erreur à la première compilation...