@itwin/core-geometry 3.0.0-dev.76 → 3.0.0-dev.80

package/lib/esm/numerics/PolarData.d.ts

@@ -0,0 +1,42 @@
+import { GeometryQuery } from "../curve/GeometryQuery";
+import { Angle } from "../geometry3d/Angle";
+/** Enumeration of how constraints have been resolved
+ * @internal
+ */
+export declare enum ConstraintState {
+    unknown = 0,
+    singlePoint = 1,
+    impossibleValues = -1,
+    onCurve = 2
+}
+/**
+ * PolarData carries (possibly incomplete) data for converting among polar and cartesian coordinates.
+ * @internal
+ */
+export declare class PolarData {
+    private static _defaultRadius;
+    static get defaultRadius(): number;
+    /** x coordinate, possibly unknown */
+    x?: number;
+    /** y coordinate, possibly unknown */
+    y?: number;
+    /** radius, possibly unknown */
+    r?: number;
+    /** angle, possibly unknown */
+    theta?: Angle;
+    /** point, line, or arc geometry, as determined by  solveFromScalars */
+    geometry?: GeometryQuery;
+    /** enumeration of resolved state validity conditions. */
+    state?: ConstraintState;
+    /** Count the number of defined values among x,y,r, theta */
+    get numberOfConstrainedScalars(): number;
+    /** Create with any combination of known and unknown scalars. */
+    static createMixedScalars(state?: ConstraintState, x?: number, y?: number, r?: number, theta?: Angle): PolarData;
+    /** Clone the scalar data, replace the state.
+     * * Geometry is NOT cloned.
+     */
+    cloneScalarsWithState(newState: ConstraintState): PolarData;
+    /** Given a possibly incomplete set of x,y,r,theta, determine the possible completions. */
+    static solveFromScalars(known: PolarData): PolarData[];
+}
+//# sourceMappingURL=PolarData.d.ts.map

package/lib/esm/numerics/PolarData.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PolarData.d.ts","sourceRoot":"","sources":["../../../src/numerics/PolarData.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAGvD,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAG5C;;GAEG;AACH,oBAAY,eAAe;IACzB,OAAO,IAAI;IACX,WAAW,IAAI;IACf,gBAAgB,KAAK;IACrB,OAAO,IAAI;CAEZ;AAED;;;GAGG;AACH,qBAAa,SAAS;IACpB,OAAO,CAAC,MAAM,CAAC,cAAc,CAAO;IACpC,WAAkB,aAAa,IAAI,MAAM,CAAqC;IAC9E,qCAAqC;IAC9B,CAAC,CAAC,EAAE,MAAM,CAAC;IAClB,qCAAqC;IAC9B,CAAC,CAAC,EAAE,MAAM,CAAC;IAClB,+BAA+B;IACxB,CAAC,CAAC,EAAE,MAAM,CAAC;IAClB,8BAA8B;IACvB,KAAK,CAAC,EAAE,KAAK,CAAC;IACrB,uEAAuE;IAChE,QAAQ,CAAC,EAAE,aAAa,CAAC;IAChC,yDAAyD;IAClD,KAAK,CAAC,EAAE,eAAe,CAAC;IAC/B,4DAA4D;IAC5D,IAAW,0BAA0B,IAAI,MAAM,CAO9C;IACD,gEAAgE;WAClD,kBAAkB,CAAC,KAAK,CAAC,EAAE,eAAe,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,CAAC,EAAE,KAAK,GAAG,SAAS;IAUvH;;OAEG;IACI,qBAAqB,CAAC,QAAQ,EAAE,eAAe,GAAG,SAAS;IAUlE,0FAA0F;WAC5E,gBAAgB,CAAC,KAAK,EAAE,SAAS,GAAG,SAAS,EAAE;CAoF9D"}

package/lib/esm/numerics/PolarData.js

