Discussion :

[driss80]

bonjour à tous;
voilà, j'ai lancé un code écrit sur fortran mais ça prend beaucoup de temps pour s'executer et ainsi donner les résultats souhaités...ma question est la suivante:comment je peux maitriser le temps de calcul pour l'obtention de meilleurs résultats avec une bonne précision au bout d'un temps réduit???..et merci à vous..

[ALT]

Tout dépend du calcul !
Il y en a qui nécessitent d'être exécutés jusqu'au bout pour être significatif (la plupart, d'ailleurs) & d'autres dont la précision s'améliore avec le temps de calcul.
La réponse à ta question est donc dans le calcul lui-même.

[Sylvain Bergeron]

Tout à fait d'accord.

En complément :

• Activer l'optimisation !!!
• Les gains les plus significatifs sont généralement obtenus en utilisant un meilleur algorithme.
• Réduire les lectures-écritures au maximum.

[FR119492]

Salut !

Pour pouvoir te répondre, il nous faudrait savoir

• ce que doit faire ton programme (quel est ton problème?)
• comment il le fait (quel algorithme?)

Jean-Marc Blanc

[driss80]

Salut !

Pour pouvoir te répondre, il nous faudrait savoir

• ce que doit faire ton programme (quel est ton problème?)
• comment il le fait (quel algorithme?)

Jean-Marc Blanc

oui pour les choses soient claires, je vous écrit mon code qui est:

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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
*********************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!             FONCTION CALCUL STATIC
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      
 
      PROGRAM pompage 
      Implicit none
      integer :: i,j,k,m,Nstep
      real (kind=8) :: Rco,Aco,Rfo,Acl,gamma_p
      real (kind=8) :: lambda_s,sigma_ap,sigma_ep,sigma_es,sigma_as
      real (kind=8) :: alpha_p,alpha_s,tau,delta_lambda_s,n
      real (kind=8) :: gamma_s,R1,R2,tau_s,h
      real (kind=8) :: c,vp,vs,A21
      real (kind=8) :: v,deltaX,Dt
      real (kind=8) :: T,deltaT,Tmin,Tmax
      real (kind=8) :: resultat,TOL
      real (kind=8) :: pi=4.D0*atan(1.D0) 
 
      Integer, parameter :: Nmax=201, Nit=200
      real (kind=8), parameter :: PsmoinsL=1.D-06,W0=30.D0, L=500.D0
      real (kind=4), Dimension(Nit) :: N1,N2,X,Pp
      real (kind=8), Dimension(Nit) :: Psplus,Psmoins
      real (kind=8), Dimension(Nit) :: R12,R21,W12,W21
      Do 10 i=1,Nit
   10   N1(i)=0.D0
!
      Do 20 i=1,Nit   
   20   N2(i)=0.D0
!
      Do 30 i=1,Nit   
   30   X(i)=0.D0
!
      Do 40 i=1,Nit       
   40   Pp(i)=0.D0
!      
      Do 50 i=1,Nit     
   50   Psplus(i)=0.D0
!     
      Do 60 i=1,Nit    
   60   Psmoins(i)=0.D0 
!     
      Do 70 i=1,Nit
   70   R12(i)=0.D0
!     
      Do 80 i=1,Nit
   80   R21(i)=0.D0
!      
      Do 90 i=1,Nit
   90   W12(i)=0.D0
!   
      Do 100 i=1,Nit
  100    W21(i)=0.D0 
 
!!!!!!!!!!!!!!!!!!!  ouverture du fichier resultats.dat
 
      open(20,file='F:\projet_these_mghaaraz\calcul_static\&
     & res_Yb1.dat',form='formatted',status='unknown')
!      open(10,file='F:\projet_these_mghaaraz\calcul_static\&
!     & données.dat',form='unformatted',status='unknown')
 
      Dt=2.D-2
      Tmin = 0.D0
      Tmax = 1000.D0
      t=Tmin
      Nstep=(Tmax-Tmin)/Dt
      write(*,*)'Nstep=',Nstep
         t=Nstep*Dt
      write(*,*) 't=',t
      TOL = 1.D-03
 
!!!!!!!!!!!!!!!!!!!  déclaration des données 
 
      data sigma_ap,sigma_ep,sigma_es,sigma_as,alpha_p,alpha_s,&
      delta_lambda_s,gamma_s,R1,R2,tau_s,h,vp,N/3.D-25,5.D-26,&
      2.5D-25,1.4D-27,3.D-03,5.D-03,2.D-08,0.75D0,0.04D0,0.4D0,3.D-07,&
      6.62D-34,3.D08,319.D12/
!
      Rco=8.D-06
      Aco=pi*Rco**2.D0
      read(*,*) Aco
      write(*,*) 'Aco=', Aco
!            
      Rfo=200.D-06     
      Acl=pi*(Rfo**2-Rco**2)
      read(*,*) Acl
      write(*,*) 'Acl=', Acl
!      
      gamma_p=Aco/Acl
      read(*,*) gamma_p
      write(*,*) 'gamma_p=', gamma_p
!     
      c=3.D08
      lambda_s=1.08D-06
      vs=c/lambda_s
      read(*,*) vs
      write(*,*) 'vs=', vs 
!
      tau=1.D-03
      A21=1.D0/tau
      read(*,*)A21
      write(*,*) 'A21=', A21
!
      n=1.45D0
      v=c/n
      read(*,*) v
      write(*,*) 'v=', v
!
      deltaX=L/float(Nit)
      read(*,*) deltaX
      write(*,*) 'deltaX=', deltaX
!     
      write(20,1000) 
      write(20,2000)     
      write(20,1000)
 
!!!!!!!!!!!!!!!!!!!  conditions initiales
 
      N=3.D19
      N2(1)=7.2D20
      N1(1)=2.28D19   ! N1(0)=N-N2(0)
      Pp(1)=W0
      Psplus(1)=0.D0
      X(1)=0.D0
      R12(1)=0.D0
      R21(1)=0.D0
      W12(1)=0.D0
      W21(1)=0.D0
      read(*,*) N2(1)
      write(*,*)'N2(1)=',N2(1)
      read(*,*) N1(1)
      write(*,*)'N1(1)=',N1(1)
      read(*,*) Pp(1)
      write(*,*)'Pp(1)=',Pp(1)
      read(*,*) Psplus(1)
      write(*,*)'Psplus(1)=',Psplus(1)
      read(*,*) X(1)
      write(*,*)'X(1)=',X(1)
      read(*,*) R12(1)
      write(*,*)'R12(1)=',R12(1)
      read(*,*) R21(1)
      write(*,*)'R21(1)=',R21(1)
      read(*,*) W12(1)
      write(*,*)'W12(1)=',W12(1)
      read(*,*) W21(1)
      write(*,*)'W21(1)=',W21(1)     
 
 
!!!!!!!!!!!!!!!!!!!  Boucles de calcul statique
 
      if (i.lt.Nmax) then 
      call init1(Psmoins, 201, 200, X)  
      read (*,*) Psmoins,X
      write (*,*) 'Psmoins(W)=',Psmoins,'X(cm)=',X
         do i=1,Nit
            N1(i+1)=((R21(i)+W21(i)+A21)/(R12(i)+W12(i)+R21(i)+W21(i)+&
            A21))*N 
            N2(i+1)=N-N1(i+1)
 
            do j=1,Nit
               R12(j+1)=(sigma_ap*Pp(j+1)*gamma_p)/(Aco*h*vp)
               R21(j+1)=(sigma_ep*Pp(j+1)*gamma_p)/(Aco*h*vp)
               W12(j+1)=(sigma_as*(Psplus(j+1)+Psmoins(j+1))*gamma_s)/&
               (Aco*h*vs)
               W21(j+1)=(sigma_es*(Psplus(j+1)+Psmoins(j+1))*gamma_s)/&
               (Aco*h*vs)
 
               do k=1,Nit
                  Psplus(k+1)=Psplus(k)+gamma_s*(sigma_es*N2(k)-&
                  sigma_as*N1(k))*Psplus(k)*deltaX+gamma_s*N2(k)*&
                  sigma_es*2.D0*int(n)*h*c*c*(delta_lambda_S/lambda_S**&
                  3)*deltaX-alpha_s*Psplus(k)*deltaX
 
                  do m=1,Nit
                     Pp(m+1)=Pp(m)+gamma_p*(sigma_ep*N2(m)-sigma_ap*&
                     N1(m))*Pp(m)*deltaX-alpha_p*Pp(m)*deltaX
                  end do
 
 
               end do
 
            end do
 
         end do
         write(20,3000)'Pp(W)=',Pp,'N2/N=',N2/N,'Psplus(W)=',Psplus
 
      else if (i.gt.1) then
         do i=Nit, 2, -1
            N1(i-1)=((R21(i)+W21(i)+A21)/(R12(i)+W12(i)+R21(i)+W21(i)+&
            A21))*N
            N2(i-1)=N-N1(i-1)
 
            do j=Nit, 2, -1
               R12(j-1)=(sigma_ap*Pp(j-1)*gamma_p)/(Aco*h*vp)
               R21(j-1)=(sigma_ep*Pp(j-1)*gamma_p)/(Aco*h*vp)
               W12(j-1)=(sigma_as*(Psplus(j-1)+Psmoins(j-1))*gamma_s)/&
               (Aco*h*vs)
               W21(j-1)=(sigma_es*(Psplus(j-1)+Psmoins(j-1))*gamma_s)/&
               (Aco*h*vs)
 
               do k=Nit, 2, -1
                  Psmoins(k-1)=Psmoins(k)+gamma_s*(sigma_es*N2(k)-&
                  sigma_as*N1(k))*Psmoins(k)*deltaX+gamma_s*N2(k)*&
                  sigma_es*2.D0*int(n)*h*c*c*(delta_lambda_S/lambda_S**&
                  3)*deltaX-alpha_s*Psmoins(k)*deltaX
               end do
 
            end do
 
         end do 
         write(20,3000) 'Psmoins(W)=',Psmoins 
      end if
 
 
 
 
!!!!!!!!!!!!!!!!!!!  test condition initiale entre Ps+ et Ps- en 0
 
      R1=R1*100.D+0
      write(*,*) 'R1=', R1
 
      resultat=Psplus(1)/Psmoins(1)
      resultat=resultat*100.D+0
      write(*,*) 'resultat=', resultat
 
      if (resultat.ge.R1-0.005D+0 .and. resultat.le.R1+0.005D+0) then
         write(*,*) 'la valeur initiale de Ps- "OK"'
      else
         write(*,*) 'la valeur initiale de Ps- "not OK"'
      endif
 
      write(20,1000)
 1000 format( 300('-') )    
 2000 format('/Pp(W)',100X,'N2/N',100X,'Psplus(W)',100X,'Psmoins(W) / ')
 3000 format(' /',2(20X,F5.3),2(20X,D10.3),' /') 
      close(20)
 
      contains 
 
!!!!!!!!!!!!!!!!!!!  subroutine initialisation de Psmoins et X
 
      subroutine init1(Psmoins, Nmax, Nit, X)
      implicit none
      real (kind=8), Dimension(Nit) :: Psmoins
      real (kind=4), Dimension(Nit) ::  X
      integer :: i,Nit,Nmax
      real (kind=8) :: PsmoinsL, deltaX, L
      L=500.D0
      PsmoinsL=1.D-06
      Psmoins(1)=PsmoinsL
      X(1)=0.
      deltaX=L/float(Nit)
      do while (i.lt.Nmax)
         do i=1, Nit
         Psmoins(i+1)=Psmoins(1)
         X(i+1)=X(i)+deltaX
         enddo
      end do
      end subroutine init1
 
      end
     *********************************

et merci à votre aide

[ALT]

En somme, tu nous demandes d'optimiser ton code, alors que tu n'as même pas cherché à savoir ce qu'il fait ?

[gregory2526]

salut.

un moyen serait de faire du profiling...
tu peux voir là : http://www.developpez.net/forums/sho...t=gprof&page=2
vers la fin de la discussion.
En utilisant un outil comme gprof, tu peux obtenir les temps processeurs des différentes instructions de ton programme, et donc de voir où le temps est consommé...
Sinon, tu peux tenter d'utiliser les fonction ETIME (au début de la discussion donnès ci-dessus)

bon courage

[FR119492]

Salut !

Pour t'aider, il nous faudrait savoir deux choses:
[liste=1][*]Quel est le temps de calcul pour exactement le même problème sous MatLab et en Fortran?[*]Quand répondras-tu enfin aux questions qu'on te pose, ce qui est nécessaire pour t'aider? Tu nous envoie ton code, qui est affreux, mais tu es débutant. Ce qu'il nous faut, c'est ton problème et ton algorithme.[/liste]

En ce qui concerne le code, quelques remarques:
[liste][*]Toujours programmer le plus simplement possible. Les normes 90 et 95 ont apporté une masse de nouvelles possibilités, certes intéressantes, mais qui compliquent considérablement le code. Pourquoi ne pas écrire

Code :

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

1
2
      Integer i,j,k,m,Nstep
real*8 Rco,Aco,Rfo,Acl,gamma_p

[*]Pourquoi utiliser l'allocation dynamique des tableaux?[*]Pourquoi définir certaines variables réelles sur 4 octets plutôt que sur 8?[*]Pourquoi utiliser 10 boucles pour initialiser les tableaux?
[liste]

Jean-Marc Blanc

[driss80]

Salut !

Pour t'aider, il nous faudrait savoir deux choses:

1. Quel est le temps de calcul pour exactement le même problème sous MatLab et en Fortran?
2. Quand répondras-tu enfin aux questions qu'on te pose, ce qui est nécessaire pour t'aider? Tu nous envoie ton code, qui est affreux, mais tu es débutant. Ce qu'il nous faut, c'est ton problème et ton algorithme.

En ce qui concerne le code, quelques remarques:
[list][*]Toujours programmer le plus simplement possible. Les normes 90 et 95 ont apporté une masse de nouvelles possibilités, certes intéressantes, mais qui compliquent considérablement le code. Pourquoi ne pas écrire

Code :

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

1
2
      Integer i,j,k,m,Nstep
real*8 Rco,Aco,Rfo,Acl,gamma_p

[*]Pourquoi utiliser l'allocation dynamique des tableaux?[*]Pourquoi définir certaines variables réelles sur 4 octets plutôt que sur 8?[*]Pourquoi utiliser 10 boucles pour initialiser les tableaux?
[list]

Jean-Marc Blanc

bonjour, mon probleme consiste à calculer les integrales ci dessous et ainsi retrouver les N2(x), Pp(x), Ps+(x), Ps-(x) d'après les équations suivantes:

0=(R12+W12)N1(x)-(R21+W21+A21)N2(x);
dPp(x)/dx=Gamma_p(sigma_epN2(x)- sigma_apN1(x))Pp(x)-alpha_pPp(x) ;
±dPs±(x)/dx=Gamma_s(sigma_esN2(x)- sigma_asN1(x)) Ps±(x)+ Gamma_s N2(x) sigma_es.2nmhc2Delta_lambda_s/lambda_s3 -alpha_s Ps±(x) ;
avec:
R12=(sigma_ap*Pp(x)*Gamma_p)/(Aco*h*vp)
R21=(sigma_ep*Pp(x)*Gamma_p)/(Aco*h*vp)
W12=(sigma_as*(Ps+(x)+Ps-(x))*Gamma_s)/(Aco*h*vs)
W12=(sigma_es*(Ps+(x)+Ps-(x))*Gamma_s)/(Aco*h*vs)
A21=1/tau

avec les conditions initiales suivantes:
Pp(0)=W0, Ps+(0)=R1*Ps-(0), Ps-(L)=0

et merci a tous...bonne journée..

[FR119492]

Salut.

Maintenant, je crois y voir un peu plus clair. Alors deux remarques:

La première est de nature algorithmique: tu as des équations différentielles ordinaires de la forme
dy/dx = f(x,y)
que tu veux intégrer. Malheureusement, tu utilises la méthode d'Euler d'ordre 1, dont on a dit: "Admirez mais n'appliquez pas!", qui consiste à remplacer l'équation différentielle par une relation de récurrence de la forme
y(k+1) = y(k) + dx * F(x(k),y(k))
Cette méthode est inefficace, imprécise et souvent instable. Tu devrais plutôt utiliser, par exemple, une méthode de Runge-Kutta d'ordre 4.

La seconde concerne directement ton programme: est-tu vraiment sur que la formule donnant pp(m+1) doit être appliquée 200*200*200*200, c'est-à-dire 1'600'000'000 fois ?

Bonne chance
Jean-Marc Blanc

[driss80]

Salut.

Maintenant, je crois y voir un peu plus clair. Alors deux remarques:

La première est de nature algorithmique: tu as des équations différentielles ordinaires de la forme
dy/dx = f(x,y)
que tu veux intégrer. Malheureusement, tu utilises la méthode d'Euler d'ordre 1, dont on a dit: "Admirez mais n'appliquez pas!", qui consiste à remplacer l'équation différentielle par une relation de récurrence de la forme
y(k+1) = y(k) + dx * F(x(k),y(k))
Cette méthode est inefficace, imprécise et souvent instable. Tu devrais plutôt utiliser, par exemple, une méthode de Runge-Kutta d'ordre 4.

La seconde concerne directement ton programme: est-tu vraiment sur que la formule donnant pp(m+1) doit être appliquée 200*200*200*200, c'est-à-dire 1'600'000'000 fois ?

Bonne chance
Jean-Marc Blanc

merci pour vos conseils, je vais les prendre en considération et on vera ce que ça donne...pour les Pp(m+1), puisque les équations sont couplées, j'ai essayé de faire des boucles imbriquées pour les calculer et je nai pas fait attention au nombre de fois; merci encore....