Есть ли стандартный способ проверки Infinite и NaN в Fortran 90/95?

Простой способ без использования ieee_arithmatic состоит в том, чтобы сделать следующее.

Бесконечность: определите свою переменную infinity = HUGE(dbl_prec_var) (или, если она у вас есть, переменную с четырехкратной точностью). Затем вы можете просто проверить, является ли ваша переменная бесконечностью на if(my_var > infinity).

НЭН: Это еще проще. По определению NAN ничему не равен, даже самому себе. Просто сравните переменную с самой собой: if(my_var /= my_var).

У меня недостаточно представителей, чтобы комментировать, поэтому я «отвечу» на предложение Рика Томпсона по тестированию бесконечности.

if (A-1 .eq. A) 

Это также будет верно, если A — очень большое число с плавающей запятой, а 1 меньше точности A.

Простой тест:

subroutine test_inf_1(A)
    real, intent(in) :: A
    print*, "Test (A-1 == A)"
    if (A-1 .eq. A) then
        print*, "    INFINITY!!!"
    else
        print*, "    NOT infinite"
    endif
end subroutine

subroutine test_inf_2(A)
    real, intent(in) :: A
    print*, "Test (A > HUGE(A))"
    if (A > HUGE(A)) then
        print*, "    INFINITY!!!"
    else
        print*, "    NOT infinite"
    endif
end subroutine


program test
    real :: A,B

    A=10
    print*, "A = ",A
    call test_inf_1(A)
    call test_inf_2(A)
    print*, ""

    A=1e20
    print*, "A = ",A
    call test_inf_1(A)
    call test_inf_2(A)
    print*, ""

    B=0.0 ! B is necessary to trick gfortran into compiling this
    A=1/B
    print*, "A = ",A
    call test_inf_1(A)
    call test_inf_2(A)
    print*, ""

end program test

выходы:

A =    10.0000000    
Test (A-1 == A)
    NOT infinite
Test (A > HUGE(A))
    NOT infinite

A =    1.00000002E+20
Test (A-1 == A)
    INFINITY!!!
Test (A > HUGE(A))
    NOT infinite

A =          Infinity
Test (A-1 == A)
    INFINITY!!!
Test (A > HUGE(A))
    INFINITY!!!

No.

Существенные части IEEE_ARITHMETIC для генерации/проверки NaN достаточно легко написать для gfortran для конкретной архитектуры.

Я использовал:

  PROGRAM MYTEST
  USE, INTRINSIC :: IEEE_ARITHMETIC, ONLY: IEEE_IS_FINITE      
  DOUBLE PRECISION :: number, test
  number = 'the expression to test'
  test = number/number
  IF (IEEE_IS_FINITE(test)) THEN
     WRITE(*,*) 'We are OK'
  ELSE
     WRITE(*,*) 'Got a problem'
  END IF         
     WRITE(*,*) number, test
  END PROGRAM MYTEST

Это напечатает «Есть проблема» для числа = 0.0D0, 1.0D0/0.0D0, 0.0D0/0.0D0, SQRT (-2.0D0), а также для переполнения и потери значимости, таких как число = EXP (1.0D800) или число = ЭКСПЕРИМЕНТ(-1.0D800). Обратите внимание, что, как правило, такие вещи, как number = EXP(1.0D-800), просто устанавливают number = 1.0 и выдают предупреждение во время компиляции, но программа печатает «We are OK», что я считаю приемлемым.

OL.

No.

В Fortran 90/95 также нет совместимого со стандартами способа проверки бесконечности или NaN, и не может быть совместимого со стандартами способа. В Фортране 90/95 не существует совместимого со стандартами способа определения любого из этих квазичисел.

Для Фортрана 1/бесконечность=0, таким образом, разделите переменную на ноль, т.е.

program test
implicit none
real :: res
integer :: i

do i=1,1000000
    res=-log((1.+(10**(-real(i))))-1.)
    print*,i,res
    if ((1./res)==0.) then
        exit
    end if
end do

end program

вот твой чек на бесконечность. Нет необходимости в усложнении.

для тестирования NaN ни одна из вещей не работала, например, если тестировать реальный s2p, чтобы увидеть, является ли это NaN, тогда

if(isnan(s2p)) 

не работал в gfortran и не

if(s2p.ne.s2p). 

Единственное, что сработало, это

if(.not.s2p<1.and..not.s2p>1)  

хотя, чтобы убедиться, что вы можете добавить

if(.not.s2p<1.and..not.s2p>1.and..not.s2p==1)    

Для Inf кажется, что если (A-1 .eq. A) верно, то A является Inf