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

package/lib/esm/geometry3d/CoincidentGeometryOps.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CoincidentGeometryOps.js","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAM/F,OAAO,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAE5F,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAW,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAExC;;;;;GAKG;AACH,MAAM,OAAO,uBAAuB;IAIlC,YAAoB,YAAoB,QAAQ,CAAC,mBAAmB;QAClE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,YAAoB,QAAQ,CAAC,mBAAmB;QACnE,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,0CAA0C,CAAC,MAA2B,EAAE,EAAU,EAAE,EAAU,EAC1G,MAAe,EAAE,MAAe,EAAE,OAAgB,KAAK;QACvD,IAAI,IAAI,EAAE;YACR,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;aAAM;YACL,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;QACD,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAE9E,CAAC;IAED;;OAEG;IACI,uBAAuB,CAAC,UAAmB,EAAE,MAAe,EAAE,MAAe;QAClF,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACvE,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,QAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAChE,OAAO,mBAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,QAAQ,EACrE,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD;;;;;;;OAOG;IACI,wBAAwB,CAAC,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,mBAA4B,IAAI;QACtI,MAAM,UAAU,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC3E,IAAI,OAAO,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACxD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,MAAM,UAAU,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC3E,IAAI,OAAO,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACxD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,UAAU,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC5C,UAAU,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC,CAAE,wCAAwC;QAChF,UAAU,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC5C,UAAU,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC;QACtC,IAAI,CAAC,gBAAgB;YACnB,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;QAEvE,MAAM,OAAO,GAAG,SAAS,CAAC,MAAM,CAAC,UAAU,CAAC,QAAQ,EAAE,UAAU,CAAC,SAAS,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,iBAAiB,EAAE,EAAE;YAC/B,OAAO,CAAC,2BAA2B,CAAC,GAAG,CAAC,CAAC;YACzC,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,UAAU,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACrD,MAAM,EAAE,GAAG,UAAU,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACrD,mCAAmC;YACnC,uBAAuB,CAAC,0CAA0C,CAAC,UAAU,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACzG,uBAAuB,CAAC,0CAA0C,CAAC,UAAU,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACzG,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;SACtE;aAAM;YACL,IAAI,OAAO,CAAC,WAAW,EAAE,GAAG,GAAG,EAAE;gBAC/B,IAAI,UAAU,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC3C,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;gBACD,IAAI,UAAU,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC5C,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;aACF;iBAAM;gBACL,IAAI,UAAU,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC3C,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;gBACD,IAAI,UAAU,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,CAAC,EAAE;oBAC5C,UAAU,CAAC,aAAa,EAAE,CAAC;oBAC3B,UAAU,CAAC,eAAe,EAAE,CAAC;oBAC7B,OAAO,uBAAuB,CAAC,aAAa,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;iBACtE;aACF;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;OAOG;IACK,gBAAgB,CAAC,GAAmB,EAAE,GAAmB,EAAE,YAAuB,EAAE,UAAkB,EAAE,UAAkB,EAAE,OAAgB;QAClJ,MAAM,MAAM,GAAG,UAAU,GAAG,UAAU,CAAC;QACvC,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;YACxC,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,YAAY,CAAC,EAAE,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC;YAChH,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtF,IAAI,OAAO,EAAE;gBACX,OAAO,CAAC,sBAAsB,EAAE,CAAC;gBACjC,OAAO,CAAC,sBAAsB,EAAE,CAAC;aAClC;YACD,OAAO,uBAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;SAChE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACO,gBAAgB,CAAC,MAA6C,EAAE,IAAyC;QAC/G,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,MAAM,CAAC;QAChB,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAW,EAAE,IAAW,EAAE,oBAA6B,IAAI;QAC5F,IAAI,MAA6C,CAAC;QAClD,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE;YAC1C,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,2BAA2B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC9E,IAAI,UAAU,EAAE;gBACd,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,GAAG,GAAG,QAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACpD,MAAM,KAAK,GAAG,QAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxD,IAAI,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACpC,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACrC,QAAQ,CAAC,mBAAmB,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;oBACzC,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAO,8CAA8C;oBAC/F,MAAM,cAAc,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAG,4DAA4D;oBAC7G,MAAM,gBAAgB,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC;oBACnF,MAAM,cAAc,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC;oBAC/E,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;oBACxG,MAAM,MAAM,GAAG,UAAU,CAAC,qBAAqB,CAAC,gBAAgB,EAAE,cAAc,CAAC,CAAC;oBAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC;oBAC1B,MAAM,eAAe,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;oBAChD,MAAM,UAAU,GAAG,MAAM,CAAC,iCAAiC,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;oBACjF,MAAM,aAAa,GAAG,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC;oBAChE,MAAM,UAAU,GAAG,UAAU,GAAG,aAAa,CAAC;oBAC9C,MAAM,cAAc,GAAG,SAAS,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;oBAChE,IAAI,cAAc,CAAC,iBAAiB,EAAE;wBACpC,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,cAAc,EAAE,UAAU,EAAE,UAAU,EAAE,sBAAsB,CAAC,CAAC,CAAC;oBAC5I,IAAI,UAAU,GAAG,eAAe,EAAE;wBAChC,MAAM,kBAAkB,GAAG,SAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,CAAC;wBACxG,IAAI,kBAAkB,CAAC,iBAAiB,EAAE;4BACxC,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,kBAAkB,EAAE,UAAU,EAAE,UAAU,EAAE,sBAAsB,CAAC,CAAC,CAAC;qBACjJ;iBACF;aACF;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF","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\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"./AngleSweep\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Segment1d } from \"./Segment1d\";\r\n\r\n/**\r\n * `CoincidentGeometryQuery` has methods useful in testing for overlapping geometry.\r\n * * Each instance carries tolerance information that can be reused over extended call sequences.\r\n * * These methods are expected to be called internally by CurveCurve intersection methods.\r\n * @internal\r\n */\r\nexport class CoincidentGeometryQuery {\r\n  private _vectorU?: Vector3d;\r\n  private _vectorV?: Vector3d;\r\n  private _tolerance: number;\r\n  private constructor(tolerance: number = Geometry.smallMetricDistance) {\r\n    this._tolerance = tolerance;\r\n  }\r\n  public static create(tolerance: number = Geometry.smallMetricDistance): CoincidentGeometryQuery {\r\n    return new CoincidentGeometryQuery(tolerance);\r\n  }\r\n  /**\r\n   * * Assign both the fraction and fraction1 values in the detail, possibly swapped.\r\n   * * reevaluate the points as simple interpolation between given points.\r\n   */\r\n  public static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number,\r\n    pointA: Point3d, pointB: Point3d, swap: boolean = false) {\r\n    if (swap) {\r\n      detail.fraction = f1;\r\n      detail.fraction1 = f0;\r\n    } else {\r\n      detail.fraction = f0;\r\n      detail.fraction1 = f1;\r\n    }\r\n    detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);\r\n    detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);\r\n\r\n  }\r\n\r\n  /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`\r\n   *\r\n   */\r\n  public projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail {\r\n    this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);\r\n    this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);\r\n    const uDotU = this._vectorU.dotProductXY(this._vectorU);\r\n    const uDotV = this._vectorU.dotProductXY(this._vectorV);\r\n    const fraction = Geometry.safeDivideFraction(uDotV, uDotU, 0.0);\r\n    return CurveLocationDetail.createCurveFractionPoint(undefined, fraction,\r\n      pointA.interpolate(fraction, pointB));\r\n  }\r\n  /**\r\n   * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`\r\n   * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   */\r\n  public coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds: boolean = true): CurveLocationDetailPair | undefined {\r\n    const detailAOnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);\r\n    if (pointA0.distanceXY(detailAOnB.point) > this._tolerance)\r\n      return undefined;\r\n    const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);\r\n    if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    const detailBOnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);\r\n    if (pointB0.distanceXY(detailBOnA.point) > this._tolerance)\r\n      return undefined;\r\n    const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);\r\n    if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    detailAOnB.fraction1 = detailA1OnB.fraction;\r\n    detailAOnB.point1 = detailA1OnB.point;  // capture -- detailB0OnA is not reused.\r\n    detailBOnA.fraction1 = detailB1OnA.fraction;\r\n    detailBOnA.point1 = detailB1OnA.point;\r\n    if (!restrictToBounds)\r\n      return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n\r\n    const segment = Segment1d.create(detailBOnA.fraction, detailBOnA.fraction1);\r\n    if (segment.clampDirectedTo01()) {\r\n      segment.reverseIfNeededForDeltaSign(1.0);\r\n      const f0 = segment.x0;\r\n      const f1 = segment.x1;\r\n      const h0 = detailBOnA.inverseInterpolateFraction(f0);\r\n      const h1 = detailBOnA.inverseInterpolateFraction(f1);\r\n      // recompute fractions and points..\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailBOnA, f0, f1, pointA0, pointA1);\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailAOnB, h0, h1, pointB0, pointB1);\r\n      return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n    } else {\r\n      if (segment.signedDelta() < 0.0) {\r\n        if (detailBOnA.point.isAlmostEqual(pointA0)) {\r\n          detailBOnA.collapseToStart();\r\n          detailAOnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n        if (detailBOnA.point1.isAlmostEqual(pointA1)) {\r\n          detailBOnA.collapseToEnd();\r\n          detailAOnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n      } else {\r\n        if (detailBOnA.point.isAlmostEqual(pointA1)) {\r\n          detailBOnA.collapseToStart();\r\n          detailAOnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n        if (detailBOnA.point1.isAlmostEqual(pointA0)) {\r\n          detailBOnA.collapseToEnd();\r\n          detailAOnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailBOnA, detailAOnB);\r\n        }\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Create a CurveLocationDetailPair from . . .\r\n   * @param cpA curveA\r\n   * @param cpB curve B\r\n   * @param fractionsOnA fractions of an active section of curveA\r\n   * @param fractionB0 fraction of an original containing B interval\r\n   * @param fractionB1 end fraction of an original containing B interval\r\n   */\r\n  private createDetailPair(cpA: CurvePrimitive, cpB: CurvePrimitive, fractionsOnA: Segment1d, fractionB0: number, fractionB1: number, reverse: boolean): CurveLocationDetailPair | undefined {\r\n    const deltaB = fractionB1 - fractionB0;\r\n    const g0 = Geometry.conditionalDivideFraction(fractionsOnA.x0 - fractionB0, deltaB);\r\n    const g1 = Geometry.conditionalDivideFraction(fractionsOnA.x1 - fractionB0, deltaB);\r\n    if (g0 !== undefined && g1 !== undefined) {\r\n      const detailA = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);\r\n      const detailB = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);\r\n      if (reverse) {\r\n        detailA.swapFractionsAndPoints();\r\n        detailB.swapFractionsAndPoints();\r\n      }\r\n      return CurveLocationDetailPair.createCapture(detailA, detailB);\r\n    }\r\n    return undefined;\r\n  }\r\n  private appendDetailPair(result: CurveLocationDetailPair[] | undefined, pair: CurveLocationDetailPair | undefined): CurveLocationDetailPair[] | undefined {\r\n    if (pair === undefined)\r\n      return result;\r\n    if (result === undefined)\r\n      return [pair];\r\n    result.push(pair);\r\n    return result;\r\n  }\r\n  /**\r\n   * Test if 2 arcs have coinciding portions.\r\n   * @param arcA\r\n   * @param arcB\r\n   * @param _restrictToBounds\r\n   * @return 0, 1, or 2 overlap intervals.\r\n   */\r\n  public coincidentArcIntersectionXY(arcA: Arc3d, arcB: Arc3d, _restrictToBounds: boolean = true): CurveLocationDetailPair[] | undefined {\r\n    let result: CurveLocationDetailPair[] | undefined;\r\n    if (arcA.center.isAlmostEqual(arcB.center)) {\r\n      const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);\r\n      if (matrixBToA) {\r\n        const ux = matrixBToA.at(0, 0); const uy = matrixBToA.at(1, 0);\r\n        const vx = matrixBToA.at(0, 1); const vy = matrixBToA.at(1, 1);\r\n        const ru = Geometry.hypotenuseXY(ux, uy);\r\n        const rv = Geometry.hypotenuseXY(vx, vy);\r\n        const dot = Geometry.dotProductXYXY(ux, uy, vx, vy);\r\n        const cross = Geometry.crossProductXYXY(ux, uy, vx, vy);\r\n        if (Geometry.isAlmostEqualNumber(ru, 1.0)\r\n          && Geometry.isAlmostEqualNumber(rv, 1.0)\r\n          && Geometry.isAlmostEqualNumber(0, dot)) {\r\n          const alphaB0Radians = Math.atan2(uy, ux);       // angular position of arcB 0 point in A sweep\r\n          const sweepDirection = cross > 0 ? 1.0 : -1.0;   // 1 if arcB's parameter space sweeps forward, -1 if reverse\r\n          const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;\r\n          const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;\r\n          const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;\r\n          const sweepB = AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);\r\n          const sweepA = arcA.sweep;\r\n          const fractionPeriodA = sweepA.fractionPeriod();\r\n          const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);\r\n          const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;\r\n          const fractionB1 = fractionB0 + fractionSweep;\r\n          const fractionSweepB = Segment1d.create(fractionB0, fractionB1);\r\n          if (fractionSweepB.clampDirectedTo01())\r\n            result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepB, fractionB0, fractionB1, fractionSpacesReversed));\r\n          if (fractionB1 > fractionPeriodA) {\r\n            const fractionSweepBWrap = Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA);\r\n            if (fractionSweepBWrap.clampDirectedTo01())\r\n              result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, fractionSweepBWrap, fractionB0, fractionB1, fractionSpacesReversed));\r\n          }\r\n        }\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry3d/Ellipsoid.d.ts

