@excalidraw/excalidraw 0.17.1-7441-0269e5a → 0.17.1-7441-4c9a6f4

package/dist/packages/excalidraw/element/collision.js

@@ -0,0 +1,727 @@
+import * as GA from "../ga";
+import * as GAPoint from "../gapoints";
+import * as GADirection from "../gadirections";
+import * as GALine from "../galines";
+import * as GATransform from "../gatransforms";
+import { distance2d, rotatePoint, isPathALoop, isPointInPolygon, rotate, } from "../math";
+import { pointsOnBezierCurves } from "points-on-curve";
+import { getElementAbsoluteCoords, getCurvePathOps, getRectangleBoxAbsoluteCoords, } from "./bounds";
+import { hasBoundTextElement, isFrameLikeElement, isIframeLikeElement, isImageElement, } from "./typeChecks";
+import { isTextElement } from ".";
+import { isTransparent } from "../utils";
+import { shouldShowBoundingBox } from "./transformHandles";
+import { getBoundTextElement } from "./textElement";
+import { ShapeCache } from "../scene/ShapeCache";
+const isElementDraggableFromInside = (element) => {
+    if (element.type === "arrow") {
+        return false;
+    }
+    if (element.type === "freedraw") {
+        return true;
+    }
+    const isDraggableFromInside = !isTransparent(element.backgroundColor) ||
+        hasBoundTextElement(element) ||
+        isIframeLikeElement(element);
+    if (element.type === "line") {
+        return isDraggableFromInside && isPathALoop(element.points);
+    }
+    return isDraggableFromInside || isImageElement(element);
+};
+export const hitTest = (element, appState, frameNameBoundsCache, x, y) => {
+    // How many pixels off the shape boundary we still consider a hit
+    const threshold = 10 / appState.zoom.value;
+    const point = [x, y];
+    if (isElementSelected(appState, element) &&
+        shouldShowBoundingBox([element], appState)) {
+        return isPointHittingElementBoundingBox(element, point, threshold, frameNameBoundsCache);
+    }
+    const boundTextElement = getBoundTextElement(element);
+    if (boundTextElement) {
+        const isHittingBoundTextElement = hitTest(boundTextElement, appState, frameNameBoundsCache, x, y);
+        if (isHittingBoundTextElement) {
+            return true;
+        }
+    }
+    return isHittingElementNotConsideringBoundingBox(element, appState, frameNameBoundsCache, point);
+};
+export const isHittingElementBoundingBoxWithoutHittingElement = (element, appState, frameNameBoundsCache, x, y) => {
+    const threshold = 10 / appState.zoom.value;
+    // So that bound text element hit is considered within bounding box of container even if its outside actual bounding box of element
+    // eg for linear elements text can be outside the element bounding box
+    const boundTextElement = getBoundTextElement(element);
+    if (boundTextElement &&
+        hitTest(boundTextElement, appState, frameNameBoundsCache, x, y)) {
+        return false;
+    }
+    return (!isHittingElementNotConsideringBoundingBox(element, appState, frameNameBoundsCache, [x, y]) &&
+        isPointHittingElementBoundingBox(element, [x, y], threshold, frameNameBoundsCache));
+};
+export const isHittingElementNotConsideringBoundingBox = (element, appState, frameNameBoundsCache, point) => {
+    const threshold = 10 / appState.zoom.value;
+    const check = isTextElement(element)
+        ? isStrictlyInside
+        : isElementDraggableFromInside(element)
+            ? isInsideCheck
+            : isNearCheck;
+    return hitTestPointAgainstElement({
+        element,
+        point,
+        threshold,
+        check,
+        frameNameBoundsCache,
+    });
+};
+const isElementSelected = (appState, element) => appState.selectedElementIds[element.id];
+export const isPointHittingElementBoundingBox = (element, [x, y], threshold, frameNameBoundsCache) => {
+    // frames needs be checked differently so as to be able to drag it
+    // by its frame, whether it has been selected or not
+    // this logic here is not ideal
+    // TODO: refactor it later...
+    if (isFrameLikeElement(element)) {
+        return hitTestPointAgainstElement({
+            element,
+            point: [x, y],
+            threshold,
+            check: isInsideCheck,
+            frameNameBoundsCache,
+        });
+    }
+    const [x1, y1, x2, y2] = getElementAbsoluteCoords(element);
+    const elementCenterX = (x1 + x2) / 2;
+    const elementCenterY = (y1 + y2) / 2;
+    // reverse rotate to take element's angle into account.
+    const [rotatedX, rotatedY] = rotate(x, y, elementCenterX, elementCenterY, -element.angle);
+    return (rotatedX > x1 - threshold &&
+        rotatedX < x2 + threshold &&
+        rotatedY > y1 - threshold &&
+        rotatedY < y2 + threshold);
+};
+export const bindingBorderTest = (element, { x, y }) => {
+    const threshold = maxBindingGap(element, element.width, element.height);
+    const check = isOutsideCheck;
+    const point = [x, y];
+    return hitTestPointAgainstElement({
+        element,
+        point,
+        threshold,
+        check,
+        frameNameBoundsCache: null,
+    });
+};
+export const maxBindingGap = (element, elementWidth, elementHeight) => {
+    // Aligns diamonds with rectangles
+    const shapeRatio = element.type === "diamond" ? 1 / Math.sqrt(2) : 1;
+    const smallerDimension = shapeRatio * Math.min(elementWidth, elementHeight);
+    // We make the bindable boundary bigger for bigger elements
+    return Math.max(16, Math.min(0.25 * smallerDimension, 32));
+};
+const hitTestPointAgainstElement = (args) => {
+    switch (args.element.type) {
+        case "rectangle":
+        case "iframe":
+        case "embeddable":
+        case "image":
+        case "text":
+        case "diamond":
+        case "ellipse":
+            const distance = distanceToBindableElement(args.element, args.point);
+            return args.check(distance, args.threshold);
+        case "freedraw": {
+            if (!args.check(distanceToRectangle(args.element, args.point), args.threshold)) {
+                return false;
+            }
+            return hitTestFreeDrawElement(args.element, args.point, args.threshold);
+        }
+        case "arrow":
+        case "line":
+            return hitTestLinear(args);
+        case "selection":
+            console.warn("This should not happen, we need to investigate why it does.");
+            return false;
+        case "frame":
+        case "magicframe": {
+            // check distance to frame element first
+            if (args.check(distanceToBindableElement(args.element, args.point), args.threshold)) {
+                return true;
+            }
+            const frameNameBounds = args.frameNameBoundsCache?.get(args.element);
+            if (frameNameBounds) {
+                return args.check(distanceToRectangleBox(frameNameBounds, args.point), args.threshold);
+            }
+            return false;
+        }
+    }
+};
+export const distanceToBindableElement = (element, point) => {
+    switch (element.type) {
+        case "rectangle":
+        case "image":
+        case "text":
+        case "iframe":
+        case "embeddable":
+        case "frame":
+        case "magicframe":
+            return distanceToRectangle(element, point);
+        case "diamond":
+            return distanceToDiamond(element, point);
+        case "ellipse":
+            return distanceToEllipse(element, point);
+    }
+};
+const isStrictlyInside = (distance, threshold) => {
+    return distance < 0;
+};
+const isInsideCheck = (distance, threshold) => {
+    return distance < threshold;
+};
+const isNearCheck = (distance, threshold) => {
+    return Math.abs(distance) < threshold;
+};
+const isOutsideCheck = (distance, threshold) => {
+    return 0 <= distance && distance < threshold;
+};
+const distanceToRectangle = (element, point) => {
+    const [, pointRel, hwidth, hheight] = pointRelativeToElement(element, point);
+    return Math.max(GAPoint.distanceToLine(pointRel, GALine.equation(0, 1, -hheight)), GAPoint.distanceToLine(pointRel, GALine.equation(1, 0, -hwidth)));
+};
+const distanceToRectangleBox = (box, point) => {
+    const [, pointRel, hwidth, hheight] = pointRelativeToDivElement(point, box);
+    return Math.max(GAPoint.distanceToLine(pointRel, GALine.equation(0, 1, -hheight)), GAPoint.distanceToLine(pointRel, GALine.equation(1, 0, -hwidth)));
+};
+const distanceToDiamond = (element, point) => {
+    const [, pointRel, hwidth, hheight] = pointRelativeToElement(element, point);
+    const side = GALine.equation(hheight, hwidth, -hheight * hwidth);
+    return GAPoint.distanceToLine(pointRel, side);
+};
+const distanceToEllipse = (element, point) => {
+    const [pointRel, tangent] = ellipseParamsForTest(element, point);
+    return -GALine.sign(tangent) * GAPoint.distanceToLine(pointRel, tangent);
+};
+const ellipseParamsForTest = (element, point) => {
+    const [, pointRel, hwidth, hheight] = pointRelativeToElement(element, point);
+    const [px, py] = GAPoint.toTuple(pointRel);
+    // We're working in positive quadrant, so start with `t = 45deg`, `tx=cos(t)`
+    let tx = 0.707;
+    let ty = 0.707;
+    const a = hwidth;
+    const b = hheight;
+    // This is a numerical method to find the params tx, ty at which
+    // the ellipse has the closest point to the given point
+    [0, 1, 2, 3].forEach((_) => {
+        const xx = a * tx;
+        const yy = b * ty;
+        const ex = ((a * a - b * b) * tx ** 3) / a;
+        const ey = ((b * b - a * a) * ty ** 3) / b;
+        const rx = xx - ex;
+        const ry = yy - ey;
+        const qx = px - ex;
+        const qy = py - ey;
+        const r = Math.hypot(ry, rx);
+        const q = Math.hypot(qy, qx);
+        tx = Math.min(1, Math.max(0, ((qx * r) / q + ex) / a));
+        ty = Math.min(1, Math.max(0, ((qy * r) / q + ey) / b));
+        const t = Math.hypot(ty, tx);
+        tx /= t;
+        ty /= t;
+    });
+    const closestPoint = GA.point(a * tx, b * ty);
+    const tangent = GALine.orthogonalThrough(pointRel, closestPoint);
+    return [pointRel, tangent];
+};
+const hitTestFreeDrawElement = (element, point, threshold) => {
+    // Check point-distance-to-line-segment for every segment in the
+    // element's points (its input points, not its outline points).
+    // This is... okay? It's plenty fast, but the GA library may
+    // have a faster option.
+    let x;
+    let y;
+    if (element.angle === 0) {
+        x = point[0] - element.x;
+        y = point[1] - element.y;
+    }
+    else {
+        // Counter-rotate the point around center before testing
+        const [minX, minY, maxX, maxY] = getElementAbsoluteCoords(element);
+        const rotatedPoint = rotatePoint(point, [minX + (maxX - minX) / 2, minY + (maxY - minY) / 2], -element.angle);
+        x = rotatedPoint[0] - element.x;
+        y = rotatedPoint[1] - element.y;
+    }
+    let [A, B] = element.points;
+    let P;
+    // For freedraw dots
+    if (distance2d(A[0], A[1], x, y) < threshold ||
+        distance2d(B[0], B[1], x, y) < threshold) {
+        return true;
+    }
+    // For freedraw lines
+    for (let i = 0; i < element.points.length; i++) {
+        const delta = [B[0] - A[0], B[1] - A[1]];
+        const length = Math.hypot(delta[1], delta[0]);
+        const U = [delta[0] / length, delta[1] / length];
+        const C = [x - A[0], y - A[1]];
+        const d = (C[0] * U[0] + C[1] * U[1]) / Math.hypot(U[1], U[0]);
+        P = [A[0] + U[0] * d, A[1] + U[1] * d];
+        const da = distance2d(P[0], P[1], A[0], A[1]);
+        const db = distance2d(P[0], P[1], B[0], B[1]);
+        P = db < da && da > length ? B : da < db && db > length ? A : P;
+        if (Math.hypot(y - P[1], x - P[0]) < threshold) {
+            return true;
+        }
+        A = B;
+        B = element.points[i + 1];
+    }
+    const shape = ShapeCache.get(element);
+    // for filled freedraw shapes, support
+    // selecting from inside
+    if (shape && shape.sets.length) {
+        return element.fillStyle === "solid"
+            ? hitTestCurveInside(shape, x, y, "round")
+            : hitTestRoughShape(shape, x, y, threshold);
+    }
+    return false;
+};
+const hitTestLinear = (args) => {
+    const { element, threshold } = args;
+    if (!ShapeCache.get(element)) {
+        return false;
+    }
+    const [point, pointAbs, hwidth, hheight] = pointRelativeToElement(args.element, args.point);
+    const side1 = GALine.equation(0, 1, -hheight);
+    const side2 = GALine.equation(1, 0, -hwidth);
+    if (!isInsideCheck(GAPoint.distanceToLine(pointAbs, side1), threshold) ||
+        !isInsideCheck(GAPoint.distanceToLine(pointAbs, side2), threshold)) {
+        return false;
+    }
+    const [relX, relY] = GAPoint.toTuple(point);
+    const shape = ShapeCache.get(element);
+    if (!shape) {
+        return false;
+    }
+    if (args.check === isInsideCheck) {
+        const hit = shape.some((subshape) => hitTestCurveInside(subshape, relX, relY, element.roundness ? "round" : "sharp"));
+        if (hit) {
+            return true;
+        }
+    }
+    // hit test all "subshapes" of the linear element
+    return shape.some((subshape) => hitTestRoughShape(subshape, relX, relY, threshold));
+};
+// Returns:
+//   1. the point relative to the elements (x, y) position
+//   2. the point relative to the element's center with positive (x, y)
+//   3. half element width
+//   4. half element height
+//
+// Note that for linear elements the (x, y) position is not at the
+// top right corner of their boundary.
+//
+// Rectangles, diamonds and ellipses are symmetrical over axes,
+// and other elements have a rectangular boundary,
+// so we only need to perform hit tests for the positive quadrant.
+const pointRelativeToElement = (element, pointTuple) => {
+    const point = GAPoint.from(pointTuple);
+    const [x1, y1, x2, y2] = getElementAbsoluteCoords(element);
+    const center = coordsCenter(x1, y1, x2, y2);
+    // GA has angle orientation opposite to `rotate`
+    const rotate = GATransform.rotation(center, element.angle);
+    const pointRotated = GATransform.apply(rotate, point);
+    const pointRelToCenter = GA.sub(pointRotated, GADirection.from(center));
+    const pointRelToCenterAbs = GAPoint.abs(pointRelToCenter);
+    const elementPos = GA.offset(element.x, element.y);
+    const pointRelToPos = GA.sub(pointRotated, elementPos);
+    const halfWidth = (x2 - x1) / 2;
+    const halfHeight = (y2 - y1) / 2;
+    return [pointRelToPos, pointRelToCenterAbs, halfWidth, halfHeight];
+};
+const pointRelativeToDivElement = (pointTuple, rectangle) => {
+    const point = GAPoint.from(pointTuple);
+    const [x1, y1, x2, y2] = getRectangleBoxAbsoluteCoords(rectangle);
+    const center = coordsCenter(x1, y1, x2, y2);
+    const rotate = GATransform.rotation(center, rectangle.angle);
+    const pointRotated = GATransform.apply(rotate, point);
+    const pointRelToCenter = GA.sub(pointRotated, GADirection.from(center));
+    const pointRelToCenterAbs = GAPoint.abs(pointRelToCenter);
+    const elementPos = GA.offset(rectangle.x, rectangle.y);
+    const pointRelToPos = GA.sub(pointRotated, elementPos);
+    const halfWidth = (x2 - x1) / 2;
+    const halfHeight = (y2 - y1) / 2;
+    return [pointRelToPos, pointRelToCenterAbs, halfWidth, halfHeight];
+};
+// Returns point in absolute coordinates
+export const pointInAbsoluteCoords = (element, 
+// Point relative to the element position
+point) => {
+    const [x, y] = point;
+    const [x1, y1, x2, y2] = getElementAbsoluteCoords(element);
+    const cx = (x2 - x1) / 2;
+    const cy = (y2 - y1) / 2;
+    const [rotatedX, rotatedY] = rotate(x, y, cx, cy, element.angle);
+    return [element.x + rotatedX, element.y + rotatedY];
+};
+const relativizationToElementCenter = (element) => {
+    const [x1, y1, x2, y2] = getElementAbsoluteCoords(element);
+    const center = coordsCenter(x1, y1, x2, y2);
+    // GA has angle orientation opposite to `rotate`
+    const rotate = GATransform.rotation(center, element.angle);
+    const translate = GA.reverse(GATransform.translation(GADirection.from(center)));
+    return GATransform.compose(rotate, translate);
+};
+const coordsCenter = (x1, y1, x2, y2) => {
+    return GA.point((x1 + x2) / 2, (y1 + y2) / 2);
+};
+// The focus distance is the oriented ratio between the size of
+// the `element` and the "focus image" of the element on which
+// all focus points lie, so it's a number between -1 and 1.
+// The line going through `a` and `b` is a tangent to the "focus image"
+// of the element.
+export const determineFocusDistance = (element, 
+// Point on the line, in absolute coordinates
+a, 
+// Another point on the line, in absolute coordinates (closer to element)
+b) => {
+    const relateToCenter = relativizationToElementCenter(element);
+    const aRel = GATransform.apply(relateToCenter, GAPoint.from(a));
+    const bRel = GATransform.apply(relateToCenter, GAPoint.from(b));
+    const line = GALine.through(aRel, bRel);
+    const q = element.height / element.width;
+    const hwidth = element.width / 2;
+    const hheight = element.height / 2;
+    const n = line[2];
+    const m = line[3];
+    const c = line[1];
+    const mabs = Math.abs(m);
+    const nabs = Math.abs(n);
+    let ret;
+    switch (element.type) {
+        case "rectangle":
+        case "image":
+        case "text":
+        case "iframe":
+        case "embeddable":
+        case "frame":
+        case "magicframe":
+            ret = c / (hwidth * (nabs + q * mabs));
+            break;
+        case "diamond":
+            ret = mabs < nabs ? c / (nabs * hwidth) : c / (mabs * hheight);
+            break;
+        case "ellipse":
+            ret = c / (hwidth * Math.sqrt(n ** 2 + q ** 2 * m ** 2));
+            break;
+    }
+    return ret || 0;
+};
+export const determineFocusPoint = (element, 
+// The oriented, relative distance from the center of `element` of the
+// returned focusPoint
+focus, adjecentPoint) => {
+    if (focus === 0) {
+        const [x1, y1, x2, y2] = getElementAbsoluteCoords(element);
+        const center = coordsCenter(x1, y1, x2, y2);
+        return GAPoint.toTuple(center);
+    }
+    const relateToCenter = relativizationToElementCenter(element);
+    const adjecentPointRel = GATransform.apply(relateToCenter, GAPoint.from(adjecentPoint));
+    const reverseRelateToCenter = GA.reverse(relateToCenter);
+    let point;
+    switch (element.type) {
+        case "rectangle":
+        case "image":
+        case "text":
+        case "diamond":
+        case "iframe":
+        case "embeddable":
+        case "frame":
+        case "magicframe":
+            point = findFocusPointForRectangulars(element, focus, adjecentPointRel);
+            break;
+        case "ellipse":
+            point = findFocusPointForEllipse(element, focus, adjecentPointRel);
+            break;
+    }
+    return GAPoint.toTuple(GATransform.apply(reverseRelateToCenter, point));
+};
+// Returns 2 or 0 intersection points between line going through `a` and `b`
+// and the `element`, in ascending order of distance from `a`.
+export const intersectElementWithLine = (element, 
+// Point on the line, in absolute coordinates
+a, 
+// Another point on the line, in absolute coordinates
+b, 
+// If given, the element is inflated by this value
+gap = 0) => {
+    const relateToCenter = relativizationToElementCenter(element);
+    const aRel = GATransform.apply(relateToCenter, GAPoint.from(a));
+    const bRel = GATransform.apply(relateToCenter, GAPoint.from(b));
+    const line = GALine.through(aRel, bRel);
+    const reverseRelateToCenter = GA.reverse(relateToCenter);
+    const intersections = getSortedElementLineIntersections(element, line, aRel, gap);
+    return intersections.map((point) => GAPoint.toTuple(GATransform.apply(reverseRelateToCenter, point)));
+};
+const getSortedElementLineIntersections = (element, 
+// Relative to element center
+line, 
+// Relative to element center
+nearPoint, gap = 0) => {
+    let intersections;
+    switch (element.type) {
+        case "rectangle":
+        case "image":
+        case "text":
+        case "diamond":
+        case "iframe":
+        case "embeddable":
+        case "frame":
+        case "magicframe":
+            const corners = getCorners(element);
+            intersections = corners
+                .flatMap((point, i) => {
+                const edge = [point, corners[(i + 1) % 4]];
+                return intersectSegment(line, offsetSegment(edge, gap));
+            })
+                .concat(corners.flatMap((point) => getCircleIntersections(point, gap, line)));
+            break;
+        case "ellipse":
+            intersections = getEllipseIntersections(element, gap, line);
+            break;
+    }
+    if (intersections.length < 2) {
+        // Ignore the "edge" case of only intersecting with a single corner
+        return [];
+    }
+    const sortedIntersections = intersections.sort((i1, i2) => GAPoint.distance(i1, nearPoint) - GAPoint.distance(i2, nearPoint));
+    return [
+        sortedIntersections[0],
+        sortedIntersections[sortedIntersections.length - 1],
+    ];
+};
+const getCorners = (element, scale = 1) => {
+    const hx = (scale * element.width) / 2;
+    const hy = (scale * element.height) / 2;
+    switch (element.type) {
+        case "rectangle":
+        case "image":
+        case "text":
+        case "iframe":
+        case "embeddable":
+        case "frame":
+        case "magicframe":
+            return [
+                GA.point(hx, hy),
+                GA.point(hx, -hy),
+                GA.point(-hx, -hy),
+                GA.point(-hx, hy),
+            ];
+        case "diamond":
+            return [
+                GA.point(0, hy),
+                GA.point(hx, 0),
+                GA.point(0, -hy),
+                GA.point(-hx, 0),
+            ];
+    }
+};
+// Returns intersection of `line` with `segment`, with `segment` moved by
+// `gap` in its polar direction.
+// If intersection coincides with second segment point returns empty array.
+const intersectSegment = (line, segment) => {
+    const [a, b] = segment;
+    const aDist = GAPoint.distanceToLine(a, line);
+    const bDist = GAPoint.distanceToLine(b, line);
+    if (aDist * bDist >= 0) {
+        // The intersection is outside segment `(a, b)`
+        return [];
+    }
+    return [GAPoint.intersect(line, GALine.through(a, b))];
+};
+const offsetSegment = (segment, distance) => {
+    const [a, b] = segment;
+    const offset = GATransform.translationOrthogonal(GADirection.fromTo(a, b), distance);
+    return [GATransform.apply(offset, a), GATransform.apply(offset, b)];
+};
+const getEllipseIntersections = (element, gap, line) => {
+    const a = element.width / 2 + gap;
+    const b = element.height / 2 + gap;
+    const m = line[2];
+    const n = line[3];
+    const c = line[1];
+    const squares = a * a * m * m + b * b * n * n;
+    const discr = squares - c * c;
+    if (squares === 0 || discr <= 0) {
+        return [];
+    }
+    const discrRoot = Math.sqrt(discr);
+    const xn = -a * a * m * c;
+    const yn = -b * b * n * c;
+    return [
+        GA.point((xn + a * b * n * discrRoot) / squares, (yn - a * b * m * discrRoot) / squares),
+        GA.point((xn - a * b * n * discrRoot) / squares, (yn + a * b * m * discrRoot) / squares),
+    ];
+};
+export const getCircleIntersections = (center, radius, line) => {
+    if (radius === 0) {
+        return GAPoint.distanceToLine(line, center) === 0 ? [center] : [];
+    }
+    const m = line[2];
+    const n = line[3];
+    const c = line[1];
+    const [a, b] = GAPoint.toTuple(center);
+    const r = radius;
+    const squares = m * m + n * n;
+    const discr = r * r * squares - (m * a + n * b + c) ** 2;
+    if (squares === 0 || discr <= 0) {
+        return [];
+    }
+    const discrRoot = Math.sqrt(discr);
+    const xn = a * n * n - b * m * n - m * c;
+    const yn = b * m * m - a * m * n - n * c;
+    return [
+        GA.point((xn + n * discrRoot) / squares, (yn - m * discrRoot) / squares),
+        GA.point((xn - n * discrRoot) / squares, (yn + m * discrRoot) / squares),
+    ];
+};
+// The focus point is the tangent point of the "focus image" of the
+// `element`, where the tangent goes through `point`.
+export const findFocusPointForEllipse = (ellipse, 
+// Between -1 and 1 (not 0) the relative size of the "focus image" of
+// the element on which the focus point lies
+relativeDistance, 
+// The point for which we're trying to find the focus point, relative
+// to the ellipse center.
+point) => {
+    const relativeDistanceAbs = Math.abs(relativeDistance);
+    const a = (ellipse.width * relativeDistanceAbs) / 2;
+    const b = (ellipse.height * relativeDistanceAbs) / 2;
+    const orientation = Math.sign(relativeDistance);
+    const [px, pyo] = GAPoint.toTuple(point);
+    // The calculation below can't handle py = 0
+    const py = pyo === 0 ? 0.0001 : pyo;
+    const squares = px ** 2 * b ** 2 + py ** 2 * a ** 2;
+    // Tangent mx + ny + 1 = 0
+    const m = (-px * b ** 2 +
+        orientation * py * Math.sqrt(Math.max(0, squares - a ** 2 * b ** 2))) /
+        squares;
+    let n = (-m * px - 1) / py;
+    if (n === 0) {
+        // if zero {-0, 0}, fall back to a same-sign value in the similar range
+        n = (Object.is(n, -0) ? -1 : 1) * 0.01;
+    }
+    const x = -(a ** 2 * m) / (n ** 2 * b ** 2 + m ** 2 * a ** 2);
+    return GA.point(x, (-m * x - 1) / n);
+};
+export const findFocusPointForRectangulars = (element, 
+// Between -1 and 1 for how far away should the focus point be relative
+// to the size of the element. Sign determines orientation.
+relativeDistance, 
+// The point for which we're trying to find the focus point, relative
+// to the element center.
+point) => {
+    const relativeDistanceAbs = Math.abs(relativeDistance);
+    const orientation = Math.sign(relativeDistance);
+    const corners = getCorners(element, relativeDistanceAbs);
+    let maxDistance = 0;
+    let tangentPoint = null;
+    corners.forEach((corner) => {
+        const distance = orientation * GALine.through(point, corner)[1];
+        if (distance > maxDistance) {
+            maxDistance = distance;
+            tangentPoint = corner;
+        }
+    });
+    return tangentPoint;
+};
+const pointInBezierEquation = (p0, p1, p2, p3, [mx, my], lineThreshold) => {
+    // B(t) = p0 * (1-t)^3 + 3p1 * t * (1-t)^2 + 3p2 * t^2 * (1-t) + p3 * t^3
+    const equation = (t, idx) => Math.pow(1 - t, 3) * p3[idx] +
+        3 * t * Math.pow(1 - t, 2) * p2[idx] +
+        3 * Math.pow(t, 2) * (1 - t) * p1[idx] +
+        p0[idx] * Math.pow(t, 3);
+    // go through t in increments of 0.01
+    let t = 0;
+    while (t <= 1.0) {
+        const tx = equation(t, 0);
+        const ty = equation(t, 1);
+        const diff = Math.sqrt(Math.pow(tx - mx, 2) + Math.pow(ty - my, 2));
+        if (diff < lineThreshold) {
+            return true;
+        }
+        t += 0.01;
+    }
+    return false;
+};
+const hitTestCurveInside = (drawable, x, y, roundness) => {
+    const ops = getCurvePathOps(drawable);
+    const points = [];
+    let odd = false; // select one line out of double lines
+    for (const operation of ops) {
+        if (operation.op === "move") {
+            odd = !odd;
+            if (odd) {
+                points.push([operation.data[0], operation.data[1]]);
+            }
+        }
+        else if (operation.op === "bcurveTo") {
+            if (odd) {
+                points.push([operation.data[0], operation.data[1]]);
+                points.push([operation.data[2], operation.data[3]]);
+                points.push([operation.data[4], operation.data[5]]);
+            }
+        }
+        else if (operation.op === "lineTo") {
+            if (odd) {
+                points.push([operation.data[0], operation.data[1]]);
+            }
+        }
+    }
+    if (points.length >= 4) {
+        if (roundness === "sharp") {
+            return isPointInPolygon(points, x, y);
+        }
+        const polygonPoints = pointsOnBezierCurves(points, 10, 5);
+        return isPointInPolygon(polygonPoints, x, y);
+    }
+    return false;
+};
+const hitTestRoughShape = (drawable, x, y, lineThreshold) => {
+    // read operations from first opSet
+    const ops = getCurvePathOps(drawable);
+    // set start position as (0,0) just in case
+    // move operation does not exist (unlikely but it is worth safekeeping it)
+    let currentP = [0, 0];
+    return ops.some(({ op, data }, idx) => {
+        // There are only four operation types:
+        // move, bcurveTo, lineTo, and curveTo
+        if (op === "move") {
+            // change starting point
+            currentP = data;
+            // move operation does not draw anything; so, it always
+            // returns false
+        }
+        else if (op === "bcurveTo") {
+            // create points from bezier curve
+            // bezier curve stores data as a flattened array of three positions
+            // [x1, y1, x2, y2, x3, y3]
+            const p1 = [data[0], data[1]];
+            const p2 = [data[2], data[3]];
+            const p3 = [data[4], data[5]];
+            const p0 = currentP;
+            currentP = p3;
+            // check if points are on the curve
+            // cubic bezier curves require four parameters
+            // the first parameter is the last stored position (p0)
+            const retVal = pointInBezierEquation(p0, p1, p2, p3, [x, y], lineThreshold);
+            // set end point of bezier curve as the new starting point for
+            // upcoming operations as each operation is based on the last drawn
+            // position of the previous operation
+            return retVal;
+        }
+        else if (op === "lineTo") {
+            return hitTestCurveInside(drawable, x, y, "sharp");
+        }
+        else if (op === "qcurveTo") {
+            // TODO: Implement this
+            console.warn("qcurveTo is not implemented yet");
+        }
+        return false;
+    });
+};