使用 dt 必须设置为固定的预定值的新信息(至少除了最后一步以外的所有值),这是我的建议:
double T0 = 0.0;
double T = 1.0;
int N = floor((T - T0)/dt);
double t = T0;
for (int step_number = 0; step_number < N; ++step_number, t += dt)
{
t = T0 + step_number * dt;
do_one_step(t, T, dt);
}
if (t < T)
{
do_one_step(t, T, T - t);
}
函数 do_one_step 为每次迭代使用 t , T 和 dt 执行必要的计算 . 必须由函数更新的数据可以作为类的成员变量,也可以作为非const引用包含在函数参数列表中 .
顺便说一句,我最后调用循环外的函数不是为了节省分支条件的可能成本,而是因为我发现代码更好地组织并且更容易理解 .
旧答案:
正如你所说,你可以在结尾轻松获得一个非常小的时间增量,因为结果 T/N 通常不精确(并且 0.03 当然不准确) .
我更喜欢像这样开发 t 和 dt :
int N = 7;
double t = 0.0;
double T = 1.0;
double dt = (T - t)/N;
for (int step_number =