Discussion :

[Marty MacFly]

Bonjour à tous,

Je souhaiterais convertir un petit programme Fortran en C.

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
real*8 function gint (g1,g2,eps,pp,Re,s1,c1,ssig,csig,along1)
implicit real*8 (a-h,o-z)
   n = 8
1  h = (g2-g1) / dfloat(n)
   hh = 0.5D0*h
   g = h*0.5773502691896D0
   y = g1 + (h-g)*0.5D0
   gint2 = eld(y,pp,Re,s1,c1,ssig,csig,along1)+
           + eld(y+g,pp,Re,s1,c1,ssig,csig,along1)
   do m = 1,n-1
         gint2 = gint2 + eld(y+h-g,pp,Re,s1,c1,ssig,csig,along1)+
                   + eld(y+h,pp,Re,s1,c1,ssig,csig,along1)
         y = y + h
   enddo
   gint2 = gint2*hh
   if (n.eq.8.or.abs(gint1-gint2).gt.eps*abs(gint1)) then
       n = n*2
       gint1 = gint2
       if (n.lt.1024) goto 1
   endif
   gint = gint2+(gint2-gint1)/15.0D0
   return
   end

Désolé pour le manque de clarté éventuelle, à ma décharge, ce n'est pas moi l'auteur de cette fonction Fortran.

J'ai pensé à la boucle do_while sur la valeur de n (inférieure à 1024)

Le problème est que ce goto est situé à l'intérieur d'une boucle IF.

Je crois pas qu'il soit possible de fermer avec le while du do_while à l'intérieur du if(n eq8 etc...) converti en C, et si je mets ce while à l'extérieur, le souci est que ce IF peut faire tourner la boucle à l'infini, pour peu que la 2ème condition ne soit pas réalisée.

Je patauge un peu, je ne demande pas à ce qu'on me serve le programme en C sur un plateau d'argent, mais si je pouvais avoir quelques pistes sur la structure à adopter pour m'affranchir de ce goto...

Je vous remercie de vos lumières, en vous souhaitant un bon week-end.

[bertry]

Salut,

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
real*8 function gint (g1,g2,eps,pp,Re,s1,c1,ssig,csig,along1)
implicit real*8 (a-h,o-z)
   n = 8
while ( 1 )  //Boucle toujours ( cette ligne serait ton label 1 )
{
     h = (g2-g1) / dfloat(n)
   hh = 0.5D0*h
   g = h*0.5773502691896D0
   y = g1 + (h-g)*0.5D0
   gint2 = eld(y,pp,Re,s1,c1,ssig,csig,along1)+
           + eld(y+g,pp,Re,s1,c1,ssig,csig,along1)
   do m = 1,n-1
         gint2 = gint2 + eld(y+h-g,pp,Re,s1,c1,ssig,csig,along1)+
                   + eld(y+h,pp,Re,s1,c1,ssig,csig,along1)
         y = y + h
   enddo
   gint2 = gint2*hh
   if (n.eq.8.or.abs(gint1-gint2).gt.eps*abs(gint1)) then
       n = n*2
       gint1 = gint2
       if ( n >= 1024 )
         break;  //Sortie de la boucle
     endif
}
   gint = gint2+(gint2-gint1)/15.0D0
   return
   end

Le code en orange est le code fortran à tranposer en c(++)

Explication :

Code :

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

1
2
3
while(1)
{
}

est une boucle avec un test toujours true ( ici 1 ), elle boucle sans fin, jusqu'a ce que tu fasse

Code :

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

break;

[Marty MacFly]

Salut,

Le souci est que dans le cas où nous n'avons ni n=8 ni la deuxième condition, alors l' "intérieur" du if du C n'est pas parcouru, il n'y a pas doublement de n ni de break, et le while se répète joyeusement à l'infini... :/

Tandis que dans le code Fortran, sauf erreur de ma part : si aucune des deux conditions du IF n'est vérifiée, alors il ne revient pas à la ligne 1, puisque l' "intérieur" du IF n'est pas parcouru. Il passe direct à la ligne "gint = gint2+(gint2-gint1)/15.0D0" et la fonction nommée gint est alors presque finie...

Peut-être que je me trompe...

Cela dit je te remercie du temps que tu as pris pour me répondre.

[souviron34]

Code C :

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
 
double gint (double g1, double g2, double eps, double pp,
             double Re, double s1, double c1, double ssig, 
             double csig, double along1)
{
   int    n=8, m ;
   double h, hh, g, y, gint1=0.0, gint2=0.0 ;
 
   while  ( (n < 1024) && ((n == 8) || (fabs(gint1-gint2) > (eps*fabs(gint1)))) )
     {
         h = (g2-g1) / (double)(n) ;
         hh = 0.5 *h ;
         g = h*0.5773502691896 ;
         y = g1 + (h-g)*0.5 ;
 
         gint2 = eld (y,pp,Re,s1,c1,ssig,csig,along1)+
                   eld(y+g,pp,Re,s1,c1,ssig,csig,along1) ;
 
         for ( m = 1 ; m <= n-1 ; m++ )
           {
              gint2 = gint2 + eld(y+h-g,pp,Re,s1,c1,ssig,csig,along1)+
                         eld(y+h,pp,Re,s1,c1,ssig,csig,along1) ;
              y = y + h ;
            }
 
         gint2 = gint2*hh ;
 
         n = n*2 ;
         gint1 = gint2 ;
     }
 
   return ( (gint2+(gint2-gint1)/15.0) );
}

Explication :

Le point central de ce qui te pose problème est :

Code Fortran :

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

1
2
3
4
5
6
   n = 8
1  ....
   if (n.eq.8.or.abs(gint1-gint2).gt.eps*abs(gint1)) then
       ...
       if (n.lt.1024) goto 1
   endif

Si on lit correctement, on voit que :

• lorsque n=8, c'est à dire au départ, on rentre dans le test, et que comme n*2 est < 1024, on va en 1.
• ensuite, si "abs(gint1-gint2).gt.eps*abs(gint1)" et que n < 1024, on va en 1

Ce qui se traduit donc par :

tant que n < 1024 (et n = 8 rentre dans ce cas) ET QUE on a la condition de > (sauf si n = 8), alors on continue

D'où un :

Code C :

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

while ( n < 1024 )

Auquel il faut ajouter la condition du test (mais qui n'est à faire QUE QUAND n n'est pas égal à 8).

D'où la deuxième condition du while :

Code C :

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

while ( (n == 8) || (fabs(gint1-gint2) > (eps*fabs(gint1)))

ici, quand n = 8 on n'évalue même pas le test, et on ne l'évalue que à partir de la première boucle..

CQFD...

[Sve@r]

Le problème est que ce goto est situé à l'intérieur d'une boucle IF.

boucle if pfff...