Die Unendlichkeit in Fortran

Was ist der sicherste Weg, um eine variable zu +Infinity in Fortran? Im moment bin ich mit:

program test
  implicit none
  print *,infinity()
contains
  real function infinity()
    implicit none
    real :: x
    x = huge(1.)
    infinity = x + x
  end function infinity
end program test

aber ich Frage mich, ob es einen besseren Weg gibt?

InformationsquelleAutor astrofrog | 2011-02-15

fortran infinity

9

Wenn Ihr compiler unterstützt ISO TR 15580 IEEE-Arithmetik, die einen Teil des sogenannten Fortran 2003 standard als Sie verwenden können, Verfahren von ieee_* Module.
```
PROGRAM main

  USE ieee_arithmetic

  IMPLICIT NONE

  REAL :: r

  IF (ieee_support_inf(r)) THEN
    r = ieee_value(r,  ieee_negative_inf)
  END IF

  PRINT *, r

END PROGRAM main
```
InformationsquelleAutor Wildcat
0

Ich würde mich nicht darauf verlassen, dass der compiler die Unterstützung der IEEE-standard, und tun ziemlich viel, was Sie Tat, mit zwei änderungen:
1. Ich würde nicht hinzufügen huge(1.)+huge(1.), da einige Compiler können Sie am Ende mit -huge(1.)+1 --- und dies kann zu einem Speicherleck (ich weiß nicht, warum, aber es ist eine experimentelle Tatsache, so zu sagen).
2. Sind Sie mit real Typen hier. Ich persönlich bevorzuge alle meine floating-point-zahlen als real*8, damit alle float-Konstanten, die qualifiziert sind, mit d0 wie diese: huge(1.d0). Dies ist nicht die Regel, natürlich; einige Menschen lieber mit beiden real-s und real*8-s.
- real*8 und double precision (1.d0) sind nicht necesarilly der gleichen realen Art. Und natürlich, ob single oder double precision ist nicht eine Frage der persönlichen Vorliebe sein, aber der mathematische Argumente und tests.
InformationsquelleAutor ev-br
0

Ich weiß nicht, über die am sichersten, aber ich kann Ihnen eine alternative Methode. Ich habe gelernt, es zu tun auf diese Weise:
```
PROGRAM infinity
  IMPLICIT NONE
  INTEGER :: inf
  REAL :: infi
  EQUIVALENCE (inf,infi) !Stores two variable at the same address
  DATA inf/z'7f800000'/!Hex for +Infinity
  WRITE(*,*)infi
END PROGRAM infinity
```
Wenn Sie mit außergewöhnlichen Werten in Ausdrücken (ich glaube nicht, dass dies generell ratsam ist) sollte man achten, wie sich Ihr compiler verarbeitet Sie, Sie könnten einige unerwartete Ergebnisse anders.

InformationsquelleAutor Usagi
0

Ich bin mir nicht sicher, ob die Lösung Balg funktioniert auf allen Compilern, aber es ist eine schöne mathematische Weise zu erreichen, infinity-als -log(0).
```
program test
  implicit none
  print *,infinity()
contains
  real function infinity()
    implicit none
    real :: x
    x = 0
    infinity=-log(x)
  end function infinity
end program test
```
Funktioniert auch schön für komplexe Variablen.
- Es wird Arbeit auf dem IEEE-konforme Maschinen, wenn Sie nicht aktivieren FPE-trapping. Aber wenn Sie arbeiten mit unendlichen, es wäre dumm, das zu tun.
InformationsquelleAutor Mikael Fremling

-1

Dies scheint für mich arbeiten.
Definieren Sie einen parameter

double precision,parameter :: inf = 1.d0/0.d0

Dann verwenden, wenn in tests.

  real :: sng
  double precision :: dbl1,dbl2

  sng = 1.0/0.0
  dbl1 = 1.d0/0.d0
  dbl2 = -log(0.d0)

  if(sng == inf) write(*,*)"sng = inf"
  if(dbl1 == inf) write(*,*)"dbl1 = inf"
  if(dbl2 == inf) write(*,*)"dbl2 = inf"
  read(*,*)

Beim kompilieren mit ifort & run, bekomme ich

sng = inf
dbl1 = inf
dbl2 = inf

Das funktioniert nicht, zumindest in gfortran. Es löst einen compiler-Fehler.

InformationsquelleAutor Ben1344

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.