GNU gfortran的expm1

Question

有什麼辦法可以從GNU Fortran中快速實現expm1嗎？ 理想情況下，最好有一個函數直接計算（exp（x）-1）/ x以避免額外檢查零參數。 expm1的元素版本會特別有用。

Answer 1

這是它是如何從libm中稱爲：

use, intrinsic :: iso_c_binding, only: c_double 

    implicit none 

    interface 
    real(c_double) function expm1(x) bind(c, name='expm1') 
     import c_double 
     real(c_double), intent(in), value :: x 
    end function expm1 
    end interface 

    print*, expm1(3.4d0) 

end program 

如果該函數的glibc的源代碼看起來並不沮喪，那麼你可能希望它以使其元素轉化成的Fortran（如果元素你的意思是Fortran關鍵字）。