关于FORTRAN的EXP函数

kyra

1. 可以确定。
2. 整形变量是无误差的，不存在双精度，只存在长度（可以用 SELECTED_INT_KIND）
3. 上代码。

kyra

EXP((ev(i,j) - fermi)/kb / te)
=exp(-17.9840108474454/8.617333262145000E-005/30)
=exp(-6956.52573015480)
=6.5957E-3022

双精度可正常表达的精度范围是 -1.79E+308 到 +1.79E+308，最小精度 2.2250738585072014e-308
即，小于 2.2250738585072014e-308 的数字认为是 0

分母为0，则 1.0d0 /(EXP((ev(i,j) - fermi)/kb / te) 为无穷大。

kyra

试试用双双精度。
dp = SELECTED_REAL_KIND(30,3000)

代码中的 1.d0 改成 1.0_dp ， 0.0d0 改成 0.0_dp  （其他常量类似）
DLOG 改成 Log