Discussion :

[felixia]

Bonjour a tous,

Je vous ecris car je me retrouve limité face a un probleme. j ai une fonction a trois indice f(k, i, it= indice sur les itérations)
la boucle k est la plus exterieure, puis viens celle sur it (les itérations) ce que je souhaite faire, c est mettre un critere d arret genre erreure relative a ce niveau, c est a dire avant je prenais une valeur pour it, maitenant j aimerai mettre une condition d arret qui compare deux valeurs de la fonction a deux iterations differentes, cest a dire si f(it+1)-f(it) inf.ou.sup a epsilon continue, sinon exit.

Je pourrais apporter plus de details si desiré.

Merci pour votre aide

[Sylvain Bergeron]

Je ne suis certain de bien comprendre, mais le code suivant résout-il ton problème ?

Code :

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

1
2
3
4
5
6
7
8
9
 
do k = ...
   f_it_prec = valeur très grande
   do it = ...
      f_it = f(k,i,it)
      if (abs(f_it - f_it_prec) < Epsilon) exit 
      f_it_prec = f_it
   enddo
enddo

Si tu veux sortir des 2 boucles sur le test, tu dois mettre une étiquette sur la boucle externe :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
boucle_externe : &
do k = ...
   f_it_prec = valeur très grande
   do it = ...
      f_it = f(k,i,it)
      if (abs(f_it - f_it_prec) < Epsilon) exit boucle_externe
      f_it_prec = f_it
   enddo
enddo boucle_externe

[felixia]

Merci Sylvain pour ta réactivité.
En fait je calcul la valeur de f a la premiere itération, si a l iteration suivante l iteration suivante la valeur differe de la premiere calculée, le calcul doit etre refait jusqu avaoir la convergence souhaitée.
f_it_prec n est pas connue. Je ne sais pas si ca donne plus d explication.

[felixia]

ci-joint un bout du code. Merci pour votre aide

Code :

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

1
2
3
4
5
6
      kk= 1              .................  boucle externe                                               
    do 99 k= 1,NT   
 
      do i=1,NU ...
    f0(i)= f1(i)
    end do

Le calcul ci dessous est réalisé it fois, ce que j'aimerais c'est définir une condition de convergence en comparant f1 de call tridag à l'iteration it de f1 de l'iteration
it+1

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
    do 88 it = 1,1000  je souhaite remplacer cette condition..............
par un calcul d'erreure relative; aui consiterqa à comparer f(it+1) à f(it)
 
       x= 0.
       y1= 0.
       small= 1e-5
 
 
      do i=1,NU
     x= x + f1(i)*usq(i)             
     y1= y1 + f1(i)*utt(i)           
     II(i)= 0.5*x*du    
     JJ1(i)= y1           
 
 
    do i=1,N
    A(i)= cc*(-II(i-1) + JJ(i-1))
    B(i)= 1. + cc*(II(i) - II(i-1) - JJ(i) - JJ(i-1))
    C(i)= cc*(II(i) + JJ(i))
      end do
 
 
      call tridag(A,B,C,f0,f1,NU)

En sortie j'obtiens f1, le calcul est effectué jusqu'à la valeur max de it, au lieu de passer par it, j'aimerai que le calcul s’arrête par exemple lorsque f(it+1)-f(it) gt.samll

Code :

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

 88   continue

[Sylvain Bergeron]

Sous réserve que j'aie manqué un détail, la solution à ce genre de problème est toujours la même : À la fin de chaque itération, tu testes avec le résultat de l'itération précédente, que tu conserves dans une variable "_prec".

Pour la première itération, tu n'as évidemment pas de valeur précédente pour tester. Il y a alors 2 solutions :

• Tu assigne à ta variable _prec une valeur qui t'assure que le test de convergence sera faux.
• Tu utilises une variable logique pour éviter le test de convergence.

Ainsi, tu as :

Code :

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

1
2
3
4
5
6
7
 
_prec = valeur inattendue ! pour que le premier test de convergence soit faux.
do
   Nouvelle valeur = ...
   if (abs(Nouvelle valeur < _prec) < seuil) exit
   _prec = Nouvelle valeur ! pour itération suivante
enddo

ou :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
logical FaireTest
FaireTest = .false.
do
   Nouvelle valeur = ...
   if (FaireTest) then
      if (abs(Nouvelle valeur < _prec) < seuil) exit
   endif
   FaireTest = .true. ! pour les itérations suivantes 
   _prec = Nouvelle valeur ! pour itération suivante
enddo

[felixia]

Merci Sylvain Bergeron, votre aide m'a été d'un grand secours. J'ai utilisé la première option que j' ai un peu modifiée.

[Hédhili Jaïdane]

Bonjour.

Juste en passant, quid des conditions d'arrêt et de sortie des boucles si des fois l'algo ne convergeait pas.

[François L.]

A titre personnel, je mets un maximum d'itérations possible, qui est une condition de sortie supplémentaire de la boucle (et je teste ensuite si on est sorti par convergence ou par un trop grand nombre d'itérations)

soit un truc du style:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
_prec = valeur inattendue ! pour que le premier test de convergence soit faux.
n=0
n_max=100 ! par exemple
do
   n=n+1   
   Nouvelle valeur = ...
   if (abs(Nouvelle valeur < _prec) < seuil) exit
   _prec = Nouvelle valeur ! pour itération suivante
   if (n>n_max) exit  
enddo
if (n>n_max) STOP 'Convergence non atteinte'