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: 详细的物理解释请恶补下范德华力之类的粒子作用力系统