Discussion :

[take_up]

Bonjour,

Je veux convertir le code suivant qui est écrit en Turbo Pascal en code Fortran90. Je ne connais rien en pascal !

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
function rse( a,b:double ) : double;
{Raises a to the b power.}
var
  sign : integer;
  c : double;
begin
  if (a > 0.0) then
    rse := sexp(b*ln(a))
  else if a = 0.0 then
    rse := 0
  else
    begin
      c := trunc(b);
      if b <> c then begin
        write('error in function rse');
        halt;
      end;
      if odd(trunc(b)) then sign := -1 else sign := 1;
        rse := sign*sexp(b*ln(abs(a)));
    end;
end; {rse}

Merci de votre aide !!

[droggo]

Lai,

Je n'en suis pas sûr, mais si ce n'est pas le record des déterrages, ça ne doit pas en être loin.

Pour le code, c'est indiqué dedans : c'est une simple fonction d'élévation du paramètre a à la puissance donnée par le paramètre b de la fonction, ça doit bien exister en Fortran, sinon c'est tout de même facile à implémenter, plutôt que chercher à traduire un langage que tu ne connais pas.

[take_up]

Voici ce que j'ai obtenu en essayant de traduire le code Turbo Pascal. Seulement le programme s'arrête quand je fais une compilation !
Auriez-vous une idée sur ce qui ne va pas dans ma traduction ?

Merci par avance.

Code Fortran :

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
FUNCTION rse(a,b)
    ! Raises a to the b power
    REAL :: a,b
    INTEGER :: mysign
    REAL :: c
    REAL :: rse
 
    IF (a > 0.0) THEN
       rse = sexp(b*LOG(a))
    ELSE IF (a == 0.0) THEN
       rse = 0.
    ELSE
        c = INT(b)
        IF (b /= c) THEN
           PRINT *,'error in function rse'
 !          PRINT *,'a = ',a,'b = ',b,',c = ',c
         STOP
        END IF
        IF (MOD(c,2.0) == 0) THEN
           mysign = 1
        ELSE
           mysign = -1
        END IF
        rse = mysign*sexp(b*LOG(ABS(a)))
    END IF 
  END FUNCTION rse

[droggo]

Qia,

Auriez-vous une idée sur ce qui ne va pas dans ma traduction.

À vrai dire, presque tout.

edit -------------

Oups, j'étais dans l'idée d'une traduction Fortran -> Pascal, au temps pour moi.

[take_up]

c'est pourquoi votre réponse m'avait un peu surpris mais bon le problème est bien posé maintenant !!

[Sylvain Bergeron]

Seulement le programme s'arrête quand je fais une compilation !
Auriez-vous une idée sur ce qui ne va pas dans ma traduction ?

Le compilateur a certainement écrit ce qui n'allait pas avant d'arrêter...

Je n'ai pas de compilateur Fortran sous la main pour tester, mais avec un survol rapide, je vois que tu utilises une fonction « sexp ». Ça n'existe pas en Fortran. Tu devrais utiliser simplement « exp », mais on ne code jamais a^b en exp(b*ln(a)) en Fortran à ma connaissance. On fait simplement a ** b.

[take_up]

sexp est une une fonction que j'avais déja défini comme suit:
FUNCTION sexp(x)
! Prevents runtime errors caused by exp underflowing.
REAL :: x
REAL :: sexp

IF (x < -2000.) THEN
sexp = 0.0
ELSE
sexp = EXP(x)
END IF
END FUNCTION sexp

Je pensais que a**b et exp(b*ln(a)) donnent le même résultat !