我正在做一个简单的n体仿真。现在,我要对其进行“强制执行”,这意味着我要计算每个对象在每一帧上对每个其他对象施加的所有力。
我现在的问题是,如果我选择了大量的对象(例如2000),在某些情况下,一开始,“行星”的对象就会消失大约2帧。当我在主循环中添加System.out.println(PlanetHandler.planets.get(0).position.x);
来检查发生了什么时,我得到了
487.0
486.99454
NaN
NaN
通过评论一些东西和反复试验,我发现问题出在这里:
private static void computeAndSetPullForce(Planet planet)
{
for(Planet otherPlanet : planets)
{
//Also here, if we are deleting the planet, don't interact with it.
if(otherPlanet != planet && !otherPlanet.delete)
{
//First we get the x,y and magnitudal distance between the two bodies.
int xDist = (int) (otherPlanet.position.x - planet.position.x);
int yDist = (int) (otherPlanet.position.y - planet.position.y);
float dist = Vector2Math.distance(planet.position, otherPlanet.position);
//Now we compute first the total and then the component forces
//Depending on choice, use r or r^2
float force = Constants.GRAVITATIONAL_CONSTANT * ((planet.mass*otherPlanet.mass)/(dist*dist));
float forceX = force * xDist/dist;
float forceY = force * yDist/dist;
//Given the component forces, we construct the force vector and apply it to the body.
Vector2 forceVec = new Vector2(forceX, forceY);
planet.force = Vector2Math.add(planet.force, forceVec);
}
}
}
“行星”列表是CopyOnWriteArray<Planets>
。
我已经花了很长时间了,但是还没有弄清楚是什么导致值(位置,速度)过大。也许曾经对此有一定经验或通常精通此类事情的人可以帮助我。
答案 0 :(得分:1)
这是JVM给您NAN的典型情况。您遇到的是零除以零( 0/0 ),在数学上这是不确定的形式。
如果float dist = Vector2Math.distance(planet.position, otherPlanet.position);
返回0。
下一条语句
float force = Constants.GRAVITATIONAL_CONSTANT * ((planet.mass*otherPlanet.mass)/(dist*dist));
其中计算力的值除以零。
此外,我建议您在需要精度时使用BigDecimal。您还可以参考答案here
之一