地球参数相关的更新函数
void CEarth::Update(const CVect3 &pos, const CVect3 &vn, int isMemsgrade)
用到位置以及速度
那么位置和速度用的哪个时刻的?
假如设计算周期为[T,2T] ;
解算时刻为2T时刻,那么地球参数用的是哪一时刻的?
看用的pos 和vn是哪个时刻.
查看调用函数的实参pos01 和vn01
可见速度是用上一个解算vn 加an*nts2 得到的;nts2是T/2周期,那么这里可以理解为是
T时刻的vn + an * T/2 其实得到是T+T/2时刻的速度,【T,--T+T/2,---2T】
位置同样是用T时刻的pos,与T+T/2的速度 得到T+T/2的pos,所以地球参数在2T时刻用的是T+T/2的地球参数,然后推到2T, [T+T/2,-----2T]
T/2用的地方有多处,当增量与影响在短时间内其引起的变化比较小的,可以用T/2处的参数等价于整个周期的参数。例如有害加速度eth.gcc 就是这样的使用。