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

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedCircle2d.js

@@ -0,0 +1,170 @@
+"use strict";
+/*---------------------------------------------------------------------------------------------
+* 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.UnboundedCircle2dByCenterAndRadius = void 0;
+const Geometry_1 = require("../../../Geometry");
+const Point2dVector2d_1 = require("../../../geometry3d/Point2dVector2d");
+const Polynomials_1 = require("../../../numerics/Polynomials");
+const ImplicitCurve2d_1 = require("./ImplicitCurve2d");
+/**
+ * Internal class for a complete circle in the xy plane, with center and radius stored.
+ * @internal
+ */
+class UnboundedCircle2dByCenterAndRadius extends ImplicitCurve2d_1.ImplicitCurve2d {
+    /** The circle's center. */
+    center;
+    /** The circle radius. */
+    radius;
+    /* Constructor - CAPTURE given center and normal */
+    constructor(center, radius) {
+        super();
+        this.center = center;
+        this.radius = radius;
+    }
+    /** Return a clone of this circle. */
+    clone() {
+        // the create method clones the inputs
+        return UnboundedCircle2dByCenterAndRadius.createPointRadius(this.center, this.radius);
+    }
+    /** Return a clone of this circle, with radius negated. */
+    cloneNegateRadius() {
+        // the create method clones the inputs
+        return UnboundedCircle2dByCenterAndRadius.createPointRadius(this.center, -this.radius);
+    }
+    /**
+     * Create an ImplicitCircle2d from XY parts of its center and its radius.
+     * @param centerX x coordinate of center
+     * @param centerY y coordinate of center
+     * @param radius circle radius
+     * @returns newly created circle object
+     */
+    static createXYRadius(centerX, centerY, radius) {
+        return new UnboundedCircle2dByCenterAndRadius(Point2dVector2d_1.Point2d.create(centerX, centerY), radius);
+    }
+    /**
+     * Create an ImplicitCircle2d from an xy object and a radius.
+     * * Zero radius is valid.
+     * * The input coordinates are copied -- the center is NOT captured.
+     * @param center xy coordinates of center
+     * @param radius circle radius
+     * @returns newly created circle object
+     */
+    static createPointRadius(center, radius) {
+        return new UnboundedCircle2dByCenterAndRadius(Point2dVector2d_1.Point2d.create(center.x, center.y), radius);
+    }
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     * @returns squared distance to center minus squared radius.
+     */
+    functionValue(xy) {
+        return Geometry_1.Geometry.distanceSquaredXYXY(xy.x, xy.y, this.center.x, this.center.y) - this.radius * this.radius;
+    }
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     */
+    gradient(xy) {
+        return Point2dVector2d_1.Vector2d.create(2 * (xy.x - this.center.x), 2 * (xy.y - this.center.y));
+    }
+    /**
+     * Emit circle points for which a vector to the space point is perpendicular to the circle.
+     * * For a non-zero radius circle, there are two perpendiculars. The one on the side of the space point is emitted first.
+     * * For a zero radius circle, the vector from center to the space point is the only perpendicular.
+     * @param spacePoint the space point.
+     * @handler the handler to receive all the points on the curve and radians where perpendicular happens.
+     */
+    emitPerpendiculars(spacePoint, handler) {
+        // Make a vector from center to space point, scaled so length equals radius.
+        // Add and subtract it to/from center to get to circumference point.
+        const radialVector = Point2dVector2d_1.Vector2d.createStartEnd(this.center, spacePoint).scaleToLength(this.radius);
+        if (radialVector !== undefined) {
+            handler(this.center.plus(radialVector), undefined);
+            handler(this.center.minus(radialVector), undefined);
+        }
+    }
+    /**
+     * Returns true if the circle radius is near zero.
+     */
+    // eslint-disable-next-line @itwin/prefer-get
+    isDegenerate() {
+        return Geometry_1.Geometry.isSameCoordinate(this.radius, 0);
+    }
+    /**
+     * Test if the centers and radii of two circles are close.
+     * @param other second circle
+     * @returns true if identical to tolerance.
+     */
+    isSameCircle(other, negatedRadiiAreEqual) {
+        if (negatedRadiiAreEqual)
+            return Geometry_1.Geometry.isSameCoordinate(Math.abs(this.radius), Math.abs(other.radius))
+                && Geometry_1.Geometry.isSamePoint2d(this.center, other.center);
+        return Geometry_1.Geometry.isSameCoordinate(this.radius, other.radius)
+            && Geometry_1.Geometry.isSamePoint2d(this.center, other.center);
+    }
+    /**
+     * Compute intersections with another circle.
+     * @param other second circle
+     * @return array of 0, 1, or 2 points of intersection
+     */
+    intersectCircle(other) {
+        const vectorAB = Point2dVector2d_1.Vector2d.createStartEnd(this.center, other.center);
+        const unitAB = vectorAB.normalize();
+        if (unitAB === undefined)
+            return [];
+        const d = vectorAB.magnitude();
+        const rA2 = Geometry_1.Geometry.square(this.radius);
+        const rB2 = Geometry_1.Geometry.square(other.radius);
+        // if intersection points are called i1 and i2 and the mid point of the line connecting i1 and i2 is called m
+        // then we can form 2 triangles:  <this.center, m, i1> and <other.center, m, i2>
+        // now if a is distance from this.center to m and h is distance from m to i1, we get rA2 = a^2 + h^2
+        // similarly we get rB2 = (d-a)^2 + h^2 and from those 2 equations we can solve for a and h
+        const a = Geometry_1.Geometry.conditionalDivideCoordinate(rA2 - rB2 + d * d, 2 * d);
+        const points = [];
+        if (a !== undefined) {
+            const h2 = rA2 - a * a;
+            if (Math.abs(h2) < Geometry_1.Geometry.smallMetricDistanceSquared)
+                points.push(this.center.plusScaled(unitAB, a));
+            else if (h2 > 0) {
+                const h = Math.sqrt(h2);
+                points.push(this.center.addForwardLeft(a, h, unitAB));
+                points.push(this.center.addForwardLeft(a, -h, unitAB));
+            }
+        }
+        return points;
+    }
+    /**
+     * Compute intersections with a line
+     * @param line the line.
+     * @return array of 0, 1, or 2 points of intersection
+     */
+    intersectLine(line) {
+        // a point on line is X = P + frac*U where P is line.point and U is the vector along line
+        // vector from X to center is X-C = (P-C) + frac*U
+        // now defined P-C = V to get X-C = V + frac*U
+        // when distance squared from X to center is equal to this.radius squared, we have intersection:
+        // (V + frac*U) dot (V + frac*U) = r^2
+        // (U.U)*frac^2 + (2U.V)*frac + (V.V-r^2) = 0
+        // now solve for frac
+        const vectorU = line.vectorAlongLine();
+        const vectorV = Point2dVector2d_1.Vector2d.createStartEnd(this.center, line.point);
+        const uDotV = vectorU.dotProduct(vectorV);
+        const vDotV = vectorV.dotProduct(vectorV);
+        const uDotU = vectorU.dotProduct(vectorU);
+        const fractions = Polynomials_1.Degree2PowerPolynomial.solveQuadratic(uDotU, 2 * uDotV, vDotV - this.radius * this.radius);
+        const points = [];
+        if (fractions !== undefined)
+            for (const f of fractions) {
+                points.push(line.point.plusScaled(vectorU, f));
+            }
+        return points;
+    }
+}
+exports.UnboundedCircle2dByCenterAndRadius = UnboundedCircle2dByCenterAndRadius;
+//# sourceMappingURL=UnboundedCircle2d.js.map

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedCircle2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedCircle2d.js","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedCircle2d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAEH,gDAA6C;AAC7C,yEAAwE;AAExE,+DAAuE;AACvE,uDAAoD;AAGpD;;;GAGG;AACH,MAAa,kCAAmC,SAAQ,iCAAe;IACrE,2BAA2B;IACpB,MAAM,CAAU;IACvB,yBAAyB;IAClB,MAAM,CAAS;IACtB,mDAAmD;IACnD,YAAoB,MAAe,EAAE,MAAc;QACjD,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;IACD,qCAAqC;IAC9B,KAAK;QACV,sCAAsC;QACtC,OAAO,kCAAkC,CAAC,iBAAiB,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IACxF,CAAC;IACD,0DAA0D;IACnD,iBAAiB;QACtB,sCAAsC;QACtC,OAAO,kCAAkC,CAAC,iBAAiB,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACzF,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,cAAc,CAAC,OAAe,EAAE,OAAe,EAAE,MAAc;QAC3E,OAAO,IAAI,kCAAkC,CAAC,yBAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,EAAE,MAAM,CAAC,CAAC;IAC1F,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,iBAAiB,CAAC,MAAa,EAAE,MAAc;QAC3D,OAAO,IAAI,kCAAkC,CAAC,yBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IAC5F,CAAC;IACD;;;;OAIG;IACa,aAAa,CAAC,EAAS;QACrC,OAAO,mBAAQ,CAAC,mBAAmB,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;IAC5G,CAAC;IACD;;;OAGG;IACa,QAAQ,CAAC,EAAS;QAChC,OAAO,0BAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IACD;;;;;;OAMG;IACa,kBAAkB,CAChC,UAAmB,EAAE,OAAkE;QAEvF,4EAA4E;QAC5E,oEAAoE;QACpE,MAAM,YAAY,GAAG,0BAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACjG,IAAI,YAAY,KAAK,SAAS,EAAE,CAAC;YAC/B,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC,EAAE,SAAS,CAAC,CAAC;YACnD,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,YAAY,CAAC,EAAE,SAAS,CAAC,CAAC;QACtD,CAAC;IACH,CAAC;IACD;;OAEG;IACH,6CAA6C;IAC7B,YAAY;QAC1B,OAAO,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,YAAY,CAAC,KAAyC,EAAE,oBAA6B;QAC1F,IAAI,oBAAoB;YACtB,OAAO,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;mBAC1E,mBAAQ,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;QACzD,OAAO,mBAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC;eACtD,mBAAQ,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;IACzD,CAAC;IACD;;;;OAIG;IACI,eAAe,CAAC,KAAyC;QAC9D,MAAM,QAAQ,GAAG,0BAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;QACpE,MAAM,MAAM,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;QACpC,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,EAAE,CAAC;QACZ,MAAM,CAAC,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,mBAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACzC,MAAM,GAAG,GAAG,mBAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAC1C,6GAA6G;QAC7G,gFAAgF;QAChF,oGAAoG;QACpG,2FAA2F;QAC3F,MAAM,CAAC,GAAG,mBAAQ,CAAC,2BAA2B,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACzE,MAAM,MAAM,GAAG,EAAE,CAAC;QAClB,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;YACpB,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;YACvB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,mBAAQ,CAAC,0BAA0B;gBACpD,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC;iBAC5C,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;gBAChB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;gBACxB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;gBACtD,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;YACzD,CAAC;QACH,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;OAIG;IACI,aAAa,CAAC,IAAqC;QACxD,yFAAyF;QACzF,kDAAkD;QAClD,8CAA8C;QAC9C,gGAAgG;QAChG,sCAAsC;QACtC,6CAA6C;QAC7C,qBAAqB;QACrB,MAAM,OAAO,GAAG,IAAI,CAAC,eAAe,EAAE,CAAC;QACvC,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACjE,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC;QAC1C,MAAM,SAAS,GAAG,oCAAsB,CAAC,cAAc,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,EAAE,KAAK,GAAG,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;QAC7G,MAAM,MAAM,GAAG,EAAE,CAAC;QAClB,IAAI,SAAS,KAAK,SAAS;YACzB,KAAK,MAAM,CAAC,IAAI,SAAS,EAAE,CAAC;gBAC1B,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC;YACjD,CAAC;QACH,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAvJD,gFAuJC","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 { XAndY } from \"../../../geometry3d/XYZProps\";\nimport { Degree2PowerPolynomial } from \"../../../numerics/Polynomials\";\nimport { ImplicitCurve2d } from \"./ImplicitCurve2d\";\nimport { UnboundedLine2dByPointAndNormal } from \"./UnboundedLine2d\";\n\n/**\n * Internal class for a complete circle in the xy plane, with center and radius stored.\n * @internal\n */\nexport class UnboundedCircle2dByCenterAndRadius extends ImplicitCurve2d {\n  /** The circle's center. */\n  public center: Point2d;\n  /** The circle radius. */\n  public radius: number;\n  /* Constructor - CAPTURE given center and normal */\n  private constructor(center: Point2d, radius: number) {\n    super();\n    this.center = center;\n    this.radius = radius;\n  }\n  /** Return a clone of this circle. */\n  public clone(): UnboundedCircle2dByCenterAndRadius {\n    // the create method clones the inputs\n    return UnboundedCircle2dByCenterAndRadius.createPointRadius(this.center, this.radius);\n  }\n  /** Return a clone of this circle, with radius negated. */\n  public cloneNegateRadius(): UnboundedCircle2dByCenterAndRadius {\n    // the create method clones the inputs\n    return UnboundedCircle2dByCenterAndRadius.createPointRadius(this.center, -this.radius);\n  }\n  /**\n   * Create an ImplicitCircle2d from XY parts of its center and its radius.\n   * @param centerX x coordinate of center\n   * @param centerY y coordinate of center\n   * @param radius circle radius\n   * @returns newly created circle object\n   */\n  public static createXYRadius(centerX: number, centerY: number, radius: number): UnboundedCircle2dByCenterAndRadius {\n    return new UnboundedCircle2dByCenterAndRadius(Point2d.create(centerX, centerY), radius);\n  }\n  /**\n   * Create an ImplicitCircle2d from an xy object and a radius.\n   * * Zero radius is valid.\n   * * The input coordinates are copied -- the center is NOT captured.\n   * @param center xy coordinates of center\n   * @param radius circle radius\n   * @returns newly created circle object\n   */\n  public static createPointRadius(center: XAndY, radius: number): UnboundedCircle2dByCenterAndRadius {\n    return new UnboundedCircle2dByCenterAndRadius(Point2d.create(center.x, center.y), radius);\n  }\n  /**\n   * Returns gradient of the implicit function.\n   * @param xy space point\n   * @returns squared distance to center minus squared radius.\n   */\n  public override functionValue(xy: XAndY): number {\n    return Geometry.distanceSquaredXYXY(xy.x, xy.y, this.center.x, this.center.y) - this.radius * this.radius;\n  }\n  /**\n   * Returns gradient of the implicit function.\n   * @param xy space point\n   */\n  public override gradient(xy: XAndY): Vector2d {\n    return Vector2d.create(2 * (xy.x - this.center.x), 2 * (xy.y - this.center.y));\n  }\n  /**\n   * Emit circle points for which a vector to the space point is perpendicular to the circle.\n   * * For a non-zero radius circle, there are two perpendiculars. The one on the side of the space point is emitted first.\n   * * For a zero radius circle, the vector from center to the space point is the only perpendicular.\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  ) {\n    // Make a vector from center to space point, scaled so length equals radius.\n    // Add and subtract it to/from center to get to circumference point.\n    const radialVector = Vector2d.createStartEnd(this.center, spacePoint).scaleToLength(this.radius);\n    if (radialVector !== undefined) {\n      handler(this.center.plus(radialVector), undefined);\n      handler(this.center.minus(radialVector), undefined);\n    }\n  }\n  /**\n   * Returns true if the circle radius is near zero.\n   */\n  // eslint-disable-next-line @itwin/prefer-get\n  public override isDegenerate(): boolean {\n    return Geometry.isSameCoordinate(this.radius, 0);\n  }\n  /**\n   * Test if the centers and radii of two circles are close.\n   * @param other second circle\n   * @returns true if identical to tolerance.\n   */\n  public isSameCircle(other: UnboundedCircle2dByCenterAndRadius, negatedRadiiAreEqual: boolean): boolean {\n    if (negatedRadiiAreEqual)\n      return Geometry.isSameCoordinate(Math.abs(this.radius), Math.abs(other.radius))\n        && Geometry.isSamePoint2d(this.center, other.center);\n    return Geometry.isSameCoordinate(this.radius, other.radius)\n      && Geometry.isSamePoint2d(this.center, other.center);\n  }\n  /**\n   * Compute intersections with another circle.\n   * @param other second circle\n   * @return array of 0, 1, or 2 points of intersection\n   */\n  public intersectCircle(other: UnboundedCircle2dByCenterAndRadius): Point2d[] {\n    const vectorAB = Vector2d.createStartEnd(this.center, other.center);\n    const unitAB = vectorAB.normalize();\n    if (unitAB === undefined)\n      return [];\n    const d = vectorAB.magnitude();\n    const rA2 = Geometry.square(this.radius);\n    const rB2 = Geometry.square(other.radius);\n    // if intersection points are called i1 and i2 and the mid point of the line connecting i1 and i2 is called m\n    // then we can form 2 triangles:  <this.center, m, i1> and <other.center, m, i2>\n    // now if a is distance from this.center to m and h is distance from m to i1, we get rA2 = a^2 + h^2\n    // similarly we get rB2 = (d-a)^2 + h^2 and from those 2 equations we can solve for a and h\n    const a = Geometry.conditionalDivideCoordinate(rA2 - rB2 + d * d, 2 * d);\n    const points = [];\n    if (a !== undefined) {\n      const h2 = rA2 - a * a;\n      if (Math.abs(h2) < Geometry.smallMetricDistanceSquared)\n        points.push(this.center.plusScaled(unitAB, a));\n      else if (h2 > 0) {\n        const h = Math.sqrt(h2);\n        points.push(this.center.addForwardLeft(a, h, unitAB));\n        points.push(this.center.addForwardLeft(a, -h, unitAB));\n      }\n    }\n    return points;\n  }\n  /**\n   * Compute intersections with a line\n   * @param line the line.\n   * @return array of 0, 1, or 2 points of intersection\n   */\n  public intersectLine(line: UnboundedLine2dByPointAndNormal): Point2d[] {\n    // a point on line is X = P + frac*U where P is line.point and U is the vector along line\n    // vector from X to center is X-C = (P-C) + frac*U\n    // now defined P-C = V to get X-C = V + frac*U\n    // when distance squared from X to center is equal to this.radius squared, we have intersection:\n    // (V + frac*U) dot (V + frac*U) = r^2\n    // (U.U)*frac^2 + (2U.V)*frac + (V.V-r^2) = 0\n    // now solve for frac\n    const vectorU = line.vectorAlongLine();\n    const vectorV = Vector2d.createStartEnd(this.center, line.point);\n    const uDotV = vectorU.dotProduct(vectorV);\n    const vDotV = vectorV.dotProduct(vectorV);\n    const uDotU = vectorU.dotProduct(vectorU);\n    const fractions = Degree2PowerPolynomial.solveQuadratic(uDotU, 2 * uDotV, vDotV - this.radius * this.radius);\n    const points = [];\n    if (fractions !== undefined)\n      for (const f of fractions) {\n        points.push(line.point.plusScaled(vectorU, f));\n      }\n    return points;\n  }\n}\n"]}

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedEllipse2d.d.ts