@@ -0,0 +1,162 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { Arc3d } from "../curve/Arc3d";
+import { CoordinateXYZ } from "../curve/CoordinateXYZ";
+import { LineSegment3d } from "../curve/LineSegment3d";
+import { Geometry } from "../Geometry";
+import { Angle } from "../geometry3d/Angle";
+import { Point3d } from "../geometry3d/Point3dVector3d";
+/** Enumeration of how constraints have been resolved
+ * @internal
+ */
+export var ConstraintState;
+(function (ConstraintState) {
+    ConstraintState[ConstraintState["unknown"] = 0] = "unknown";
+    ConstraintState[ConstraintState["singlePoint"] = 1] = "singlePoint";
+    ConstraintState[ConstraintState["impossibleValues"] = -1] = "impossibleValues";
+    ConstraintState[ConstraintState["onCurve"] = 2] = "onCurve";
+})(ConstraintState || (ConstraintState = {}));
+/**
+ * PolarData carries (possibly incomplete) data for converting among polar and cartesian coordinates.
+ * @internal
+ */
+export class PolarData {
+    static get defaultRadius() { return PolarData._defaultRadius; }
+    /** Count the number of defined values among x,y,r, theta */
+    get numberOfConstrainedScalars() {
+        let n = 0;
+        if (this.x !== undefined)
+            n++;
+        if (this.y !== undefined)
+            n++;
+        if (this.r !== undefined)
+            n++;
+        if (this.theta !== undefined)
+            n++;
+        return n;
+    }
+    /** Create with any combination of known and unknown scalars. */
+    static createMixedScalars(state, x, y, r, theta) {
+        const result = new PolarData();
+        result.x = x;
+        result.y = y;
+        result.r = r;
+        result.theta = theta !== undefined ? theta.clone() : undefined;
+        result.state = state;
+        return result;
+    }
+    /** Clone the scalar data, replace the state.
+     * * Geometry is NOT cloned.
+     */
+    cloneScalarsWithState(newState) {
+        const result = new PolarData();
+        result.x = this.x;
+        result.y = this.y;
+        result.r = this.r;
+        result.theta = this.theta !== undefined ? this.theta.clone() : undefined;
+        result.state = newState;
+        return result;
+    }
+    /** Given a possibly incomplete set of x,y,r,theta, determine the possible completions. */
+    static solveFromScalars(known) {
+        const result = [];
+        // .. fill out known combinations (x,y) (x,r) etc.
+        // .. any "singlePoints" gets created with no geometry, but geometry fills in at end.
+        // .. "free on curve" has curve filled in, with finite curve range controlled by PolarData.defaultRadius.
+        if (known.x !== undefined) {
+            if (known.y !== undefined) {
+                result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, known.y, Geometry.hypotenuseXY(known.x, known.y), Angle.createAtan2(known.y, known.x)));
+            }
+            else if (known.r !== undefined) {
+                const rr = known.r * known.r;
+                const xx = known.x * known.x;
+                const yy = rr - xx;
+                if (yy < 0.0) {
+                    result.push(known.cloneScalarsWithState(ConstraintState.impossibleValues));
+                }
+                else {
+                    const y = Math.sqrt(yy);
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, y, known.r, Angle.createAtan2(y, known.x)));
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, -y, known.r, Angle.createAtan2(-y, known.x)));
+                }
+            }
+            else if (known.theta !== undefined) {
+                const r = Geometry.conditionalDivideFraction(known.x, known.theta.cos());
+                // test 4*theta to isolate either positive or negative PI/2 case ...
+                if (r === undefined) {
+                    // Anywhere on y axis does has x=0, this angle . . .
+                    const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+                    q.geometry = LineSegment3d.createXYXY(0, PolarData._defaultRadius, 0, PolarData._defaultRadius);
+                }
+                else {
+                    const y = r * known.theta.sin();
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, y, r, known.theta));
+                }
+            }
+            else { // only x known --- fill out a vertical line
+                const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+                q.geometry = LineSegment3d.createXYZXYZ(known.x, PolarData._defaultRadius, 0, known.x, PolarData._defaultRadius, 0);
+            }
+        }
+        else if (known.y !== undefined) { // and we already know x is undefined .....
+            if (known.r !== undefined) {
+                const rr = known.r * known.r;
+                const yy = known.y * known.y;
+                const xx = rr - yy;
+                if (xx < 0.0) {
+                    result.push(known.cloneScalarsWithState(ConstraintState.impossibleValues));
+                }
+                else {
+                    const x = Math.sqrt(xx);
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, x, known.y, known.r, Angle.createAtan2(known.y, x)));
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, -x, known.y, known.r, Angle.createAtan2(known.y, -x)));
+                }
+            }
+            else if (known.theta !== undefined) {
+                const r = Geometry.conditionalDivideFraction(known.y, known.theta.sin());
+                // test 4*theta to isolate either positive or negative PI/2 case ...
+                if (r === undefined) {
+                    // Anywhere on y axis does has x=0, this angle . . .
+                    const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+                    q.geometry = LineSegment3d.createXYZXYZ(0, PolarData._defaultRadius, 0, 0, PolarData._defaultRadius, 0);
+                }
+                else {
+                    const x = r * known.theta.cos();
+                    result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, x, known.y, r, known.theta));
+                }
+            }
+            else { // only x known --- fill out a horizontal line
+                const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+                q.geometry = LineSegment3d.createXYZXYZ(PolarData._defaultRadius, known.y, 0, PolarData._defaultRadius, known.y, 0);
+            }
+        }
+        else if (known.r !== undefined) {
+            if (known.theta !== undefined) {
+                result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.r * known.theta.cos(), known.r * known.theta.sin(), known.r, known.theta));
+            }
+            else {
+                const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+                q.geometry = Arc3d.createXY(Point3d.create(0, 0, 0), known.r);
+                result.push(q);
+            }
+        }
+        else if (known.theta !== undefined) {
+            const q = known.cloneScalarsWithState(ConstraintState.onCurve);
+            const x = PolarData._defaultRadius * known.theta.cos();
+            const y = PolarData._defaultRadius * known.theta.sin();
+            q.geometry = LineSegment3d.createXYXY(-x, -y, x, y);
+            result.push(q);
+        }
+        // ----------------------------------------------------------------------------------------
+        // add tangible geometry to single points ...
+        for (const r of result) {
+            if (r.state === ConstraintState.singlePoint && r.geometry === undefined)
+                r.geometry = CoordinateXYZ.createXYZ(r.x, r.y, 0.0);
+        }
+        return result;
+    }
+}
+PolarData._defaultRadius = 1.0;
+//# sourceMappingURL=PolarData.js.map

