[Fortran 90] Le temps de calcul d'un code

Version imprimable

07/11/2007, 12h44
driss80

[Fortran 90] Le temps de calcul d'un code

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..
07/11/2007, 17h54
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.
07/11/2007, 21h53
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.
07/11/2007, 23h26
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

Citation:

Envoyé par 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

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

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

08/11/2007, 10h50
ALT

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

Citation:

Envoyé par ALT

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

c mon code je l'ai fait moi meme, il calcule des tableaux avec des boucles imbriquées, le probleme c'est que mon programme prend beaucoup de temps à les calculer, l'execussion est tres lente...s'il vous plait les gens de critiques s'abstenir...on est là pour discuter mais non critiquer les uns les autres et merci...
08/11/2007, 11h27
ALT

Alors si c'est toi qui l'as écrit, tu dois pouvoir nous dire où l'exécution est lente, mettre éventuellement des commentaires à chaque boucle, expliquer (en début de programme, si ,si, ça se fait !) à quoi sert le programme, comment il doit le faire, etc., utiliser des noms de variables plus explicites (W11, i, j... il y a peut-être plus clair).

Autre chose, il y a au moins quatre boucles imbriquées. N'y a-t-il pas moyen de faire autrement ?
Car il est plus facile d'introduire une erreur (un bug) dans une telle structure & c'est assez lent à l'exécution...

Citation:

Envoyé par ALT

Alors si c'est toi qui l'as écrit, tu dois pouvoir nous dire où l'exécution est lente, mettre éventuellement des commentaires à chaque boucle, expliquer (en début de programme, si ,si, ça se fait !) à quoi sert le programme, comment il doit le faire, etc., utiliser des noms de variables plus explicites (W11, i, j... il y a peut-être plus clair).

Autre chose, il y a au moins quatre boucles imbriquées. N'y a-t-il pas moyen de faire autrement ?
Car il est plus facile d'introduire une erreur (un bug) dans une telle structure & c'est assez lent à l'exécution...

**************************************************

premierement vous devez savoir que je ne suis pas un informaticien mais je suis plutot physicien, je travaille sur les lasers et j'ai un programme en matlab que je dois faire en fortran et c'est la premiere fois que j'y travaille, ce qui fait je trouve des difficultés...

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
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
263
264
265
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!             FONCTION CALCUL STATIC
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!      
 
      PROGRAM pompage 
      Implicit none
!!!!!!!!!!!!!! declaration des reels et des tableaux
      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
!!!!!!!!!!!!!!! initialisation des tableaux
      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
!!!!!!!!!!!!!!!!!!!  Boucle de calcul aller
 
      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
!!!!!!!!!!!!!!  Boucle de calcul retour
 
      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

08/11/2007, 16h01
ALT

Ouh là !
Faut pas se fâcher !

Mais qu'on soit informaticien ou physicien ou n'importe quoi d'autre, les règles de la programmation s'appliquent à tous les programmeurs.
Ainsi ça permet de faire des codes propres, facilement lisibles par tous & y compris par le programmeur lui-même quand il reprend son programme deux ans plus tard.

Mes remarques n'étaient que des conseils en ce sens.

Pour le reste, comme tu n'avais rien précisé au départ, personne ne pouvait deviner le contexte, le problème, &c.
Navré, mais je ne suis pas devin.
Et effectivement, j'avoue que je t'avais pris pour un étudiant (& il y en a un grand nombre !) qui essayait de se faire faire son exo gratos.
Toutes mes excuses si ce n'est pas le cas.
08/11/2007, 17h49
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
10/11/2007, 18h57
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:

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
12/11/2007, 11h15
driss80
Citation:
Envoyé par FR119492

Salut !

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

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.

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:

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..
12/11/2007, 18h13
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
13/11/2007, 09h41
driss80

Citation:

Envoyé par 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

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....
18/11/2007, 13h56
qwilldune

Matlab vs Fortran

Salut,

Je suis à peu près dans le même cas de figure que toi.
J'ai un code Matlab qui marche mais il est un peu lent.
Je veux le passer sur Fortran mais cela fait une décennie que je n'en ai pas fait.
Il faut que je m'exerce...;)

Je suis partant pour m'essayer à ton exercice si ton code Matlab n'est pas confidentiel. Je veux bien essayer d'en faire une version sous fortran

En passant, tu utilises quel compilateur?
19/11/2007, 10h52
driss80

Citation:

Envoyé par qwilldune

Salut,

Je suis à peu près dans le même cas de figure que toi.
J'ai un code Matlab qui marche mais il est un peu lent.
Je veux le passer sur Fortran mais cela fait une décennie que je n'en ai pas fait.
Il faut que je m'exerce...;)

Je suis partant pour m'essayer à ton exercice si ton code Matlab n'est pas confidentiel. Je veux bien essayer d'en faire une version sous fortran

En passant, tu utilises quel compilateur?

je suis vraiment desole mais je ne peux pas vous passer mon code en matlab, parce que je travaille dans un milieu professionnel et ça ne m'autorise pas à disperser nos programmes car c confidentiel, vous comprenez j'espere...
19/11/2007, 11h03
Jerome Briot

Citation:

Envoyé par qwilldune

J'ai un code Matlab qui marche mais il est un peu lent.

Si ce code n'a pas déjà été optimisé, je t'invite à venir en discuter sur le forum MATLAB afin d'essayer de l'accélérer un peu (au moins ;))