@@ -0,0 +1,96 @@
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { Point2d, Vector2d } from "../../../geometry3d/Point2dVector2d";
+import { XAndY } from "../../../geometry3d/XYZProps";
+import { ImplicitCurve2d } from "./ImplicitCurve2d";
+/**
+ * Class for an ellipse in the xy plane.The ellipse equation in angular parameterization is
+ * * X = A + U * cos(theta) + V * sin(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
+ * * Note that the U and V vectors DO NOT need to be perpendicular or have any particular length relationship.
+ * * If U and V ARE perpendicular, their lengths are commonly called the major and minor axis lengths,
+ * and the major and minor axis points are in the U and V directions.
+ * * If U and V are NOT perpendicular, the major and minor axis points are at other directions.
+ * @internal
+ */
+export declare class UnboundedEllipse2d extends ImplicitCurve2d {
+    /** The Cartesian coordinates of any center on the line. */
+    center: Point2d;
+    /** The local u axis direction. */
+    vectorU: Vector2d;
+    /** The local v axis direction. */
+    vectorV: Vector2d;
+    private constructor();
+    /** Return a clone of this circle. */
+    clone(): UnboundedEllipse2d;
+    /**
+     * Create an axis-aligned UnboundedEllipse2d with axis lengths.
+     * * The implicit equation is (x/scaleU)^2 + (y/scaleV)^2 = 1
+     * @param centerX x coordinate of center
+     * @param centerY y coordinate of center
+     * @param scaleX x axis radius
+     * @param scaleY y axis radius
+     */
+    static createFromCenterAndScales(centerX: number, centerY: number, scaleX: number, scaleY: number): UnboundedEllipse2d;
+    /**
+     * Create an UnboundedEllipse2d from xy coordinates and axis vectors.
+     * @param center xy coordinates of center
+     * @param vectorU vector from center to to theta=0 point
+     * @param vectorV vector from center to the theta=90 point
+     */
+    static createCenterAndAxisVectors(center: XAndY, vectorU: XAndY, vectorV: XAndY): UnboundedEllipse2d;
+    /**
+     * Convert the coordinates of spacePoint to local coordinates relative to the ellipse vectors.
+     * @param spacePoint point for coordinate conversion
+     */
+    globalToLocal(spacePoint: XAndY): Vector2d | undefined;
+    /**
+     * Return the implicit function value at xy.
+     * @param xy space point
+     * @returns squared local (uv) coordinates minus 1
+     */
+    functionValue(xy: XAndY): number;
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     */
+    gradient(xy: XAndY): Vector2d;
+    /**
+     * Find points that are on the ellipse at foot of perpendiculars to the space point.
+     * Pass each point to the handler.
+     * * 0 to 4 handler calls 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: Point2d, handler: (curvePoint: Point2d, radians: number | undefined) => void): void;
+    /**
+     * Evaluate the curve point at parametric angle given in radians.
+     * * Note that the radians is as it appears in sin(radians) and cos( radians),
+     * and this is NOT a geometric angle in the xy plane.
+     * * The radians value is only geometric if the U and V are perpendicular and of equal length,
+     * i.e., the ellipse is a circle.
+     * @param radians angular coordinate in the ellipse.
+     */
+    radiansToPoint2d(radians: number): Point2d | undefined;
+    /**
+     * Returns the tangent at given radians value.
+     * @param radians parametric angle on the ellipse
+     */
+    radiansToTangentVector2d(radians: number): Vector2d | undefined;
+    /** Returns true if the ellipse is collapsed to a point or line.*/
+    isDegenerate(): boolean;
+    /**
+     * 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.
+     */
+    isSameEllipse(other: UnboundedEllipse2d, negatedAndExchangedAxesAreEqual?: boolean): boolean;
+}
+//# sourceMappingURL=UnboundedEllipse2d.d.ts.map

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedEllipse2d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedEllipse2d.d.ts","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedEllipse2d.ts"],"names":[],"mappings":"AAIA;;GAEG;AAIH,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qCAAqC,CAAC;AACxE,OAAO,EAAE,KAAK,EAAE,MAAM,8BAA8B,CAAC;AAErD,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEpD;;;;;;;;;;GAUG;AACH,qBAAa,kBAAmB,SAAQ,eAAe;IACrD,2DAA2D;IACpD,MAAM,EAAE,OAAO,CAAC;IACvB,kCAAkC;IAC3B,OAAO,EAAE,QAAQ,CAAC;IACzB,kCAAkC;IAC3B,OAAO,EAAE,QAAQ,CAAC;IAEzB,OAAO;IAMP,qCAAqC;IAC9B,KAAK,IAAI,kBAAkB;IAIlC;;;;;;;OAOG;WACW,yBAAyB,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,kBAAkB;IAI7H;;;;;OAKG;WACW,0BAA0B,CAAC,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,GAAG,kBAAkB;IAO3G;;;OAGG;IACI,aAAa,CAAC,UAAU,EAAE,KAAK,GAAG,QAAQ,GAAG,SAAS;IAU7D;;;;OAIG;IACa,aAAa,CAAC,EAAE,EAAE,KAAK,GAAG,MAAM;IAMhD;;;OAGG;IACa,QAAQ,CAAC,EAAE,EAAE,KAAK,GAAG,QAAQ;IAM7C;;;;;;OAMG;IACa,kBAAkB,CAChC,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS,KAAK,IAAI,GACvF,IAAI;IAwBP;;;;;;;OAOG;IACa,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,GAAG,SAAS;IAGtE;;;OAGG;IACa,wBAAwB,CAAC,OAAO,EAAE,MAAM,GAAG,QAAQ,GAAG,SAAS;IAK/E,kEAAkE;IAElD,YAAY,IAAI,OAAO;IAGvC;;;;;;;;;OASG;IACI,aAAa,CAAC,KAAK,EAAE,kBAAkB,EAAE,+BAA+B,GAAE,OAAe,GAAG,OAAO;CAsB3G"}

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedEllipse2d.js

