@itwin/core-geometry 5.8.0-dev.9 → 5.8.0

package/lib/esm/curve/internalContexts/geometry2d/UnboundedHyperbola2d.js

@@ -0,0 +1,193 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { SmallSystem } from "../../../numerics/SmallSystem";
+import { Geometry } from "../../../Geometry";
+import { Point2d, Vector2d } from "../../../geometry3d/Point2dVector2d";
+import { TrigPolynomial } from "../../../numerics/Polynomials";
+import { ImplicitCurve2d } from "./ImplicitCurve2d";
+/**
+ * Internal class for hyperbola in the xy plane. The hyperbola equation in angular parameterization is
+ * * X = A + U * sec(theta) + V * tan(theta)
+ * which means that in the (skewed and scaled) coordinate system with origin at A and local u and v as
+ * multiples of U and V the implicit equation is u^2 - v^2 = 1
+ * which means the hyperbola opens along the positive and negative U axis with asymptotes at 45 degrees,
+ * i.e., along the directions (U+V) and (U-V)
+ * @internal
+ */
+export class UnboundedHyperbola2d extends ImplicitCurve2d {
+    /** The Cartesian coordinates of any center on the line. */
+    center;
+    /** The local u axis direction. */
+    vectorU;
+    /** The local v axis direction. */
+    vectorV;
+    /* Constructor - CAPTURE given center and axis vectors */
+    constructor(center, vectorU, vectorV) {
+        super();
+        this.center = center;
+        this.vectorU = vectorU;
+        this.vectorV = vectorV;
+    }
+    /** Return a clone of this hyperbola. */
+    clone() {
+        // the create method clones the inputs
+        return UnboundedHyperbola2d.createCenterAndAxisVectors(this.center, this.vectorU, this.vectorV);
+    }
+    /**
+     * Create UnboundedHyperbola2d with scale factors.
+     * * The implicit equation is (x/a)^2 - (y/b)^2 = 1
+     * @param centerX x coordinate of center
+     * @param centerY y coordinate of center
+     * @param scaleU scale factor along the U axis
+     * @param scaleV scale factor along the V axis
+     * @returns newly created hyperbola object
+     */
+    static createFromCenterAndScales(centerX, centerY, scaleU, scaleV) {
+        return new UnboundedHyperbola2d(Point2d.create(centerX, centerY), Vector2d.create(scaleU, 0), Vector2d.create(0, scaleV));
+    }
+    /**
+     * Create an UnboundedHyperbola2d from an xy object and a radius.
+     * @param center xy coordinates of center
+     * @param vectorU vector from center to to theta=0 vertex, i.e. along axis on the "inside" of the curve
+     * @param vectorV vector from center to the transverse direction
+     */
+    static createCenterAndAxisVectors(center, vectorU, vectorV) {
+        return new UnboundedHyperbola2d(Point2d.create(center.x, center.y), Vector2d.create(vectorU.x, vectorU.y), Vector2d.create(vectorV.x, vectorV.y));
+    }
+    /**
+     * Return the local (uv) coordinate of a global point.
+     * @param spacePoint point for coordinate conversion
+     */
+    globalToLocal(spacePoint) {
+        const result = Vector2d.create();
+        if (SmallSystem.linearSystem2d(this.vectorU.x, this.vectorV.x, this.vectorU.y, this.vectorV.y, spacePoint.x - this.center.x, spacePoint.y - this.center.y, result))
+            return result;
+        return undefined;
+    }
+    /**
+     * Returns xy coordinate on the hyperbola.
+     * @param radians parametric angle on the hyperbola
+     */
+    radiansToPoint2d(radians) {
+        const c = Math.cos(radians);
+        const s = Math.sin(radians);
+        if (Geometry.isSmallMetricDistance(Math.abs(c)))
+            return undefined;
+        return this.center.plus2Scaled(this.vectorU, 1.0 / c, this.vectorV, s / c);
+    }
+    /**
+     * Returns the tangent at given radians value.
+     * @param radians parametric angle on the hyperbola
+     */
+    radiansToTangentVector2d(radians) {
+        // The curve is parameterized as     X = A + U * sec(theta) + V * tan(theta)
+        // Its tangent vector (derivative) is X' = U * sec(theta) * tan(theta + V * tan^2(theta))
+        const c = Math.cos(radians);
+        const s = Math.sin(radians);
+        const cc = c * c;
+        const sec2 = Geometry.conditionalDivideCoordinate(1.0, cc);
+        const secTan = Geometry.conditionalDivideCoordinate(s, cc);
+        if (sec2 === undefined || secTan === undefined)
+            return undefined;
+        return Vector2d.createAdd2Scaled(this.vectorU, secTan, this.vectorV, sec2);
+    }
+    /**
+     * Return the implicit function value at xy.
+     * @param xy space point
+     * @returns squared local (uv) coordinates minus 1.
+     */
+    functionValue(xy) {
+        const vectorUV = this.globalToLocal(xy);
+        if (vectorUV === undefined)
+            return 0;
+        return vectorUV.x * vectorUV.x - vectorUV.y * vectorUV.y - 1.0;
+    }
+    /**
+     * Evaluate the gradient (steepest descent) of the implicit function.
+     * @param xy space point
+     * @returns gradient vector for the implicit function
+     */
+    gradient(xy) {
+        const vectorUV = this.globalToLocal(xy);
+        const result = Vector2d.create();
+        // use INVERSE of TRANSPOSE of [UV] matrix to map gradient terms
+        if (vectorUV !== undefined && SmallSystem.linearSystem2d(this.vectorU.x, this.vectorU.y, this.vectorV.x, this.vectorV.y, 2 * vectorUV.x, -2 * vectorUV.y, result))
+            return result;
+        return Vector2d.create(0, 0);
+    }
+    /**
+     * Solve for parametric radians at hyperbola points which are perpendicular projections of spacePoint.
+     * * Up to 4 solutions are possible.
+     * @param spacePoint the space point.
+     * @handler the handler to receive all the points on the curve and radians where perpendicular happens.
+     */
+    emitPerpendiculars(spacePoint, handler) {
+        const centerToSpacePoint = Vector2d.createStartEnd(spacePoint, this.center);
+        // coefficients of C and S where C and S are the unit circle points parameterized by theta.
+        // Hyperbola point at theta is X = A + U sec + V tan
+        // tangent is T = U sec tan + V sec^2
+        // vector W = spacePOint - this.center
+        // perpendicular condition for vector V from space point is
+        //     0 = T dot (spacePoint - X) = T dot (W - U sec - V tan)
+        // multiply through by cos^4 and there are various appearances of sine and cosine
+        //    in the numerator.   There are constant term, terms with one trig, and terms with product
+        //    of 2 trigs.   This fits the call list for solving "intersection" of a quadric with unit circle.
+        //
+        const coffSC = centerToSpacePoint.dotProduct(this.vectorU);
+        const coffC = centerToSpacePoint.dotProduct(this.vectorV);
+        const coffS = this.vectorU.dotProduct(this.vectorU) + this.vectorV.dotProduct(this.vectorV);
+        const coff1 = this.vectorU.dotProduct(this.vectorV);
+        const coffSS = this.vectorV.dotProduct(this.vectorU);
+        const radiansSolutions = [];
+        TrigPolynomial.solveUnitCircleImplicitQuadricIntersection(0, coffSC, coffSS, coffC, coffS, coff1, radiansSolutions);
+        for (const radians of radiansSolutions) {
+            const s = Math.sin(radians);
+            const c = Math.cos(radians);
+            const curvePoint = this.center.plus2Scaled(this.vectorU, 1.0 / c, this.vectorV, s / c);
+            handler(curvePoint, radians);
+        }
+    }
+    /** Returns true if the circle radius is near zero. */
+    // eslint-disable-next-line @itwin/prefer-get
+    isDegenerate() {
+        return undefined === this.globalToLocal(Point2d.create(0, 0));
+    }
+    /**
+     * Test if the centers and axes of the vectors are close
+     * @param other second hyperbola
+     * @param negatedAndExchangedAxesAreEqual
+     *   * if false, a strong equality test requiring both U and V vectors to match
+     *      * This strong test is a test for identical parameterization
+     *   * if true, a weak equality test that allows U and V to be negated and or exchanged.
+     *      * The weak test is a test for the same implicit set
+     * @returns true if identical to tolerance.
+     */
+    isSameHyperbola(other, negatedAndExchangedAxesAreEqual = false) {
+        const almostEqualOrNegated = (vectorU, vectorV) => {
+            if (Geometry.isSameCoordinate(vectorU.x, vectorV.x) && Geometry.isSameCoordinate(vectorU.y, vectorV.y))
+                return 1;
+            if (Geometry.isSameCoordinate(vectorU.x, -vectorV.x) && Geometry.isSameCoordinate(vectorU.y, -vectorV.y))
+                return -1;
+            return 0;
+        };
+        if (Geometry.isSamePoint2d(this.center, other.center)) {
+            if (!negatedAndExchangedAxesAreEqual)
+                return this.vectorU.isAlmostEqualMetric(other.vectorU) && this.vectorV.isAlmostEqualMetric(other.vectorV);
+            const uuCode = almostEqualOrNegated(this.vectorU, other.vectorU);
+            const vvCode = almostEqualOrNegated(this.vectorV, other.vectorV);
+            if (uuCode !== 0 && vvCode !== 0)
+                return true;
+            const uvCode = almostEqualOrNegated(this.vectorU, other.vectorV);
+            const vuCode = almostEqualOrNegated(this.vectorV, other.vectorU);
+            if (uvCode !== 0 && vuCode !== 0)
+                return true;
+        }
+        return false;
+    }
+}
+//# sourceMappingURL=UnboundedHyperbola2d.js.map

