Fortran中Heun方法迭代满足条件却未按预期进入if循环的问题

你的核心问题在于代码的结构是单次迭代而非循环迭代——当前的嵌套if只会在误差第一次超标时执行一次修正操作，但之后不会重新检查误差是否仍然超标，也不会继续迭代直到满足终止条件（误差达标或达到最大迭代次数）。

具体问题拆解 #

看你当前的代码逻辑：

1. 计算初始的y_refined和误差
2. 如果误差超过容差，先把iteration_index加1
3. 然后检查迭代次数是否小于最大值，执行一次修正计算，得到新的误差
4. 但这之后没有回到误差判断的步骤，直接退出整个if块，导致哪怕新误差（比如你输出里的1.66%）仍然远大于0.01%的容差，程序也不会继续迭代，直接跳到下一个x_next

修复方案：改用循环结构 #

你需要用循环（比如do while）来持续检查终止条件，直到误差达标或者迭代次数用尽。修改后的代码如下：

iteration_index = 1
y_refined = y_current + (0.5d0*dx*(dydx_predictor + dydx_refined_corrector))
relative_percent_error = 100.d0*(abs((y_refined - y_next)/y_refined))

! 循环迭代：直到误差达标 或 达到最大迭代次数
do while (relative_percent_error > heun_percent_tolerance .and. iteration_index < max_heun_number)
    iteration_index = iteration_index + 1
    print*, 'Iteration ', iteration_index, ' | x_next: ', x_next, ' | Current error: ', relative_percent_error
    
    ! 更新迭代变量
    y_next = y_refined
    call dydx(x_next, y_next, dydx_refined_corrector)
    y_refined = y_current + (0.5d0*dx*(dydx_predictor + dydx_refined_corrector))
    
    ! 重新计算误差，供下一次循环判断
    relative_percent_error = 100.d0*(abs((y_refined - y_next)/y_refined))
end do

! 可选：打印最终迭代结果
print*, '--- Done with x_next = ', x_next, ' | Total iterations: ', iteration_index, ' | Final error: ', relative_percent_error

修复后的逻辑说明 #

• 用do while循环替代原来的嵌套if，每次循环都会重新检查误差是否超标和迭代次数是否未达上限两个条件
• 每次循环内完成迭代更新：更新y_next、计算新的导数、生成新的y_refined、重新计算误差
• 只有当误差≤0.01% 或者 迭代次数达到15次时，才会退出循环，继续处理下一个x_next

这样就能实现你期望的“持续迭代直到达到容差或次数上限”的效果了。

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