@@ -0,0 +1,187 @@
+"use strict";
+/*---------------------------------------------------------------------------------------------
+* 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.UnboundedEllipse2d = void 0;
+const SmallSystem_1 = require("../../../numerics/SmallSystem");
+const Geometry_1 = require("../../../Geometry");
+const Point2dVector2d_1 = require("../../../geometry3d/Point2dVector2d");
+const Polynomials_1 = require("../../../numerics/Polynomials");
+const ImplicitCurve2d_1 = require("./ImplicitCurve2d");
+/**
+ * Class for an ellipse in the xy plane.The ellipse equation in angular parameterization is
+ * * X = A + U * cos(theta) + V * sin(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
+ * * Note that the U and V vectors DO NOT need to be perpendicular or have any particular length relationship.
+ * * If U and V ARE perpendicular, their lengths are commonly called the major and minor axis lengths,
+ * and the major and minor axis points are in the U and V directions.
+ * * If U and V are NOT perpendicular, the major and minor axis points are at other directions.
+ * @internal
+ */
+class UnboundedEllipse2d extends ImplicitCurve2d_1.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 circle. */
+    clone() {
+        // the create method clones the inputs
+        return UnboundedEllipse2d.createCenterAndAxisVectors(this.center, this.vectorU, this.vectorV);
+    }
+    /**
+     * Create an axis-aligned UnboundedEllipse2d with axis lengths.
+     * * The implicit equation is (x/scaleU)^2 + (y/scaleV)^2 = 1
+     * @param centerX x coordinate of center
+     * @param centerY y coordinate of center
+     * @param scaleX x axis radius
+     * @param scaleY y axis radius
+     */
+    static createFromCenterAndScales(centerX, centerY, scaleX, scaleY) {
+        return new UnboundedEllipse2d(Point2dVector2d_1.Point2d.create(centerX, centerY), Point2dVector2d_1.Vector2d.create(scaleX, 0), Point2dVector2d_1.Vector2d.create(0, scaleY));
+    }
+    /**
+     * Create an UnboundedEllipse2d from xy coordinates and axis vectors.
+     * @param center xy coordinates of center
+     * @param vectorU vector from center to to theta=0 point
+     * @param vectorV vector from center to the theta=90 point
+     */
+    static createCenterAndAxisVectors(center, vectorU, vectorV) {
+        return new UnboundedEllipse2d(Point2dVector2d_1.Point2d.create(center.x, center.y), Point2dVector2d_1.Vector2d.create(vectorU.x, vectorU.y), Point2dVector2d_1.Vector2d.create(vectorV.x, vectorV.y));
+    }
+    /**
+     * Convert the coordinates of spacePoint to local coordinates relative to the ellipse vectors.
+     * @param spacePoint point for coordinate conversion
+     */
+    globalToLocal(spacePoint) {
+        const result = Point2dVector2d_1.Vector2d.create();
+        if (SmallSystem_1.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;
+    }
+    /**
+     * 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;
+    }
+    /**
+     * Returns gradient of the implicit function.
+     * @param xy space point
+     */
+    gradient(xy) {
+        const vectorUV = this.globalToLocal(xy);
+        if (vectorUV === undefined)
+            return Point2dVector2d_1.Vector2d.create(0, 0);
+        return ImplicitCurve2d_1.ImplicitCurve2d.gradientLocalToGlobal(2 * vectorUV.x, 2 * vectorUV.y, this.vectorU, this.vectorV);
+    }
+    /**
+     * Find points that are on the ellipse at foot of perpendiculars to the space point.
+     * Pass each point to the handler.
+     * * 0 to 4 handler calls 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 vectorW = Point2dVector2d_1.Vector2d.createStartEnd(spacePoint, this.center);
+        // vector from space point to point on ellipse is
+        //    R = center + U*c + V * s - spacePoint
+        //      = W + Uc * Vs           where W = center - spacePOInt
+        // curve tangent is X' = (-U*s + V * c)
+        // dot product R and X' is
+        //    0 = (W + Uc + Vs ) dot (-Us + Vc)
+        // coefficients of c and s and combinations are:
+        const coffC = vectorW.dotProduct(this.vectorV);
+        const coffS = -vectorW.dotProduct(this.vectorU);
+        const coff1 = 0;
+        const dotUV = this.vectorU.dotProduct(this.vectorV);
+        const coffCC = dotUV;
+        const coffSC = this.vectorV.dotProduct(this.vectorV) - this.vectorU.dotProduct(this.vectorU);
+        const coffSS = -dotUV;
+        const radiansSolutions = [];
+        Polynomials_1.TrigPolynomial.solveUnitCircleImplicitQuadricIntersection(coffCC, coffSC, coffSS, coffC, coffS, coff1, radiansSolutions);
+        for (const radians of radiansSolutions) {
+            const curvePoint = this.radiansToPoint2d(radians);
+            if (curvePoint !== undefined)
+                handler(curvePoint, radians);
+        }
+    }
+    /**
+     * Evaluate the curve point at parametric angle given in radians.
+     * * Note that the radians is as it appears in sin(radians) and cos( radians),
+     * and this is NOT a geometric angle in the xy plane.
+     * * The radians value is only geometric if the U and V are perpendicular and of equal length,
+     * i.e., the ellipse is a circle.
+     * @param radians angular coordinate in the ellipse.
+     */
+    radiansToPoint2d(radians) {
+        return this.center.plus2Scaled(this.vectorU, Math.cos(radians), this.vectorV, Math.sin(radians));
+    }
+    /**
+     * Returns the tangent at given radians value.
+     * @param radians parametric angle on the ellipse
+     */
+    radiansToTangentVector2d(radians) {
+        const c = Math.cos(radians);
+        const s = Math.sin(radians);
+        return Point2dVector2d_1.Vector2d.createAdd2Scaled(this.vectorU, -s, this.vectorV, c);
+    }
+    /** Returns true if the ellipse is collapsed to a point or line.*/
+    // eslint-disable-next-line @itwin/prefer-get
+    isDegenerate() {
+        return undefined === this.globalToLocal(Point2dVector2d_1.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.
+     */
+    isSameEllipse(other, negatedAndExchangedAxesAreEqual = false) {
+        const almostEqualOrNegated = (vectorU, vectorV) => {
+            if (Geometry_1.Geometry.isSameCoordinate(vectorU.x, vectorV.x) && Geometry_1.Geometry.isSameCoordinate(vectorU.y, vectorV.y))
+                return 1;
+            if (Geometry_1.Geometry.isSameCoordinate(vectorU.x, -vectorV.x) && Geometry_1.Geometry.isSameCoordinate(vectorU.y, -vectorV.y))
+                return -1;
+            return 0;
+        };
+        if (Geometry_1.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;
+    }
+}
+exports.UnboundedEllipse2d = UnboundedEllipse2d;
+//# sourceMappingURL=UnboundedEllipse2d.js.map

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedEllipse2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedEllipse2d.js","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedEllipse2d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAEH,+DAA4D;AAC5D,gDAA6C;AAC7C,yEAAwE;AAExE,+DAA+D;AAC/D,uDAAoD;AAEpD;;;;;;;;;;GAUG;AACH,MAAa,kBAAmB,SAAQ,iCAAe;IACrD,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,qCAAqC;IAC9B,KAAK;QACV,sCAAsC;QACtC,OAAO,kBAAkB,CAAC,0BAA0B,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAChG,CAAC;IACD;;;;;;;OAOG;IACI,MAAM,CAAC,yBAAyB,CAAC,OAAe,EAAE,OAAe,EAAE,MAAc,EAAE,MAAc;QACtG,OAAO,IAAI,kBAAkB,CAAC,yBAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,CAAC,EAC5D,0BAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,0BAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;IAC5D,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,0BAA0B,CAAC,MAAa,EAAE,OAAc,EAAE,OAAc;QACpF,OAAO,IAAI,kBAAkB,CAC3B,yBAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAClC,0BAAQ,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EACrC,0BAAQ,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,0BAAQ,CAAC,MAAM,EAAE,CAAC;QACjC,IAAI,yBAAW,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;;;;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;;;OAGG;IACa,QAAQ,CAAC,EAAS;QAChC,MAAM,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,EAAE,CAAC,CAAC;QACxC,IAAI,QAAQ,KAAK,SAAS;YACxB,OAAO,0BAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/B,OAAO,iCAAe,CAAC,qBAAqB,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAC3G,CAAC;IACD;;;;;;OAMG;IACa,kBAAkB,CAChC,UAAmB,EAAE,OAAmE;QAExF,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,UAAU,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACjE,iDAAiD;QACjD,2CAA2C;QAC3C,6DAA6D;QAC7D,uCAAuC;QACvC,0BAA0B;QAC1B,uCAAuC;QACvC,gDAAgD;QAChD,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC/C,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAChD,MAAM,KAAK,GAAG,CAAC,CAAC;QAChB,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,KAAK,CAAC;QACrB,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC7F,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC;QACtB,MAAM,gBAAgB,GAAa,EAAE,CAAC;QACtC,4BAAc,CAAC,0CAA0C,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,gBAAgB,CAAC,CAAC;QACzH,KAAK,MAAM,OAAO,IAAI,gBAAgB,EAAE,CAAC;YACvC,MAAM,UAAU,GAAG,IAAI,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;YAClD,IAAI,UAAU,KAAK,SAAS;gBAC1B,OAAO,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD;;;;;;;OAOG;IACa,gBAAgB,CAAC,OAAe;QAC9C,OAAO,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC;IACnG,CAAC;IACD;;;OAGG;IACa,wBAAwB,CAAC,OAAe;QACtD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QAC5B,OAAO,0BAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACtE,CAAC;IACD,kEAAkE;IAClE,6CAA6C;IAC7B,YAAY;QAC1B,OAAO,SAAS,KAAK,IAAI,CAAC,aAAa,CAAC,yBAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAChE,CAAC;IACD;;;;;;;;;OASG;IACI,aAAa,CAAC,KAAyB,EAAE,kCAA2C,KAAK;QAC9F,MAAM,oBAAoB,GAAG,CAAC,OAAiB,EAAE,OAAiB,EAAc,EAAE;YAChF,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC;gBACpG,OAAO,CAAC,CAAC;YACX,IAAI,mBAAQ,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,mBAAQ,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,mBAAQ,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;AAzKD,gDAyKC","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 * Class for an ellipse in the xy plane.The ellipse equation in angular parameterization is\n * * X = A + U * cos(theta) + V * sin(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 * * Note that the U and V vectors DO NOT need to be perpendicular or have any particular length relationship.\n * * If U and V ARE perpendicular, their lengths are commonly called the major and minor axis lengths,\n * and the major and minor axis points are in the U and V directions.\n * * If U and V are NOT perpendicular, the major and minor axis points are at other directions.\n * @internal\n */\nexport class UnboundedEllipse2d 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 circle. */\n  public clone(): UnboundedEllipse2d {\n    // the create method clones the inputs\n    return UnboundedEllipse2d.createCenterAndAxisVectors(this.center, this.vectorU, this.vectorV);\n  }\n  /**\n   * Create an axis-aligned UnboundedEllipse2d with axis lengths.\n   * * The implicit equation is (x/scaleU)^2 + (y/scaleV)^2 = 1\n   * @param centerX x coordinate of center\n   * @param centerY y coordinate of center\n   * @param scaleX x axis radius\n   * @param scaleY y axis radius\n   */\n  public static createFromCenterAndScales(centerX: number, centerY: number, scaleX: number, scaleY: number): UnboundedEllipse2d {\n    return new UnboundedEllipse2d(Point2d.create(centerX, centerY),\n      Vector2d.create(scaleX, 0), Vector2d.create(0, scaleY));\n  }\n  /**\n   * Create an UnboundedEllipse2d from xy coordinates and axis vectors.\n   * @param center xy coordinates of center\n   * @param vectorU vector from center to to theta=0 point\n   * @param vectorV vector from center to the theta=90 point\n   */\n  public static createCenterAndAxisVectors(center: XAndY, vectorU: XAndY, vectorV: XAndY): UnboundedEllipse2d {\n    return new UnboundedEllipse2d(\n      Point2d.create(center.x, center.y),\n      Vector2d.create(vectorU.x, vectorU.y),\n      Vector2d.create(vectorV.x, vectorV.y),\n    );\n  }\n  /**\n   * Convert the coordinates of spacePoint to local coordinates relative to the ellipse vectors.\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   * 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   * Returns gradient of the implicit function.\n   * @param xy space point\n   */\n  public override gradient(xy: XAndY): Vector2d {\n    const vectorUV = this.globalToLocal(xy);\n    if (vectorUV === undefined)\n      return Vector2d.create(0, 0);\n    return ImplicitCurve2d.gradientLocalToGlobal(2 * vectorUV.x, 2 * vectorUV.y, this.vectorU, this.vectorV);\n  }\n  /**\n   * Find points that are on the ellipse at foot of perpendiculars to the space point.\n   * Pass each point to the handler.\n   * * 0 to 4 handler calls 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) => void,\n  ): void {\n    const vectorW = Vector2d.createStartEnd(spacePoint, this.center);\n    // vector from space point to point on ellipse is\n    //    R = center + U*c + V * s - spacePoint\n    //      = W + Uc * Vs           where W = center - spacePOInt\n    // curve tangent is X' = (-U*s + V * c)\n    // dot product R and X' is\n    //    0 = (W + Uc + Vs ) dot (-Us + Vc)\n    // coefficients of c and s and combinations are:\n    const coffC = vectorW.dotProduct(this.vectorV);\n    const coffS = -vectorW.dotProduct(this.vectorU);\n    const coff1 = 0;\n    const dotUV = this.vectorU.dotProduct(this.vectorV);\n    const coffCC = dotUV;\n    const coffSC = this.vectorV.dotProduct(this.vectorV) - this.vectorU.dotProduct(this.vectorU);\n    const coffSS = -dotUV;\n    const radiansSolutions: number[] = [];\n    TrigPolynomial.solveUnitCircleImplicitQuadricIntersection(coffCC, coffSC, coffSS, coffC, coffS, coff1, radiansSolutions);\n    for (const radians of radiansSolutions) {\n      const curvePoint = this.radiansToPoint2d(radians);\n      if (curvePoint !== undefined)\n        handler(curvePoint, radians);\n    }\n  }\n  /**\n   * Evaluate the curve point at parametric angle given in radians.\n   * * Note that the radians is as it appears in sin(radians) and cos( radians),\n   * and this is NOT a geometric angle in the xy plane.\n   * * The radians value is only geometric if the U and V are perpendicular and of equal length,\n   * i.e., the ellipse is a circle.\n   * @param radians angular coordinate in the ellipse.\n   */\n  public override radiansToPoint2d(radians: number): Point2d | undefined {\n    return this.center.plus2Scaled(this.vectorU, Math.cos(radians), this.vectorV, Math.sin(radians));\n  }\n  /**\n   * Returns the tangent at given radians value.\n   * @param radians parametric angle on the ellipse\n   */\n  public override radiansToTangentVector2d(radians: number): Vector2d | undefined {\n    const c = Math.cos(radians);\n    const s = Math.sin(radians);\n    return Vector2d.createAdd2Scaled(this.vectorU, -s, this.vectorV, c);\n  }\n  /** Returns true if the ellipse is collapsed to a point or line.*/\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 isSameEllipse(other: UnboundedEllipse2d, 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/cjs/curve/internalContexts/geometry2d/UnboundedHyperbola2d.d.ts

@@ -0,0 +1,91 @@
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { Point2d, Vector2d } from "../../../geometry3d/Point2dVector2d";
+import { XAndY } from "../../../geometry3d/XYZProps";
+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 declare class UnboundedHyperbola2d extends ImplicitCurve2d {
+    /** The Cartesian coordinates of any center on the line. */
+    center: Point2d;
+    /** The local u axis direction. */
+    vectorU: Vector2d;
+    /** The local v axis direction. */
+    vectorV: Vector2d;
+    private constructor();
+    /** Return a clone of this hyperbola. */
+    clone(): UnboundedHyperbola2d;
+    /**
+     * 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: number, centerY: number, scaleU: number, scaleV: number): UnboundedHyperbola2d;
+    /**
+     * 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: XAndY, vectorU: XAndY, vectorV: XAndY): UnboundedHyperbola2d;
+    /**
+     * Return the local (uv) coordinate of a global point.
+     * @param spacePoint point for coordinate conversion
+     */
+    globalToLocal(spacePoint: XAndY): Vector2d | undefined;
+    /**
+     * Returns xy coordinate on the hyperbola.
+     * @param radians parametric angle on the hyperbola
+     */
+    radiansToPoint2d(radians: number): Point2d | undefined;
+    /**
+     * Returns the tangent at given radians value.
+     * @param radians parametric angle on the hyperbola
+     */
+    radiansToTangentVector2d(radians: number): Vector2d | undefined;
+    /**
+     * Return the implicit function value at xy.
+     * @param xy space point
+     * @returns squared local (uv) coordinates minus 1.
+     */
+    functionValue(xy: XAndY): number;
+    /**
+     * Evaluate the gradient (steepest descent) of the implicit function.
+     * @param xy space point
+     * @returns gradient vector for the implicit function
+     */
+    gradient(xy: XAndY): Vector2d;
+    /**
+     * 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: Point2d, handler: (curvePoint: Point2d, radians: number | undefined) => any): any;
+    /** Returns true if the circle radius is near zero. */
+    isDegenerate(): boolean;
+    /**
+     * 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: UnboundedHyperbola2d, negatedAndExchangedAxesAreEqual?: boolean): boolean;
+}
+//# sourceMappingURL=UnboundedHyperbola2d.d.ts.map

package/lib/cjs/curve/internalContexts/geometry2d/UnboundedHyperbola2d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"UnboundedHyperbola2d.d.ts","sourceRoot":"","sources":["../../../../../src/curve/internalContexts/geometry2d/UnboundedHyperbola2d.ts"],"names":[],"mappings":"AAIA;;GAEG;AAIH,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qCAAqC,CAAC;AACxE,OAAO,EAAE,KAAK,EAAE,MAAM,8BAA8B,CAAC;AAErD,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAEpD;;;;;;;;GAQG;AACH,qBAAa,oBAAqB,SAAQ,eAAe;IACvD,2DAA2D;IACpD,MAAM,EAAE,OAAO,CAAC;IACvB,kCAAkC;IAC3B,OAAO,EAAE,QAAQ,CAAC;IACzB,kCAAkC;IAC3B,OAAO,EAAE,QAAQ,CAAC;IAEzB,OAAO;IAMP,wCAAwC;IACjC,KAAK,IAAI,oBAAoB;IAIpC;;;;;;;;OAQG;WACW,yBAAyB,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,oBAAoB;IAK/H;;;;;OAKG;WACW,0BAA0B,CAAC,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,EAAE,OAAO,EAAE,KAAK,GAAG,oBAAoB;IAM7G;;;OAGG;IACI,aAAa,CAAC,UAAU,EAAE,KAAK,GAAG,QAAQ,GAAG,SAAS;IAU7D;;;OAGG;IACa,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,GAAG,SAAS;IAOtE;;;OAGG;IACa,wBAAwB,CAAC,OAAO,EAAE,MAAM,GAAG,QAAQ,GAAG,SAAS;IAY/E;;;;OAIG;IACa,aAAa,CAAC,EAAE,EAAE,KAAK,GAAG,MAAM;IAMhD;;;;OAIG;IACa,QAAQ,CAAC,EAAE,EAAE,KAAK,GAAG,QAAQ;IAY7C;;;;;OAKG;IACa,kBAAkB,CAChC,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,UAAU,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS,KAAK,GAAG,GACtF,GAAG;IA0BN,sDAAsD;IAEtC,YAAY,IAAI,OAAO;IAGvC;;;;;;;;;OASG;IACI,eAAe,CAAC,KAAK,EAAE,oBAAoB,EAAE,+BAA+B,GAAE,OAAe,GAAG,OAAO;CAsB/G"}