package/lib/esm/curve/internalContexts/geometry2d/UnboundedHyperbola2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedHyperbola2d.js","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedHyperbola2d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,WAAW,EAAE,MAAM,+BAA+B,CAAC;AAC5D,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AAC7C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qCAAqC,CAAC;AAExE,OAAO,EAAE,cAAc,EAAE,MAAM,+BAA+B,CAAC;AAC/D,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEpD;;;;;;;;GAQG;AACH,MAAM,OAAO,oBAAqB,SAAQ,eAAe;IACvD,2DAA2D;IACpD,MAAM,CAAU;IACvB,kCAAkC;IAC3B,OAAO,CAAW;IACzB,kCAAkC;IAC3B,OAAO,CAAW;IACzB,yDAAyD;IACzD,YAAoB,MAAe,EAAE,OAAiB,EAAE,OAAiB;QACvE,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;IACzB,CAAC;IACD,wCAAwC;IACjC,KAAK;QACV,sCAAsC;QACtC,OAAO,oBAAoB,CAAC,0BAA0B,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAClG,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,yBAAyB,CAAC,OAAe,EAAE,OAAe,EAAE,MAAc,EAAE,MAAc;QACtG,OAAO,IAAI,oBAAoB,CAC7B,OAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CACzF,CAAC;IACJ,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,0BAA0B,CAAC,MAAa,EAAE,OAAc,EAAE,OAAc;QACpF,OAAO,IAAI,oBAAoB,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAChE,QAAQ,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EACrC,QAAQ,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CACtC,CAAC;IACJ,CAAC;IACD;;;OAGG;IACI,aAAa,CAAC,UAAiB;QACpC,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QACjC,IAAI,WAAW,CAAC,cAAc,CAC5B,IAAI,CAAC,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,EAC9B,IAAI,CAAC,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,EAC9B,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,EAC1D,MAAM,CAAC;YACP,OAAO,MAAM,CAAC;QAChB,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;OAGG;IACa,gBAAgB,CAAC,OAAe;QAC9C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,IAAI,QAAQ,CAAC,qBAAqB,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7C,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,CAAC;IACD;;;OAGG;IACa,wBAAwB,CAAC,OAAe;QACtD,4EAA4E;QAC5E,yFAAyF;QACzF,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,IAAI,GAAG,QAAQ,CAAC,2BAA2B,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAC3D,MAAM,MAAM,GAAG,QAAQ,CAAC,2BAA2B,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC3D,IAAI,IAAI,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS;YAC5C,OAAO,SAAS,CAAC;QACnB,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC;IAC7E,CAAC;IACD;;;;OAIG;IACa,aAAa,CAAC,EAAS;QACrC,MAAM,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,EAAE,CAAC,CAAC;QACxC,IAAI,QAAQ,KAAK,SAAS;YACxB,OAAO,CAAC,CAAC;QACX,OAAO,QAAQ,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,GAAG,CAAC;IACjE,CAAC;IACD;;;;OAIG;IACa,QAAQ,CAAC,EAAS;QAChC,MAAM,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,EAAE,CAAC,CAAC;QACxC,MAAM,MAAM,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QACjC,gEAAgE;QAChE,IAAI,QAAQ,KAAK,SAAS,IAAI,WAAW,CAAC,cAAc,CACtD,IAAI,CAAC,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,EAC9B,IAAI,CAAC,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC,EAC9B,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,EAC/B,MAAM,CAAC;YACP,OAAO,MAAM,CAAC;QAChB,OAAO,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC;IACD;;;;;OAKG;IACa,kBAAkB,CAChC,UAAmB,EAAE,OAAkE;QAEvF,MAAM,kBAAkB,GAAG,QAAQ,CAAC,cAAc,CAAC,UAAU,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAC5E,2FAA2F;QAC3F,oDAAoD;QACpD,qCAAqC;QACrC,sCAAsC;QACtC,2DAA2D;QAC3D,6DAA6D;QAC7D,iFAAiF;QACjF,8FAA8F;QAC9F,qGAAqG;QACrG,EAAE;QACF,MAAM,MAAM,GAAG,kBAAkB,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC3D,MAAM,KAAK,GAAG,kBAAkB,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC1D,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC5F,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACrD,MAAM,gBAAgB,GAAa,EAAE,CAAC;QACtC,cAAc,CAAC,0CAA0C,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,gBAAgB,CAAC,CAAC;QACpH,KAAK,MAAM,OAAO,IAAI,gBAAgB,EAAE,CAAC;YACvC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YAC5B,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACvF,OAAO,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,sDAAsD;IACtD,6CAA6C;IAC7B,YAAY;QAC1B,OAAO,SAAS,KAAK,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAChE,CAAC;IACD;;;;;;;;;OASG;IACI,eAAe,CAAC,KAA2B,EAAE,kCAA2C,KAAK;QAClG,MAAM,oBAAoB,GAAG,CAAC,OAAiB,EAAE,OAAiB,EAAc,EAAE;YAChF,IAAI,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;gBACpG,OAAO,CAAC,CAAC;YACX,IAAI,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC;gBACtG,OAAO,CAAC,CAAC,CAAC;YACZ,OAAO,CAAC,CAAC;QACX,CAAC,CAAA;QACD,IAAI,QAAQ,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,EAAE,CAAC;YACtD,IAAI,CAAC,+BAA+B;gBAClC,OAAO,IAAI,CAAC,OAAO,CAAC,mBAAmB,CAAC,KAAK,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,OAAO,CAAC,mBAAmB,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;YAC5G,MAAM,MAAM,GAAG,oBAAoB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC,CAAC;YACjE,MAAM,MAAM,GAAG,oBAAoB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC,CAAC;YACjE,IAAI,MAAM,KAAK,CAAC,IAAI,MAAM,KAAK,CAAC;gBAC9B,OAAO,IAAI,CAAA;YACb,MAAM,MAAM,GAAG,oBAAoB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC,CAAC;YACjE,MAAM,MAAM,GAAG,oBAAoB,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,OAAO,CAAC,CAAC;YACjE,IAAI,MAAM,KAAK,CAAC,IAAI,MAAM,KAAK,CAAC;gBAC9B,OAAO,IAAI,CAAC;QAChB,CAAC;QACD,OAAO,KAAK,CAAC;IACf,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module CartesianGeometry\n */\n\nimport { SmallSystem } from \"../../../numerics/SmallSystem\";\nimport { Geometry } from \"../../../Geometry\";\nimport { Point2d, Vector2d } from \"../../../geometry3d/Point2dVector2d\";\nimport { XAndY } from \"../../../geometry3d/XYZProps\";\nimport { TrigPolynomial } from \"../../../numerics/Polynomials\";\nimport { ImplicitCurve2d } from \"./ImplicitCurve2d\";\n\n/**\n * Internal class for hyperbola in the xy plane. The hyperbola equation in angular parameterization is\n * * X = A + U * sec(theta) + V * tan(theta)\n * which means that in the (skewed and scaled) coordinate system with origin at A and local u and v as\n * multiples of U and V the implicit equation is u^2 - v^2 = 1\n * which means the hyperbola opens along the positive and negative U axis with asymptotes at 45 degrees,\n * i.e., along the directions (U+V) and (U-V)\n * @internal\n */\nexport class UnboundedHyperbola2d extends ImplicitCurve2d {\n  /** The Cartesian coordinates of any center on the line. */\n  public center: Point2d;\n  /** The local u axis direction. */\n  public vectorU: Vector2d;\n  /** The local v axis direction. */\n  public vectorV: Vector2d;\n  /* Constructor - CAPTURE given center and axis vectors */\n  private constructor(center: Point2d, vectorU: Vector2d, vectorV: Vector2d) {\n    super();\n    this.center = center;\n    this.vectorU = vectorU;\n    this.vectorV = vectorV;\n  }\n  /** Return a clone of this hyperbola. */\n  public clone(): UnboundedHyperbola2d {\n    // the create method clones the inputs\n    return UnboundedHyperbola2d.createCenterAndAxisVectors(this.center, this.vectorU, this.vectorV);\n  }\n  /**\n   * Create UnboundedHyperbola2d with scale factors.\n   * * The implicit equation is (x/a)^2 - (y/b)^2 = 1\n   * @param centerX x coordinate of center\n   * @param centerY y coordinate of center\n   * @param scaleU scale factor along the U axis\n   * @param scaleV scale factor along the V axis\n   * @returns newly created hyperbola object\n   */\n  public static createFromCenterAndScales(centerX: number, centerY: number, scaleU: number, scaleV: number): UnboundedHyperbola2d {\n    return new UnboundedHyperbola2d(\n      Point2d.create(centerX, centerY), Vector2d.create(scaleU, 0), Vector2d.create(0, scaleV),\n    );\n  }\n  /**\n   * Create an UnboundedHyperbola2d from an xy object and a radius.\n   * @param center xy coordinates of center\n   * @param vectorU vector from center to to theta=0 vertex, i.e. along axis on the \"inside\" of the curve\n   * @param vectorV vector from center to the transverse direction\n   */\n  public static createCenterAndAxisVectors(center: XAndY, vectorU: XAndY, vectorV: XAndY): UnboundedHyperbola2d {\n    return new UnboundedHyperbola2d(Point2d.create(center.x, center.y),\n      Vector2d.create(vectorU.x, vectorU.y),\n      Vector2d.create(vectorV.x, vectorV.y),\n    );\n  }\n  /**\n   * Return the local (uv) coordinate of a global point.\n   * @param spacePoint point for coordinate conversion\n   */\n  public globalToLocal(spacePoint: XAndY): Vector2d | undefined {\n    const result = Vector2d.create();\n    if (SmallSystem.linearSystem2d(\n      this.vectorU.x, this.vectorV.x,\n      this.vectorU.y, this.vectorV.y,\n      spacePoint.x - this.center.x, spacePoint.y - this.center.y,\n      result))\n      return result;\n    return undefined;\n  }\n  /**\n   * Returns xy coordinate on the hyperbola.\n   * @param radians parametric angle on the hyperbola\n   */\n  public override radiansToPoint2d(radians: number): Point2d | undefined {\n    const c = Math.cos(radians);\n    const s = Math.sin(radians);\n    if (Geometry.isSmallMetricDistance(Math.abs(c)))\n      return undefined;\n    return this.center.plus2Scaled(this.vectorU, 1.0 / c, this.vectorV, s / c);\n  }\n  /**\n   * Returns the tangent at given radians value.\n   * @param radians parametric angle on the hyperbola\n   */\n  public override radiansToTangentVector2d(radians: number): Vector2d | undefined {\n    // The curve is parameterized as     X = A + U * sec(theta) + V * tan(theta)\n    // Its tangent vector (derivative) is X' = U * sec(theta) * tan(theta + V * tan^2(theta))\n    const c = Math.cos(radians);\n    const s = Math.sin(radians);\n    const cc = c * c;\n    const sec2 = Geometry.conditionalDivideCoordinate(1.0, cc);\n    const secTan = Geometry.conditionalDivideCoordinate(s, cc);\n    if (sec2 === undefined || secTan === undefined)\n      return undefined;\n    return Vector2d.createAdd2Scaled(this.vectorU, secTan, this.vectorV, sec2);\n  }\n  /**\n   * Return the implicit function value at xy.\n   * @param xy space point\n   * @returns squared local (uv) coordinates minus 1.\n   */\n  public override functionValue(xy: XAndY): number {\n    const vectorUV = this.globalToLocal(xy);\n    if (vectorUV === undefined)\n      return 0;\n    return vectorUV.x * vectorUV.x - vectorUV.y * vectorUV.y - 1.0;\n  }\n  /**\n   * Evaluate the gradient (steepest descent) of the implicit function.\n   * @param xy space point\n   * @returns gradient vector for the implicit function\n   */\n  public override gradient(xy: XAndY): Vector2d {\n    const vectorUV = this.globalToLocal(xy);\n    const result = Vector2d.create();\n    // use INVERSE of TRANSPOSE of [UV] matrix to map gradient terms\n    if (vectorUV !== undefined && SmallSystem.linearSystem2d(\n      this.vectorU.x, this.vectorU.y,\n      this.vectorV.x, this.vectorV.y,\n      2 * vectorUV.x, -2 * vectorUV.y,\n      result))\n      return result;\n    return Vector2d.create(0, 0);\n  }\n  /**\n   * Solve for parametric radians at hyperbola points which are perpendicular projections of spacePoint.\n   * * Up to 4 solutions are possible.\n   * @param spacePoint the space point.\n   * @handler the handler to receive all the points on the curve and radians where perpendicular happens.\n   */\n  public override emitPerpendiculars(\n    spacePoint: Point2d, handler: (curvePoint: Point2d, radians: number | undefined) => any,\n  ): any {\n    const centerToSpacePoint = Vector2d.createStartEnd(spacePoint, this.center);\n    // coefficients of C and S where C and S are the unit circle points parameterized by theta.\n    // Hyperbola point at theta is X = A + U sec + V tan\n    // tangent is T = U sec tan + V sec^2\n    // vector W = spacePOint - this.center\n    // perpendicular condition for vector V from space point is\n    //     0 = T dot (spacePoint - X) = T dot (W - U sec - V tan)\n    // multiply through by cos^4 and there are various appearances of sine and cosine\n    //    in the numerator.   There are constant term, terms with one trig, and terms with product\n    //    of 2 trigs.   This fits the call list for solving \"intersection\" of a quadric with unit circle.\n    //\n    const coffSC = centerToSpacePoint.dotProduct(this.vectorU);\n    const coffC = centerToSpacePoint.dotProduct(this.vectorV);\n    const coffS = this.vectorU.dotProduct(this.vectorU) + this.vectorV.dotProduct(this.vectorV);\n    const coff1 = this.vectorU.dotProduct(this.vectorV);\n    const coffSS = this.vectorV.dotProduct(this.vectorU);\n    const radiansSolutions: number[] = [];\n    TrigPolynomial.solveUnitCircleImplicitQuadricIntersection(0, coffSC, coffSS, coffC, coffS, coff1, radiansSolutions);\n    for (const radians of radiansSolutions) {\n      const s = Math.sin(radians);\n      const c = Math.cos(radians);\n      const curvePoint = this.center.plus2Scaled(this.vectorU, 1.0 / c, this.vectorV, s / c);\n      handler(curvePoint, radians);\n    }\n  }\n  /** Returns true if the circle radius is near zero. */\n  // eslint-disable-next-line @itwin/prefer-get\n  public override isDegenerate(): boolean {\n    return undefined === this.globalToLocal(Point2d.create(0, 0));\n  }\n  /**\n   * Test if the centers and axes of the vectors are close\n   * @param other second hyperbola\n   * @param negatedAndExchangedAxesAreEqual\n   *   * if false, a strong equality test requiring both U and V vectors to match\n   *      * This strong test is a test for identical parameterization\n   *   * if true, a weak equality test that allows U and V to be negated and or exchanged.\n   *      * The weak test is a test for the same implicit set\n   * @returns true if identical to tolerance.\n   */\n  public isSameHyperbola(other: UnboundedHyperbola2d, negatedAndExchangedAxesAreEqual: boolean = false): boolean {\n    const almostEqualOrNegated = (vectorU: Vector2d, vectorV: Vector2d): -1 | 0 | 1 => {\n      if (Geometry.isSameCoordinate(vectorU.x, vectorV.x) && Geometry.isSameCoordinate(vectorU.y, vectorV.y))\n        return 1;\n      if (Geometry.isSameCoordinate(vectorU.x, -vectorV.x) && Geometry.isSameCoordinate(vectorU.y, -vectorV.y))\n        return -1;\n      return 0;\n    }\n    if (Geometry.isSamePoint2d(this.center, other.center)) {\n      if (!negatedAndExchangedAxesAreEqual)\n        return this.vectorU.isAlmostEqualMetric(other.vectorU) && this.vectorV.isAlmostEqualMetric(other.vectorV);\n      const uuCode = almostEqualOrNegated(this.vectorU, other.vectorU);\n      const vvCode = almostEqualOrNegated(this.vectorV, other.vectorV);\n      if (uuCode !== 0 && vvCode !== 0)\n        return true\n      const uvCode = almostEqualOrNegated(this.vectorU, other.vectorV);\n      const vuCode = almostEqualOrNegated(this.vectorV, other.vectorU);\n      if (uvCode !== 0 && vuCode !== 0)\n        return true;\n    }\n    return false;\n  }\n}\n"]}

