Traduction de code Pascal en Fortran90

Version imprimable

Bonjour,

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

Code:

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

08/09/2011, 20h28
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.

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:

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

09/09/2011, 20h02
droggo

Qia,

Citation:

Envoyé par take_up

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.
09/09/2011, 20h21
take_up

c'est pourquoi votre réponse m'avait un peu surpris mais bon le problème est bien posé maintenant !!
10/09/2011, 14h23
Sylvain Bergeron

Citation:

Envoyé par take_up

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.
19/09/2011, 17h17
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 !