Modern Fortran中如何确定数值计算的时间步总数

你遇到的差1问题本质是浮点精度误差导致的常规计算问题，当(tstop - tstart)刚好是dt整数倍的场景下，浮点运算可能得到略小于理论值的结果，直接用FLOOR就会出现少算1步的情况，下面是两类可落地的实现方案：

方案1：固定步数DO循环（优先推荐） #

固定步数DO循环的编译器优化效率更高，性能远优于动态判断的DO WHILE循环，适合dt固定的绝大多数数值计算场景，只需要在计算nsteps时增加浮点容差校正即可避免差1问题：

1. 引入极小的容差阈值，通常取dt的1e-12倍即可（远小于常规数值计算的精度要求）
2. 校正后的nsteps计算逻辑：

real :: eps = 1e-12 * dt
nsteps = FLOOR( (tstop - tstart + eps) / dt )
! 极端场景额外校验
if ( tstart + (nsteps + 1) * dt <= tstop + eps ) then
    nsteps = nsteps + 1
end if

3. 后续继续用你原来的DO i=0,nsteps循环即可，计算t时优先用t = tstart + i*dt而非累加dt，可避免浮点误差累积。

方案2：DO WHILE循环（适合特殊精度要求/可变dt场景） #

如果你的计算场景中dt是动态变化的，或者对时间边界的准确性要求极高，可以用DO WHILE循环，完全通过当前时间判断循环是否终止，避免提前计数的误差问题，实现逻辑：

integer :: i = 0
real :: t, eps = 1e-12 * dt
t = tstart
DO WHILE (t <= tstop + eps)
    ! 放入你的时间步计算逻辑
    i = i + 1
    t = tstart + i * dt
END DO

该方案的缺点是编译器无法提前获取迭代次数，优化空间小，计算性能会略低于固定步数循环。

内容的提问来源于stack exchange，提问作者user106306