package/lib/esm/curve/internalContexts/geometry2d/UnboundedLine2d.d.ts

@@ -0,0 +1,129 @@
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { Point2d, Vector2d } from "../../../geometry3d/Point2dVector2d";
+import { ImplicitCurve2d } from "./ImplicitCurve2d";
+import { XAndY } from "../../../geometry3d/XYZProps";
+/**
+ * An UnboundedLine2dByPointAndNormal represents an infinite line by a single point on the line and a normal vector.
+ * * The normal vector is NOT required to be a unit (normalized) vector.
+ * * Use method `cloneNormalizedFromOrigin` to create a line with unit normal.
+ * @internal
+ */
+export declare class UnboundedLine2dByPointAndNormal extends ImplicitCurve2d {
+    /** The Cartesian coordinates of any point on the line. */
+    point: Point2d;
+    /** The vector perpendicular to the line */
+    normal: Vector2d;
+    private constructor();
+    /** Return a clone of this line. */
+    clone(): UnboundedLine2dByPointAndNormal;
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and the normal vector.
+     * @param pointX x coordinate of a reference point on the line
+     * @param pointY y coordinate of the reference point
+     * @param normalX x component of normal vector
+     * @param normalY y component of normal vector
+     * @returns new line object.
+     */
+    static createPointXYNormalXY(pointX: number, pointY: number, normalX: number, normalY: number): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from a point on the line and a normal vector.
+     * * The xy data from the inputs is copied to the line.
+     * * i.e. the inputs are NOT captured.
+     * @param point any point on the line
+     * @param normal the normal vector
+     * @returns new line object.
+     */
+    static createPointNormal(point: Point2d, normal: Vector2d): UnboundedLine2dByPointAndNormal;
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and a vector along the line
+     * @param pointX x coordinate of a reference point on the line
+     * @param pointY y coordinate of the reference point
+     * @param directionX x component of the vector along the line
+     * @param directionY y component of the vector along the line
+     * @returns new line object or undefined.
+     */
+    static createPointXYDirectionXY(pointX: number, pointY: number, directionX: number, directionY: number): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of two points on the line.
+     * * The nominal origin is at pointA.
+     * * The line direction is from pointA to pointB, converted to a normal by createPointXYDirectionXY.
+     * @param pointAX x coordinate of first point on the line
+     * @param pointAY y coordinate of first point on the line
+     * @param pointBX x coordinate of second point on the line
+     * @param pointBY y component of second point on the line
+     * @returns new line object or undefined.
+     */
+    static createPointXYPointXY(pointAX: number, pointAY: number, pointBX: number, pointBY: number): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from two points on the line.
+     * @param pointA first point on the line
+     * @param pointB second point on the line
+     * @returns new line object or undefined.
+     */
+    static createPointPoint(pointA: XAndY, pointB: XAndY): UnboundedLine2dByPointAndNormal | undefined;
+    /** Return true if the normal vector has zero length. */
+    isDegenerate(): boolean;
+    /**
+     * Return the implicit function value at xy.
+     * @param xy space point for evaluation
+     * @returns dot product of the line normal with the vector from the line's point to the space point.
+     */
+    functionValue(xy: XAndY): number;
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     * @returns unit normal of the line.
+     */
+    gradient(_xy: XAndY): Vector2d;
+    /** Returns a vector along the line, i.e., rotated 90 degrees from the normal vector, with normal to the left. */
+    vectorAlongLine(): Vector2d;
+    /**
+     * Rotate the line normal by 90 degrees counterclockwise and normalize it
+     * @returns the unit vector, or undefined if the normal is zero length.
+     */
+    unitVectorAlongLine(): Vector2d | undefined;
+    /**
+     * Return unit vector for the line normal.
+     * @returns the unit vector, or undefined if the normal is zero length.
+     */
+    unitNormal(): Vector2d | undefined;
+    /**
+     * Drop a perpendicular from spacePoint to the line. Emit that point to the handler.
+     * @param spacePoint the space point to be projected.
+     * @handler the handler to receive the projection point.
+     */
+    emitPerpendiculars(spacePoint: Point2d, handler: (curvePoint: Point2d, radians: number | undefined) => any): any;
+    /**
+     * Return a new implicit line with its reference point given relative to newOrigin, and its perpendicular vector normalized.
+     * @return undefined if perpendicular vector has zero length.
+     * @param newOrigin origin coordinates to subtract from existing origin.
+    */
+    cloneNormalizedFromOrigin(newOrigin?: Point2d): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Return a new implicit line with given origin and unit normal based on this.
+     * * WARNING: Beware of confusing related function cloneNormalizedFromOrigin, which has different origin logic.
+     * @return undefined if perpendicular vector has zero length.
+    */
+    cloneNormalizedWithOrigin(newOrigin?: Point2d): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Return a new implicit line with its reference point shifted by given multiple of its normal vector.
+     * * The shift is applied to the point, and the normal copied unchanged-- no normalization or test for zero.
+     * @param shiftFactor multiplier for normal.
+     * @return shifted line
+    */
+    cloneShifted(shiftFactor: number): UnboundedLine2dByPointAndNormal | undefined;
+    /**
+     * Compute the intersection of two lines.
+     * * Each can be offset.
+     * * Offsets are in units of the dot products with normal vectors.
+     *   * If normal vectors are unit vectors, the offsets are simple distance.
+     *   * Otherwise the offset scaling is scaled --- use carefully.
+     * @param other the other line.
+     * @param thisOffset target function value for dot products with this line
+     * @param otherOffset target function value for dot products with other line
+     */
+    intersectUnboundedLine2dByPointAndNormalWithOffsets(other: UnboundedLine2dByPointAndNormal, thisOffset?: number, otherOffset?: number): Point2d | undefined;
+}
+//# sourceMappingURL=UnboundedLine2d.d.ts.map

