return;
xdelta = ND_pos(q)[0] - ND_pos(p)[0];
ydelta = ND_pos(q)[1] - ND_pos(p)[1];
- dist = sqrt(xdelta * xdelta + ydelta * ydelta);
+ dist = hypot(xdelta, ydelta);
force = (dout * dout) / (X_K * dist);
#elif defined(ALT)
xdelta = ND_pos(q)[0] - ND_pos(p)[0];
ydelta = ND_pos(q)[1] - ND_pos(p)[1];
- dist = sqrt(xdelta * xdelta + ydelta * ydelta);
+ dist = hypot(xdelta, ydelta);
din = RAD(p) + RAD(q);
if (dist < X_K + din)
return;
}
xdelta = ND_pos(q)[0] - ND_pos(p)[0];
ydelta = ND_pos(q)[1] - ND_pos(p)[1];
- dist = sqrt(xdelta * xdelta + ydelta * ydelta);
+ dist = hypot(xdelta, ydelta);
din = RAD(p) + RAD(q);
dout = dist - din;
force = dout * dout / ((X_K + din) * dist);