javascript - box2dweb中的这段代码求解释，关于万有引力的实现

Question

b2GravityController.prototype.Step = function (step) {
      var i = null;
      var body1 = null;
      var p1 = null;
      var mass1 = 0;
      var j = null;
      var body2 = null;
      var p2 = null;
      var dx = 0;
      var dy = 0;
      var r2 = 0;
      var f = null;
      if (this.invSqr) {
         for (i = this.m_bodyList;
         i; i = i.nextBody) {
            body1 = i.body;
            p1 = body1.GetWorldCenter();
            mass1 = body1.GetMass();
            for (j = this.m_bodyList;
            j != i; j = j.nextBody) {
               body2 = j.body;
               p2 = body2.GetWorldCenter();
               dx = p2.x - p1.x;
               dy = p2.y - p1.y;
               r2 = dx * dx + dy * dy;
               if (r2 < Number.MIN_VALUE) continue;
               f = new b2Vec2(dx, dy);
               f.Multiply(this.G / r2 / Math.sqrt(r2) * mass1 * body2.GetMass());
               if (body1.IsAwake()) body1.ApplyForce(f, p1);
               f.Multiply((-1));
               if (body2.IsAwake()) body2.ApplyForce(f, p2);
            }
         }
      }
      else {
         for (i = this.m_bodyList;
         i; i = i.nextBody) {
            body1 = i.body;
            p1 = body1.GetWorldCenter();
            mass1 = body1.GetMass();
            for (j = this.m_bodyList;
            j != i; j = j.nextBody) {
               body2 = j.body;
               p2 = body2.GetWorldCenter();
               dx = p2.x - p1.x;
               dy = p2.y - p1.y;
               r2 = dx * dx + dy * dy;
               if (r2 < Number.MIN_VALUE) continue;
               f = new b2Vec2(dx, dy);
               f.Multiply(this.G / r2 * mass1 * body2.GetMass());
               if (body1.IsAwake()) body1.ApplyForce(f, p1);
               f.Multiply((-1));
               if (body2.IsAwake()) body2.ApplyForce(f, p2);
            }
         }
      }
   }

这个if语句是为什么区分什么？

Answer 1

if (this.invSqr) {
//blablabla
}else{
//blabla
}

勘误： invSqr 应该是inverse-square 的简写，平方反比(1/r²)，但据下文f.Multiply(this.G / r2 / Math.sqrt(r2) * mass1 * body2.GetMass())来看明显拼写错了，作者应该是想写inverse-cube立方反比(1/r³)，但却误写成invSqr(inverse-square)。

万有引力只适用于宏观尺度，随着粒子直径不断缩小，粒子间的万有引力不断缩水，计算出来的误差越来越大。为了修正这个偏差，GmM/r² =>GmM/r³、甚至GmM/r² =>GmM/r⁴。。。。直至粒子变成了分子原子进入电磁力的主导范围时，这种误差修正就被完全抛弃了。

因此第一个if分支，立方反比(1/r³)是模拟微粒(0.01m~)的运动情况；第二个if分支，平方反比(1/r²)就是正常的粒子(0.1m~10m)运动。

PS: 详细的物理解释请恶补下范德华力之类的粒子作用力系统