package/lib/esm/curve/internalContexts/geometry2d/UnboundedLine2d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedLine2d.d.ts","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedLine2d.ts"],"names":[],"mappings":"AAIA;;GAEG;AAGH,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qCAAqC,CAAC;AACxE,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,KAAK,EAAE,MAAM,8BAA8B,CAAC;AAGrD;;;;;GAKG;AACH,qBAAa,+BAAgC,SAAQ,eAAe;IAClE,0DAA0D;IACnD,KAAK,EAAE,OAAO,CAAC;IACtB,2CAA2C;IACpC,MAAM,EAAE,QAAQ,CAAC;IAExB,OAAO;IAKP,mCAAmC;IAC5B,KAAK,IAAI,+BAA+B;IAI/C;;;;;;;OAOG;WACW,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,+BAA+B,GAAG,SAAS;IAMlJ;;;;;;;OAOG;WACW,iBAAiB,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,EAAE,QAAQ,GAAG,+BAA+B;IAGlG;;;;;;;OAOG;WACW,wBAAwB,CACpC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GACrE,+BAA+B,GAAG,SAAS;IAM9C;;;;;;;;;OASG;WACW,oBAAoB,CAChC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GACjE,+BAA+B,GAAG,SAAS;IAG9C;;;;;OAKG;WACW,gBAAgB,CAAC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,GAAG,+BAA+B,GAAG,SAAS;IAGzG,wDAAwD;IAExC,YAAY,IAAI,OAAO;IAGvC;;;;OAIG;IACa,aAAa,CAAC,EAAE,EAAE,KAAK,GAAG,MAAM;IAGhD;;;;OAIG;IACa,QAAQ,CAAC,GAAG,EAAE,KAAK,GAAG,QAAQ;IAM9C,iHAAiH;IAC1G,eAAe,IAAI,QAAQ;IAGlC;;;OAGG;IACI,mBAAmB,IAAI,QAAQ,GAAG,SAAS;IAGlD;;;OAGG;IACI,UAAU,IAAI,QAAQ,GAAG,SAAS;IAGzC;;;;OAIG;IACa,kBAAkB,CAChC,UAAU,EAAE,OAAO,EACnB,OAAO,EAAE,CAAC,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS,KAAK,GAAG,GACjE,GAAG;IAMN;;;;MAIE;IACK,yBAAyB,CAAC,SAAS,CAAC,EAAE,OAAO,GAAG,+BAA+B,GAAG,SAAS;IAUlG;;;;MAIE;IACK,yBAAyB,CAAC,SAAS,CAAC,EAAE,OAAO,GAAG,+BAA+B,GAAG,SAAS;IAUlG;;;;;MAKE;IACK,YAAY,CAAC,WAAW,EAAE,MAAM,GAAG,+BAA+B,GAAG,SAAS;IAKrF;;;;;;;;;OASG;IACI,mDAAmD,CACxD,KAAK,EAAE,+BAA+B,EAAE,UAAU,GAAE,MAAU,EAAE,WAAW,GAAE,MAAU,GACtF,OAAO,GAAG,SAAS;CAgBvB"}

