/**
* Determine closest point on hyperbolic arm to a point pt
*
* @param pt : point
* @return t value for closest point; not necessarily within clipped range
*/
public double closestPointTo(FPoint2 pt) {
final boolean db = false;
pt = toCurveSpace(pt, null);
double U = B + 1;
double V = -pt.y;
double W = B * pt.x;
Polyn p =
new Polyn( //
B * U * U, //
2 * B * U * V, //
B * V * V + A * U * U - W * W, //
2 * A * U * V, //
A * V * V //
);
if (db) Streams.out.println("closestPointTo, pt=" + pt + "\n" + p);
double ret = 0;
try {
DArray r = new DArray();
if (Math.abs(p.c(0)) < 1e-5) r.addDouble(0);
else p.solve(r);
if (r.isEmpty()) {
throw new FPError("can't find closest point, poly=\n" + p);
}
double bestDist = 0;
for (int i = 0; i < r.size(); i++) {
double t = r.getDouble(i);
FPoint2 apt = calcPoint(t);
double dist = apt.distance(pt);
if (i == 0 || dist < bestDist) {
bestDist = dist;
ret = t;
}
}
} catch (FPError e) {
Tools.warn("caught FPError");
// Streams.out.println("caught:\n" + e);
ret = (this.minParameter() + this.maxParameter()) * .5;
}
return ret;
}
