public synchronized double move(MappableGraph map) {
boolean doitall = false;
double x = 0;
double y = 0;
double xx = 0;
double yy = 0;
double r = 0;
double s = 0;
double t = 0;
double dx = 0;
double dy = 0;
double d2x = 0;
double d2y = 0;
double d2 = 0;
double d2xy = 0;
Mappable a = null;
Mappable b = null;
CartesianPoint[] n = null;
Mappable[] o = (Mappable[]) map.getShownVertices().toArray(new Mappable[0]);
RadixPlaneSorter<Mappable> p = null;
double delta = 0;
if (alpha > Double.MIN_VALUE) {
p = new RadixPlaneSorter<Mappable>(Math.sqrt(gamma), Math.sqrt(gamma), 40);
for (int i = 0; i < o.length; i++) p.add(o[i]);
doitall = true;
}
for (int i = 0; i < o.length; i++) {
if (!(a = o[i]).isMobile()) continue;
dx = dy = d2x = d2y = d2xy = 0;
if (doitall) {
n = p.getLocal(a, gamma);
p.remove(a);
for (int j = 0; j < n.length; j++) {
b = (Mappable) n[j];
x = a.getX() - b.getX();
xx = x * x;
y = a.getY() - b.getY();
yy = y * y;
r = xx + yy;
if (r < Double.MIN_VALUE) continue;
s = r - gamma;
t = 3 * gamma - r;
r = 1 / Math.sqrt(r);
s *= r * r * r;
t *= r * r * r * r * r;
/*
r = Math.sqrt(r);
s = (r*r -gamma) /r/r/r;
t = (3*gamma - r*r)/r/r/r/r/r;
*/
dx += x * s;
dy += y * s;
d2x += s + xx * t;
d2y += s + yy * t;
}
dx *= alpha;
dy *= alpha;
d2x *= alpha;
d2y *= alpha;
}
x = y = r = 0;
n = (Mappable[]) map.getShownNeighbours(a).toArray(new Mappable[0]);
for (int j = 0; j < n.length; j++) {
if ((b = (Mappable) n[j]) == a) continue;
x += b.getX();
y += b.getY();
r += 1;
}
r *= 2;
dx += r * a.getX() - 2 * x;
dy += r * a.getY() - 2 * y;
d2x += r;
d2y += r;
d2 = Math.abs(d2x + d2y);
if (d2 > Double.MIN_VALUE) {
a.setX(a.getX() - dx / d2);
a.setY(a.getY() - dy / d2);
delta += ((dx * dx) + (dy * dy)) / d2 / d2;
}
if (doitall) p.add(a);
}
return delta;
}
public synchronized double move(MappableGraph map) {
boolean doitall = false;
double x = 0;
double y = 0;
double xx = 0;
double yy = 0;
double r = 0;
double s = 0;
double t = 0;
double dx = 0;
double dy = 0;
double d2x = 0;
double d2y = 0;
double d2xy = 0;
double d2 = 0;
Mappable a = null;
Mappable b = null;
CartesianPoint[] n = null;
Mappable[] o = (Mappable[]) map.getShownVertices().toArray(new Mappable[0]);
RadixPlaneSorter<Mappable> p = null;
double change = 0;
if (gamma > Double.MIN_VALUE && alpha > Double.MIN_VALUE) {
// p = new RadixPlaneSorter(Math.sqrt(gamma)/4,Math.sqrt(gamma)/4,40);
p = new RadixPlaneSorter<Mappable>(Math.sqrt(gamma), Math.sqrt(gamma), 40);
for (int i = 0; i < o.length; i++) p.add(o[i]);
doitall = true;
}
for (int i = 0; i < o.length; i++) {
if (!(a = o[i]).isMobile()) continue;
dx = dy = d2x = d2y = 0;
if (doitall) {
n = p.getLocal(a, gamma);
p.remove(a);
for (int j = 0; j < n.length; j++) {
b = (Mappable) n[j];
x = a.getX() - b.getX();
xx = x * x;
y = a.getY() - b.getY();
yy = y * y;
r = xx + yy;
if (r < Double.MIN_VALUE) continue;
// s = df(r);
// t = 2*d2f(r);
r = 1 / r;
s = gamma * r;
s = s * s;
t = 4 * s * r;
s = 1 - s;
dx += s * x;
dy += s * y;
d2x += s + t * xx;
d2y += s + t * yy;
// d2xy += t*x*y;
}
dx *= alpha;
dy *= alpha;
d2x *= alpha;
d2y *= alpha;
// d2xy *= alpha;
}
x = y = r = 0;
n = (Mappable[]) map.getShownNeighbours(a).toArray(new Mappable[0]);
for (int j = 0; j < n.length; j++) {
if ((b = (Mappable) n[j]) == a) continue;
x += b.getX();
y += b.getY();
r += 1;
}
dx += r * a.getX() - x;
dy += r * a.getY() - y;
d2x += r;
d2y += r;
n = (CartesianPoint[]) map.getShownOutNeighbours(a).toArray(new Mappable[0]);
s = 0;
t = 0;
for (int j = 0; j < n.length; j++) {
if ((b = (Mappable) n[j]) == a) continue;
s += a.getY() - b.getY() + delta;
t += 1;
}
n = (CartesianPoint[]) map.getShownInNeighbours(a).toArray(new Mappable[0]);
for (int j = 0; j < n.length; j++) {
if ((b = (Mappable) n[j]) == a) continue;
s += a.getY() - b.getY() - delta;
t += 1;
}
dy += beta * s;
d2y += beta * t;
d2 = Math.abs(d2x + d2y);
// d2 = Math.abs(d2x+d2y+2*d2xy);
if (d2 > Double.MIN_VALUE) {
a.setX(a.getX() - dx / d2);
a.setY(a.getY() - dy / d2);
change += ((dx * dx) + (dy * dy)) / d2 / d2;
}
if (doitall) p.add(a);
}
return change;
}