package/lib/esm/curve/internalContexts/geometry2d/UnboundedLine2d.js

@@ -0,0 +1,204 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { Geometry } from "../../../Geometry";
+import { Point2d, Vector2d } from "../../../geometry3d/Point2dVector2d";
+import { ImplicitCurve2d } from "./ImplicitCurve2d";
+import { SmallSystem } from "../../../numerics/SmallSystem";
+/**
+ * An UnboundedLine2dByPointAndNormal represents an infinite line by a single point on the line and a normal vector.
+ * * The normal vector is NOT required to be a unit (normalized) vector.
+ * * Use method `cloneNormalizedFromOrigin` to create a line with unit normal.
+ * @internal
+ */
+export class UnboundedLine2dByPointAndNormal extends ImplicitCurve2d {
+    /** The Cartesian coordinates of any point on the line. */
+    point;
+    /** The vector perpendicular to the line */
+    normal;
+    /* Constructor - CAPTURE given point and normal */
+    constructor(point, normal) {
+        super();
+        this.point = point;
+        this.normal = normal;
+    }
+    /** Return a clone of this line. */
+    clone() {
+        // the create method clones the inputs
+        return UnboundedLine2dByPointAndNormal.createPointNormal(this.point, this.normal);
+    }
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and the normal vector.
+     * @param pointX x coordinate of a reference point on the line
+     * @param pointY y coordinate of the reference point
+     * @param normalX x component of normal vector
+     * @param normalY y component of normal vector
+     * @returns new line object.
+     */
+    static createPointXYNormalXY(pointX, pointY, normalX, normalY) {
+        const unitVector = Vector2d.create(normalX, normalY).normalize();
+        if (unitVector === undefined)
+            return undefined;
+        return new UnboundedLine2dByPointAndNormal(Point2d.create(pointX, pointY), unitVector);
+    }
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from a point on the line and a normal vector.
+     * * The xy data from the inputs is copied to the line.
+     * * i.e. the inputs are NOT captured.
+     * @param point any point on the line
+     * @param normal the normal vector
+     * @returns new line object.
+     */
+    static createPointNormal(point, normal) {
+        return new UnboundedLine2dByPointAndNormal(Point2d.create(point.x, point.y), Vector2d.create(normal.x, normal.y));
+    }
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and a vector along the line
+     * @param pointX x coordinate of a reference point on the line
+     * @param pointY y coordinate of the reference point
+     * @param directionX x component of the vector along the line
+     * @param directionY y component of the vector along the line
+     * @returns new line object or undefined.
+     */
+    static createPointXYDirectionXY(pointX, pointY, directionX, directionY) {
+        const unitVector = Vector2d.create(directionY, -directionX).normalize();
+        if (unitVector === undefined)
+            return undefined;
+        return new UnboundedLine2dByPointAndNormal(Point2d.create(pointX, pointY), unitVector);
+    }
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from XY parts of two points on the line.
+     * * The nominal origin is at pointA.
+     * * The line direction is from pointA to pointB, converted to a normal by createPointXYDirectionXY.
+     * @param pointAX x coordinate of first point on the line
+     * @param pointAY y coordinate of first point on the line
+     * @param pointBX x coordinate of second point on the line
+     * @param pointBY y component of second point on the line
+     * @returns new line object or undefined.
+     */
+    static createPointXYPointXY(pointAX, pointAY, pointBX, pointBY) {
+        return this.createPointXYDirectionXY(pointAX, pointAY, pointBX - pointAX, pointBY - pointAY);
+    }
+    /**
+     * Create an UnboundedLine2dByPointAndNormal from two points on the line.
+     * @param pointA first point on the line
+     * @param pointB second point on the line
+     * @returns new line object or undefined.
+     */
+    static createPointPoint(pointA, pointB) {
+        return this.createPointXYDirectionXY(pointA.x, pointA.y, pointB.x - pointA.x, pointB.y - pointA.y);
+    }
+    /** Return true if the normal vector has zero length. */
+    // eslint-disable-next-line @itwin/prefer-get
+    isDegenerate() {
+        return this.normal.isAlmostZero;
+    }
+    /**
+     * Return the implicit function value at xy.
+     * @param xy space point for evaluation
+     * @returns dot product of the line normal with the vector from the line's point to the space point.
+     */
+    functionValue(xy) {
+        return this.normal.dotProductStartEnd(this.point, xy);
+    }
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     * @returns unit normal of the line.
+     */
+    gradient(_xy) {
+        const unit = this.normal.normalize();
+        if (unit !== undefined)
+            return unit;
+        return Vector2d.create(0, 0);
+    }
+    /** Returns a vector along the line, i.e., rotated 90 degrees from the normal vector, with normal to the left. */
+    vectorAlongLine() {
+        return this.normal.rotate90CWXY();
+    }
+    /**
+     * Rotate the line normal by 90 degrees counterclockwise and normalize it
+     * @returns the unit vector, or undefined if the normal is zero length.
+     */
+    unitVectorAlongLine() {
+        return this.normal.rotate90CCWXY().normalize();
+    }
+    /**
+     * Return unit vector for the line normal.
+     * @returns the unit vector, or undefined if the normal is zero length.
+     */
+    unitNormal() {
+        return this.normal.normalize();
+    }
+    /**
+     * Drop a perpendicular from spacePoint to the line. Emit that point to the handler.
+     * @param spacePoint the space point to be projected.
+     * @handler the handler to receive the projection point.
+     */
+    emitPerpendiculars(spacePoint, handler) {
+        const fraction = Geometry.fractionOfProjectionToVectorXYXY(spacePoint.x - this.point.x, spacePoint.y - this.point.y, this.normal.x, this.normal.y);
+        handler(spacePoint.plusScaled(this.normal, -fraction), undefined);
+    }
+    /**
+     * Return a new implicit line with its reference point given relative to newOrigin, and its perpendicular vector normalized.
+     * @return undefined if perpendicular vector has zero length.
+     * @param newOrigin origin coordinates to subtract from existing origin.
+    */
+    cloneNormalizedFromOrigin(newOrigin) {
+        const unitNormal = this.normal.normalize();
+        if (unitNormal === undefined)
+            return undefined;
+        if (newOrigin === undefined)
+            return new UnboundedLine2dByPointAndNormal(this.point.clone(), Vector2d.create(unitNormal.x, unitNormal.y));
+        return new UnboundedLine2dByPointAndNormal(Point2d.create(this.point.x - newOrigin.x, this.point.y - newOrigin.y), Vector2d.create(unitNormal.x, unitNormal.y));
+    }
+    /**
+     * Return a new implicit line with given origin and unit normal based on this.
+     * * WARNING: Beware of confusing related function cloneNormalizedFromOrigin, which has different origin logic.
+     * @return undefined if perpendicular vector has zero length.
+    */
+    cloneNormalizedWithOrigin(newOrigin) {
+        const unitNormal = this.normal.normalize();
+        if (unitNormal === undefined)
+            return undefined;
+        if (newOrigin === undefined)
+            return new UnboundedLine2dByPointAndNormal(this.point, unitNormal);
+        return new UnboundedLine2dByPointAndNormal(Point2d.create(newOrigin.x, newOrigin.y), Vector2d.create(unitNormal.x, unitNormal.y));
+    }
+    /**
+     * Return a new implicit line with its reference point shifted by given multiple of its normal vector.
+     * * The shift is applied to the point, and the normal copied unchanged-- no normalization or test for zero.
+     * @param shiftFactor multiplier for normal.
+     * @return shifted line
+    */
+    cloneShifted(shiftFactor) {
+        return UnboundedLine2dByPointAndNormal.createPointNormal(this.point.plusScaled(this.normal, shiftFactor), this.normal);
+    }
+    /**
+     * Compute the intersection of two lines.
+     * * Each can be offset.
+     * * Offsets are in units of the dot products with normal vectors.
+     *   * If normal vectors are unit vectors, the offsets are simple distance.
+     *   * Otherwise the offset scaling is scaled --- use carefully.
+     * @param other the other line.
+     * @param thisOffset target function value for dot products with this line
+     * @param otherOffset target function value for dot products with other line
+     */
+    intersectUnboundedLine2dByPointAndNormalWithOffsets(other, thisOffset = 0, otherOffset = 0) {
+        // points p = (x,y) on the line satisfies equation "normal dot (p−point​) = 0" and
+        // points p = (x,y) on the  offset line satisfies equation "normal dot (p−point​) = offset"
+        // therefore:
+        // normal.x ​(x−point.x​) + normal.y ​(y−point.y​) = thisOffset
+        // normal.x ​x + normal.y ​y = thisOffset + normal dot point.x​
+        const xy = Vector2d.create();
+        if (SmallSystem.linearSystem2d(this.normal.x, this.normal.y, other.normal.x, other.normal.y, thisOffset + this.normal.dotProduct(this.point), otherOffset + other.normal.dotProduct(other.point), xy)) {
+            return Point2d.create(xy.x, xy.y);
+        }
+        return undefined;
+    }
+}
+//# sourceMappingURL=UnboundedLine2d.js.map