package/lib/esm/numerics/PolarData.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"PolarData.js","sourceRoot":"","sources":["../../../src/numerics/PolarData.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AACvC,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAEvD,OAAO,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AACvD,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AAExD;;GAEG;AACH,MAAM,CAAN,IAAY,eAMX;AAND,WAAY,eAAe;IACzB,2DAAW,CAAA;IACX,mEAAe,CAAA;IACf,8EAAqB,CAAA;IACrB,2DAAW,CAAA;AAEb,CAAC,EANW,eAAe,KAAf,eAAe,QAM1B;AAED;;;GAGG;AACH,MAAM,OAAO,SAAS;IAEb,MAAM,KAAK,aAAa,KAAa,OAAO,SAAS,CAAC,cAAc,CAAC,CAAC,CAAC;IAa9E,4DAA4D;IAC5D,IAAW,0BAA0B;QACnC,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,IAAI,CAAC,CAAC,KAAK,SAAS;YAAE,CAAC,EAAE,CAAC;QAC9B,IAAI,IAAI,CAAC,CAAC,KAAK,SAAS;YAAE,CAAC,EAAE,CAAC;QAC9B,IAAI,IAAI,CAAC,CAAC,KAAK,SAAS;YAAE,CAAC,EAAE,CAAC;QAC9B,IAAI,IAAI,CAAC,KAAK,KAAK,SAAS;YAAE,CAAC,EAAE,CAAC;QAClC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,gEAAgE;IACzD,MAAM,CAAC,kBAAkB,CAAC,KAAuB,EAAE,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,KAAa;QACzG,MAAM,MAAM,GAAG,IAAI,SAAS,EAAE,CAAC;QAC/B,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,KAAK,GAAG,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC;QAC/D,MAAM,CAAC,KAAK,GAAG,KAAK,CAAC;QACrB,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;OAEG;IACI,qBAAqB,CAAC,QAAyB;QACpD,MAAM,MAAM,GAAG,IAAI,SAAS,EAAE,CAAC;QAC/B,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAClB,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAClB,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAClB,MAAM,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC;QACzE,MAAM,CAAC,KAAK,GAAG,QAAQ,CAAC;QACxB,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,0FAA0F;IACnF,MAAM,CAAC,gBAAgB,CAAC,KAAgB;QAC7C,MAAM,MAAM,GAAG,EAAE,CAAC;QAClB,kDAAkD;QAClD,qFAAqF;QACrF,yGAAyG;QACzG,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE;YACzB,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE;gBACzB,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,QAAQ,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACxK;iBAAM,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE;gBAChC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;gBAC7B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;gBAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBACnB,IAAI,EAAE,GAAG,GAAG,EAAE;oBACZ,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;iBAC5E;qBAAM;oBACL,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;oBACxB,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3H,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;iBAC9H;aACF;iBAAM,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE;gBACpC,MAAM,CAAC,GAAG,QAAQ,CAAC,yBAAyB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;gBACzE,oEAAoE;gBACpE,IAAI,CAAC,KAAK,SAAS,EAAE;oBACnB,oDAAoD;oBACpD,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;oBAC/D,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,UAAU,CAAC,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,CAAC,EAAE,SAAS,CAAC,cAAc,CAAC,CAAC;iBACjG;qBAAM;oBACL,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;oBAChC,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;iBACpG;aACF;iBAAM,EAAG,4CAA4C;gBACpD,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;gBAC/D,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,CAAC,CAAC,CAAC;aACrH;SACF;aAAM,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE,EAAI,2CAA2C;YAC/E,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE;gBACzB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;gBAC7B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;gBAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBACnB,IAAI,EAAE,GAAG,GAAG,EAAE;oBACZ,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,gBAAgB,CAAC,CAAC,CAAC;iBAC5E;qBAAM;oBACL,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;oBACxB,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3H,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;iBAC9H;aACF;iBAAM,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE;gBACpC,MAAM,CAAC,GAAG,QAAQ,CAAC,yBAAyB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;gBACzE,oEAAoE;gBACpE,IAAI,CAAC,KAAK,SAAS,EAAE;oBACnB,oDAAoD;oBACpD,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;oBAC/D,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,YAAY,CAAC,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,CAAC,CAAC,CAAC;iBACzG;qBAAM;oBACL,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;oBAChC,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;iBACpG;aACF;iBAAM,EAAG,8CAA8C;gBACtD,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;gBAC/D,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,YAAY,CAAC,SAAS,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,SAAS,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;aACrH;SACF;aAAM,IAAI,KAAK,CAAC,CAAC,KAAK,SAAS,EAAE;YAChC,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE;gBAC7B,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,eAAe,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,EAAE,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;aACxJ;iBAAM;gBACL,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;gBAC/D,CAAC,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC9D,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aAChB;SACF;aAAM,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE;YACpC,MAAM,CAAC,GAAG,KAAK,CAAC,qBAAqB,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;YAC/D,MAAM,CAAC,GAAG,SAAS,CAAC,cAAc,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;YACvD,MAAM,CAAC,GAAG,SAAS,CAAC,cAAc,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;YACvD,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACpD,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;SAChB;QACD,2FAA2F;QAC3F,6CAA6C;QAC7C,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE;YACtB,IAAI,CAAC,CAAC,KAAK,KAAK,eAAe,CAAC,WAAW,IAAI,CAAC,CAAC,QAAQ,KAAK,SAAS;gBACrE,CAAC,CAAC,QAAQ,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;SACvD;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;;AAnIc,wBAAc,GAAG,GAAG,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CoordinateXYZ } from \"../curve/CoordinateXYZ\";\r\nimport { GeometryQuery } from \"../curve/GeometryQuery\";\r\nimport { LineSegment3d } from \"../curve/LineSegment3d\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { Angle } from \"../geometry3d/Angle\";\r\nimport { Point3d } from \"../geometry3d/Point3dVector3d\";\r\n\r\n/** Enumeration of how constraints have been resolved\r\n * @internal\r\n */\r\nexport enum ConstraintState {\r\n  unknown = 0,\r\n  singlePoint = 1,\r\n  impossibleValues = -1,\r\n  onCurve = 2,\r\n\r\n}\r\n\r\n/**\r\n * PolarData carries (possibly incomplete) data for converting among polar and cartesian coordinates.\r\n * @internal\r\n */\r\nexport class PolarData {\r\n  private static _defaultRadius = 1.0;\r\n  public static get defaultRadius(): number { return PolarData._defaultRadius; }\r\n  /** x coordinate, possibly unknown */\r\n  public x?: number;\r\n  /** y coordinate, possibly unknown */\r\n  public y?: number;\r\n  /** radius, possibly unknown */\r\n  public r?: number;\r\n  /** angle, possibly unknown */\r\n  public theta?: Angle;\r\n  /** point, line, or arc geometry, as determined by  solveFromScalars */\r\n  public geometry?: GeometryQuery;\r\n  /** enumeration of resolved state validity conditions. */\r\n  public state?: ConstraintState;\r\n  /** Count the number of defined values among x,y,r, theta */\r\n  public get numberOfConstrainedScalars(): number {\r\n    let n = 0;\r\n    if (this.x !== undefined) n++;\r\n    if (this.y !== undefined) n++;\r\n    if (this.r !== undefined) n++;\r\n    if (this.theta !== undefined) n++;\r\n    return n;\r\n  }\r\n  /** Create with any combination of known and unknown scalars. */\r\n  public static createMixedScalars(state?: ConstraintState, x?: number, y?: number, r?: number, theta?: Angle): PolarData {\r\n    const result = new PolarData();\r\n    result.x = x;\r\n    result.y = y;\r\n    result.r = r;\r\n    result.theta = theta !== undefined ? theta.clone() : undefined;\r\n    result.state = state;\r\n    return result;\r\n  }\r\n\r\n  /** Clone the scalar data, replace the state.\r\n   * * Geometry is NOT cloned.\r\n   */\r\n  public cloneScalarsWithState(newState: ConstraintState): PolarData {\r\n    const result = new PolarData();\r\n    result.x = this.x;\r\n    result.y = this.y;\r\n    result.r = this.r;\r\n    result.theta = this.theta !== undefined ? this.theta.clone() : undefined;\r\n    result.state = newState;\r\n    return result;\r\n  }\r\n\r\n  /** Given a possibly incomplete set of x,y,r,theta, determine the possible completions. */\r\n  public static solveFromScalars(known: PolarData): PolarData[] {\r\n    const result = [];\r\n    // .. fill out known combinations (x,y) (x,r) etc.\r\n    // .. any \"singlePoints\" gets created with no geometry, but geometry fills in at end.\r\n    // .. \"free on curve\" has curve filled in, with finite curve range controlled by PolarData.defaultRadius.\r\n    if (known.x !== undefined) {\r\n      if (known.y !== undefined) {\r\n        result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, known.y, Geometry.hypotenuseXY(known.x, known.y), Angle.createAtan2(known.y, known.x)));\r\n      } else if (known.r !== undefined) {\r\n        const rr = known.r * known.r;\r\n        const xx = known.x * known.x;\r\n        const yy = rr - xx;\r\n        if (yy < 0.0) {\r\n          result.push(known.cloneScalarsWithState(ConstraintState.impossibleValues));\r\n        } else {\r\n          const y = Math.sqrt(yy);\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, y, known.r, Angle.createAtan2(y, known.x)));\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, -y, known.r, Angle.createAtan2(-y, known.x)));\r\n        }\r\n      } else if (known.theta !== undefined) {\r\n        const r = Geometry.conditionalDivideFraction(known.x, known.theta.cos());\r\n        // test 4*theta to isolate either positive or negative PI/2 case ...\r\n        if (r === undefined) {\r\n          // Anywhere on y axis does has x=0, this angle . . .\r\n          const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n          q.geometry = LineSegment3d.createXYXY(0, PolarData._defaultRadius, 0, PolarData._defaultRadius);\r\n        } else {\r\n          const y = r * known.theta.sin();\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.x, y, r, known.theta));\r\n        }\r\n      } else {  // only x known --- fill out a vertical line\r\n        const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n        q.geometry = LineSegment3d.createXYZXYZ(known.x, PolarData._defaultRadius, 0, known.x, PolarData._defaultRadius, 0);\r\n      }\r\n    } else if (known.y !== undefined) {   // and we already know x is undefined .....\r\n      if (known.r !== undefined) {\r\n        const rr = known.r * known.r;\r\n        const yy = known.y * known.y;\r\n        const xx = rr - yy;\r\n        if (xx < 0.0) {\r\n          result.push(known.cloneScalarsWithState(ConstraintState.impossibleValues));\r\n        } else {\r\n          const x = Math.sqrt(xx);\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, x, known.y, known.r, Angle.createAtan2(known.y, x)));\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, -x, known.y, known.r, Angle.createAtan2(known.y, -x)));\r\n        }\r\n      } else if (known.theta !== undefined) {\r\n        const r = Geometry.conditionalDivideFraction(known.y, known.theta.sin());\r\n        // test 4*theta to isolate either positive or negative PI/2 case ...\r\n        if (r === undefined) {\r\n          // Anywhere on y axis does has x=0, this angle . . .\r\n          const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n          q.geometry = LineSegment3d.createXYZXYZ(0, PolarData._defaultRadius, 0, 0, PolarData._defaultRadius, 0);\r\n        } else {\r\n          const x = r * known.theta.cos();\r\n          result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, x, known.y, r, known.theta));\r\n        }\r\n      } else {  // only x known --- fill out a horizontal line\r\n        const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n        q.geometry = LineSegment3d.createXYZXYZ(PolarData._defaultRadius, known.y, 0, PolarData._defaultRadius, known.y, 0);\r\n      }\r\n    } else if (known.r !== undefined) {\r\n      if (known.theta !== undefined) {\r\n        result.push(PolarData.createMixedScalars(ConstraintState.singlePoint, known.r * known.theta.cos(), known.r * known.theta.sin(), known.r, known.theta));\r\n      } else {\r\n        const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n        q.geometry = Arc3d.createXY(Point3d.create(0, 0, 0), known.r);\r\n        result.push(q);\r\n      }\r\n    } else if (known.theta !== undefined) {\r\n      const q = known.cloneScalarsWithState(ConstraintState.onCurve);\r\n      const x = PolarData._defaultRadius * known.theta.cos();\r\n      const y = PolarData._defaultRadius * known.theta.sin();\r\n      q.geometry = LineSegment3d.createXYXY(-x, -y, x, y);\r\n      result.push(q);\r\n    }\r\n    // ----------------------------------------------------------------------------------------\r\n    // add tangible geometry to single points ...\r\n    for (const r of result) {\r\n      if (r.state === ConstraintState.singlePoint && r.geometry === undefined)\r\n        r.geometry = CoordinateXYZ.createXYZ(r.x, r.y, 0.0);\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}