关于FORTRAN的EXP函数

Kieran

大家好，

我想咨询下关于精度的问题。

1.如果我想把变量a设置成双精度浮点数，是否可以如下设置？
INTEGER, PARAMETER :: dp=SELECTED_REAL_KIND(15,14)
REAL (KIND=dp) :: a
我是想问，15位有效小数位，取值在10的-14次方到10的14次方的范围。这样设置是否可以，确定变量a为双精度浮点数呢？

2.如果想把整型变量也设置成双精度的话，SELECTED_REAL_KIND中的两个参数要如何设置呢？

3.我按上面的方式设置了双精度的浮点数和整数，并将它们放在了EXP函数里运算，但这个函数却返回Infinity类型的数据。请问要如何解决这个问题呢？这是不是由双精度的浮点数和整数设置不合理导致的呢？

谢谢大家啦。

Kieran

谢谢你的回复。

下面是我出问题的代码部分。

[Fortran] 纯文本查看 复制代码

INTEGER, PARAMETER                      :: dp = SELECTED_REAL_KIND(15,14)
INTEGER                                                  :: i, j
REAL (KIND=dp), ALLOCATABLE   :: ev(:,:)
REAL (KIND=dp)                                   :: fermi
REAL (KIND=dp)                                   :: kb
REAL (KIND=dp)                                   :: te
ei = (ev(i,j) - fermi) * (1.0d0 /(EXP((ev(i,j) - fermi)/kb / te) + 1.0d0))&
         + kb * te * DLOG(1.0d0 + EXP(0.0d0 - (ev(i,j) - fermi) /kb /te))
write (unit=*, fmt=*) 'kb', kb
write (unit=*, fmt=*) 'fermi', fermi
write (unit=*, fmt=*) 'te', te
write (unit=*, fmt=*) 'ev', ev(i,j)
write (unit=*, fmt=*) '(ev(i,j) - fermi)', (ev(i,j) - fermi)
write (unit=*, fmt=*) 'f_function', (1.0d0 /(EXP((ev(i,j) - fermi)/kb / te) + 1.0d0))
write (unit=*, fmt=*) 'EXP(0.0d0 - (ev(i,j) - fermi) /kb /te)', EXP(0.0d0 - (ev(i,j) - fermi) /kb /te)
write (unit=*, fmt=*) 'DLOG_function', DLOG(1.0d0 + EXP(0.0d0 - (ev(i,j) - fermi) /kb /te))
write (unit=*, fmt=*) 'ei=', ei

其中，ev(i,j)是从主程序传递进来的一个二维双精度数组。

编译并运行程序后，出现下面的结果。
kb  8.617333262145000E-005
fermi -0.531700000000000
te   30.0000000000000
ev  -20.4013773387809
(ev(i,j) - fermi)  -17.9840108474454
f_function   1.00000000000000
EXP(0.0d0 - (ev(i,j) - fermi) /kb /te)                Infinity
DLOG_function                Infinity
ei=                Infinity

因为EXP函数的结果是Infinity类型的数据，导致后面的DLOG函数结果和ei变量也都变成了Infinity类型的数据。

请问，1. 我对于双精度浮点数据类型的设定是否有误呢？
2. 要如何调整代码，让EXP函数的输出结果合理呢？

谢谢啦。

Kieran

kyra 发表于 2021-5-8 14:31
试试用双双精度。
dp = SELECTED_REAL_KIND(30,3000)

非常感谢你的回复。

我还想请教一个很基础的问题。selected_real_kind（[p，r]）函数中的p和r分别指代的是小数精度p和十进制指数范围r。既然FORTRAN中的双精度浮点数范围是从2.2250738585072008D-308到4.94065645841246544D-324。那是不是说，双精度的浮点数据类型的设定应该是selected_real_kind（16，324），而非selected_real_kind（15，14）呢？我是指最大范围地包括所有的双精度浮点数据。

谢谢啦。

Kieran

vvt 发表于 2021-5-12 08:56
编译器为你提供几种类型的 real，比如 单精度，双精度，四精度（也叫双双精度，有的编译器不提供四精度）
...

谢谢你的指点。

也就是说，SELECTED_REAL_KIND函数是把，由p和r指定的浮点数，归类成单精度，双精度或者四精度的类型；而其归类的方式，主要是根据p确定的有效数字，来确定的。至于r这个参数的作用不大。因为（p=16,r=324）和(p=15,r=14)的区别主要是由p这个变量来确定的，r的作用似乎可以忽略不计。

我是否可以这样来理解呢？

Kieran

vvt 发表于 2021-5-18 08:30
不是的，p 和 r 都起作用。
SELECTED_REAL_KIND 的作用是：“选择一个精度类型，使其至少满足 p，也同时满 ...

非常感谢你的解释，我现在懂了。