package/lib/esm/curve/internalContexts/geometry2d/UnboundedLine2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedLine2d.js","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedLine2d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AAC7C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qCAAqC,CAAC;AACxE,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEpD,OAAO,EAAE,WAAW,EAAE,MAAM,+BAA+B,CAAC;AAE5D;;;;;GAKG;AACH,MAAM,OAAO,+BAAgC,SAAQ,eAAe;IAClE,0DAA0D;IACnD,KAAK,CAAU;IACtB,2CAA2C;IACpC,MAAM,CAAW;IACxB,kDAAkD;IAClD,YAAoB,KAAc,EAAE,MAAgB;QAClD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;IACD,mCAAmC;IAC5B,KAAK;QACV,sCAAsC;QACtC,OAAO,+BAA+B,CAAC,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IACpF,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,qBAAqB,CAAC,MAAc,EAAE,MAAc,EAAE,OAAe,EAAE,OAAe;QAClG,MAAM,UAAU,GAAG,QAAQ,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,SAAS,EAAE,CAAC;QACjE,IAAI,UAAU,KAAK,SAAS;YAC1B,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,+BAA+B,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,UAAU,CAAC,CAAC;IACzF,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,iBAAiB,CAAC,KAAc,EAAE,MAAgB;QAC9D,OAAO,IAAI,+BAA+B,CAAC,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACpH,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,wBAAwB,CACpC,MAAc,EAAE,MAAc,EAAE,UAAkB,EAAE,UAAkB;QAEtE,MAAM,UAAU,GAAG,QAAQ,CAAC,MAAM,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC,CAAC,SAAS,EAAE,CAAC;QACxE,IAAI,UAAU,KAAK,SAAS;YAC1B,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,+BAA+B,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,UAAU,CAAC,CAAC;IACzF,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,oBAAoB,CAChC,OAAe,EAAE,OAAe,EAAE,OAAe,EAAE,OAAe;QAElE,OAAO,IAAI,CAAC,wBAAwB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,GAAG,OAAO,EAAE,OAAO,GAAG,OAAO,CAAC,CAAC;IAC/F,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,gBAAgB,CAAC,MAAa,EAAE,MAAa;QACzD,OAAO,IAAI,CAAC,wBAAwB,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACrG,CAAC;IACD,wDAAwD;IACxD,6CAA6C;IAC7B,YAAY;QAC1B,OAAO,IAAI,CAAC,MAAM,CAAC,YAAY,CAAC;IAClC,CAAC;IACD;;;;OAIG;IACa,aAAa,CAAC,EAAS;QACrC,OAAO,IAAI,CAAC,MAAM,CAAC,kBAAkB,CAAC,IAAI,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;IACxD,CAAC;IACD;;;;OAIG;IACa,QAAQ,CAAC,GAAU;QACjC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC;QACrC,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,IAAI,CAAC;QACd,OAAO,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC;IACD,iHAAiH;IAC1G,eAAe;QACpB,OAAO,IAAI,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC;IACpC,CAAC;IACD;;;OAGG;IACI,mBAAmB;QACxB,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,EAAE,CAAC,SAAS,EAAE,CAAC;IACjD,CAAC;IACD;;;OAGG;IACI,UAAU;QACf,OAAO,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC;IACjC,CAAC;IACD;;;;OAIG;IACa,kBAAkB,CAChC,UAAmB,EACnB,OAAkE;QAElE,MAAM,QAAQ,GAAG,QAAQ,CAAC,gCAAgC,CACxD,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CACvF,CAAC;QACF,OAAO,CAAC,UAAU,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,EAAE,SAAS,CAAC,CAAC;IACpE,CAAC;IACD;;;;MAIE;IACK,yBAAyB,CAAC,SAAmB;QAClD,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC;QAC3C,IAAI,UAAU,KAAK,SAAS;YAC1B,OAAO,SAAS,CAAC;QACnB,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,IAAI,+BAA+B,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,EAAE,QAAQ,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9G,OAAO,IAAI,+BAA+B,CACxC,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EACtE,QAAQ,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC;IACD;;;;MAIE;IACK,yBAAyB,CAAC,SAAmB;QAClD,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,CAAC;QAC3C,IAAI,UAAU,KAAK,SAAS;YAC1B,OAAO,SAAS,CAAC;QACnB,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,IAAI,+BAA+B,CAAC,IAAI,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;QACrE,OAAO,IAAI,+BAA+B,CACxC,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CACtF,CAAC;IACJ,CAAC;IACD;;;;;MAKE;IACK,YAAY,CAAC,WAAmB;QACrC,OAAO,+BAA+B,CAAC,iBAAiB,CACtD,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,EAC/C,IAAI,CAAC,MAAM,CAAC,CAAC;IACjB,CAAC;IACD;;;;;;;;;OASG;IACI,mDAAmD,CACxD,KAAsC,EAAE,aAAqB,CAAC,EAAE,cAAsB,CAAC;QAEvF,kFAAkF;QAClF,2FAA2F;QAC3F,aAAa;QACb,+DAA+D;QAC/D,+DAA+D;QAC/D,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;QAC7B,IAAI,WAAW,CAAC,cAAc,CAC5B,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAC5B,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC,EAC9B,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,WAAW,GAAG,KAAK,CAAC,MAAM,CAAC,UAAU,CAAC,KAAK,CAAC,KAAK,CAAC,EACnG,EAAE,CAAC,EAAE,CAAC;YACN,OAAO,OAAO,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACpC,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module CartesianGeometry\n */\n\nimport { Geometry } from \"../../../Geometry\";\nimport { Point2d, Vector2d } from \"../../../geometry3d/Point2dVector2d\";\nimport { ImplicitCurve2d } from \"./ImplicitCurve2d\";\nimport { XAndY } from \"../../../geometry3d/XYZProps\";\nimport { SmallSystem } from \"../../../numerics/SmallSystem\";\n\n/**\n * An UnboundedLine2dByPointAndNormal represents an infinite line by a single point on the line and a normal vector.\n * * The normal vector is NOT required to be a unit (normalized) vector.\n * * Use method `cloneNormalizedFromOrigin` to create a line with unit normal.\n * @internal\n */\nexport class UnboundedLine2dByPointAndNormal extends ImplicitCurve2d {\n  /** The Cartesian coordinates of any point on the line. */\n  public point: Point2d;\n  /** The vector perpendicular to the line */\n  public normal: Vector2d;\n  /* Constructor - CAPTURE given point and normal */\n  private constructor(point: Point2d, normal: Vector2d) {\n    super();\n    this.point = point;\n    this.normal = normal;\n  }\n  /** Return a clone of this line. */\n  public clone(): UnboundedLine2dByPointAndNormal {\n    // the create method clones the inputs\n    return UnboundedLine2dByPointAndNormal.createPointNormal(this.point, this.normal);\n  }\n  /**\n   * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and the normal vector.\n   * @param pointX x coordinate of a reference point on the line\n   * @param pointY y coordinate of the reference point\n   * @param normalX x component of normal vector\n   * @param normalY y component of normal vector\n   * @returns new line object.\n   */\n  public static createPointXYNormalXY(pointX: number, pointY: number, normalX: number, normalY: number): UnboundedLine2dByPointAndNormal | undefined {\n    const unitVector = Vector2d.create(normalX, normalY).normalize();\n    if (unitVector === undefined)\n      return undefined;\n    return new UnboundedLine2dByPointAndNormal(Point2d.create(pointX, pointY), unitVector);\n  }\n  /**\n   * Create an UnboundedLine2dByPointAndNormal from a point on the line and a normal vector.\n   * * The xy data from the inputs is copied to the line.\n   * * i.e. the inputs are NOT captured.\n   * @param point any point on the line\n   * @param normal the normal vector\n   * @returns new line object.\n   */\n  public static createPointNormal(point: Point2d, normal: Vector2d): UnboundedLine2dByPointAndNormal {\n    return new UnboundedLine2dByPointAndNormal(Point2d.create(point.x, point.y), Vector2d.create(normal.x, normal.y));\n  }\n  /**\n   * Create an UnboundedLine2dByPointAndNormal from XY parts of a point on the line and a vector along the line\n   * @param pointX x coordinate of a reference point on the line\n   * @param pointY y coordinate of the reference point\n   * @param directionX x component of the vector along the line\n   * @param directionY y component of the vector along the line\n   * @returns new line object or undefined.\n   */\n  public static createPointXYDirectionXY(\n    pointX: number, pointY: number, directionX: number, directionY: number\n  ): UnboundedLine2dByPointAndNormal | undefined {\n    const unitVector = Vector2d.create(directionY, -directionX).normalize();\n    if (unitVector === undefined)\n      return undefined;\n    return new UnboundedLine2dByPointAndNormal(Point2d.create(pointX, pointY), unitVector);\n  }\n  /**\n   * Create an UnboundedLine2dByPointAndNormal from XY parts of two points on the line.\n   * * The nominal origin is at pointA.\n   * * The line direction is from pointA to pointB, converted to a normal by createPointXYDirectionXY.\n   * @param pointAX x coordinate of first point on the line\n   * @param pointAY y coordinate of first point on the line\n   * @param pointBX x coordinate of second point on the line\n   * @param pointBY y component of second point on the line\n   * @returns new line object or undefined.\n   */\n  public static createPointXYPointXY(\n    pointAX: number, pointAY: number, pointBX: number, pointBY: number\n  ): UnboundedLine2dByPointAndNormal | undefined {\n    return this.createPointXYDirectionXY(pointAX, pointAY, pointBX - pointAX, pointBY - pointAY);\n  }\n  /**\n   * Create an UnboundedLine2dByPointAndNormal from two points on the line.\n   * @param pointA first point on the line\n   * @param pointB second point on the line\n   * @returns new line object or undefined.\n   */\n  public static createPointPoint(pointA: XAndY, pointB: XAndY): UnboundedLine2dByPointAndNormal | undefined {\n    return this.createPointXYDirectionXY(pointA.x, pointA.y, pointB.x - pointA.x, pointB.y - pointA.y);\n  }\n  /** Return true if the normal vector has zero length. */\n  // eslint-disable-next-line @itwin/prefer-get\n  public override isDegenerate(): boolean {\n    return this.normal.isAlmostZero;\n  }\n  /**\n   * Return the implicit function value at xy.\n   * @param xy space point for evaluation\n   * @returns dot product of the line normal with the vector from the line's point to the space point.\n   */\n  public override functionValue(xy: XAndY): number {\n    return this.normal.dotProductStartEnd(this.point, xy);\n  }\n  /**\n   * Returns gradient of the implicit function.\n   * @param xy space point\n   * @returns unit normal of the line.\n   */\n  public override gradient(_xy: XAndY): Vector2d {\n    const unit = this.normal.normalize();\n    if (unit !== undefined)\n      return unit;\n    return Vector2d.create(0, 0);\n  }\n  /** Returns a vector along the line, i.e., rotated 90 degrees from the normal vector, with normal to the left. */\n  public vectorAlongLine(): Vector2d {\n    return this.normal.rotate90CWXY();\n  }\n  /**\n   * Rotate the line normal by 90 degrees counterclockwise and normalize it\n   * @returns the unit vector, or undefined if the normal is zero length.\n   */\n  public unitVectorAlongLine(): Vector2d | undefined {\n    return this.normal.rotate90CCWXY().normalize();\n  }\n  /**\n   * Return unit vector for the line normal.\n   * @returns the unit vector, or undefined if the normal is zero length.\n   */\n  public unitNormal(): Vector2d | undefined {\n    return this.normal.normalize();\n  }\n  /**\n   * Drop a perpendicular from spacePoint to the line. Emit that point to the handler.\n   * @param spacePoint the space point to be projected.\n   * @handler the handler to receive the projection point.\n   */\n  public override emitPerpendiculars(\n    spacePoint: Point2d,\n    handler: (curvePoint: Point2d, radians: number | undefined) => any,\n  ): any {\n    const fraction = Geometry.fractionOfProjectionToVectorXYXY(\n      spacePoint.x - this.point.x, spacePoint.y - this.point.y, this.normal.x, this.normal.y,\n    );\n    handler(spacePoint.plusScaled(this.normal, -fraction), undefined);\n  }\n  /**\n   * Return a new implicit line with its reference point given relative to newOrigin, and its perpendicular vector normalized.\n   * @return undefined if perpendicular vector has zero length.\n   * @param newOrigin origin coordinates to subtract from existing origin.\n  */\n  public cloneNormalizedFromOrigin(newOrigin?: Point2d): UnboundedLine2dByPointAndNormal | undefined {\n    const unitNormal = this.normal.normalize();\n    if (unitNormal === undefined)\n      return undefined;\n    if (newOrigin === undefined)\n      return new UnboundedLine2dByPointAndNormal(this.point.clone(), Vector2d.create(unitNormal.x, unitNormal.y));\n    return new UnboundedLine2dByPointAndNormal(\n      Point2d.create(this.point.x - newOrigin.x, this.point.y - newOrigin.y),\n      Vector2d.create(unitNormal.x, unitNormal.y));\n  }\n  /**\n   * Return a new implicit line with given origin and unit normal based on this.\n   * * WARNING: Beware of confusing related function cloneNormalizedFromOrigin, which has different origin logic.\n   * @return undefined if perpendicular vector has zero length.\n  */\n  public cloneNormalizedWithOrigin(newOrigin?: Point2d): UnboundedLine2dByPointAndNormal | undefined {\n    const unitNormal = this.normal.normalize();\n    if (unitNormal === undefined)\n      return undefined;\n    if (newOrigin === undefined)\n      return new UnboundedLine2dByPointAndNormal(this.point, unitNormal);\n    return new UnboundedLine2dByPointAndNormal(\n      Point2d.create(newOrigin.x, newOrigin.y), Vector2d.create(unitNormal.x, unitNormal.y),\n    );\n  }\n  /**\n   * Return a new implicit line with its reference point shifted by given multiple of its normal vector.\n   * * The shift is applied to the point, and the normal copied unchanged-- no normalization or test for zero.\n   * @param shiftFactor multiplier for normal.\n   * @return shifted line\n  */\n  public cloneShifted(shiftFactor: number): UnboundedLine2dByPointAndNormal | undefined {\n    return UnboundedLine2dByPointAndNormal.createPointNormal(\n      this.point.plusScaled(this.normal, shiftFactor),\n      this.normal);\n  }\n  /**\n   * Compute the intersection of two lines.\n   * * Each can be offset.\n   * * Offsets are in units of the dot products with normal vectors.\n   *   * If normal vectors are unit vectors, the offsets are simple distance.\n   *   * Otherwise the offset scaling is scaled --- use carefully.\n   * @param other the other line.\n   * @param thisOffset target function value for dot products with this line\n   * @param otherOffset target function value for dot products with other line\n   */\n  public intersectUnboundedLine2dByPointAndNormalWithOffsets(\n    other: UnboundedLine2dByPointAndNormal, thisOffset: number = 0, otherOffset: number = 0,\n  ): Point2d | undefined {\n    // points p = (x,y) on the line satisfies equation \"normal dot (p−point​) = 0\" and\n    // points p = (x,y) on the  offset line satisfies equation \"normal dot (p−point​) = offset\"\n    // therefore:\n    // normal.x ​(x−point.x​) + normal.y ​(y−point.y​) = thisOffset\n    // normal.x ​x + normal.y ​y = thisOffset + normal dot point.x​\n    const xy = Vector2d.create();\n    if (SmallSystem.linearSystem2d(\n      this.normal.x, this.normal.y,\n      other.normal.x, other.normal.y,\n      thisOffset + this.normal.dotProduct(this.point), otherOffset + other.normal.dotProduct(other.point),\n      xy)) {\n      return Point2d.create(xy.x, xy.y);\n    }\n    return undefined;\n  }\n}\n"]}