@@ -0,0 +1,388 @@
+/** @packageDocumentation
+ * @module CartesianGeometry
+ */
+import { CurveAndSurfaceLocationDetail } from "../bspline/SurfaceLocationDetail";
+import { Clipper } from "../clipping/ClipUtils";
+import { Arc3d } from "../curve/Arc3d";
+import { AnnounceNumberNumber, AnnounceNumberNumberCurvePrimitive } from "../curve/CurvePrimitive";
+import { Point4d } from "../geometry4d/Point4d";
+import { Angle } from "./Angle";
+import { AngleSweep } from "./AngleSweep";
+import { UVSurface } from "./GeometryHandler";
+import { LongitudeLatitudeNumber } from "./LongitudeLatitudeAltitude";
+import { Matrix3d } from "./Matrix3d";
+import { Plane3dByOriginAndUnitNormal } from "./Plane3dByOriginAndUnitNormal";
+import { Plane3dByOriginAndVectors } from "./Plane3dByOriginAndVectors";
+import { Point3d, Vector3d } from "./Point3dVector3d";
+import { Range3d } from "./Range";
+import { Ray3d } from "./Ray3d";
+import { Transform } from "./Transform";
+import { XYAndZ } from "./XYZProps";
+/**
+ * * An Ellipsoid is a (complete) unit sphere with an arbitrary (possibly skewed) `Transform` to 3d.
+ * * The (unit) sphere parameterization with respect to longitude `theta` and latitude `phi` is
+ *    * `u = cos(theta) * cos (phi)`
+ *    * `v = sin(theta) * cos(phi)`
+ *    * `w = sin(phi)`
+ *  * The sphere (u,v,w) multiply the x,y,z columns of the Ellipsoid transform.
+ * @public
+ */
+export declare class Ellipsoid implements Clipper {
+    private _transform;
+    private _unitVectorA;
+    private _unitVectorB;
+    private _workPointA;
+    private _workPointB;
+    private constructor();
+    /** Create with a clone (not capture) with given transform.
+     * * If transform is undefined, create a unit sphere.
+     */
+    static create(matrixOrTransform?: Transform | Matrix3d): Ellipsoid;
+    /**
+     * Create a transform with given center and directions, applying the radii as multipliers for the respective columns of the axes.
+     * @param center center of ellipsoid
+     * @param axes x,y,z directions are columns of this matrix
+     * @param radiusX multiplier to be applied to the x direction
+     * @param radiusY multiplier to be applied to the y direction
+     * @param radiusZ  multiplier to be applied to the z direction
+     */
+    static createCenterMatrixRadii(center: Point3d, axes: Matrix3d | undefined, radiusX: number, radiusY: number, radiusZ: number): Ellipsoid;
+    /** Return a (REFERENCE TO) the transform from world space to the mapped sphere space.
+     * * This maps coordinates "relative to the sphere" to world.
+     * * Its inverse maps world coordinates into the sphere space.
+     *   * In the sphere space, an xyz (vector from origin) with magnitude equal to 1 is ON the sphere (hence its world image is ON the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude less than 1 is INSIDE the sphere (hence its world image is INSIDE the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude greater than 1 is OUTSIDE the sphere (hence its world image is OUTSIDE the ellipsoid)
+     */
+    get transformRef(): Transform;
+    /**
+     * * Convert a world point to point within the underlying mapped sphere space.
+     *   * In the sphere space, an xyz (vector from origin) with magnitude equal to 1 is ON the sphere (hence its world image is ON the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude less than 1 is INSIDE the sphere (hence its world image is INSIDE the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude greater than 1 is OUTSIDE the sphere (hence its world image is OUTSIDE the ellipsoid)
+     * * This is undefined in the highly unusual case that the ellipsoid frame is singular.
+     */
+    worldToLocal(worldPoint: XYAndZ, result?: Point3d): Point3d | undefined;
+    /**
+     * * Convert a point within the underlying mapped sphere space to world coordinates.
+     *   * In the sphere space, an xyz (vector from origin) with magnitude equal to 1 is ON the sphere (hence its world image is ON the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude less than 1 is INSIDE the sphere (hence its world image is INSIDE the ellipsoid)
+     *   * In the sphere space, an xyz (vector from origin) with magnitude greater than 1 is OUTSIDE the sphere (hence its world image is OUTSIDE the ellipsoid)
+     */
+    localToWorld(localPoint: XYAndZ, result?: Point3d): Point3d;
+    /** return a clone with same coordinates */
+    clone(): Ellipsoid;
+    /** test equality of the 4 points */
+    isAlmostEqual(other: Ellipsoid): boolean;
+    /** Apply the transform to each point */
+    tryTransformInPlace(transform: Transform): boolean;
+    /**
+     * return a cloned and transformed ellipsoid.
+     * @param transform
+     */
+    cloneTransformed(transform: Transform): Ellipsoid | undefined;
+    /** Find the closest point of the (patch of the) ellipsoid.
+     * * In general there are multiple points where a space point projects onto an ellipse.
+     * * This searches for only one point, using heuristics which are reliable for points close to the surface but not for points distant from highly skewed ellipsoid
+     */
+    projectPointToSurface(spacePoint: Point3d): LongitudeLatitudeNumber | undefined;
+    /** Find the silhouette of the ellipsoid as viewed from a homogeneous eyepoint.
+     * * Returns undefined if the eyepoint is inside the ellipsoid
+     */
+    silhouetteArc(eyePoint: Point4d): Arc3d | undefined;
+    /** Compute intersections with a ray.
+     * * Return the number of intersections
+     * * Fill any combinations of arrays of
+     *    * rayFractions = fractions along the ray
+     *    * xyz = xyz intersection coordinates points in space
+     *    * thetaPhiRadians = sphere longitude and latitude in radians.
+     * * For each optional array, caller must of course initialize an array (usually empty)
+     * * return 0 if ray length is too small.
+     */
+    intersectRay(ray: Ray3d, rayFractions: number[] | undefined, xyz: Point3d[] | undefined, thetaPhiRadians: LongitudeLatitudeNumber[] | undefined): number;
+    /** Return the range of a uv-aligned patch of the sphere. */
+    patchRangeStartEndRadians(theta0Radians: number, theta1Radians: number, phi0Radians: number, phi1Radians: number, result?: Range3d): Range3d;
+    /**
+     * Evaluate a point on the ellipsoid at angles give in radians.
+     * @param thetaRadians longitude, in radians
+     * @param phiRadians latitude, in radians
+     * @param result optional point result
+     */
+    radiansToPoint(thetaRadians: number, phiRadians: number, result?: Point3d): Point3d;
+    /**
+     * * For a given pair of points on an ellipsoid, construct an arc (possibly elliptical) which
+     *   * passes through both points
+     *   * is completely within the ellipsoid surface
+     *   * has its centerEvaluate a point on the ellipsoid at angles give in radians.
+     * * If the ellipsoid is a sphere, this is the shortest great-circle arc between the two points.
+     * * If the ellipsoid is not a sphere, this is close to but not precisely the shortest path.
+     * @param thetaARadians longitude, in radians, for pointA
+     * @param phiARadians latitude, in radians, for pointA
+     * @param thetaBRadians longitude, in radians, for pointB
+     * @param phiBRadians latitude, in radians, for pointB
+     * @param result optional preallocated result
+     */
+    radiansPairToGreatArc(thetaARadians: number, phiARadians: number, thetaBRadians: number, phiBRadians: number, result?: Arc3d): Arc3d | undefined;
+    /**
+     * See radiansPairToGreatArc, which does this computation with positions from `angleA` and `angleB` directly as radians
+     */
+    anglePairToGreatArc(angleA: LongitudeLatitudeNumber, angleB: LongitudeLatitudeNumber, result?: Arc3d): Arc3d | undefined;
+    /**
+     * Construct an arc for the section cut of a plane with the ellipsoid.
+     * * this is undefined if the plane does not intersect the ellipsoid.
+     */
+    createPlaneSection(plane: Plane3dByOriginAndUnitNormal): Arc3d | undefined;
+    /**
+     * Construct an arc which
+     *  * start at pointA (defined by its angle position)
+     *  * ends at pointB (defined by its angle position)
+     *  * contains the 3rd vector as an in-plane point.
+     */
+    createSectionArcPointPointVectorInPlane(pointAnglesA: LongitudeLatitudeNumber, pointAnglesB: LongitudeLatitudeNumber, inPlaneVector: Vector3d, result?: Arc3d): Arc3d | undefined;
+    /**
+     * * For a given pair of points on an ellipsoid, construct another ellipsoid
+     *   * touches the same xyz points in space
+     *   * has transformation modified so that the original two points are on the equator.
+     * * Note that except for true sphere inputs, the result axes can be both non-perpendicular axes and of different lengths.
+     * @param thetaARadians longitude, in radians, for pointA
+     * @param phiARadians latitude, in radians, for pointA
+     * @param thetaBRadians longitude, in radians, for pointB
+     * @param phiBRadians latitude, in radians, for pointB
+     * @param result optional preallocated result
+     */
+    radiansPairToEquatorialEllipsoid(thetaARadians: number, phiARadians: number, thetaBRadians: number, phiBRadians: number, result?: Ellipsoid): Ellipsoid | undefined;
+    /**
+     * Return an arc (circular or elliptical) at constant longitude
+     * @param longitude (strongly typed) longitude
+     * @param latitude latitude sweep angles
+     * @param result
+     */
+    constantLongitudeArc(longitude: Angle, latitudeSweep: AngleSweep, result?: Arc3d): Arc3d | undefined;
+    /**
+     * Return an arc (circular or elliptical) at constant longitude
+     * @param latitude sweep angles
+     * @param latitude (strongly typed) latitude
+     * @param result
+     */
+    constantLatitudeArc(longitudeSweep: AngleSweep, latitude: Angle, result?: Arc3d): Arc3d | undefined;
+    /**
+     * * create a section arc with and end at positions A and B, and in plane with the normal at a fractional
+     *    interpolation between.
+     * @param angleA start point of arc (given as angles on this ellipsoid)
+     * @param intermediateNormalFraction
+     * @param angleB end point of arc (given as angles on this ellipsoid)
+     */
+    sectionArcWithIntermediateNormal(angleA: LongitudeLatitudeNumber, intermediateNormalFraction: number, angleB: LongitudeLatitudeNumber): Arc3d;
+    /**
+     * Evaluate a point and derivatives with respect to angle on the ellipsoid at angles give in radians.
+     * * "u direction" vector of the returned plane is derivative with respect to longitude.
+     * * "v direction" vector fo the returned plane is derivative with respect ot latitude.
+     * @param thetaRadians longitude, in radians
+     * @param phiRadians latitude, in radians
+     * @param applyCosPhiFactor selector for handling of theta (around equator derivative)
+     *   * if true, compute the properly scaled derivative, which goes to zero at the poles.
+     *    * If false, omit he cos(phi) factor on the derivative wrt theta.  This ensures it is always nonzero and can be safely used in cross product for surface normal.
+     * @param result optional plane result
+     */
+    radiansToPointAndDerivatives(thetaRadians: number, phiRadians: number, applyCosPhiFactor?: boolean, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors;
+    /**
+     * Evaluate a point and derivatives wrt to theta, phi, thetaTheta, phiPhi, and thetaPhi.
+     * All outputs are to caller-allocated points and vectors.
+     * @param thetaRadians longitude, in radians
+     * @param phiRadians latitude, in radians
+     * @param point (returned) surface point
+     * @param d1Theta (returned) derivative wrt theta
+     * @param d1Phi (returned) derivative wrt phi
+     * @param d2ThetaTheta (returned) second derivative wrt theta twice
+     * @param d2PhiPhi (returned) second derivative wrt phi twice
+     * @param d2ThetaPhi (returned) second derivative wrt theta and phi
+     * @param result optional plane result
+     */
+    radiansToPointAnd2Derivatives(thetaRadians: number, phiRadians: number, point: Point3d, d1Theta: Vector3d, d1Phi: Vector3d, d2ThetaTheta: Vector3d, d2PhiPhi: Vector3d, d2ThetaPhi: Vector3d): void;
+    /**
+     * Evaluate a point and rigid local coordinate frame the ellipsoid at angles give in radians.
+     * * The undefined return is only possible if the placement transform is singular (and even then only at critical angles)
+     * @param thetaRadians longitude, in radians
+     * @param phiRadians latitude, in radians
+     * @param result optional transform result
+     *
+     */
+    radiansToFrenetFrame(thetaRadians: number, phiRadians: number, result?: Transform): Transform | undefined;
+    /**
+     * Evaluate a point and unit normal at given angles.
+     * @param thetaRadians longitude, in radians
+     * @param phiRadians latitude, in radians
+     * @param result optional transform result
+     *
+     */
+    radiansToUnitNormalRay(thetaRadians: number, phiRadians: number, result?: Ray3d): Ray3d | undefined;
+    /**
+     * Find the (unique) extreme point for a given true surface perpendicular vector (outward)
+     */
+    surfaceNormalToAngles(normal: Vector3d, result?: LongitudeLatitudeNumber): LongitudeLatitudeNumber;
+    /**
+     * * Evaluate the surface normal on `other` ellipsoid at given angles
+     *    * If `other` is undefined, default to unit sphere.
+     * * Find the angles for the same normal on `this` ellipsoid
+     */
+    otherEllipsoidAnglesToThisEllipsoidAngles(otherEllipsoid: Ellipsoid | undefined, otherAngles: LongitudeLatitudeNumber, result?: LongitudeLatitudeNumber): LongitudeLatitudeNumber | undefined;
+    /**
+     * * if ellipsoid is given, return its surface point and unit normal as a Ray3d.
+     * * if not given, return surface point and unit normal for unit sphere.
+     */
+    static radiansToUnitNormalRay(ellipsoid: Ellipsoid | undefined, thetaRadians: number, phiRadians: number, result?: Ray3d): Ray3d | undefined;
+    /** Implement the `isPointInOnOrOutside` test fom the `interface` */
+    isPointOnOrInside(point: Point3d): boolean;
+    /** Announce "in" portions of a line segment.  See `Clipper.announceClippedSegmentIntervals` */
+    announceClippedSegmentIntervals(f0: number, f1: number, pointA: Point3d, pointB: Point3d, announce?: AnnounceNumberNumber): boolean;
+    /** Announce "in" portions of a line segment.  See `Clipper.announceClippedSegmentIntervals` */
+    announceClippedArcIntervals(arc: Arc3d, announce?: AnnounceNumberNumberCurvePrimitive): boolean;
+}
+/**
+ * * An `EllipsoidPatch` is
+ *   * An underlying (full) `Ellipsoid` object
+ *   * an angular range (`AngleSweep`) of longitudes around the equator
+ *   * an angular range (`AngleSweep`) of latitudes, with 0 at the equator, +90 degrees at north pole.
+ * * The `EllipsoidPatch` implements `UVSurface` methods, so a `PolyfaceBuilder` can generate facets in its method `addUVGridBody`
+ * @public
+ */
+export declare class EllipsoidPatch implements UVSurface {
+    ellipsoid: Ellipsoid;
+    longitudeSweep: AngleSweep;
+    latitudeSweep: AngleSweep;
+    /**
+     * CAPTURE ellipsoid and sweeps as an EllipsoidPatch.
+     * @param ellipsoid
+     * @param longitudeSweep
+     * @param latitudeSweep
+     */
+    private constructor();
+    /**
+     * Create a new EllipsoidPatch, capturing (not cloning) all input object references.
+     * @param ellipsoid  full ellipsoid
+     * @param longitudeSweep sweep of longitudes in the active patch
+     * @param latitudeSweep sweep of latitudes in the active patch.
+     */
+    static createCapture(ellipsoid: Ellipsoid, longitudeSweep: AngleSweep, latitudeSweep: AngleSweep): EllipsoidPatch;
+    /** Return the point on the ellipsoid at fractional positions in the angular ranges. */
+    uvFractionToPoint(longitudeFraction: number, latitudeFraction: number, result?: Point3d): Point3d;
+    /** Return the point and derivative vectors on the ellipsoid at fractional positions in the angular ranges.
+     * * Derivatives are with respect to fractional position.
+     */
+    uvFractionToPointAndTangents(longitudeFraction: number, latitudeFraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors;
+    /** Return the range of the patch, considering both boundary and internal extrema. */
+    range(result?: Range3d): Range3d;
+    /** Return intersections of the ray and surface.
+     * * uv values in the intersections are in radians unless `convertIntersectionRadiansToFractions` is true requesting conversion to patch fractions.
+     */
+    intersectRay(ray: Ray3d, restrictToPatch: boolean, convertIntersectionRadiansToFractions?: boolean): CurveAndSurfaceLocationDetail[];
+    /**
+     * test if the angles of the `LongitudeLatitudeNumber` are within the sweep ranges.
+     * @param position longitude and latitude to test.
+     * @param `allowPeriodicLongitude` true to allow the longitude to be in when shifted by a multiple of 2 PI
+     *    (latitude is never periodic for patches)
+     */
+    containsAngles(position: LongitudeLatitudeNumber, allowPeriodicLongitude?: boolean): boolean;
+    /**
+     * Compute point (with altitude) at given angles and altitude.
+     * * Never fails for non-singular ellipsoid.
+     * * In the returned ray,
+     *    * ray.origin is the point at requested altitude.
+     *    * ray.direction is an outward-directed unit vector
+     * @param position longitude, latitude, and height
+     *
+     */
+    anglesToUnitNormalRay(position: LongitudeLatitudeNumber, result?: Ray3d): Ray3d | undefined;
+    /**
+     * Return simple angles of a fractional position in the patch.
+     * @param thetaFraction fractional position in longitude (theta) interval
+     * @param phiFraction fractional position in latitude (phi) interval
+     * @param h optional altitude
+     * @param result optional preallocated result.
+     */
+    uvFractionToAngles(longitudeFraction: number, phiFraction: number, h?: number, result?: LongitudeLatitudeNumber): LongitudeLatitudeNumber;
+    /** Find the closest point of the (patch of the) ellipsoid. */
+    projectPointToSurface(spacePoint: Point3d): LongitudeLatitudeNumber | undefined;
+}
+/**
+ * Detailed data for a point on a 2-angle parameter space.
+ * @public
+ */
+export declare class GeodesicPathPoint {
+    /** First angle, in radians */
+    thetaRadians: number;
+    /** Second angle, in radians */
+    phiRadians: number;
+    point: Point3d;
+    dTheta: Vector3d;
+    dPhi: Vector3d;
+    d2Theta: Vector3d;
+    d2Phi: Vector3d;
+    d2ThetaPhi: Vector3d;
+    d1Cross: Vector3d;
+    constructor();
+    /** Fill all evaluations at given theta and phi. */
+    evaluateDerivativesAtCurrentAngles(ellipsoid: Ellipsoid): void;
+    private static _vectorAB?;
+    private static _vectorCB?;
+    private static _vectorCross?;
+    /** Evaluate the newton function and derivatives:
+     *          `(UAB cross UCB) dot d1cross`
+     * with as the central data, UAB = vector from pointA to pointB, UCB = vector from pointC to pointA.
+     * * Return order is:
+     *   * values[0] = the function
+     *   * values[1] = derivative wrt pointA.phi
+     *   * values[2] = derivative wrt pointB.phi
+     *   * values[3] = derivative wrt pointC.phi
+     */
+    static evaluateNewtonFunction(pointA: GeodesicPathPoint, pointB: GeodesicPathPoint, pointC: GeodesicPathPoint, values: Float64Array): void;
+    /**
+     * Extract the two angles form this structure to a LongitudeLatitudeNumber structure.
+     */
+    toAngles(): LongitudeLatitudeNumber;
+}
+/**
+ * Algorithm implementation class for computing approximate optimal (shortest) path points.
+ * * Call the static method `createGeodesicPath` to compute path points.
+ * @public
+ */
+export declare class GeodesicPathSolver {
+    private _defaultArc;
+    private _pathPoints;
+    private _tridiagonalSolver;
+    private constructor();
+    /**
+     *
+     * @param originalEllipsoid Given start and endpoints on an ellipsoid, compute points along a near-optimal shortest path.
+     * * The points are located so that at each point the local surface normal is contained in the plane of the point and its two neighbors.
+     * @param startAngles angles for the start of the path
+     * @param endAngles angles for the end of the path
+     * @param density If this is a number, it is the requested edge count.  If this is an angle, it ias an angular spacing measured in the great arc through the two points.
+     */
+    static createGeodesicPath(originalEllipsoid: Ellipsoid, startAngles: LongitudeLatitudeNumber, endAngles: LongitudeLatitudeNumber, density: number | Angle): GeodesicPathPoint[] | undefined;
+    private createInitialPointsAndTridiagonalSystem;
+    private applyUpdate;
+    /**
+     * Set up a step with specified ellipsoid.
+     * * ASSUME angles in _pathPoints are valid on given ellipsoid.
+     * @param ellipsoid
+     */
+    private setupStep;
+    private solve;
+    /**
+     * Construct various section arcs (on the ellipsoid), using planes that (a) pass through the two given points and (b) have in-plane vector sampled between the normals of the two points.
+     * * Each candidate ellipse has is in a plane with ellipsoid normal at vector constructed "between" the endpoint normals.
+     * * The intermediate construction is by interpolation between stated fractions (which maybe outside 0 to 1)
+     * @param ellipsoid
+     * @param angleA start point of all candidates
+     * @param angleB end point of all candidates
+     * @param numSample number of ellipses to construct as candidates.
+     * @param normalInterpolationFraction0
+     * @param normalInterpolationFraction1
+     */
+    static approximateMinimumLengthSectionArc(ellipsoid: Ellipsoid, angleA: LongitudeLatitudeNumber, angleB: LongitudeLatitudeNumber, numSample: number, normalInterpolationFraction0: number, normalInterpolationFraction1: number): {
+        minLengthArc: Arc3d;
+        minLengthNormalInterpolationFraction: number;
+    } | undefined;
+}
+//# sourceMappingURL=Ellipsoid.d.ts.map

package/lib/esm/geometry3d/Ellipsoid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Ellipsoid.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Ellipsoid.ts"],"names":[],"mappings":"AAKA;;GAEG;AAEH,OAAO,EAAE,6BAA6B,EAA2B,MAAM,kCAAkC,CAAC;AAC1G,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AAEvC,OAAO,EAAE,oBAAoB,EAAE,kCAAkC,EAAE,MAAM,yBAAyB,CAAC;AAEnG,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAKhD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,SAAS,EAAE,MAAM,mBAAmB,CAAC;AAC9C,OAAO,EAAE,uBAAuB,EAAE,MAAM,6BAA6B,CAAC;AACtE,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,4BAA4B,EAAE,MAAM,gCAAgC,CAAC;AAC9E,OAAO,EAAE,yBAAyB,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAW,OAAO,EAAE,MAAM,SAAS,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AA4GpC;;;;;;;;GAQG;AACH,qBAAa,SAAU,YAAW,OAAO;IACvC,OAAO,CAAC,UAAU,CAAY;IAC9B,OAAO,CAAC,YAAY,CAAW;IAC/B,OAAO,CAAC,YAAY,CAAW;IAC/B,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO,CAAC,WAAW,CAAU;IAC7B,OAAO;IAOP;;OAEG;WACW,MAAM,CAAC,iBAAiB,CAAC,EAAE,SAAS,GAAG,QAAQ,GAAG,SAAS;IAQzE;;;;;;;OAOG;WACW,uBAAuB,CAAC,MAAM,EAAE,OAAO,EAAE,IAAI,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,SAAS;IAQhJ;;;;;;OAMG;IACH,IAAW,YAAY,IAAI,SAAS,CAA4B;IAChE;;;;;;OAMG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAG9E;;;;;OAKG;IACI,YAAY,CAAC,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIlE,2CAA2C;IACpC,KAAK,IAAI,SAAS;IAGzB,oCAAoC;IAC7B,aAAa,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAG/C,wCAAwC;IACjC,mBAAmB,CAAC,SAAS,EAAE,SAAS,GAAG,OAAO;IAIzD;;;OAGG;IACI,gBAAgB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAKpE;;;OAGG;IACI,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;IAKtF;;OAEG;IACI,aAAa,CAAC,QAAQ,EAAE,OAAO,GAAG,KAAK,GAAG,SAAS;IAoB1D;;;;;;;;OAQG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,YAAY,EAAE,MAAM,EAAE,GAAG,SAAS,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,SAAS,EAAE,eAAe,EAAE,uBAAuB,EAAE,GAAG,SAAS,GAAG,MAAM;IA+B/J,4DAA4D;IACrD,yBAAyB,CAAC,aAAa,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAqCnJ;;;;;OAKG;IACI,cAAc,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ1F;;;;;;;;;;;;OAYG;IACI,qBAAqB,CAC1B,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAapC;;OAEG;IACI,mBAAmB,CAAC,MAAM,EAAE,uBAAuB,EAAE,MAAM,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAI/H;;;OAGG;IACI,kBAAkB,CAAC,KAAK,EAAE,4BAA4B,GAAG,KAAK,GAAG,SAAS;IAsBjF;;;;;OAKG;IACI,uCAAuC,CAAC,YAAY,EAAE,uBAAuB,EAAE,YAAY,EAAE,uBAAuB,EAAE,aAAa,EAAE,QAAQ,EAClJ,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IA6BpC;;;;;;;;;;OAUG;IACI,gCAAgC,CACrC,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,aAAa,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAC1C,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAc5C;;;;;OAKG;IACI,oBAAoB,CAAC,SAAS,EAAE,KAAK,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAU3G;;;;;OAKG;IACI,mBAAmB,CAAC,cAAc,EAAE,UAAU,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAY1G;;;;;;OAMG;IACI,gCAAgC,CACrC,MAAM,EAAE,uBAAuB,EAC/B,0BAA0B,EAAE,MAAM,EAClC,MAAM,EAAE,uBAAuB,GAAG,KAAK;IAQzC;;;;;;;;;;OAUG;IACI,4BAA4B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,iBAAiB,UAAO,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAkBtK;;;;;;;;;;;;OAYG;IACI,6BAA6B,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAC3E,KAAK,EAAE,OAAO,EACd,OAAO,EAAE,QAAQ,EACjB,KAAK,EAAE,QAAQ,EACf,YAAY,EAAE,QAAQ,EACtB,QAAQ,EAAE,QAAQ,EAClB,UAAU,EAAE,QAAQ;IAmBtB;;;;;;;OAOG;IACI,oBAAoB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,GAAG,SAAS;IAIhH;;;;;;OAMG;IACI,sBAAsB,CAAC,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAK1G;;OAEG;IACI,qBAAqB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAWzG;;;;OAIG;IACI,yCAAyC,CAAC,cAAc,EAAE,SAAS,GAAG,SAAS,EAAE,WAAW,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB,GAAG,SAAS;IAMpM;;;OAGG;WACW,sBAAsB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,EAAE,YAAY,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAWnJ,oEAAoE;IAC7D,iBAAiB,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAMjD,+FAA+F;IACxF,+BAA+B,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,CAAC,EAAE,oBAAoB,GAAG,OAAO;IAqC1I,+FAA+F;IACxF,2BAA2B,CAAC,GAAG,EAAE,KAAK,EAAE,QAAQ,CAAC,EAAE,kCAAkC,GAAG,OAAO;CA6CvG;AACD;;;;;;;GAOG;AACH,qBAAa,cAAe,YAAW,SAAS;IACvC,SAAS,EAAE,SAAS,CAAC;IACrB,cAAc,EAAE,UAAU,CAAC;IAC3B,aAAa,EAAE,UAAU,CAAC;IACjC;;;;;OAKG;IACH,OAAO;IAKP;;;;;OAKG;WACW,aAAa,CAAC,SAAS,EAAE,SAAS,EAAE,cAAc,EAAE,UAAU,EAAE,aAAa,EAAE,UAAU,GAAG,cAAc;IAGxH,uFAAuF;IAChF,iBAAiB,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGxG;;OAEG;IACI,4BAA4B,CAAC,iBAAiB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,yBAAyB,GAAG,yBAAyB;IAUvJ,qFAAqF;IAC9E,KAAK,CAAC,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAGvC;;OAEG;IACI,YAAY,CAAC,GAAG,EAAE,KAAK,EAAE,eAAe,EAAE,OAAO,EAAE,qCAAqC,GAAE,OAAe,GAAG,6BAA6B,EAAE;IA4BlJ;;;;;OAKG;IACI,cAAc,CAAC,QAAQ,EAAE,uBAAuB,EAAE,sBAAsB,GAAE,OAAc,GAAG,OAAO;IAKzG;;;;;;;;OAQG;IACI,qBAAqB,CAAC,QAAQ,EAAE,uBAAuB,EAAE,MAAM,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,SAAS;IAOlG;;;;;;OAMG;IACI,kBAAkB,CAAC,iBAAiB,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,CAAC,GAAE,MAAU,EAAE,MAAM,CAAC,EAAE,uBAAuB,GAAG,uBAAuB;IAGnJ,8DAA8D;IACvD,qBAAqB,CAAC,UAAU,EAAE,OAAO,GAAG,uBAAuB,GAAG,SAAS;CAGvF;AA4DD;;;GAGG;AACH,qBAAa,iBAAiB;IAC5B,8BAA8B;IACvB,YAAY,EAAE,MAAM,CAAC;IAC5B,+BAA+B;IACxB,UAAU,EAAE,MAAM,CAAC;IACnB,KAAK,EAAE,OAAO,CAAC;IACf,MAAM,EAAE,QAAQ,CAAC;IACjB,IAAI,EAAE,QAAQ,CAAC;IACf,OAAO,EAAE,QAAQ,CAAC;IAClB,KAAK,EAAE,QAAQ,CAAC;IAChB,UAAU,EAAE,QAAQ,CAAC;IACrB,OAAO,EAAE,QAAQ,CAAC;;IAYzB,mDAAmD;IAC5C,kCAAkC,CAAC,SAAS,EAAE,SAAS;IAI9D,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,SAAS,CAAC,CAAW;IACpC,OAAO,CAAC,MAAM,CAAC,YAAY,CAAC,CAAW;IACvC;;;;;;;;OAQG;WACW,sBAAsB,CAAC,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,EAAE,YAAY;IAkB1I;;OAEG;IACI,QAAQ,IAAI,uBAAuB;CAG3C;AACD;;;;GAIG;AACH,qBAAa,kBAAkB;IAC7B,OAAO,CAAC,WAAW,CAAQ;IAC3B,OAAO,CAAC,WAAW,CAAsB;IACzC,OAAO,CAAC,kBAAkB,CAAqB;IAC/C,OAAO;IAIP;;;;;;;OAOG;WACW,kBAAkB,CAAC,iBAAiB,EAAE,SAAS,EAC3D,WAAW,EAAE,uBAAuB,EAAE,SAAS,EAAE,uBAAuB,EAAE,OAAO,EAAE,MAAM,GAAG,KAAK,GAAG,iBAAiB,EAAE,GAAG,SAAS;IAsDrI,OAAO,CAAC,uCAAuC;IAc/C,OAAO,CAAC,WAAW;IAUnB;;;;OAIG;IACH,OAAO,CAAC,SAAS;IAqBjB,OAAO,CAAC,KAAK;IAGb;;;;;;;;;;OAUG;WACW,kCAAkC,CAAC,SAAS,EAAE,SAAS,EACnE,MAAM,EAAE,uBAAuB,EAC/B,MAAM,EAAE,uBAAuB,EAAE,SAAS,EAAE,MAAM,EAClD,4BAA4B,EAAE,MAAM,EACpC,4BAA4B,EAAE,MAAM,GAAG;QAAE,YAAY,EAAE,KAAK,CAAC;QAAC,oCAAoC,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;CA6B3H"}