计算机模拟中的高阶数

Question

我写了一个 Fortran 代码来求解 PDE（比如连续性方程）但是初始 未知值的顺序是 1.0e20 这意味着我的代码给出无穷大的 NANE（不是数字），因为它乘以或除以大数 运行 模拟这么大的数字怎么办？
方程是：泊松方程和连续性方程

Answer 1

您可以使用 extended precision、real*8 或 double precision（它们是 64 位浮点表示）作为类型，而不是 real（这是32 位）。这将给出至少 308 的指数范围，而不是 38 的较小范围。