@itwin/core-geometry 4.0.0-dev.54 → 4.0.0-dev.56

package/lib/cjs/curve/spiral/TransitionSpiral3d.js.map

@@ -1 +1 @@
-{"version":3,"file":"TransitionSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionSpiral3d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,6CAA0C;AAE1C,kDAA0D;AAC1D,0DAAuD;AACvD,0DAAuD;AACvD,sDAAmD;AACnD,mFAAgF;AAChF,6FAA0F;AAC1F,mFAAyE;AAwCzE;;;;;;;;;;;;GAYG;AACH,MAAsB,kBAAmB,SAAQ,+BAAc;IAU7D,wEAAwE;IACxE,IAAW,sBAAsB,KAAgB,OAAO,IAAI,CAAC,uBAAuB,CAAC,CAAC,CAAC;IAKvF,0CAA0C;IAC1C,IAAW,YAAY,KAAgB,OAAO,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IAEnE,YAAsB,UAA8B,EAAE,YAAuB,EAAE,sBAA6C,EAAE,gBAA6D;QACzL,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACjE,IAAI,CAAC,iBAAiB,GAAG,gBAAgB,CAAC;QAC1C,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,uBAAuB,GAAG,sBAAsB,CAAC,CAAC,CAAC,sBAAsB,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1G,CAAC;IAED,IAAW,UAAU,KAAa,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAC5D,iIAAiI;IAC1H,MAAM,CAAC,iBAAiB,CAAC,MAAc,IAAY,OAAO,CAAC,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC;IAEzG,+GAA+G;IACxG,MAAM,CAAC,iBAAiB,CAAC,SAAiB;QAC/C,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,GAAG,mBAAQ,CAAC,iBAAiB;YAClD,OAAO,GAAG,CAAC;QACb,OAAO,GAAG,GAAG,SAAS,CAAC;IACzB,CAAC;IAQD,sEAAsE;IACtD,iBAAiB,CAAC,SAAiB,EAAE,SAAiB;QACpE,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC7B,MAAM,eAAe,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;QAChF,MAAM,eAAe,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;QAChF,OAAO,CAAC,uBAAuB,CAAC,GAAG,CAAC,eAAe,EAAE,eAAe,CAAC,CAAC;QACtE,OAAO,CAAC,yBAAyB,EAAE,CAAC;QACpC,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,sCAAsC;IACtB,gBAAgB,CAAC,SAAoB;QACnD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC,CAAC,CAAC,2CAA2C;QAClF,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,0DAA0D;IACnD,MAAM,CAAC,gBAAgB,CAAC,YAAuB;QACpD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/H,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,oBAAoB,CAAC,EAAU,EAAE,EAAU;QACvD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IACrG,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,wBAAwB,CAAC,EAAU,EAAE,QAAgB,EAAE,EAAU;QAC7E,OAAO,mBAAQ,CAAC,WAAW,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5H,CAAC;IAED,qFAAqF;IAC9E,MAAM,CAAC,mCAAmC,CAAC,OAAe,EAAE,OAAe,EAAE,YAAoB;QACtG,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;IAC5F,CAAC;IAED,yEAAyE;IAClE,MAAM,CAAC,gCAAgC,CAAC,OAAe,EAAE,OAAe,EAAE,SAAiB;QAChG,OAAO,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,GAAG,SAAS,CAAC;IAC/E,CAAC;IAED,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,iFAAiF;IACjF,IAAW,gBAAgB,KAAkD,OAAO,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC;IAC7G;;;;;OAKG;IACO,yBAAyB,CAAC,UAAqB;QACvD,MAAM,SAAS,GAAG,UAAU,CAAC,MAAM,CAAC,0BAA0B,EAAE,CAAC;QACjE,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,6BAA6B;YAC7B,8DAA8D;YAC9D,MAAM,WAAW,GAAG,UAAU,CAAC,0BAA0B,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;YAC7E,gDAAgD;YAChD,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,CAAC,oBAAoB,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC;YACnF,IAAI,CAAC,aAAa,GAAG,qBAAS,CAAC,qBAAqB,CAAC,WAAW,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;YAClF,IAAI,IAAI,CAAC,gBAAgB;gBACvB,IAAI,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;YAE1D,OAAO,SAAS,CAAC;SAClB;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;;;OAKG;IACa,iBAAiB,CAAC,uBAA+C;QAC/E,MAAM,OAAO,GAAG,oCAAa,CAAC,MAAM,CAAC,uBAAuB,CAAC,CAAC;QAC9D,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,IAAI,EAAE,OAAO,CAAC,kBAAkB,CAAC,CAAC;QAC3E,IAAI,CAAC,kBAAkB,CAAC,OAAO,EAAE,OAAO,CAAC,aAAa,CAAC,CAAC;QACxD,OAAO,OAAO,CAAC,WAAW,EAAE,CAAC;IAC/B,CAAC;IACD,oDAAoD;IACpC,WAAW,CAAC,aAAsB,EAAE,SAAqB;QACvE,MAAM,OAAO,GAAG,IAAI,CAAC,qBAAqB,CAAE,GAAG,EAAE,GAAG,EAAE,SAAS,CAAC,CAAC;QACjE,aAAa,CAAC,WAAW,CAAE,OAAO,CAAC,CAAC;IACtC,CAAC;IAED;;OAEG;IACc,qBAAqB,CAAC,SAAiB,EAAE,SAAiB,EAAE,SAAqB;QAChG,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC;QACnC,IAAI,SAAS,KAAK,OAAO;YACvB,OAAO,eAAO,CAAC,UAAU,EAAG,CAAC;QAC/B,IAAI,KAAK,GAAG,IAAI,CAAC,IAAI,CAAE,OAAO,CAAC,SAAS,EAAE,GAAG,IAAI,CAAC,GAAG,CAAE,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC;QAC/E,KAAK,GAAG,mBAAQ,CAAC,KAAK,CAAE,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;QACtC,OAAO,IAAI,CAAC,4BAA4B,CAAE,SAAS,EAAE,SAAS,EAAE,KAAK,EAAE,SAAS,EAAE,GAAG,CAAC,CAAC;IACzF,CAAC;IACD;;;;OAIG;IACa,uBAAuB,CAAC,GAAqB,EAAE,OAAiB;QAC9E,OAAO,qDAAyB,CAAC,kCAAkC,CAAC,IAAI,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAC1F,CAAC;CACF;AAlKD,gDAkKC","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\n/** @packageDocumentation\n * @module Curve\n */\n\nimport { Geometry } from \"../../Geometry\";\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\nimport { Range1d, Range3d } from \"../../geometry3d/Range\";\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\nimport { Transform } from \"../../geometry3d/Transform\";\nimport { CurvePrimitive } from \"../CurvePrimitive\";\nimport { CurveOffsetXYHandler } from \"../internalContexts/CurveOffsetXYHandler\";\nimport { PlaneAltitudeRangeContext } from \"../internalContexts/PlaneAltitudeRangeContext\";\nimport { OffsetOptions } from \"../internalContexts/PolygonOffsetContext\";\nimport { LineString3d } from \"../LineString3d\";\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\n\nimport type { Vector3d } from \"../../geometry3d/Point3dVector3d\";\nimport type { Ray3d } from \"../../geometry3d/Ray3d\";\n\n/**\n * This is the set of valid type names for \"integrated\" spirals\n * * Behavior is expressed by a `NormalizedTransition` snap function.\n * * The snap function varies smoothly from f(0)=0 to f(1)=1\n * * The various snap functions are:\n *   * clothoid: linear\n *   * biquadratic: 2 quadratics pieced together, joining with 1st derivative continuity at f(0.) = 0.5, with zero slope f'(0)=0 and f'(1)= 0\n *   * bloss: A single cubic with zero slope at 0 and 1\n *   * cosine: half of a cosine wave, centered around 0.5\n *   * sine: full period of a sine wave added to the line f(u)=u\n * *\n * @public\n */\nexport type IntegratedSpiralTypeName = \"clothoid\" | \"bloss\" | \"biquadratic\" | \"cosine\" | \"sine\";\n\n/**\n * This is the set of valid type names for \"direct\" spirals.\n * \"Direct\" spirals can evaluate fractionToPoint by direct equations, i.e. not requiring the numeric integrations in \"Integrated\" spiral types.\n * @public\n */\nexport type DirectSpiralTypeName =\n   \"JapaneseCubic\"  // 1 term from each of the X,Y clothoid series expansions:  y = x^3 / (6RL)\n  | \"Arema\"       // 2 terms from each of the X,Y clothoid series expansions.  Identical to ChineseCubic!\n  | \"ChineseCubic\"  // Identical to Arema!\n  | \"HalfCosine\"  // high continuity cosine variation from quadratic.\n  | \"AustralianRailCorp\" // cubic with high accuracy distance series\n  | \"WesternAustralian\"  // simple cubic -- 2 terms of x series, 1 term of y series.\n  | \"Czech\"  // simple cubic with two term distance approximation\n  | \"MXCubicAlongArc\"  // x obtained from fractional distance via 2-terms from series, y = x^3/ (6RL)\n  | \"Polish\"\n  | \"Italian\"\n  ;\n\n/**\n * TransitionSpiral3d is a base class for multiple variants of spirals.\n * * The menagerie of spiral types have 2 broad categories:\n *   * IntegratedSpiral3d -- a spiral whose direct function for curvature versus distance must be integrated to determine x,y\n *     * The IntegratedSpiral3d types are enumerated in `IntegratedSpiralTypes`\n *   * DirectSpiral3d -- a spiral implemented with direct calculation of x,y from fractional position along the spiral.\n *     * The direct spiral types are enumerated in the `DirectSpiralType`\n * * The method set for CurvePrimitive support includes a `handleTransitionSpiral(g: TransitionSpiral3d)` which receives all the spiral types.\n * * The spiral class may impose expectations that its inflection is at the origin, with tangent along the x axis.\n *   * This is generally necessary for direct spirals.\n *   * This is not necessary for integrated spirals.\n * @public\n */\nexport abstract class TransitionSpiral3d extends CurvePrimitive {\n  /** string name of spiral type */\n  protected _spiralType: string;\n  /** Original defining properties. */\n  protected _designProperties: TransitionConditionalProperties | undefined;\n\n  /** Fractional interval for the \"active\" part of a containing spiral.\n   * (The radius, angle, and length conditions define a complete spiral, and some portion of it is \"active\")\n   */\n  protected _activeFractionInterval: Segment1d;\n  /** Return (reference to) the active portion of the reference spiral. */\n  public get activeFractionInterval(): Segment1d { return this._activeFractionInterval; }\n  /** strokes in the active portion */\n  public abstract get activeStrokes(): LineString3d;\n  /** Placement transform */\n  protected _localToWorld: Transform;\n  /** (reference to) placement transform. */\n  public get localToWorld(): Transform { return this._localToWorld; }\n\n  protected constructor(spiralType: string | undefined, localToWorld: Transform, activeFractionInterval: Segment1d | undefined, designProperties: TransitionConditionalProperties | undefined) {\n    super();\n    this._spiralType = spiralType ? spiralType : \"unknownSpiralType\";\n    this._designProperties = designProperties;\n    this._localToWorld = localToWorld;\n    this._activeFractionInterval = activeFractionInterval ? activeFractionInterval : Segment1d.create(0, 1);\n  }\n\n  public get spiralType(): string { return this._spiralType; }\n  /** Return 1/r with convention that if true zero is given as radius it represents infinite radius (0 curvature, straight line) */\n  public static radiusToCurvature(radius: number): number { return (radius === 0.0) ? 0.0 : 1.0 / radius; }\n\n  /** Return 1/k with convention that if near-zero is given as curvature, its infinite radius is returned as 0 */\n  public static curvatureToRadius(curvature: number): number {\n    if (Math.abs(curvature) < Geometry.smallAngleRadians)\n      return 0.0;\n    return 1.0 / curvature;\n  }\n\n  /** Return a deep clone. */\n  public abstract override clone(): TransitionSpiral3d;\n\n  /** Recompute strokes */\n  public abstract refreshComputedProperties(): void;\n\n  /** Return (if possible) a spiral which is a portion of this curve. */\n  public override clonePartialCurve(fractionA: number, fractionB: number): TransitionSpiral3d {\n    const spiralB = this.clone();\n    const globalFractionA = this._activeFractionInterval.fractionToPoint(fractionA);\n    const globalFractionB = this._activeFractionInterval.fractionToPoint(fractionB);\n    spiralB._activeFractionInterval.set(globalFractionA, globalFractionB);\n    spiralB.refreshComputedProperties();\n    return spiralB;\n  }\n\n  /** Clone with a transform applied  */\n  public override cloneTransformed(transform: Transform): TransitionSpiral3d {\n    const result = this.clone();\n    result.tryTransformInPlace(transform); // ok, we're confident it will always work.\n    return result;\n  }\n\n  /** Return the average of the start and end curvatures. */\n  public static averageCurvature(radiusLimits: Segment1d): number {\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(radiusLimits.x0) + TransitionSpiral3d.radiusToCurvature(radiusLimits.x1));\n  }\n  /**\n   * Given two radii (or zeros for 0 curvature) return the average curvature\n   * @param r0 start radius, or 0 for line\n   * @param r1 end radius, or 0 for line\n   */\n  public static averageCurvatureR0R1(r0: number, r1: number): number {\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(r0) + TransitionSpiral3d.radiusToCurvature(r1));\n  }\n  /**\n   * Given two radii (or zeros for 0 curvature) return the average curvature\n   * @param r0 start radius, or 0 for line\n   * @param r1 end radius, or 0 for line\n   */\n  public static interpolateCurvatureR0R1(r0: number, fraction: number, r1: number): number {\n    return Geometry.interpolate(TransitionSpiral3d.radiusToCurvature(r0), fraction, TransitionSpiral3d.radiusToCurvature(r1));\n  }\n\n  /** Return the arc length of a transition spiral with given sweep and radius pair. */\n  public static radiusRadiusSweepRadiansToArcLength(radius0: number, radius1: number, sweepRadians: number): number {\n    return Math.abs(sweepRadians / TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1));\n  }\n\n  /** Return the turn angle for spiral of given length between two radii */\n  public static radiusRadiusLengthToSweepRadians(radius0: number, radius1: number, arcLength: number): number {\n    return TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1) * arcLength;\n  }\n\n  /** Return the end radius for spiral of given start radius, length, and turn angle. */\n  public static radius0LengthSweepRadiansToRadius1(radius0: number, arcLength: number, sweepRadians: number) {\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius0));\n  }\n  /** Return the start radius for spiral of given end radius, length, and turn angle. */\n  public static radius1LengthSweepRadiansToRadius0(radius1: number, arcLength: number, sweepRadians: number) {\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius1));\n  }\n  /** Return the original defining properties (if any) saved by the constructor. */\n  public get designProperties(): TransitionConditionalProperties | undefined { return this._designProperties; }\n  /**\n   * * If transformA is rigid with uniform scale, apply the rigid part of transformA to the localToWorld transform and return the scale and rigid separation.\n   * * If not rigid, do nothing and return undefined.\n   * * Also apply the scale factor to the designProperties.\n   * @param transformA\n   */\n  protected applyRigidPartOfTransform(transformA: Transform): { rigidAxes: Matrix3d, scale: number } | undefined {\n    const rigidData = transformA.matrix.factorRigidWithSignedScale();\n    if (rigidData !== undefined) {\n      // [sQ a][R b] = [sQ*R sQb+a]\n      // but we save it as [Q*R sQb+a] with spiral data scaled by s.\n      const transformC0 = transformA.multiplyTransformTransform(this.localToWorld);\n      // BUT pull the scale part out of the matrix ...\n      const matrixC = rigidData.rigidAxes.multiplyMatrixMatrix(this.localToWorld.matrix);\n      this._localToWorld = Transform.createOriginAndMatrix(transformC0.origin, matrixC);\n      if (this.designProperties)\n        this.designProperties.applyScaleFactor(rigidData.scale);\n\n      return rigidData;\n    }\n    return undefined;\n  }\n\n  /**\n   * Construct an offset of the instance curve as viewed in the xy-plane (ignoring z).\n   * * No attempt is made to join the offsets of smaller constituent primitives. To construct a fully joined offset\n   *   for an aggregate instance (e.g., LineString3d, CurveChainWithDistanceIndex), use RegionOps.constructCurveXYOffset() instead.\n   * @param offsetDistanceOrOptions offset distance (positive to left of the instance curve), or options object\n   */\n  public override constructOffsetXY(offsetDistanceOrOptions: number | OffsetOptions): CurvePrimitive | CurvePrimitive[] | undefined {\n    const options = OffsetOptions.create(offsetDistanceOrOptions);\n    const handler = new CurveOffsetXYHandler(this, options.leftOffsetDistance);\n    this.emitStrokableParts(handler, options.strokeOptions);\n    return handler.claimResult();\n  }\n  /** extend the range by the strokes of the spiral */\n  public override extendRange(rangeToExtend: Range3d, transform?: Transform): void {\n    const myRange = this.rangeBetweenFractions (0.0, 1.0, transform);\n    rangeToExtend.extendRange (myRange);\n  }\n\n  /** return the range of spiral between fractions of the activeStrokes.\n   * * Use activeStrokes point count times interval factor for initial evaluation count, but do at least 5\n   */\n   public override rangeBetweenFractions(fractionA: number, fractionB: number, transform?: Transform): Range3d {\n    const strokes = this.activeStrokes;\n    if (undefined === strokes)\n      return Range3d.createNull ();\n    let count = Math.ceil (strokes.numPoints() * Math.abs (fractionB - fractionA));\n    count = Geometry.clamp (5, count, 30);\n    return this.rangeBetweenFractionsByCount (fractionA, fractionB, count, transform, 0.5);\n  }\n  /** Project instance geometry (via dispatch) onto the given ray, and return the extreme fractional parameters of projection.\n   * @param ray ray onto which the instance is projected. A `Vector3d` is treated as a `Ray3d` with zero origin.\n   * @param lowHigh optional receiver for output\n   * @returns range of fractional projection parameters onto the ray, where 0.0 is start of the ray and 1.0 is the end of the ray.\n   */\n  public override projectedParameterRange(ray: Vector3d | Ray3d, lowHigh?: Range1d): Range1d | undefined {\n    return PlaneAltitudeRangeContext.findExtremeFractionsAlongDirection(this, ray, lowHigh);\n  }\n}\n"]}
+{"version":3,"file":"TransitionSpiral3d.js","sourceRoot":"","sources":["../../../../src/curve/spiral/TransitionSpiral3d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AAEH,6CAA0C;AAE1C,kDAA0D;AAC1D,0DAAuD;AACvD,0DAAuD;AACvD,sDAAmD;AACnD,mFAAgF;AAChF,6FAA0F;AAC1F,mFAAyE;AAwCzE;;;;;;;;;;;;GAYG;AACH,MAAsB,kBAAmB,SAAQ,+BAAc;IAU7D,wEAAwE;IACxE,IAAW,sBAAsB,KAAgB,OAAO,IAAI,CAAC,uBAAuB,CAAC,CAAC,CAAC;IAKvF,0CAA0C;IAC1C,IAAW,YAAY,KAAgB,OAAO,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC;IAEnE,YAAsB,UAA8B,EAAE,YAAuB,EAAE,sBAA6C,EAAE,gBAA6D;QACzL,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACjE,IAAI,CAAC,iBAAiB,GAAG,gBAAgB,CAAC;QAC1C,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,uBAAuB,GAAG,sBAAsB,CAAC,CAAC,CAAC,sBAAsB,CAAC,CAAC,CAAC,qBAAS,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1G,CAAC;IAED,IAAW,UAAU,KAAa,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC;IAC5D,iIAAiI;IAC1H,MAAM,CAAC,iBAAiB,CAAC,MAAc,IAAY,OAAO,CAAC,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC;IAEzG,+GAA+G;IACxG,MAAM,CAAC,iBAAiB,CAAC,SAAiB;QAC/C,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,GAAG,mBAAQ,CAAC,iBAAiB;YAClD,OAAO,GAAG,CAAC;QACb,OAAO,GAAG,GAAG,SAAS,CAAC;IACzB,CAAC;IAQD,sEAAsE;IACtD,iBAAiB,CAAC,SAAiB,EAAE,SAAiB;QACpE,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC7B,MAAM,eAAe,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;QAChF,MAAM,eAAe,GAAG,IAAI,CAAC,uBAAuB,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;QAChF,OAAO,CAAC,uBAAuB,CAAC,GAAG,CAAC,eAAe,EAAE,eAAe,CAAC,CAAC;QACtE,OAAO,CAAC,yBAAyB,EAAE,CAAC;QACpC,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,sCAAsC;IACtB,gBAAgB,CAAC,SAAoB;QACnD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QAC5B,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC,CAAC,CAAC,2CAA2C;QAClF,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,0DAA0D;IACnD,MAAM,CAAC,gBAAgB,CAAC,YAAuB;QACpD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/H,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,oBAAoB,CAAC,EAAU,EAAE,EAAU;QACvD,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IACrG,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,wBAAwB,CAAC,EAAU,EAAE,QAAgB,EAAE,EAAU;QAC7E,OAAO,mBAAQ,CAAC,WAAW,CAAC,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,kBAAkB,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5H,CAAC;IAED,qFAAqF;IAC9E,MAAM,CAAC,mCAAmC,CAAC,OAAe,EAAE,OAAe,EAAE,YAAoB;QACtG,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,GAAG,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;IAC5F,CAAC;IAED,yEAAyE;IAClE,MAAM,CAAC,gCAAgC,CAAC,OAAe,EAAE,OAAe,EAAE,SAAiB;QAChG,OAAO,kBAAkB,CAAC,oBAAoB,CAAC,OAAO,EAAE,OAAO,CAAC,GAAG,SAAS,CAAC;IAC/E,CAAC;IAED,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,sFAAsF;IAC/E,MAAM,CAAC,kCAAkC,CAAC,OAAe,EAAE,SAAiB,EAAE,YAAoB;QACvG,OAAO,kBAAkB,CAAC,iBAAiB,CAAC,CAAC,GAAG,GAAG,YAAY,GAAG,SAAS,CAAC,GAAG,kBAAkB,CAAC,iBAAiB,CAAC,OAAO,CAAC,CAAC,CAAC;IAChI,CAAC;IACD,iFAAiF;IACjF,IAAW,gBAAgB,KAAkD,OAAO,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC;IAC7G;;;;;OAKG;IACO,yBAAyB,CAAC,UAAqB;QACvD,MAAM,SAAS,GAAG,UAAU,CAAC,MAAM,CAAC,0BAA0B,EAAE,CAAC;QACjE,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,6BAA6B;YAC7B,8DAA8D;YAC9D,MAAM,WAAW,GAAG,UAAU,CAAC,0BAA0B,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;YAC7E,gDAAgD;YAChD,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,CAAC,oBAAoB,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC;YACnF,IAAI,CAAC,aAAa,GAAG,qBAAS,CAAC,qBAAqB,CAAC,WAAW,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;YAClF,IAAI,IAAI,CAAC,gBAAgB;gBACvB,IAAI,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;YAE1D,OAAO,SAAS,CAAC;SAClB;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;;;OAKG;IACa,iBAAiB,CAAC,uBAA+C;QAC/E,MAAM,OAAO,GAAG,oCAAa,CAAC,MAAM,CAAC,uBAAuB,CAAC,CAAC;QAC9D,MAAM,OAAO,GAAG,IAAI,2CAAoB,CAAC,IAAI,EAAE,OAAO,CAAC,kBAAkB,CAAC,CAAC;QAC3E,IAAI,CAAC,kBAAkB,CAAC,OAAO,EAAE,OAAO,CAAC,aAAa,CAAC,CAAC;QACxD,OAAO,OAAO,CAAC,WAAW,EAAE,CAAC;IAC/B,CAAC;IACD,oDAAoD;IACpC,WAAW,CAAC,aAAsB,EAAE,SAAqB;QACvE,MAAM,OAAO,GAAG,IAAI,CAAC,qBAAqB,CAAE,GAAG,EAAE,GAAG,EAAE,SAAS,CAAC,CAAC;QACjE,aAAa,CAAC,WAAW,CAAE,OAAO,CAAC,CAAC;IACtC,CAAC;IAED;;OAEG;IACc,qBAAqB,CAAC,SAAiB,EAAE,SAAiB,EAAE,SAAqB;QAChG,MAAM,OAAO,GAAG,IAAI,CAAC,aAAa,CAAC;QACnC,IAAI,SAAS,KAAK,OAAO;YACvB,OAAO,eAAO,CAAC,UAAU,EAAG,CAAC;QAC/B,IAAI,KAAK,GAAG,IAAI,CAAC,IAAI,CAAE,OAAO,CAAC,SAAS,EAAE,GAAG,IAAI,CAAC,GAAG,CAAE,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC;QAC/E,KAAK,GAAG,mBAAQ,CAAC,KAAK,CAAE,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;QACtC,OAAO,IAAI,CAAC,4BAA4B,CAAE,SAAS,EAAE,SAAS,EAAE,KAAK,EAAE,SAAS,EAAE,GAAG,CAAC,CAAC;IACzF,CAAC;IACD;;;;OAIG;IACa,uBAAuB,CAAC,GAAqB,EAAE,OAAiB;QAC9E,OAAO,qDAAyB,CAAC,kCAAkC,CAAC,IAAI,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;IAC1F,CAAC;CACF;AAlKD,gDAkKC","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 Curve\r\n */\r\n\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Matrix3d } from \"../../geometry3d/Matrix3d\";\r\nimport { Range1d, Range3d } from \"../../geometry3d/Range\";\r\nimport { Segment1d } from \"../../geometry3d/Segment1d\";\r\nimport { Transform } from \"../../geometry3d/Transform\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { CurveOffsetXYHandler } from \"../internalContexts/CurveOffsetXYHandler\";\r\nimport { PlaneAltitudeRangeContext } from \"../internalContexts/PlaneAltitudeRangeContext\";\r\nimport { OffsetOptions } from \"../internalContexts/PolygonOffsetContext\";\r\nimport { LineString3d } from \"../LineString3d\";\r\nimport { TransitionConditionalProperties } from \"./TransitionConditionalProperties\";\r\n\r\nimport type { Vector3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport type { Ray3d } from \"../../geometry3d/Ray3d\";\r\n\r\n/**\r\n * This is the set of valid type names for \"integrated\" spirals\r\n * * Behavior is expressed by a `NormalizedTransition` snap function.\r\n * * The snap function varies smoothly from f(0)=0 to f(1)=1\r\n * * The various snap functions are:\r\n *   * clothoid: linear\r\n *   * biquadratic: 2 quadratics pieced together, joining with 1st derivative continuity at f(0.) = 0.5, with zero slope f'(0)=0 and f'(1)= 0\r\n *   * bloss: A single cubic with zero slope at 0 and 1\r\n *   * cosine: half of a cosine wave, centered around 0.5\r\n *   * sine: full period of a sine wave added to the line f(u)=u\r\n * *\r\n * @public\r\n */\r\nexport type IntegratedSpiralTypeName = \"clothoid\" | \"bloss\" | \"biquadratic\" | \"cosine\" | \"sine\";\r\n\r\n/**\r\n * This is the set of valid type names for \"direct\" spirals.\r\n * \"Direct\" spirals can evaluate fractionToPoint by direct equations, i.e. not requiring the numeric integrations in \"Integrated\" spiral types.\r\n * @public\r\n */\r\nexport type DirectSpiralTypeName =\r\n   \"JapaneseCubic\"  // 1 term from each of the X,Y clothoid series expansions:  y = x^3 / (6RL)\r\n  | \"Arema\"       // 2 terms from each of the X,Y clothoid series expansions.  Identical to ChineseCubic!\r\n  | \"ChineseCubic\"  // Identical to Arema!\r\n  | \"HalfCosine\"  // high continuity cosine variation from quadratic.\r\n  | \"AustralianRailCorp\" // cubic with high accuracy distance series\r\n  | \"WesternAustralian\"  // simple cubic -- 2 terms of x series, 1 term of y series.\r\n  | \"Czech\"  // simple cubic with two term distance approximation\r\n  | \"MXCubicAlongArc\"  // x obtained from fractional distance via 2-terms from series, y = x^3/ (6RL)\r\n  | \"Polish\"\r\n  | \"Italian\"\r\n  ;\r\n\r\n/**\r\n * TransitionSpiral3d is a base class for multiple variants of spirals.\r\n * * The menagerie of spiral types have 2 broad categories:\r\n *   * IntegratedSpiral3d -- a spiral whose direct function for curvature versus distance must be integrated to determine x,y\r\n *     * The IntegratedSpiral3d types are enumerated in `IntegratedSpiralTypes`\r\n *   * DirectSpiral3d -- a spiral implemented with direct calculation of x,y from fractional position along the spiral.\r\n *     * The direct spiral types are enumerated in the `DirectSpiralType`\r\n * * The method set for CurvePrimitive support includes a `handleTransitionSpiral(g: TransitionSpiral3d)` which receives all the spiral types.\r\n * * The spiral class may impose expectations that its inflection is at the origin, with tangent along the x axis.\r\n *   * This is generally necessary for direct spirals.\r\n *   * This is not necessary for integrated spirals.\r\n * @public\r\n */\r\nexport abstract class TransitionSpiral3d extends CurvePrimitive {\r\n  /** string name of spiral type */\r\n  protected _spiralType: string;\r\n  /** Original defining properties. */\r\n  protected _designProperties: TransitionConditionalProperties | undefined;\r\n\r\n  /** Fractional interval for the \"active\" part of a containing spiral.\r\n   * (The radius, angle, and length conditions define a complete spiral, and some portion of it is \"active\")\r\n   */\r\n  protected _activeFractionInterval: Segment1d;\r\n  /** Return (reference to) the active portion of the reference spiral. */\r\n  public get activeFractionInterval(): Segment1d { return this._activeFractionInterval; }\r\n  /** strokes in the active portion */\r\n  public abstract get activeStrokes(): LineString3d;\r\n  /** Placement transform */\r\n  protected _localToWorld: Transform;\r\n  /** (reference to) placement transform. */\r\n  public get localToWorld(): Transform { return this._localToWorld; }\r\n\r\n  protected constructor(spiralType: string | undefined, localToWorld: Transform, activeFractionInterval: Segment1d | undefined, designProperties: TransitionConditionalProperties | undefined) {\r\n    super();\r\n    this._spiralType = spiralType ? spiralType : \"unknownSpiralType\";\r\n    this._designProperties = designProperties;\r\n    this._localToWorld = localToWorld;\r\n    this._activeFractionInterval = activeFractionInterval ? activeFractionInterval : Segment1d.create(0, 1);\r\n  }\r\n\r\n  public get spiralType(): string { return this._spiralType; }\r\n  /** Return 1/r with convention that if true zero is given as radius it represents infinite radius (0 curvature, straight line) */\r\n  public static radiusToCurvature(radius: number): number { return (radius === 0.0) ? 0.0 : 1.0 / radius; }\r\n\r\n  /** Return 1/k with convention that if near-zero is given as curvature, its infinite radius is returned as 0 */\r\n  public static curvatureToRadius(curvature: number): number {\r\n    if (Math.abs(curvature) < Geometry.smallAngleRadians)\r\n      return 0.0;\r\n    return 1.0 / curvature;\r\n  }\r\n\r\n  /** Return a deep clone. */\r\n  public abstract override clone(): TransitionSpiral3d;\r\n\r\n  /** Recompute strokes */\r\n  public abstract refreshComputedProperties(): void;\r\n\r\n  /** Return (if possible) a spiral which is a portion of this curve. */\r\n  public override clonePartialCurve(fractionA: number, fractionB: number): TransitionSpiral3d {\r\n    const spiralB = this.clone();\r\n    const globalFractionA = this._activeFractionInterval.fractionToPoint(fractionA);\r\n    const globalFractionB = this._activeFractionInterval.fractionToPoint(fractionB);\r\n    spiralB._activeFractionInterval.set(globalFractionA, globalFractionB);\r\n    spiralB.refreshComputedProperties();\r\n    return spiralB;\r\n  }\r\n\r\n  /** Clone with a transform applied  */\r\n  public override cloneTransformed(transform: Transform): TransitionSpiral3d {\r\n    const result = this.clone();\r\n    result.tryTransformInPlace(transform); // ok, we're confident it will always work.\r\n    return result;\r\n  }\r\n\r\n  /** Return the average of the start and end curvatures. */\r\n  public static averageCurvature(radiusLimits: Segment1d): number {\r\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(radiusLimits.x0) + TransitionSpiral3d.radiusToCurvature(radiusLimits.x1));\r\n  }\r\n  /**\r\n   * Given two radii (or zeros for 0 curvature) return the average curvature\r\n   * @param r0 start radius, or 0 for line\r\n   * @param r1 end radius, or 0 for line\r\n   */\r\n  public static averageCurvatureR0R1(r0: number, r1: number): number {\r\n    return 0.5 * (TransitionSpiral3d.radiusToCurvature(r0) + TransitionSpiral3d.radiusToCurvature(r1));\r\n  }\r\n  /**\r\n   * Given two radii (or zeros for 0 curvature) return the average curvature\r\n   * @param r0 start radius, or 0 for line\r\n   * @param r1 end radius, or 0 for line\r\n   */\r\n  public static interpolateCurvatureR0R1(r0: number, fraction: number, r1: number): number {\r\n    return Geometry.interpolate(TransitionSpiral3d.radiusToCurvature(r0), fraction, TransitionSpiral3d.radiusToCurvature(r1));\r\n  }\r\n\r\n  /** Return the arc length of a transition spiral with given sweep and radius pair. */\r\n  public static radiusRadiusSweepRadiansToArcLength(radius0: number, radius1: number, sweepRadians: number): number {\r\n    return Math.abs(sweepRadians / TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1));\r\n  }\r\n\r\n  /** Return the turn angle for spiral of given length between two radii */\r\n  public static radiusRadiusLengthToSweepRadians(radius0: number, radius1: number, arcLength: number): number {\r\n    return TransitionSpiral3d.averageCurvatureR0R1(radius0, radius1) * arcLength;\r\n  }\r\n\r\n  /** Return the end radius for spiral of given start radius, length, and turn angle. */\r\n  public static radius0LengthSweepRadiansToRadius1(radius0: number, arcLength: number, sweepRadians: number) {\r\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius0));\r\n  }\r\n  /** Return the start radius for spiral of given end radius, length, and turn angle. */\r\n  public static radius1LengthSweepRadiansToRadius0(radius1: number, arcLength: number, sweepRadians: number) {\r\n    return TransitionSpiral3d.curvatureToRadius((2.0 * sweepRadians / arcLength) - TransitionSpiral3d.radiusToCurvature(radius1));\r\n  }\r\n  /** Return the original defining properties (if any) saved by the constructor. */\r\n  public get designProperties(): TransitionConditionalProperties | undefined { return this._designProperties; }\r\n  /**\r\n   * * If transformA is rigid with uniform scale, apply the rigid part of transformA to the localToWorld transform and return the scale and rigid separation.\r\n   * * If not rigid, do nothing and return undefined.\r\n   * * Also apply the scale factor to the designProperties.\r\n   * @param transformA\r\n   */\r\n  protected applyRigidPartOfTransform(transformA: Transform): { rigidAxes: Matrix3d, scale: number } | undefined {\r\n    const rigidData = transformA.matrix.factorRigidWithSignedScale();\r\n    if (rigidData !== undefined) {\r\n      // [sQ a][R b] = [sQ*R sQb+a]\r\n      // but we save it as [Q*R sQb+a] with spiral data scaled by s.\r\n      const transformC0 = transformA.multiplyTransformTransform(this.localToWorld);\r\n      // BUT pull the scale part out of the matrix ...\r\n      const matrixC = rigidData.rigidAxes.multiplyMatrixMatrix(this.localToWorld.matrix);\r\n      this._localToWorld = Transform.createOriginAndMatrix(transformC0.origin, matrixC);\r\n      if (this.designProperties)\r\n        this.designProperties.applyScaleFactor(rigidData.scale);\r\n\r\n      return rigidData;\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /**\r\n   * Construct an offset of the instance curve as viewed in the xy-plane (ignoring z).\r\n   * * No attempt is made to join the offsets of smaller constituent primitives. To construct a fully joined offset\r\n   *   for an aggregate instance (e.g., LineString3d, CurveChainWithDistanceIndex), use RegionOps.constructCurveXYOffset() instead.\r\n   * @param offsetDistanceOrOptions offset distance (positive to left of the instance curve), or options object\r\n   */\r\n  public override constructOffsetXY(offsetDistanceOrOptions: number | OffsetOptions): CurvePrimitive | CurvePrimitive[] | undefined {\r\n    const options = OffsetOptions.create(offsetDistanceOrOptions);\r\n    const handler = new CurveOffsetXYHandler(this, options.leftOffsetDistance);\r\n    this.emitStrokableParts(handler, options.strokeOptions);\r\n    return handler.claimResult();\r\n  }\r\n  /** extend the range by the strokes of the spiral */\r\n  public override extendRange(rangeToExtend: Range3d, transform?: Transform): void {\r\n    const myRange = this.rangeBetweenFractions (0.0, 1.0, transform);\r\n    rangeToExtend.extendRange (myRange);\r\n  }\r\n\r\n  /** return the range of spiral between fractions of the activeStrokes.\r\n   * * Use activeStrokes point count times interval factor for initial evaluation count, but do at least 5\r\n   */\r\n   public override rangeBetweenFractions(fractionA: number, fractionB: number, transform?: Transform): Range3d {\r\n    const strokes = this.activeStrokes;\r\n    if (undefined === strokes)\r\n      return Range3d.createNull ();\r\n    let count = Math.ceil (strokes.numPoints() * Math.abs (fractionB - fractionA));\r\n    count = Geometry.clamp (5, count, 30);\r\n    return this.rangeBetweenFractionsByCount (fractionA, fractionB, count, transform, 0.5);\r\n  }\r\n  /** Project instance geometry (via dispatch) onto the given ray, and return the extreme fractional parameters of projection.\r\n   * @param ray ray onto which the instance is projected. A `Vector3d` is treated as a `Ray3d` with zero origin.\r\n   * @param lowHigh optional receiver for output\r\n   * @returns range of fractional projection parameters onto the ray, where 0.0 is start of the ray and 1.0 is the end of the ray.\r\n   */\r\n  public override projectedParameterRange(ray: Vector3d | Ray3d, lowHigh?: Range1d): Range1d | undefined {\r\n    return PlaneAltitudeRangeContext.findExtremeFractionsAlongDirection(this, ray, lowHigh);\r\n  }\r\n}\r\n"]}

package/lib/cjs/curve/spiral/XYCurveEvaluator.js.map

@@ -1 +1 @@
-{"version":3,"file":"XYCurveEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/XYCurveEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,6CAA0C;AAC1C,0FAAuF;AACvF,sEAAqE;AACrE,kDAA+C;AAC/C,0DAAuD;AACvD,kDAAqD;AACrD;;;;;;;GAOG;AACH,MAAsB,gBAAgB;IAsBpC,iFAAiF;IAC1E,eAAe,CAAC,QAAgB,EAAE,MAAgB;QACvD,OAAO,yBAAO,CAAC,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;IAC7F,CAAC;IACD,4IAA4I;IACrI,4BAA4B,CAAC,QAAgB,EAAE,MAAc;QAClE,OAAO,aAAK,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EACnF,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC3D,MAAM,CAAC,CAAC;IACZ,CAAC;IACD,+IAA+I;IACxI,8BAA8B,CAAC,QAAgB,EAAE,MAAkC;QACxF,OAAO,qDAAyB,CAAC,yBAAyB,CACxD,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EAC3D,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC3D,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC7D,MAAM,CAAC,CAAC;IACZ,CAAC;IACD;;;OAGG;IACI,0BAA0B,CAAC,QAAgB;QAChD,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACtC,OAAO,mBAAQ,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,CAAC;IAWM,MAAM,CAAC,aAAa;QACzB,gBAAgB,CAAC,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/C,gBAAgB,CAAC,YAAY,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACpD,gBAAgB,CAAC,YAAY,GAAG,uBAAU,CAAC,WAAW,CAAC;IACzD,CAAC;IACD;;;;;OAKG;IACI,iCAAiC,CAAC,SAAiB,EAAE,SAAiB;QAC3E,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,CAAC;QACxC,MAAM,WAAW,GAAG,gBAAgB,CAAC,YAAY,CAAC;QAClD,MAAM,OAAO,GAAG,gBAAgB,CAAC,YAAY,CAAC,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,WAAW,CAAC,CAAC;QACzF,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE;YAChC,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,0BAA0B,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;SACpE;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IACD;;;;;;;OAOG;IACI,uBAAuB,CAAC,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EAAE,cAAsB;QAC/H,MAAM,aAAa,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;QAC9G,IAAI,aAAa,KAAK,SAAS,EAAE;YAC/B,OAAO,qBAAY,CAAC,WAAW,CAAC,aAAa,EAC3C,CAAC,QAAgB,EAAE,EAAE;gBACnB,MAAM,CAAC,GAAG,IAAI,CAAC,iCAAiC,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;gBACtE,OAAO,SAAS,GAAG,CAAC,GAAG,cAAc,CAAC;YACxC,CAAC,EACD,CAAC,QAAgB,EAAE,EAAE,CAAC,IAAI,CAAC,0BAA0B,CAAC,QAAQ,CAAC,CAAC,CAAC;SACpE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;;;;;OAOG;IACI,8BAA8B,CAAC,QAAgB,EAAE,EAAW,EAAE,IAAe,EAAE,IAAe,EAAE,IAAe;QACpH,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QAClE,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QACxE,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1E,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5E,CAAC;CAGF;AAxHD,4CAwHC;AACD,sDAAsD;AACtD,gBAAgB,CAAC,aAAa,EAAE,CAAC","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\n/** @packageDocumentation\n * @module Curve\n */\nimport { Geometry } from \"../../Geometry\";\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\nimport { Point3d, Vector3d } from \"../../geometry3d/Point3dVector3d\";\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\nimport { Quadrature } from \"../../numerics/Quadrature\";\nimport { SimpleNewton } from \"../../numerics/Newton\";\n/**\n * XYCurveEvaluator is an abstract with methods for evaluating X and Y parts of a curve parameterized by a fraction.\n * * The required methods call for independent X and Y evaluation.\n * * Base class methods package those (multiple) calls into point, ray, and plane structures.\n *    * A implementation that has evaluation substantial cost that can be shared among x,y parts or between\n *       primary functions and derivatives might choose to implement the point and derivative methods directly.\n * @internal\n */\nexport abstract class XYCurveEvaluator {\n  /** return a deep copy of the evaluator */\n  public abstract clone(): XYCurveEvaluator;\n  /** test for near identical evaluator */\n  public abstract isAlmostEqual(other: any): boolean;\n  /** Evaluate X at fractional position. */\n  public abstract fractionToX(fraction: number): number;\n  /** Evaluate Y at fractional position. */\n  public abstract fractionToY(fraction: number): number;\n  /** Evaluate derivative of X with respect to fraction at fractional position. */\n  public abstract fractionToDX(fraction: number): number;\n  /** Evaluate derivative of Y with respect to fraction at fractional position. */\n  public abstract fractionToDY(fraction: number): number;\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\n  public abstract fractionToDDX(fraction: number): number;\n  /** Evaluate second derivative of Y with respect to fraction at fractional position. */\n  public abstract fractionToDDY(fraction: number): number;\n  /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\n  public abstract fractionToD3X(fraction: number): number;\n  /** Evaluate second derivative of Y with respect to fraction at fractional position. */\n  public abstract fractionToD3Y(fraction: number): number;\n  /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */\n  public fractionToPoint(fraction: number, result?: Point3d): Point3d {\n    return Point3d.create(this.fractionToX(fraction), this.fractionToY(fraction), 0.0, result);\n  }\n  /** Evaluate both X and Y and their first derivatives at fractional coordinate, return bundled as origin and (non-unit) direction vector. */\n  public fractionToPointAndDerivative(fraction: number, result?: Ray3d): Ray3d {\n    return Ray3d.createXYZUVW(this.fractionToX(fraction), this.fractionToY(fraction), 0.0,\n      this.fractionToDX(fraction), this.fractionToDY(fraction), 0,\n      result);\n  }\n  /** Evaluate both X and Y and their second derivatives at fractional coordinate, return bundled as origin and (non-unit) vectorU an vectorV. */\n  public fractionToPointAnd2Derivatives(fraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors {\n    return Plane3dByOriginAndVectors.createOriginAndVectorsXYZ(\n      this.fractionToX(fraction), this.fractionToY(fraction), 0.0,\n      this.fractionToDX(fraction), this.fractionToDY(fraction), 0,\n      this.fractionToDDX(fraction), this.fractionToDDY(fraction), 0,\n      result);\n  }\n  /**\n   * Return the magnitude of the tangent vector at fraction.\n   * @param fraction fractional position along the curve\n   */\n  public fractionToTangentMagnitude(fraction: number): number {\n    const u = this.fractionToDX(fraction);\n    const v = this.fractionToDY(fraction);\n    return Geometry.hypotenuseXY(u, v);\n  }\n\n  /** Invert the fractionToX function for given X. */\n\n  public abstract xToFraction(x: number): number | undefined;\n  /** Initialize class level work arrays for 5 point Gaussian Quadrature. */\n  // Class resources for integration . . .\n  // These static variables are reused on calls to integrateFromStartFraction\n  protected static _gaussX: Float64Array;\n  protected static _gaussWeight: Float64Array;\n  protected static _gaussMapper: (xA: number, xB: number, arrayX: Float64Array, arrayW: Float64Array) => number;\n  public static initWorkSpace() {\n    XYCurveEvaluator._gaussX = new Float64Array(5);\n    XYCurveEvaluator._gaussWeight = new Float64Array(5);\n    XYCurveEvaluator._gaussMapper = Quadrature.setupGauss5;\n  }\n  /**\n   * Integrate between nominal fractions with default gauss rule.\n   * * The caller is expected to choose nearby fractions so that the single gauss interval accuracy is good.\n   * @param fraction0\n   * @param fraction1\n   */\n  public integrateDistanceBetweenFractions(fraction0: number, fraction1: number): number {\n    const gaussX = XYCurveEvaluator._gaussX;\n    const gaussWeight = XYCurveEvaluator._gaussWeight;\n    const numEval = XYCurveEvaluator._gaussMapper(fraction0, fraction1, gaussX, gaussWeight);\n    let sum = 0;\n    for (let k = 0; k < numEval; k++) {\n      sum += gaussWeight[k] * this.fractionToTangentMagnitude(gaussX[k]);\n    }\n    return sum;\n  }\n  /**\n   * Inverse integrated distance\n   * @param fraction0 start of fraction interval\n   * @param fraction1 end of fraction interval\n   * @param distance0 distance at start\n   * @param distance1 distance at end\n   * @param targetDistance intermediate distance.\n   */\n  public inverseDistanceFraction(fraction0: number, fraction1: number, distance0: number, distance1: number, targetDistance: number): number | undefined {\n    const startFraction = Geometry.inverseInterpolate(fraction0, distance0, fraction1, distance1, targetDistance);\n    if (startFraction !== undefined) {\n      return SimpleNewton.runNewton1D(startFraction,\n        (fraction: number) => {\n          const d = this.integrateDistanceBetweenFractions(fraction0, fraction);\n          return distance0 + d - targetDistance;\n        },\n        (fraction: number) => this.fractionToTangentMagnitude(fraction));\n    }\n    return undefined;\n  }\n\n  /**\n   *\n   * @param fraction fractional position along x axis\n   * @param xy xy coordinates of point on the curve\n   * @param d1xy\n   * @param d2xy\n   * @param d3xy\n   */\n  public fractionToPointAnd3Derivatives(fraction: number, xy: Point3d, d1xy?: Vector3d, d2xy?: Vector3d, d3xy?: Vector3d) {\n    xy.set(this.fractionToX(fraction), this.fractionToY(fraction), 0);\n    if (d1xy)\n      d1xy.set(this.fractionToDX(fraction), this.fractionToDY(fraction), 0);\n    if (d2xy)\n      d2xy.set(this.fractionToDDX(fraction), this.fractionToDDY(fraction), 0);\n    if (d3xy)\n      d3xy.set(this.fractionToD3X(fraction), this.fractionToD3Y(fraction), 0);\n  }\n  /** Apply a uniform scale around the origin. */\n  public abstract scaleInPlace(scaleFactor: number): void;\n}\n// at load time, initialize gauss quadrature workspace\nXYCurveEvaluator.initWorkSpace();\n"]}
+{"version":3,"file":"XYCurveEvaluator.js","sourceRoot":"","sources":["../../../../src/curve/spiral/XYCurveEvaluator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,6CAA0C;AAC1C,0FAAuF;AACvF,sEAAqE;AACrE,kDAA+C;AAC/C,0DAAuD;AACvD,kDAAqD;AACrD;;;;;;;GAOG;AACH,MAAsB,gBAAgB;IAsBpC,iFAAiF;IAC1E,eAAe,CAAC,QAAgB,EAAE,MAAgB;QACvD,OAAO,yBAAO,CAAC,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;IAC7F,CAAC;IACD,4IAA4I;IACrI,4BAA4B,CAAC,QAAgB,EAAE,MAAc;QAClE,OAAO,aAAK,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EACnF,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC3D,MAAM,CAAC,CAAC;IACZ,CAAC;IACD,+IAA+I;IACxI,8BAA8B,CAAC,QAAgB,EAAE,MAAkC;QACxF,OAAO,qDAAyB,CAAC,yBAAyB,CACxD,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,GAAG,EAC3D,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC3D,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC7D,MAAM,CAAC,CAAC;IACZ,CAAC;IACD;;;OAGG;IACI,0BAA0B,CAAC,QAAgB;QAChD,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACtC,OAAO,mBAAQ,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,CAAC;IAWM,MAAM,CAAC,aAAa;QACzB,gBAAgB,CAAC,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/C,gBAAgB,CAAC,YAAY,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACpD,gBAAgB,CAAC,YAAY,GAAG,uBAAU,CAAC,WAAW,CAAC;IACzD,CAAC;IACD;;;;;OAKG;IACI,iCAAiC,CAAC,SAAiB,EAAE,SAAiB;QAC3E,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,CAAC;QACxC,MAAM,WAAW,GAAG,gBAAgB,CAAC,YAAY,CAAC;QAClD,MAAM,OAAO,GAAG,gBAAgB,CAAC,YAAY,CAAC,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,WAAW,CAAC,CAAC;QACzF,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE;YAChC,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,0BAA0B,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;SACpE;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IACD;;;;;;;OAOG;IACI,uBAAuB,CAAC,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EAAE,cAAsB;QAC/H,MAAM,aAAa,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,SAAS,EAAE,cAAc,CAAC,CAAC;QAC9G,IAAI,aAAa,KAAK,SAAS,EAAE;YAC/B,OAAO,qBAAY,CAAC,WAAW,CAAC,aAAa,EAC3C,CAAC,QAAgB,EAAE,EAAE;gBACnB,MAAM,CAAC,GAAG,IAAI,CAAC,iCAAiC,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;gBACtE,OAAO,SAAS,GAAG,CAAC,GAAG,cAAc,CAAC;YACxC,CAAC,EACD,CAAC,QAAgB,EAAE,EAAE,CAAC,IAAI,CAAC,0BAA0B,CAAC,QAAQ,CAAC,CAAC,CAAC;SACpE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;;;;;OAOG;IACI,8BAA8B,CAAC,QAAgB,EAAE,EAAW,EAAE,IAAe,EAAE,IAAe,EAAE,IAAe;QACpH,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QAClE,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QACxE,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1E,IAAI,IAAI;YACN,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5E,CAAC;CAGF;AAxHD,4CAwHC;AACD,sDAAsD;AACtD,gBAAgB,CAAC,aAAa,EAAE,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\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Plane3dByOriginAndVectors } from \"../../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { Point3d, Vector3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { Ray3d } from \"../../geometry3d/Ray3d\";\r\nimport { Quadrature } from \"../../numerics/Quadrature\";\r\nimport { SimpleNewton } from \"../../numerics/Newton\";\r\n/**\r\n * XYCurveEvaluator is an abstract with methods for evaluating X and Y parts of a curve parameterized by a fraction.\r\n * * The required methods call for independent X and Y evaluation.\r\n * * Base class methods package those (multiple) calls into point, ray, and plane structures.\r\n *    * A implementation that has evaluation substantial cost that can be shared among x,y parts or between\r\n *       primary functions and derivatives might choose to implement the point and derivative methods directly.\r\n * @internal\r\n */\r\nexport abstract class XYCurveEvaluator {\r\n  /** return a deep copy of the evaluator */\r\n  public abstract clone(): XYCurveEvaluator;\r\n  /** test for near identical evaluator */\r\n  public abstract isAlmostEqual(other: any): boolean;\r\n  /** Evaluate X at fractional position. */\r\n  public abstract fractionToX(fraction: number): number;\r\n  /** Evaluate Y at fractional position. */\r\n  public abstract fractionToY(fraction: number): number;\r\n  /** Evaluate derivative of X with respect to fraction at fractional position. */\r\n  public abstract fractionToDX(fraction: number): number;\r\n  /** Evaluate derivative of Y with respect to fraction at fractional position. */\r\n  public abstract fractionToDY(fraction: number): number;\r\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\r\n  public abstract fractionToDDX(fraction: number): number;\r\n  /** Evaluate second derivative of Y with respect to fraction at fractional position. */\r\n  public abstract fractionToDDY(fraction: number): number;\r\n  /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */\r\n  /** Evaluate second derivative of X with respect to fraction at fractional position. */\r\n  public abstract fractionToD3X(fraction: number): number;\r\n  /** Evaluate second derivative of Y with respect to fraction at fractional position. */\r\n  public abstract fractionToD3Y(fraction: number): number;\r\n  /** Evaluate both X and Y at fractional coordinate, return bundled as a point. */\r\n  public fractionToPoint(fraction: number, result?: Point3d): Point3d {\r\n    return Point3d.create(this.fractionToX(fraction), this.fractionToY(fraction), 0.0, result);\r\n  }\r\n  /** Evaluate both X and Y and their first derivatives at fractional coordinate, return bundled as origin and (non-unit) direction vector. */\r\n  public fractionToPointAndDerivative(fraction: number, result?: Ray3d): Ray3d {\r\n    return Ray3d.createXYZUVW(this.fractionToX(fraction), this.fractionToY(fraction), 0.0,\r\n      this.fractionToDX(fraction), this.fractionToDY(fraction), 0,\r\n      result);\r\n  }\r\n  /** Evaluate both X and Y and their second derivatives at fractional coordinate, return bundled as origin and (non-unit) vectorU an vectorV. */\r\n  public fractionToPointAnd2Derivatives(fraction: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors {\r\n    return Plane3dByOriginAndVectors.createOriginAndVectorsXYZ(\r\n      this.fractionToX(fraction), this.fractionToY(fraction), 0.0,\r\n      this.fractionToDX(fraction), this.fractionToDY(fraction), 0,\r\n      this.fractionToDDX(fraction), this.fractionToDDY(fraction), 0,\r\n      result);\r\n  }\r\n  /**\r\n   * Return the magnitude of the tangent vector at fraction.\r\n   * @param fraction fractional position along the curve\r\n   */\r\n  public fractionToTangentMagnitude(fraction: number): number {\r\n    const u = this.fractionToDX(fraction);\r\n    const v = this.fractionToDY(fraction);\r\n    return Geometry.hypotenuseXY(u, v);\r\n  }\r\n\r\n  /** Invert the fractionToX function for given X. */\r\n\r\n  public abstract xToFraction(x: number): number | undefined;\r\n  /** Initialize class level work arrays for 5 point Gaussian Quadrature. */\r\n  // Class resources for integration . . .\r\n  // These static variables are reused on calls to integrateFromStartFraction\r\n  protected static _gaussX: Float64Array;\r\n  protected static _gaussWeight: Float64Array;\r\n  protected static _gaussMapper: (xA: number, xB: number, arrayX: Float64Array, arrayW: Float64Array) => number;\r\n  public static initWorkSpace() {\r\n    XYCurveEvaluator._gaussX = new Float64Array(5);\r\n    XYCurveEvaluator._gaussWeight = new Float64Array(5);\r\n    XYCurveEvaluator._gaussMapper = Quadrature.setupGauss5;\r\n  }\r\n  /**\r\n   * Integrate between nominal fractions with default gauss rule.\r\n   * * The caller is expected to choose nearby fractions so that the single gauss interval accuracy is good.\r\n   * @param fraction0\r\n   * @param fraction1\r\n   */\r\n  public integrateDistanceBetweenFractions(fraction0: number, fraction1: number): number {\r\n    const gaussX = XYCurveEvaluator._gaussX;\r\n    const gaussWeight = XYCurveEvaluator._gaussWeight;\r\n    const numEval = XYCurveEvaluator._gaussMapper(fraction0, fraction1, gaussX, gaussWeight);\r\n    let sum = 0;\r\n    for (let k = 0; k < numEval; k++) {\r\n      sum += gaussWeight[k] * this.fractionToTangentMagnitude(gaussX[k]);\r\n    }\r\n    return sum;\r\n  }\r\n  /**\r\n   * Inverse integrated distance\r\n   * @param fraction0 start of fraction interval\r\n   * @param fraction1 end of fraction interval\r\n   * @param distance0 distance at start\r\n   * @param distance1 distance at end\r\n   * @param targetDistance intermediate distance.\r\n   */\r\n  public inverseDistanceFraction(fraction0: number, fraction1: number, distance0: number, distance1: number, targetDistance: number): number | undefined {\r\n    const startFraction = Geometry.inverseInterpolate(fraction0, distance0, fraction1, distance1, targetDistance);\r\n    if (startFraction !== undefined) {\r\n      return SimpleNewton.runNewton1D(startFraction,\r\n        (fraction: number) => {\r\n          const d = this.integrateDistanceBetweenFractions(fraction0, fraction);\r\n          return distance0 + d - targetDistance;\r\n        },\r\n        (fraction: number) => this.fractionToTangentMagnitude(fraction));\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /**\r\n   *\r\n   * @param fraction fractional position along x axis\r\n   * @param xy xy coordinates of point on the curve\r\n   * @param d1xy\r\n   * @param d2xy\r\n   * @param d3xy\r\n   */\r\n  public fractionToPointAnd3Derivatives(fraction: number, xy: Point3d, d1xy?: Vector3d, d2xy?: Vector3d, d3xy?: Vector3d) {\r\n    xy.set(this.fractionToX(fraction), this.fractionToY(fraction), 0);\r\n    if (d1xy)\r\n      d1xy.set(this.fractionToDX(fraction), this.fractionToDY(fraction), 0);\r\n    if (d2xy)\r\n      d2xy.set(this.fractionToDDX(fraction), this.fractionToDDY(fraction), 0);\r\n    if (d3xy)\r\n      d3xy.set(this.fractionToD3X(fraction), this.fractionToD3Y(fraction), 0);\r\n  }\r\n  /** Apply a uniform scale around the origin. */\r\n  public abstract scaleInPlace(scaleFactor: number): void;\r\n}\r\n// at load time, initialize gauss quadrature workspace\r\nXYCurveEvaluator.initWorkSpace();\r\n"]}

package/lib/cjs/geometry3d/Angle.d.ts

@@ -90,11 +90,11 @@ export declare class Angle implements BeJSONFunctions {
      */
     setFromJSON(json?: AngleProps, defaultValRadians?: number): void;
     /**
-   * Create an Angle from a JSON object
-   * @param json object from JSON.parse. If a number, value is in *DEGREES*
-   * @param defaultValRadians if json is undefined, default value in radians.
-   * @return a new Angle
-   */
+     * Create an Angle from a JSON object
+     * @param json object from JSON.parse. If a number, value is in *DEGREES*
+     * @param defaultValRadians if json is undefined, default value in radians.
+     * @return a new Angle
+     */
     static fromJSON(json?: AngleProps, defaultValRadians?: number): Angle;
     /** Convert an Angle to a JSON object as a number in degrees */
     toJSON(): AngleProps;
@@ -114,25 +114,19 @@ export declare class Angle implements BeJSONFunctions {
      * @param degrees angle in radians
      */
     static radiansToDegrees(radians: number): number;
-    /**
-     * Return the cosine of this Angle object's angle.
-     */
+    /** Return the cosine of this Angle object's angle */
     cos(): number;
-    /**
-     * Return the sine of this Angle object's angle.
-     */
+    /** Return the sine of this Angle object's angle */
     sin(): number;
-    /**
-     * Return the tangent of this Angle object's angle.
-     */
+    /** Return the tangent of this Angle object's angle */
     tan(): number;
     /** Test if a radians (absolute) value is nearly 2PI or larger! */
     static isFullCircleRadians(radians: number): boolean;
     /** Test if the radians value is a half circle */
     static isHalfCircleRadians(radians: number): boolean;
-    /** test if the angle is a full circle */
+    /** Test if the angle is a full circle */
     get isFullCircle(): boolean;
-    /** test if the angle is a half circle (in either direction) */
+    /** Test if the angle is a half circle (in either direction) */
     get isHalfCircle(): boolean;
     /** Adjust a radians value so it is positive in 0..360 */
     static adjustDegrees0To360(degrees: number): number;
@@ -142,7 +136,7 @@ export declare class Angle implements BeJSONFunctions {
     static adjustRadians0To2Pi(radians: number): number;
     /** Adjust a radians value so it is positive in -PI..PI */
     static adjustRadiansMinusPiPlusPi(radians: number): number;
-    /** return a (newly allocated) Angle object with value 0 radians */
+    /** Return a (newly allocated) Angle object with value 0 radians */
     static zero(): Angle;
     /** Test if the angle is exactly zero. */
     get isExactZero(): boolean;
@@ -205,24 +199,23 @@ export declare class Angle implements BeJSONFunctions {
     /** If value is close to -1, -0.5, 0, 0.5, 1, adjust it to the exact value. */
     static cleanupTrigValue(value: number, tolerance?: number): number;
     /**
-       * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define
-       * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The
-       * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.
-       * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.
-       * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is
-       * perpendicular to the ellipse.
-       * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)
-       * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.
-       * math details can be found at docs/learning/geometry/Angle.md
-       * @param dotUU dot product of vectorU with itself
-       * @param dotVV dot product of vectorV with itself
-       * @param dotUV dot product of vectorU with vectorV
-       */
+     * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define
+     * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The
+     * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.
+     * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.
+     * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is
+     * perpendicular to the ellipse.
+     * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)
+     * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.
+     * math details can be found at docs/learning/geometry/Angle.md
+     * @param dotUU dot product of vectorU with itself
+     * @param dotVV dot product of vectorV with itself
+     * @param dotUV dot product of vectorU with vectorV
+     */
     static dotProductsToHalfAngleTrigValues(dotUU: number, dotVV: number, dotUV: number, favorZero?: boolean): TrigValues;
     /**
-     * * Returns the angle between two vectors, with the vectors given as xyz components
+     * Returns the angle between two vectors, with the vectors given as xyz components
      * * The returned angle is between 0 and PI
-     *
      * @param ux x component of vector u
      * @param uy y component of vector u
      * @param uz z component of vector u
@@ -232,7 +225,8 @@ export declare class Angle implements BeJSONFunctions {
      */
     static radiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number): number;
     /**
-     * * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve angle to a full 2PI range.
+     * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve
+     * angle to a full 2PI range.
      * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)
      * * The angle is in the plane of the U and V vectors.
      * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.
@@ -246,7 +240,8 @@ export declare class Angle implements BeJSONFunctions {
      * @param upVectorX x component of vector to positive side of plane.
      * @param upVectorY y component of vector to positive side of plane.
      * @param upVectorZ z component of vector to positive side of plane.
-     * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI).  If false, return signed (-PI < radians <= PI)
+     * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI). If false, return
+     * signed (-PI < radians <= PI)
      */
     static orientedRadiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number, upVectorX: number, upVectorY: number, upVectorZ: number, adjustToPositive?: boolean): number;
     /**

package/lib/cjs/geometry3d/Angle.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Angle.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":"AAIA;;GAEG;AAEH,OAAO,EAAE,UAAU,EAAE,eAAe,EAAY,UAAU,EAAE,MAAM,aAAa,CAAC;AAEhF;;;;;;;;;GASG;AACH,qBAAa,KAAM,YAAW,eAAe;IAC3C,+DAA+D;IAC/D,gBAAuB,eAAe,uBAAuB;IAC7D,8DAA8D;IAC9D,gBAAuB,cAAc,sBAA0B;IAC/D,8DAA8D;IAC9D,gBAAuB,cAAc,sBAA2B;IAChE,6DAA6D;IAC7D,gBAAuB,SAAS,qBAA0B;IAC1D,+DAA+D;IAC/D,gBAAuB,UAAU,qBAA0B;IAC3D,qDAAqD;IACrD,gBAAuB,gBAAgB,SAAiC;IACxE,qDAAqD;IACrD,gBAAuB,gBAAgB,SAAiC;IACxE,OAAO,CAAC,QAAQ,CAAS;IACzB,OAAO,CAAC,QAAQ,CAAC,CAAS;IAC1B,OAAO;IAIP,gDAAgD;IACzC,KAAK,IAAI,KAAK;IAGrB,8CAA8C;IACvC,MAAM,IAAI,QAAQ,CAAC,IAAI,CAAC;IAG/B;;;OAGG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGnD;;;OAGG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGnD;;;;;OAKG;WACW,iBAAiB,CAAC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,GAAG,KAAK;IAGtF;;;OAGG;IACI,WAAW,CAAC,KAAK,EAAE,MAAM,GAAG,KAAK;IAGxC;;;OAGG;IACI,UAAU,CAAC,OAAO,EAAE,MAAM;IAIjC;;;OAGG;IACI,UAAU,CAAC,OAAO,EAAE,MAAM;IAIjC,yCAAyC;WAC3B,SAAS,IAAI,KAAK;IAGhC;;;;;OAKG;WACW,WAAW,CAAC,SAAS,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,GAAG,KAAK;IAGxE;;;OAGG;IACI,OAAO,CAAC,KAAK,EAAE,KAAK;IAI3B;;;;;;;OAOG;IACI,WAAW,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,iBAAiB,CAAC,EAAE,MAAM;IAgBhE;;;;;KAKC;WACa,QAAQ,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,iBAAiB,CAAC,EAAE,MAAM,GAAG,KAAK;IAK5E,+DAA+D;IACxD,MAAM,IAAI,UAAU;IAG3B,yEAAyE;IAClE,aAAa,IAAI,UAAU;IAGlC,4CAA4C;IAC5C,IAAW,OAAO,IAAI,MAAM,CAE3B;IACD,4CAA4C;IAC5C,IAAW,OAAO,IAAI,MAAM,CAE3B;IACD;;;OAGG;WACW,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAGvD;;;OAGG;WACW,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAsBvD;;OAEG;IACI,GAAG,IAAI,MAAM;IAGpB;;OAEG;IACI,GAAG,IAAI,MAAM;IAGpB;;OAEG;IACI,GAAG,IAAI,MAAM;IAGpB,kEAAkE;WACpD,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO;IAG3D,iDAAiD;WACnC,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO;IAG3D,yCAAyC;IACzC,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,+DAA+D;IAC/D,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,yDAAyD;WAC3C,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAe1D,mDAAmD;WACrC,sBAAsB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAc7D,yDAAyD;WAC3C,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAc1D,0DAA0D;WAC5C,0BAA0B,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAcjE,mEAAmE;WACrD,IAAI,IAAI,KAAK;IAG3B,yCAAyC;IACzC,IAAW,WAAW,IAAI,OAAO,CAEhC;IACD,uFAAuF;IACvF,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,wGAAwG;IACxG,IAAW,wBAAwB,IAAI,OAAO,CAE7C;IACD,gEAAgE;WAClD,2BAA2B,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGjE,mEAAmE;WACrD,4BAA4B,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGlE;;;;;OAKG;WACW,oCAAoC,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EACnF,SAAS,GAAE,MAAmC,GAAG,OAAO;IAW1D;;;;OAIG;IACI,6BAA6B,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG3G;;;;;OAKG;WACW,iCAAiC,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAChF,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG1D;;;;OAIG;IACI,0BAA0B,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAGxG;;;;;OAKG;IACI,aAAa,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG3F;;;;;OAKG;WACW,qBAAqB,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,OAAO;IAKzF;;;;;;OAMG;WACW,+BAA+B,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,UAAU;IA4CzF,8EAA8E;WAChE,gBAAgB,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,GAAE,MAAgB,GAAG,MAAM;IAYlF;;;;;;;;;;;;;SAaK;WACS,gCAAgC,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,SAAS,GAAE,OAAc,GAAG,UAAU;IAQlI;;;;;;;;;;OAUG;WACW,wBAAwB,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAItH;;;;;;;;;;;;;;;;OAgBG;WACW,gCAAgC,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EACnH,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EACvD,gBAAgB,GAAE,OAAe,GAAG,MAAM;IAoB5C;;;OAGG;IACI,uBAAuB,CAAC,QAAQ,EAAE,MAAM;CAQhD"}
+{"version":3,"file":"Angle.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":"AAIA;;GAEG;AAEH,OAAO,EAAE,UAAU,EAAE,eAAe,EAAY,UAAU,EAAE,MAAM,aAAa,CAAC;AAEhF;;;;;;;;;GASG;AACH,qBAAa,KAAM,YAAW,eAAe;IAC3C,+DAA+D;IAC/D,gBAAuB,eAAe,uBAAuB;IAC7D,8DAA8D;IAC9D,gBAAuB,cAAc,sBAA0B;IAC/D,8DAA8D;IAC9D,gBAAuB,cAAc,sBAA2B;IAChE,6DAA6D;IAC7D,gBAAuB,SAAS,qBAA0B;IAC1D,+DAA+D;IAC/D,gBAAuB,UAAU,qBAA0B;IAC3D,qDAAqD;IACrD,gBAAuB,gBAAgB,SAAiC;IACxE,qDAAqD;IACrD,gBAAuB,gBAAgB,SAAiC;IACxE,OAAO,CAAC,QAAQ,CAAS;IACzB,OAAO,CAAC,QAAQ,CAAC,CAAS;IAC1B,OAAO;IAIP,gDAAgD;IACzC,KAAK,IAAI,KAAK;IAGrB,8CAA8C;IACvC,MAAM,IAAI,QAAQ,CAAC,IAAI,CAAC;IAG/B;;;OAGG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGnD;;;OAGG;WACW,aAAa,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGnD;;;;;OAKG;WACW,iBAAiB,CAAC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,GAAG,KAAK;IAGtF;;;OAGG;IACI,WAAW,CAAC,KAAK,EAAE,MAAM,GAAG,KAAK;IAGxC;;;OAGG;IACI,UAAU,CAAC,OAAO,EAAE,MAAM;IAIjC;;;OAGG;IACI,UAAU,CAAC,OAAO,EAAE,MAAM;IAIjC,yCAAyC;WAC3B,SAAS,IAAI,KAAK;IAGhC;;;;;OAKG;WACW,WAAW,CAAC,SAAS,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,GAAG,KAAK;IAGxE;;;OAGG;IACI,OAAO,CAAC,KAAK,EAAE,KAAK;IAI3B;;;;;;;OAOG;IACI,WAAW,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,iBAAiB,CAAC,EAAE,MAAM;IAgBhE;;;;;OAKG;WACW,QAAQ,CAAC,IAAI,CAAC,EAAE,UAAU,EAAE,iBAAiB,CAAC,EAAE,MAAM,GAAG,KAAK;IAK5E,+DAA+D;IACxD,MAAM,IAAI,UAAU;IAG3B,yEAAyE;IAClE,aAAa,IAAI,UAAU;IAGlC,4CAA4C;IAC5C,IAAW,OAAO,IAAI,MAAM,CAE3B;IACD,4CAA4C;IAC5C,IAAW,OAAO,IAAI,MAAM,CAE3B;IACD;;;OAGG;WACW,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAGvD;;;OAGG;WACW,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAsBvD,qDAAqD;IAC9C,GAAG,IAAI,MAAM;IAGpB,mDAAmD;IAC5C,GAAG,IAAI,MAAM;IAGpB,sDAAsD;IAC/C,GAAG,IAAI,MAAM;IAGpB,kEAAkE;WACpD,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO;IAG3D,iDAAiD;WACnC,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO;IAG3D,yCAAyC;IACzC,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,+DAA+D;IAC/D,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,yDAAyD;WAC3C,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAe1D,mDAAmD;WACrC,sBAAsB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAc7D,yDAAyD;WAC3C,mBAAmB,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAc1D,0DAA0D;WAC5C,0BAA0B,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM;IAcjE,mEAAmE;WACrD,IAAI,IAAI,KAAK;IAG3B,yCAAyC;IACzC,IAAW,WAAW,IAAI,OAAO,CAEhC;IACD,uFAAuF;IACvF,IAAW,YAAY,IAAI,OAAO,CAEjC;IACD,wGAAwG;IACxG,IAAW,wBAAwB,IAAI,OAAO,CAE7C;IACD,gEAAgE;WAClD,2BAA2B,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGjE,mEAAmE;WACrD,4BAA4B,CAAC,OAAO,EAAE,MAAM,GAAG,KAAK;IAGlE;;;;;OAKG;WACW,oCAAoC,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EACnF,SAAS,GAAE,MAAmC,GAAG,OAAO;IAW1D;;;;OAIG;IACI,6BAA6B,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG3G;;;;;OAKG;WACW,iCAAiC,CAAC,QAAQ,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAChF,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG1D;;;;OAIG;IACI,0BAA0B,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAGxG;;;;;OAKG;IACI,aAAa,CAAC,KAAK,EAAE,KAAK,EAAE,SAAS,GAAE,MAAmC,GAAG,OAAO;IAG3F;;;;;OAKG;WACW,qBAAqB,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,OAAO;IAKzF;;;;;;OAMG;WACW,+BAA+B,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,UAAU;IA4CzF,8EAA8E;WAChE,gBAAgB,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,GAAE,MAAgB,GAAG,MAAM;IAYlF;;;;;;;;;;;;;OAaG;WACW,gCAAgC,CAC5C,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,SAAS,GAAE,OAAc,GACrE,UAAU;IAQb;;;;;;;;;OASG;WACW,wBAAwB,CACpC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GACrE,MAAM;IAIT;;;;;;;;;;;;;;;;;;OAkBG;WACW,gCAAgC,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EACnH,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EACvD,gBAAgB,GAAE,OAAe,GAAG,MAAM;IAmB5C;;;OAGG;IACI,uBAAuB,CAAC,QAAQ,EAAE,MAAM;CAQhD"}

package/lib/cjs/geometry3d/Angle.js

@@ -128,11 +128,11 @@ class Angle {
         }
     }
     /**
-   * Create an Angle from a JSON object
-   * @param json object from JSON.parse. If a number, value is in *DEGREES*
-   * @param defaultValRadians if json is undefined, default value in radians.
-   * @return a new Angle
-   */
+     * Create an Angle from a JSON object
+     * @param json object from JSON.parse. If a number, value is in *DEGREES*
+     * @param defaultValRadians if json is undefined, default value in radians.
+     * @return a new Angle
+     */
     static fromJSON(json, defaultValRadians) {
         const val = new Angle();
         val.setFromJSON(json, defaultValRadians);
@@ -187,21 +187,15 @@ class Angle {
         // all larger radians reference from 360 degrees (2PI)
         return 360.0 + 180 * ((radians - 2.0 * pi) / pi);
     }
-    /**
-     * Return the cosine of this Angle object's angle.
-     */
+    /** Return the cosine of this Angle object's angle */
     cos() {
         return Math.cos(this._radians);
     }
-    /**
-     * Return the sine of this Angle object's angle.
-     */
+    /** Return the sine of this Angle object's angle */
     sin() {
         return Math.sin(this._radians);
     }
-    /**
-     * Return the tangent of this Angle object's angle.
-     */
+    /** Return the tangent of this Angle object's angle */
     tan() {
         return Math.tan(this._radians);
     }
@@ -213,11 +207,11 @@ class Angle {
     static isHalfCircleRadians(radians) {
         return Math.abs(Math.abs(radians) - Math.PI) <= Geometry_1.Geometry.smallAngleRadians;
     }
-    /** test if the angle is a full circle */
+    /** Test if the angle is a full circle */
     get isFullCircle() {
         return Angle.isFullCircleRadians(this._radians);
     }
-    /** test if the angle is a half circle (in either direction) */
+    /** Test if the angle is a half circle (in either direction) */
     get isHalfCircle() {
         return Angle.isHalfCircleRadians(this._radians);
     }
@@ -286,7 +280,7 @@ class Angle {
         // fall through for NaN disaster.
         return 0;
     }
-    /** return a (newly allocated) Angle object with value 0 radians */
+    /** Return a (newly allocated) Angle object with value 0 radians */
     static zero() {
         return new Angle(0);
     }
@@ -440,19 +434,19 @@ class Angle {
         return value;
     }
     /**
-       * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define
-       * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The
-       * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.
-       * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.
-       * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is
-       * perpendicular to the ellipse.
-       * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)
-       * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.
-       * math details can be found at docs/learning/geometry/Angle.md
-       * @param dotUU dot product of vectorU with itself
-       * @param dotVV dot product of vectorV with itself
-       * @param dotUV dot product of vectorU with vectorV
-       */
+     * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define
+     * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The
+     * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.
+     * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.
+     * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is
+     * perpendicular to the ellipse.
+     * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)
+     * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.
+     * math details can be found at docs/learning/geometry/Angle.md
+     * @param dotUU dot product of vectorU with itself
+     * @param dotVV dot product of vectorV with itself
+     * @param dotUV dot product of vectorU with vectorV
+     */
     static dotProductsToHalfAngleTrigValues(dotUU, dotVV, dotUV, favorZero = true) {
         const cos2t0 = dotUU - dotVV;
         const sin2t0 = 2.0 * dotUV;
@@ -461,9 +455,8 @@ class Angle {
         return Angle.trigValuesToHalfAngleTrigValues(cos2t0, sin2t0);
     }
     /**
-     * * Returns the angle between two vectors, with the vectors given as xyz components
+     * Returns the angle between two vectors, with the vectors given as xyz components
      * * The returned angle is between 0 and PI
-     *
      * @param ux x component of vector u
      * @param uy y component of vector u
      * @param uz z component of vector u
@@ -476,7 +469,8 @@ class Angle {
         return Math.atan2(Geometry_1.Geometry.crossProductMagnitude(ux, uy, uz, vx, vy, vz), uDotV);
     }
     /**
-     * * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve angle to a full 2PI range.
+     * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve
+     * angle to a full 2PI range.
      * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)
      * * The angle is in the plane of the U and V vectors.
      * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.
@@ -490,7 +484,8 @@ class Angle {
      * @param upVectorX x component of vector to positive side of plane.
      * @param upVectorY y component of vector to positive side of plane.
      * @param upVectorZ z component of vector to positive side of plane.
-     * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI).  If false, return signed (-PI < radians <= PI)
+     * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI). If false, return
+     * signed (-PI < radians <= PI)
      */
     static orientedRadiansBetweenVectorsXYZ(ux, uy, uz, vx, vy, vz, upVectorX, upVectorY, upVectorZ, adjustToPositive = false) {
         const uDotV = ux * vx + uy * vy + uz * vz;
@@ -501,8 +496,8 @@ class Angle {
         const crossMagnitude = Geometry_1.Geometry.hypotenuseXYZ(wx, wy, wz);
         if (upDotW < 0.0) {
             if (adjustToPositive) {
-                // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as addition to PI.
-                // This gives the smoothest numerical transition passing PI.
+                // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as
+                // addition to PI. This gives the smoothest numerical transition passing PI.
                 return Math.PI + Math.atan2(crossMagnitude, -uDotV);
             }
             else {

package/lib/cjs/geometry3d/Angle.js.map

@@ -1 +1 @@
-{"version":3,"file":"Angle.js","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAEH,0CAAgF;AAEhF;;;;;;;;;GASG;AACH,MAAa,KAAK;IAiBhB,YAAoB,OAAO,GAAG,CAAC,EAAE,OAAgB;QAC/C,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,gDAAgD;IACzC,KAAK;QACV,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjD,CAAC;IACD,8CAA8C;IACvC,MAAM;QACX,OAAO,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7D,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC;IAC5B,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iBAAiB,CAAC,MAAa,EAAE,QAAgB,EAAE,MAAa;QAC5E,OAAO,IAAI,KAAK,CAAC,mBAAQ,CAAC,WAAW,CAAC,MAAM,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC;IACnF,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,KAAa;QAC9B,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,SAAS,CAAC;IAC5B,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,yCAAyC;IAClC,MAAM,CAAC,SAAS;QACrB,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,WAAW,CAAC,SAAiB,EAAE,WAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,SAAS,EAAE,WAAW,CAAC,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,OAAO,CAAC,KAAY;QACzB,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;IACjC,CAAC;IACD;;;;;;;OAOG;IACI,WAAW,CAAC,IAAiB,EAAE,iBAA0B;QAC9D,IAAI,CAAC,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1D,IAAI,CAAC,IAAI;YACP,OAAO;QACT,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE;YAC5B,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACvB;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;IACH,CAAC;IACD;;;;;KAKC;IACM,MAAM,CAAC,QAAQ,CAAC,IAAiB,EAAE,iBAA0B;QAClE,MAAM,GAAG,GAAG,IAAI,KAAK,EAAE,CAAC;QACxB,GAAG,CAAC,WAAW,CAAC,IAAI,EAAE,iBAAiB,CAAC,CAAC;QACzC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,+DAA+D;IACxD,MAAM;QACX,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IACD,yEAAyE;IAClE,aAAa;QAClB,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7F,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,OAAO,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IACjC,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,IAAI,OAAO,GAAG,CAAC;YACb,OAAO,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,OAAO,CAAC,CAAC;QAC3C,8BAA8B;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,CAAC;QAC1B;;;;YAII;QACJ,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,MAAM,GAAG,OAAO,CAAC;QAC1B,IAAI,OAAO,GAAG,IAAI,GAAG,EAAE;YACrB,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAClD,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QACnD,sDAAsD;QACtD,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACnD,CAAC;IACD;;OAEG;IACI,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD;;OAEG;IACI,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD;;OAEG;IACI,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,kEAAkE;IAC3D,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,mBAAQ,CAAC,gCAAgC,CAAC;IACxE,CAAC;IACD,iDAAiD;IAC1C,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,mBAAQ,CAAC,iBAAiB,CAAC;IAC7E,CAAC;IACD,yCAAyC;IACzC,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,+DAA+D;IAC/D,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,MAAM,OAAO,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;YACpD,OAAO,KAAK,GAAG,OAAO,CAAC;SACxB;QACD,sCAAsC;QACtC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,sBAAsB,CAAC,OAAe;QAClD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,KAAK;YAC5B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;YAC9D,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC,OAAO,CAAC,CAAC;SAChD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;SAC5D;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,0DAA0D;IACnD,MAAM,CAAC,0BAA0B,CAAC,OAAe;QACtD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,EAAE;YAC9B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;YAChE,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,0BAA0B,CAAC,CAAC,OAAO,CAAC,CAAC;SACpD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,IAAI;QAChB,OAAO,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,yCAAyC;IACzC,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,OAAO,KAAK,CAAC,CAAC;IAC5B,CAAC;IACD,uFAAuF;IACvF,IAAW,YAAY;QACrB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,mBAAQ,CAAC,iBAAiB,CAAC;IAC7D,CAAC;IACD,wGAAwG;IACxG,IAAW,wBAAwB;QACjC,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,MAAM,CAAC,2BAA2B,CAAC,OAAe;QACvD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,mBAAmB,CAAC,OAAO,CAAC,CAAC,CAAC;IACjE,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,4BAA4B,CAAC,OAAe;QACxD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC,OAAO,CAAC,CAAC,CAAC;IACpE,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,oCAAoC,CAAC,QAAgB,EAAE,QAAgB,EACnF,YAAoB,mBAAQ,CAAC,iBAAiB;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC;QAC5C,IAAI,KAAK,IAAI,SAAS;YACpB,OAAO,IAAI,CAAC;QACd,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC7B,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,SAAS;YACvC,OAAO,IAAI,CAAC;QACd,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,KAAK,GAAG,SAAS,GAAG,MAAM,CAAC;QAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,SAAS,CAAC;IACvC,CAAC;IACD;;;;OAIG;IACI,6BAA6B,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC/F,OAAO,KAAK,CAAC,oCAAoC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC9F,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iCAAiC,CAAC,QAAgB,EAAE,QAAgB,EAChF,YAAoB,mBAAQ,CAAC,iBAAiB;QAC9C,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,0BAA0B,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC5F,OAAO,KAAK,CAAC,iCAAiC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC3F,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC/E,OAAO,IAAI,CAAC,0BAA0B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,qBAAqB,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa;QAC7E,OAAO,KAAK,GAAG,mBAAQ,CAAC,0BAA0B;eAC7C,KAAK,GAAG,mBAAQ,CAAC,0BAA0B;eAC3C,KAAK,GAAG,KAAK,IAAI,mBAAQ,CAAC,wBAAwB,GAAG,KAAK,GAAG,KAAK,CAAC;IAC1E,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,+BAA+B,CAAC,MAAc,EAAE,MAAc;QAC1E,MAAM,CAAC,GAAG,mBAAQ,CAAC,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,mBAAQ,CAAC,mBAAmB,EAAE;YACpC,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC,YAAY;SACtD;aAAM;YACL;;;;eAIG;YACH,IAAI,IAAI,CAAC;YACT,IAAI,IAAI,GAAG,GAAG,CAAC;YACf,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,qEAAqE;YACrE,IAAI,KAAK,IAAI,GAAG,EAAE;gBAChB;;;;mBAIG;gBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;gBACtC,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;aAC7B;iBAAM;gBACL,2DAA2D;gBAC3D,IAAI,KAAK,GAAG,GAAG,EAAE;oBACf;;;;uBAIG;oBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBACtC,2DAA2D;iBAC5D;qBAAM;oBACL,kDAAkD;oBAClD,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;iBACxC;gBACD,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB;aAChD;YACD,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;SAC9D;IACH,CAAC;IACD,8EAA8E;IACvE,MAAM,CAAC,gBAAgB,CAAC,KAAa,EAAE,YAAoB,OAAO;QACvE,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,IAAI,SAAS;YACvB,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;;;;SAaK;IACE,MAAM,CAAC,gCAAgC,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa,EAAE,YAAqB,IAAI;QAEnH,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;QAC7B,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QAC3B,IAAI,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,mBAAQ,CAAC,iBAAiB,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAClG,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,KAAK,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;;OAUG;IACI,MAAM,CAAC,wBAAwB,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAC3G,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,IAAI,CAAC,KAAK,CAAC,mBAAQ,CAAC,qBAAqB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACnF,CAAC;IACD;;;;;;;;;;;;;;;;OAgBG;IACI,MAAM,CAAC,gCAAgC,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EACnH,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,mBAA4B,KAAK;QACjC,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,CAAC;QAChE,MAAM,cAAc,GAAG,mBAAQ,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC1D,IAAI,MAAM,GAAG,GAAG,EAAE;YAChB,IAAI,gBAAgB,EAAE;gBACpB,wHAAwH;gBACxH,4DAA4D;gBAC5D,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,CAAC,KAAK,CAAC,CAAC;aACrD;iBAAM;gBACL,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;aAC3C;SACF;aAAM;YACL,OAAO,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IAED;;;OAGG;IACI,uBAAuB,CAAC,QAAgB;QAC7C,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,EAAE;YAC/B,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC;YAClC,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SACvD;aAAM;YACL,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;SAC9C;IACH,CAAC;;AAjgBD,+DAA+D;AACxC,qBAAe,GAAG,mBAAmB,CAAC;AAC7D,8DAA8D;AACvC,oBAAc,GAAG,sBAAsB,CAAC;AAC/D,8DAA8D;AACvC,oBAAc,GAAG,uBAAuB,CAAC;AAChE,6DAA6D;AACtC,eAAS,GAAG,sBAAsB,CAAC;AAC1D,+DAA+D;AACxC,gBAAU,GAAG,sBAAsB,CAAC;AAC3D,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,cAAc,CAAC,CAAC;AACxE,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,KAAK,CAAC,cAAc,GAAG,IAAI,CAAC,CAAC;AAd7D,sBAAK","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 { AngleProps, BeJSONFunctions, Geometry, TrigValues } from \"../Geometry\";\n\n/**\n * An `Angle` carries the numeric value of an angle, with methods to allow (require!) callers to\n * be clear about whether their angle is degrees or radians.\n * * After the Angle object is created, the callers should not know or care whether it is stored in\n * `degrees` or `radians` because both are available if requested by caller.\n * * The various access method are named so that callers can specify whether untyped numbers passed in or\n * out are degrees or radians.\n * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/AngleSweep\n * @public\n */\nexport class Angle implements BeJSONFunctions {\n  /** maximal accuracy value of pi/12 (15 degrees), in radians */\n  public static readonly piOver12Radians = 0.26179938779914946;\n  /** maximal accuracy value of pi/4 (45 degrees), in radians */\n  public static readonly piOver4Radians = 7.853981633974483e-001;\n  /** maximal accuracy value of pi/2 (90 degrees), in radians */\n  public static readonly piOver2Radians = 1.5707963267948966e+000;\n  /** maximal accuracy value of pi (180 degrees), in radians */\n  public static readonly piRadians = 3.141592653589793e+000;\n  /** maximal accuracy value of 2*pi (360 degrees), in radians */\n  public static readonly pi2Radians = 6.283185307179586e+000;\n  /** scale factor for converting radians to degrees */\n  public static readonly degreesPerRadian = (45.0 / Angle.piOver4Radians);\n  /** scale factor for converting degrees to radians */\n  public static readonly radiansPerDegree = (Angle.piOver4Radians / 45.0);\n  private _radians: number;\n  private _degrees?: number;\n  private constructor(radians = 0, degrees?: number) {\n    this._radians = radians;\n    this._degrees = degrees;\n  }\n  /** Return a new angle with the same content. */\n  public clone(): Angle {\n    return new Angle(this._radians, this._degrees);\n  }\n  /** Freeze this instance so it is read-only */\n  public freeze(): Readonly<this> {\n    return Object.freeze(this);\n  }\n  /**\n   * Return a new Angle object for angle given in degrees.\n   * @param degrees angle in degrees\n   */\n  public static createDegrees(degrees: number): Angle {\n    return new Angle(Angle.degreesToRadians(degrees), degrees);\n  }\n  /**\n   * Return a (new) Angle object for a value given in radians.\n   * @param radians angle in radians\n   */\n  public static createRadians(radians: number): Angle {\n    return new Angle(radians);\n  }\n  /**\n   * Return a (new) Angle object that is interpolated between two inputs (based on a fraction)\n   * @param angle0 first angle in radians\n   * @param fraction the interpolation fraction\n   * @param angle1 second angle in radians\n   */\n  public static createInterpolate(angle0: Angle, fraction: number, angle1: Angle): Angle {\n    return new Angle(Geometry.interpolate(angle0.radians, fraction, angle1.radians));\n  }\n  /**\n   * Return a (new) Angle object, with angle scaled from existing angle.\n   * @param scale scale factor to apply to angle.\n   */\n  public cloneScaled(scale: number): Angle {\n    return new Angle(this.radians * scale);\n  }\n  /**\n   * Set this angle to a value given in radians.\n   * @param radians angle given in radians\n   */\n  public setRadians(radians: number) {\n    this._radians = radians;\n    this._degrees = undefined;\n  }\n  /**\n   * Set this angle to a value given in degrees.\n   * @param degrees angle given in degrees.\n   */\n  public setDegrees(degrees: number) {\n    this._radians = Angle.degreesToRadians(degrees);\n    this._degrees = degrees;\n  }\n  /** Create an angle for a full circle. */\n  public static create360(): Angle {\n    return new Angle(Math.PI * 2.0, 360.0);\n  }\n  /**\n   * Create a (strongly typed) Angle whose tangent is `numerator/denominator`, using the signs of both in\n   * determining the (otherwise ambiguous) quadrant.\n   * @param numerator numerator for tangent\n   * @param denominator denominator for tangent\n   */\n  public static createAtan2(numerator: number, denominator: number): Angle {\n    return new Angle(Math.atan2(numerator, denominator));\n  }\n  /**\n   * Copy all contents of `other` to this Angle.\n   * @param other source data\n   */\n  public setFrom(other: Angle) {\n    this._radians = other._radians;\n    this._degrees = other._degrees;\n  }\n  /**\n   * Set an Angle from a JSON object\n   * * A simple number is considered as degrees.\n   * * specified `json.degrees` or `json._degrees` is degree value.\n   * * specified `son.radians` or `json._radians` is radians value.\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\n   * @param defaultValRadians if json is undefined, default value in radians.\n   */\n  public setFromJSON(json?: AngleProps, defaultValRadians?: number) {\n    this._radians = defaultValRadians ? defaultValRadians : 0;\n    if (!json)\n      return;\n    if (typeof json === \"number\") {\n      this.setDegrees(json);\n    } else if (typeof (json as any).degrees === \"number\") {\n      this.setDegrees((json as any).degrees);\n    } else if (typeof (json as any)._degrees === \"number\") {\n      this.setDegrees((json as any)._degrees);\n    } else if (typeof (json as any).radians === \"number\") {\n      this.setRadians((json as any).radians);\n    } else if (typeof (json as any)._radians === \"number\") {\n      this.setRadians((json as any)._radians);\n    }\n  }\n  /**\n * Create an Angle from a JSON object\n * @param json object from JSON.parse. If a number, value is in *DEGREES*\n * @param defaultValRadians if json is undefined, default value in radians.\n * @return a new Angle\n */\n  public static fromJSON(json?: AngleProps, defaultValRadians?: number): Angle {\n    const val = new Angle();\n    val.setFromJSON(json, defaultValRadians);\n    return val;\n  }\n  /** Convert an Angle to a JSON object as a number in degrees */\n  public toJSON(): AngleProps {\n    return this.degrees;\n  }\n  /** Return a json object with radians keyword, e.g. `{ radians: 0.10}` */\n  public toJSONRadians(): AngleProps {\n    return { radians: this.radians };\n  }\n  /** Return the angle measured in radians. */\n  public get radians(): number {\n    return this._radians;\n  }\n  /** Return the angle measured in degrees. */\n  public get degrees(): number {\n    return this._degrees !== undefined ? this._degrees : Angle.radiansToDegrees(this._radians);\n  }\n  /**\n   * Convert an angle in degrees to radians.\n   * @param degrees angle in degrees\n   */\n  public static degreesToRadians(degrees: number): number {\n    return degrees * Math.PI / 180;\n  }\n  /**\n   * Convert an angle in radians to degrees.\n   * @param degrees angle in radians\n   */\n  public static radiansToDegrees(radians: number): number {\n    if (radians < 0)\n      return -Angle.radiansToDegrees(-radians);\n    // Now radians is positive ...\n    const pi = Math.PI;\n    const factor = 180.0 / pi;\n    /* the following if statements are for round-off reasons. The problem is that no IEEE number is\n      * an exact hit for any primary multiple of pi (90, 180, etc). The following is supposed to have\n      * a better chance that if the input was computed by direct assignment from 90, 180, etc degrees\n      * it will return exactly 90,180 etc.\n      */\n    if (radians <= 0.25 * pi)\n      return factor * radians;\n    if (radians < 0.75 * pi)\n      return 90.0 + 180 * ((radians - 0.5 * pi) / pi);\n    if (radians <= 1.25 * pi)\n      return 180.0 + 180 * ((radians - pi) / pi);\n    if (radians <= 1.75 * pi)\n      return 270.0 + 180 * ((radians - 1.5 * pi) / pi);\n    // all larger radians reference from 360 degrees (2PI)\n    return 360.0 + 180 * ((radians - 2.0 * pi) / pi);\n  }\n  /**\n   * Return the cosine of this Angle object's angle.\n   */\n  public cos(): number {\n    return Math.cos(this._radians);\n  }\n  /**\n   * Return the sine of this Angle object's angle.\n   */\n  public sin(): number {\n    return Math.sin(this._radians);\n  }\n  /**\n   * Return the tangent of this Angle object's angle.\n   */\n  public tan(): number {\n    return Math.tan(this._radians);\n  }\n  /** Test if a radians (absolute) value is nearly 2PI or larger! */\n  public static isFullCircleRadians(radians: number): boolean {\n    return Math.abs(radians) >= Geometry.fullCircleRadiansMinusSmallAngle;\n  }\n  /** Test if the radians value is a half circle */\n  public static isHalfCircleRadians(radians: number): boolean {\n    return Math.abs(Math.abs(radians) - Math.PI) <= Geometry.smallAngleRadians;\n  }\n  /** test if the angle is a full circle */\n  public get isFullCircle(): boolean {\n    return Angle.isFullCircleRadians(this._radians);\n  }\n  /** test if the angle is a half circle (in either direction) */\n  public get isHalfCircle(): boolean {\n    return Angle.isHalfCircleRadians(this._radians);\n  }\n  /** Adjust a radians value so it is positive in 0..360 */\n  public static adjustDegrees0To360(degrees: number): number {\n    if (degrees >= 0) {\n      const period = 360.0;\n      if (degrees < period)\n        return degrees;\n      const numPeriods = Math.floor(degrees / period);\n      return degrees - numPeriods * period;\n    } else if (degrees < 0) {\n      // negative angle ...\n      const radians = Angle.adjustDegrees0To360(-degrees);\n      return 360.0 - radians;\n    }\n    // fall through for Nan (disaster) !!!\n    return 0;\n  }\n  /** Adjust a radians value so it is in -180..180 */\n  public static adjustDegreesSigned180(degrees: number): number {\n    if (Math.abs(degrees) <= 180.0)\n      return degrees;\n    if (degrees >= 0) {\n      const period = 360.0;\n      const numPeriods = 1 + Math.floor((degrees - 180.0) / period);\n      return degrees - numPeriods * period;\n    } else if (degrees < 0) {\n      // negative angle ...\n      return -Angle.adjustDegreesSigned180(-degrees);\n    }\n    // fall through for NaN disaster.\n    return 0;\n  }\n  /** Adjust a radians value so it is positive in 0..2Pi */\n  public static adjustRadians0To2Pi(radians: number): number {\n    if (radians >= 0) {\n      const period = Math.PI * 2.0;\n      if (radians < period)\n        return radians;\n      const numPeriods = Math.floor(radians / period);\n      return radians - numPeriods * period;\n    } else if (radians < 0) {\n      // negative angle ...\n      return Math.PI * 2.0 - Angle.adjustRadians0To2Pi(-radians);\n    }\n    // fall through for NaN disaster.\n    return 0;\n  }\n  /** Adjust a radians value so it is positive in -PI..PI */\n  public static adjustRadiansMinusPiPlusPi(radians: number): number {\n    if (Math.abs(radians) <= Math.PI)\n      return radians;\n    if (radians >= 0) {\n      const period = Math.PI * 2.0;\n      const numPeriods = 1 + Math.floor((radians - Math.PI) / period);\n      return radians - numPeriods * period;\n    } else if (radians < 0) {\n      // negative angle ...\n      return -Angle.adjustRadiansMinusPiPlusPi(-radians);\n    }\n    // fall through for NaN disaster.\n    return 0;\n  }\n  /** return a (newly allocated) Angle object with value 0 radians */\n  public static zero(): Angle {\n    return new Angle(0);\n  }\n  /** Test if the angle is exactly zero. */\n  public get isExactZero(): boolean {\n    return this.radians === 0;\n  }\n  /** Test if the angle is almost zero (within tolerance `Geometry.smallAngleRadians`) */\n  public get isAlmostZero(): boolean {\n    return Math.abs(this.radians) < Geometry.smallAngleRadians;\n  }\n  /** Test if the angle is almost a north or south pole (within tolerance `Geometry.smallAngleRadians`) */\n  public get isAlmostNorthOrSouthPole(): boolean {\n    return Angle.isHalfCircleRadians(this.radians * 2.0);\n  }\n  /** Create an angle object with degrees adjusted into 0..360. */\n  public static createDegreesAdjustPositive(degrees: number): Angle {\n    return Angle.createDegrees(Angle.adjustDegrees0To360(degrees));\n  }\n  /** Create an angle object with degrees adjusted into -180..180. */\n  public static createDegreesAdjustSigned180(degrees: number): Angle {\n    return Angle.createDegrees(Angle.adjustDegreesSigned180(degrees));\n  }\n  /**\n   * Test if two radians values are equivalent, allowing shift by full circle (i.e. by a multiple of `2*PI`)\n   * @param radiansA first radians value\n   * @param radiansB second radians value\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\n   */\n  public static isAlmostEqualRadiansAllowPeriodShift(radiansA: number, radiansB: number,\n    radianTol: number = Geometry.smallAngleRadians): boolean {\n    const delta = Math.abs(radiansA - radiansB);\n    if (delta <= radianTol)\n      return true;\n    const period = Math.PI * 2.0;\n    if (Math.abs(delta - period) <= radianTol)\n      return true;\n    const numPeriod = Math.round(delta / period);\n    const delta1 = delta - numPeriod * period;\n    return Math.abs(delta1) <= radianTol;\n  }\n  /**\n   * Test if this angle and other are equivalent, allowing shift by full circle (i.e., multiples of `2 * PI`).\n   * @param other the other angle\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\n   */\n  public isAlmostEqualAllowPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\n    return Angle.isAlmostEqualRadiansAllowPeriodShift(this._radians, other._radians, radianTol);\n  }\n  /**\n   * Test if two angle (in radians)  almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\n   * @param radiansA first radians value\n   * @param radiansB second radians value\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\n   */\n  public static isAlmostEqualRadiansNoPeriodShift(radiansA: number, radiansB: number,\n    radianTol: number = Geometry.smallAngleRadians): boolean {\n    return Math.abs(radiansA - radiansB) < radianTol;\n  }\n  /**\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\n   * @param other the other angle\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\n   */\n  public isAlmostEqualNoPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\n    return Angle.isAlmostEqualRadiansNoPeriodShift(this._radians, other._radians, radianTol);\n  }\n  /**\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\n   * * This function is same as isAlmostEqualRadiansNoPeriodShift. Please use isAlmostEqualRadiansNoPeriodShift.\n   * @param other the other angle\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\n   */\n  public isAlmostEqual(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\n    return this.isAlmostEqualNoPeriodShift(other, radianTol);\n  }\n  /**\n   * Test if dot product values indicate non-zero length perpendicular vectors.\n   * @param dotUU dot product of vectorU with itself\n   * @param dotVV dot product of vectorV with itself\n   * @param dotUV dot product of vectorU with vectorV\n   */\n  public static isPerpendicularDotSet(dotUU: number, dotVV: number, dotUV: number): boolean {\n    return dotUU > Geometry.smallMetricDistanceSquared\n      && dotVV > Geometry.smallMetricDistanceSquared\n      && dotUV * dotUV <= Geometry.smallAngleRadiansSquared * dotUU * dotVV;\n  }\n  /**\n   * Return cosine, sine, and radians for the half angle of a \"cosine,sine\" pair.\n   * * This function assumes the input arguments are related to an angle between -PI and PI\n   * * This function returns an angle between -PI and PI\n   * @param rCos2A cosine value (scaled by radius) for initial angle.\n   * @param rSin2A sine value (scaled by radius) for final angle.\n   */\n  public static trigValuesToHalfAngleTrigValues(rCos2A: number, rSin2A: number): TrigValues {\n    const r = Geometry.hypotenuseXY(rCos2A, rSin2A);\n    if (r < Geometry.smallMetricDistance) {\n      return { c: 1.0, s: 0.0, radians: 0.0 }; // angle = 0\n    } else {\n      /* If the caller really gave you sine and cosine values, r should be 1.  However,\n       * to allow scaled values -- e.g. the x and y components of any vector -- we normalize\n       * right here. This adds an extra sqrt and 2 divides to the whole process, but improves\n       * both the usefulness and robustness of the computation.\n       */\n      let cosA;\n      let sinA = 0.0;\n      const cos2A = rCos2A / r;\n      const sin2A = rSin2A / r;\n      // Original angle in NE and SE quadrants. Half angle in same quadrant\n      if (cos2A >= 0.0) {\n        /*\n         * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\n         * so 1 + cos2A = 2(cosA)^2 and therefore, cosA = sqrt((1+cos2A)/2)\n         * cosine is positive in NE and SE quadrants so we use +sqrt\n         */\n        cosA = Math.sqrt(0.5 * (1.0 + cos2A));\n        // We know sin2A = 2 sinA cosA so sinA = sin2A/(2*cosA)\n        sinA = sin2A / (2.0 * cosA);\n      } else {\n        // Original angle in NW quadrant. Half angle in NE quadrant\n        if (sin2A > 0.0) {\n          /*\n           * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\n           * so 1 - cos2A = 2(sinA)^2 and therefore, sinA = sqrt((1-cos2A)/2)\n           * sine is positive in NE quadrant so we use +sqrt\n           */\n          sinA = Math.sqrt(0.5 * (1.0 - cos2A));\n          // Original angle in SW quadrant. Half angle in SE quadrant\n        } else {\n          // sine is negative in SE quadrant so we use -sqrt\n          sinA = -Math.sqrt(0.5 * (1.0 - cos2A));\n        }\n        // We know sin2A = 2 sinA cosA so cosA = sin2A/(2*sinA)\n        cosA = sin2A / (2.0 * sinA); // always positive\n      }\n      return { c: cosA, s: sinA, radians: Math.atan2(sinA, cosA) };\n    }\n  }\n  /** If value is close to -1, -0.5, 0, 0.5, 1, adjust it to the exact value. */\n  public static cleanupTrigValue(value: number, tolerance: number = 1.0e-15): number {\n    const absValue = Math.abs(value);\n    if (absValue <= tolerance)\n      return 0;\n    let a = Math.abs(absValue - 0.5);\n    if (a <= tolerance)\n      return value < 0.0 ? -0.5 : 0.5;\n    a = Math.abs(absValue - 1.0);\n    if (a <= tolerance)\n      return value < 0.0 ? -1.0 : 1.0;\n    return value;\n  }\n  /**\n     * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define\n     * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The\n     * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.\n     * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.\n     * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is\n     * perpendicular to the ellipse.\n     * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)\n     * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.\n     * math details can be found at docs/learning/geometry/Angle.md\n     * @param dotUU dot product of vectorU with itself\n     * @param dotVV dot product of vectorV with itself\n     * @param dotUV dot product of vectorU with vectorV\n     */\n  public static dotProductsToHalfAngleTrigValues(dotUU: number, dotVV: number, dotUV: number, favorZero: boolean = true): TrigValues {\n\n    const cos2t0 = dotUU - dotVV;\n    const sin2t0 = 2.0 * dotUV;\n    if (favorZero && Math.abs(sin2t0) < Geometry.smallAngleRadians * (Math.abs(dotUU) + Math.abs(dotVV)))\n      return { c: 1.0, s: 0.0, radians: 0.0 };\n    return Angle.trigValuesToHalfAngleTrigValues(cos2t0, sin2t0);\n  }\n  /**\n   * * Returns the angle between two vectors, with the vectors given as xyz components\n   * * The returned angle is between 0 and PI\n   *\n   * @param ux x component of vector u\n   * @param uy y component of vector u\n   * @param uz z component of vector u\n   * @param vx x component of vector v\n   * @param vy y component of vector v\n   * @param vz z component of vector v\n   */\n  public static radiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number): number {\n    const uDotV = ux * vx + uy * vy + uz * vz;\n    return Math.atan2(Geometry.crossProductMagnitude(ux, uy, uz, vx, vy, vz), uDotV);\n  }\n  /**\n   * * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve angle to a full 2PI range.\n   * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)\n   * * The angle is in the plane of the U and V vectors.\n   * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.\n   *\n   * @param ux x component of vector u\n   * @param uy y component of vector u\n   * @param uz z component of vector u\n   * @param vx x component of vector v\n   * @param vy y component of vector v\n   * @param vz z component of vector v\n   * @param upVectorX x component of vector to positive side of plane.\n   * @param upVectorY y component of vector to positive side of plane.\n   * @param upVectorZ z component of vector to positive side of plane.\n   * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI).  If false, return signed (-PI < radians <= PI)\n   */\n  public static orientedRadiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number,\n    upVectorX: number, upVectorY: number, upVectorZ: number,\n    adjustToPositive: boolean = false): number {\n    const uDotV = ux * vx + uy * vy + uz * vz;\n    const wx = uy * vz - uz * vy;\n    const wy = uz * vx - ux * vz;\n    const wz = ux * vy - uy * vx;\n    const upDotW = upVectorX * wx + upVectorY * wy + upVectorZ * wz;\n    const crossMagnitude = Geometry.hypotenuseXYZ(wx, wy, wz);\n    if (upDotW < 0.0) {\n      if (adjustToPositive) {\n        // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as addition to PI.\n        // This gives the smoothest numerical transition passing PI.\n        return Math.PI + Math.atan2(crossMagnitude, -uDotV);\n      } else {\n        return -Math.atan2(crossMagnitude, uDotV);\n      }\n    } else {\n      return Math.atan2(crossMagnitude, uDotV);\n    }\n  }\n\n  /**\n   * Add a multiple of a full circle angle (360 degrees, 2PI) in place.\n   * @param multiple multiplier factor\n   */\n  public addMultipleOf2PiInPlace(multiple: number) {\n    if (this._degrees !== undefined) {\n      this._degrees += multiple * 360.0;\n      this._radians = Angle.degreesToRadians(this._degrees);\n    } else {\n      this._radians += multiple * Angle.pi2Radians;\n    }\n  }\n}\n"]}
+{"version":3,"file":"Angle.js","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAEH,0CAAgF;AAEhF;;;;;;;;;GASG;AACH,MAAa,KAAK;IAiBhB,YAAoB,OAAO,GAAG,CAAC,EAAE,OAAgB;QAC/C,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,gDAAgD;IACzC,KAAK;QACV,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjD,CAAC;IACD,8CAA8C;IACvC,MAAM;QACX,OAAO,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7D,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC;IAC5B,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iBAAiB,CAAC,MAAa,EAAE,QAAgB,EAAE,MAAa;QAC5E,OAAO,IAAI,KAAK,CAAC,mBAAQ,CAAC,WAAW,CAAC,MAAM,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC;IACnF,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,KAAa;QAC9B,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,SAAS,CAAC;IAC5B,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,yCAAyC;IAClC,MAAM,CAAC,SAAS;QACrB,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,WAAW,CAAC,SAAiB,EAAE,WAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,SAAS,EAAE,WAAW,CAAC,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,OAAO,CAAC,KAAY;QACzB,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;IACjC,CAAC;IACD;;;;;;;OAOG;IACI,WAAW,CAAC,IAAiB,EAAE,iBAA0B;QAC9D,IAAI,CAAC,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1D,IAAI,CAAC,IAAI;YACP,OAAO;QACT,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE;YAC5B,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACvB;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;IACH,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,QAAQ,CAAC,IAAiB,EAAE,iBAA0B;QAClE,MAAM,GAAG,GAAG,IAAI,KAAK,EAAE,CAAC;QACxB,GAAG,CAAC,WAAW,CAAC,IAAI,EAAE,iBAAiB,CAAC,CAAC;QACzC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,+DAA+D;IACxD,MAAM;QACX,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IACD,yEAAyE;IAClE,aAAa;QAClB,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7F,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,OAAO,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IACjC,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,IAAI,OAAO,GAAG,CAAC;YACb,OAAO,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,OAAO,CAAC,CAAC;QAC3C,8BAA8B;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,CAAC;QAC1B;;;;YAII;QACJ,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,MAAM,GAAG,OAAO,CAAC;QAC1B,IAAI,OAAO,GAAG,IAAI,GAAG,EAAE;YACrB,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAClD,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QACnD,sDAAsD;QACtD,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACnD,CAAC;IACD,qDAAqD;IAC9C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,mDAAmD;IAC5C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,sDAAsD;IAC/C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,kEAAkE;IAC3D,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,mBAAQ,CAAC,gCAAgC,CAAC;IACxE,CAAC;IACD,iDAAiD;IAC1C,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,mBAAQ,CAAC,iBAAiB,CAAC;IAC7E,CAAC;IACD,yCAAyC;IACzC,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,+DAA+D;IAC/D,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,MAAM,OAAO,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;YACpD,OAAO,KAAK,GAAG,OAAO,CAAC;SACxB;QACD,sCAAsC;QACtC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,sBAAsB,CAAC,OAAe;QAClD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,KAAK;YAC5B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;YAC9D,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC,OAAO,CAAC,CAAC;SAChD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;SAC5D;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,0DAA0D;IACnD,MAAM,CAAC,0BAA0B,CAAC,OAAe;QACtD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,EAAE;YAC9B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;YAChE,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,0BAA0B,CAAC,CAAC,OAAO,CAAC,CAAC;SACpD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,IAAI;QAChB,OAAO,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,yCAAyC;IACzC,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,OAAO,KAAK,CAAC,CAAC;IAC5B,CAAC;IACD,uFAAuF;IACvF,IAAW,YAAY;QACrB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,mBAAQ,CAAC,iBAAiB,CAAC;IAC7D,CAAC;IACD,wGAAwG;IACxG,IAAW,wBAAwB;QACjC,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,MAAM,CAAC,2BAA2B,CAAC,OAAe;QACvD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,mBAAmB,CAAC,OAAO,CAAC,CAAC,CAAC;IACjE,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,4BAA4B,CAAC,OAAe;QACxD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC,OAAO,CAAC,CAAC,CAAC;IACpE,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,oCAAoC,CAAC,QAAgB,EAAE,QAAgB,EACnF,YAAoB,mBAAQ,CAAC,iBAAiB;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC;QAC5C,IAAI,KAAK,IAAI,SAAS;YACpB,OAAO,IAAI,CAAC;QACd,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC7B,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,SAAS;YACvC,OAAO,IAAI,CAAC;QACd,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,KAAK,GAAG,SAAS,GAAG,MAAM,CAAC;QAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,SAAS,CAAC;IACvC,CAAC;IACD;;;;OAIG;IACI,6BAA6B,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC/F,OAAO,KAAK,CAAC,oCAAoC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC9F,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iCAAiC,CAAC,QAAgB,EAAE,QAAgB,EAChF,YAAoB,mBAAQ,CAAC,iBAAiB;QAC9C,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,0BAA0B,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC5F,OAAO,KAAK,CAAC,iCAAiC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC3F,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY,EAAE,YAAoB,mBAAQ,CAAC,iBAAiB;QAC/E,OAAO,IAAI,CAAC,0BAA0B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,qBAAqB,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa;QAC7E,OAAO,KAAK,GAAG,mBAAQ,CAAC,0BAA0B;eAC7C,KAAK,GAAG,mBAAQ,CAAC,0BAA0B;eAC3C,KAAK,GAAG,KAAK,IAAI,mBAAQ,CAAC,wBAAwB,GAAG,KAAK,GAAG,KAAK,CAAC;IAC1E,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,+BAA+B,CAAC,MAAc,EAAE,MAAc;QAC1E,MAAM,CAAC,GAAG,mBAAQ,CAAC,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,mBAAQ,CAAC,mBAAmB,EAAE;YACpC,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC,YAAY;SACtD;aAAM;YACL;;;;eAIG;YACH,IAAI,IAAI,CAAC;YACT,IAAI,IAAI,GAAG,GAAG,CAAC;YACf,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,qEAAqE;YACrE,IAAI,KAAK,IAAI,GAAG,EAAE;gBAChB;;;;mBAIG;gBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;gBACtC,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;aAC7B;iBAAM;gBACL,2DAA2D;gBAC3D,IAAI,KAAK,GAAG,GAAG,EAAE;oBACf;;;;uBAIG;oBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBACtC,2DAA2D;iBAC5D;qBAAM;oBACL,kDAAkD;oBAClD,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;iBACxC;gBACD,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB;aAChD;YACD,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;SAC9D;IACH,CAAC;IACD,8EAA8E;IACvE,MAAM,CAAC,gBAAgB,CAAC,KAAa,EAAE,YAAoB,OAAO;QACvE,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,IAAI,SAAS;YACvB,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;;;;OAaG;IACI,MAAM,CAAC,gCAAgC,CAC5C,KAAa,EAAE,KAAa,EAAE,KAAa,EAAE,YAAqB,IAAI;QAGtE,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;QAC7B,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QAC3B,IAAI,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,mBAAQ,CAAC,iBAAiB,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAClG,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,KAAK,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,wBAAwB,CACpC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAEtE,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,IAAI,CAAC,KAAK,CAAC,mBAAQ,CAAC,qBAAqB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACnF,CAAC;IACD;;;;;;;;;;;;;;;;;;OAkBG;IACI,MAAM,CAAC,gCAAgC,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EACnH,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,mBAA4B,KAAK;QACjC,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,CAAC;QAChE,MAAM,cAAc,GAAG,mBAAQ,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC1D,IAAI,MAAM,GAAG,GAAG,EAAE;YAChB,IAAI,gBAAgB,EAAE;gBACpB,wGAAwG;gBACxG,4EAA4E;gBAC5E,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,CAAC,KAAK,CAAC,CAAC;aACrD;iBAAM;gBACL,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;aAC3C;SACF;aAAM;YACL,OAAO,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IACD;;;OAGG;IACI,uBAAuB,CAAC,QAAgB;QAC7C,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,EAAE;YAC/B,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC;YAClC,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SACvD;aAAM;YACL,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;SAC9C;IACH,CAAC;;AA/fD,+DAA+D;AACxC,qBAAe,GAAG,mBAAmB,CAAC;AAC7D,8DAA8D;AACvC,oBAAc,GAAG,sBAAsB,CAAC;AAC/D,8DAA8D;AACvC,oBAAc,GAAG,uBAAuB,CAAC;AAChE,6DAA6D;AACtC,eAAS,GAAG,sBAAsB,CAAC;AAC1D,+DAA+D;AACxC,gBAAU,GAAG,sBAAsB,CAAC;AAC3D,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,cAAc,CAAC,CAAC;AACxE,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,KAAK,CAAC,cAAc,GAAG,IAAI,CAAC,CAAC;AAd7D,sBAAK","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/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\n\r\nimport { AngleProps, BeJSONFunctions, Geometry, TrigValues } from \"../Geometry\";\r\n\r\n/**\r\n * An `Angle` carries the numeric value of an angle, with methods to allow (require!) callers to\r\n * be clear about whether their angle is degrees or radians.\r\n * * After the Angle object is created, the callers should not know or care whether it is stored in\r\n * `degrees` or `radians` because both are available if requested by caller.\r\n * * The various access method are named so that callers can specify whether untyped numbers passed in or\r\n * out are degrees or radians.\r\n * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/AngleSweep\r\n * @public\r\n */\r\nexport class Angle implements BeJSONFunctions {\r\n  /** maximal accuracy value of pi/12 (15 degrees), in radians */\r\n  public static readonly piOver12Radians = 0.26179938779914946;\r\n  /** maximal accuracy value of pi/4 (45 degrees), in radians */\r\n  public static readonly piOver4Radians = 7.853981633974483e-001;\r\n  /** maximal accuracy value of pi/2 (90 degrees), in radians */\r\n  public static readonly piOver2Radians = 1.5707963267948966e+000;\r\n  /** maximal accuracy value of pi (180 degrees), in radians */\r\n  public static readonly piRadians = 3.141592653589793e+000;\r\n  /** maximal accuracy value of 2*pi (360 degrees), in radians */\r\n  public static readonly pi2Radians = 6.283185307179586e+000;\r\n  /** scale factor for converting radians to degrees */\r\n  public static readonly degreesPerRadian = (45.0 / Angle.piOver4Radians);\r\n  /** scale factor for converting degrees to radians */\r\n  public static readonly radiansPerDegree = (Angle.piOver4Radians / 45.0);\r\n  private _radians: number;\r\n  private _degrees?: number;\r\n  private constructor(radians = 0, degrees?: number) {\r\n    this._radians = radians;\r\n    this._degrees = degrees;\r\n  }\r\n  /** Return a new angle with the same content. */\r\n  public clone(): Angle {\r\n    return new Angle(this._radians, this._degrees);\r\n  }\r\n  /** Freeze this instance so it is read-only */\r\n  public freeze(): Readonly<this> {\r\n    return Object.freeze(this);\r\n  }\r\n  /**\r\n   * Return a new Angle object for angle given in degrees.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static createDegrees(degrees: number): Angle {\r\n    return new Angle(Angle.degreesToRadians(degrees), degrees);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object for a value given in radians.\r\n   * @param radians angle in radians\r\n   */\r\n  public static createRadians(radians: number): Angle {\r\n    return new Angle(radians);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object that is interpolated between two inputs (based on a fraction)\r\n   * @param angle0 first angle in radians\r\n   * @param fraction the interpolation fraction\r\n   * @param angle1 second angle in radians\r\n   */\r\n  public static createInterpolate(angle0: Angle, fraction: number, angle1: Angle): Angle {\r\n    return new Angle(Geometry.interpolate(angle0.radians, fraction, angle1.radians));\r\n  }\r\n  /**\r\n   * Return a (new) Angle object, with angle scaled from existing angle.\r\n   * @param scale scale factor to apply to angle.\r\n   */\r\n  public cloneScaled(scale: number): Angle {\r\n    return new Angle(this.radians * scale);\r\n  }\r\n  /**\r\n   * Set this angle to a value given in radians.\r\n   * @param radians angle given in radians\r\n   */\r\n  public setRadians(radians: number) {\r\n    this._radians = radians;\r\n    this._degrees = undefined;\r\n  }\r\n  /**\r\n   * Set this angle to a value given in degrees.\r\n   * @param degrees angle given in degrees.\r\n   */\r\n  public setDegrees(degrees: number) {\r\n    this._radians = Angle.degreesToRadians(degrees);\r\n    this._degrees = degrees;\r\n  }\r\n  /** Create an angle for a full circle. */\r\n  public static create360(): Angle {\r\n    return new Angle(Math.PI * 2.0, 360.0);\r\n  }\r\n  /**\r\n   * Create a (strongly typed) Angle whose tangent is `numerator/denominator`, using the signs of both in\r\n   * determining the (otherwise ambiguous) quadrant.\r\n   * @param numerator numerator for tangent\r\n   * @param denominator denominator for tangent\r\n   */\r\n  public static createAtan2(numerator: number, denominator: number): Angle {\r\n    return new Angle(Math.atan2(numerator, denominator));\r\n  }\r\n  /**\r\n   * Copy all contents of `other` to this Angle.\r\n   * @param other source data\r\n   */\r\n  public setFrom(other: Angle) {\r\n    this._radians = other._radians;\r\n    this._degrees = other._degrees;\r\n  }\r\n  /**\r\n   * Set an Angle from a JSON object\r\n   * * A simple number is considered as degrees.\r\n   * * specified `json.degrees` or `json._degrees` is degree value.\r\n   * * specified `son.radians` or `json._radians` is radians value.\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   */\r\n  public setFromJSON(json?: AngleProps, defaultValRadians?: number) {\r\n    this._radians = defaultValRadians ? defaultValRadians : 0;\r\n    if (!json)\r\n      return;\r\n    if (typeof json === \"number\") {\r\n      this.setDegrees(json);\r\n    } else if (typeof (json as any).degrees === \"number\") {\r\n      this.setDegrees((json as any).degrees);\r\n    } else if (typeof (json as any)._degrees === \"number\") {\r\n      this.setDegrees((json as any)._degrees);\r\n    } else if (typeof (json as any).radians === \"number\") {\r\n      this.setRadians((json as any).radians);\r\n    } else if (typeof (json as any)._radians === \"number\") {\r\n      this.setRadians((json as any)._radians);\r\n    }\r\n  }\r\n  /**\r\n   * Create an Angle from a JSON object\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   * @return a new Angle\r\n   */\r\n  public static fromJSON(json?: AngleProps, defaultValRadians?: number): Angle {\r\n    const val = new Angle();\r\n    val.setFromJSON(json, defaultValRadians);\r\n    return val;\r\n  }\r\n  /** Convert an Angle to a JSON object as a number in degrees */\r\n  public toJSON(): AngleProps {\r\n    return this.degrees;\r\n  }\r\n  /** Return a json object with radians keyword, e.g. `{ radians: 0.10}` */\r\n  public toJSONRadians(): AngleProps {\r\n    return { radians: this.radians };\r\n  }\r\n  /** Return the angle measured in radians. */\r\n  public get radians(): number {\r\n    return this._radians;\r\n  }\r\n  /** Return the angle measured in degrees. */\r\n  public get degrees(): number {\r\n    return this._degrees !== undefined ? this._degrees : Angle.radiansToDegrees(this._radians);\r\n  }\r\n  /**\r\n   * Convert an angle in degrees to radians.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static degreesToRadians(degrees: number): number {\r\n    return degrees * Math.PI / 180;\r\n  }\r\n  /**\r\n   * Convert an angle in radians to degrees.\r\n   * @param degrees angle in radians\r\n   */\r\n  public static radiansToDegrees(radians: number): number {\r\n    if (radians < 0)\r\n      return -Angle.radiansToDegrees(-radians);\r\n    // Now radians is positive ...\r\n    const pi = Math.PI;\r\n    const factor = 180.0 / pi;\r\n    /* the following if statements are for round-off reasons. The problem is that no IEEE number is\r\n      * an exact hit for any primary multiple of pi (90, 180, etc). The following is supposed to have\r\n      * a better chance that if the input was computed by direct assignment from 90, 180, etc degrees\r\n      * it will return exactly 90,180 etc.\r\n      */\r\n    if (radians <= 0.25 * pi)\r\n      return factor * radians;\r\n    if (radians < 0.75 * pi)\r\n      return 90.0 + 180 * ((radians - 0.5 * pi) / pi);\r\n    if (radians <= 1.25 * pi)\r\n      return 180.0 + 180 * ((radians - pi) / pi);\r\n    if (radians <= 1.75 * pi)\r\n      return 270.0 + 180 * ((radians - 1.5 * pi) / pi);\r\n    // all larger radians reference from 360 degrees (2PI)\r\n    return 360.0 + 180 * ((radians - 2.0 * pi) / pi);\r\n  }\r\n  /** Return the cosine of this Angle object's angle */\r\n  public cos(): number {\r\n    return Math.cos(this._radians);\r\n  }\r\n  /** Return the sine of this Angle object's angle */\r\n  public sin(): number {\r\n    return Math.sin(this._radians);\r\n  }\r\n  /** Return the tangent of this Angle object's angle */\r\n  public tan(): number {\r\n    return Math.tan(this._radians);\r\n  }\r\n  /** Test if a radians (absolute) value is nearly 2PI or larger! */\r\n  public static isFullCircleRadians(radians: number): boolean {\r\n    return Math.abs(radians) >= Geometry.fullCircleRadiansMinusSmallAngle;\r\n  }\r\n  /** Test if the radians value is a half circle */\r\n  public static isHalfCircleRadians(radians: number): boolean {\r\n    return Math.abs(Math.abs(radians) - Math.PI) <= Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is a full circle */\r\n  public get isFullCircle(): boolean {\r\n    return Angle.isFullCircleRadians(this._radians);\r\n  }\r\n  /** Test if the angle is a half circle (in either direction) */\r\n  public get isHalfCircle(): boolean {\r\n    return Angle.isHalfCircleRadians(this._radians);\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..360 */\r\n  public static adjustDegrees0To360(degrees: number): number {\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      if (degrees < period)\r\n        return degrees;\r\n      const numPeriods = Math.floor(degrees / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      const radians = Angle.adjustDegrees0To360(-degrees);\r\n      return 360.0 - radians;\r\n    }\r\n    // fall through for Nan (disaster) !!!\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is in -180..180 */\r\n  public static adjustDegreesSigned180(degrees: number): number {\r\n    if (Math.abs(degrees) <= 180.0)\r\n      return degrees;\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      const numPeriods = 1 + Math.floor((degrees - 180.0) / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustDegreesSigned180(-degrees);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..2Pi */\r\n  public static adjustRadians0To2Pi(radians: number): number {\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      if (radians < period)\r\n        return radians;\r\n      const numPeriods = Math.floor(radians / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return Math.PI * 2.0 - Angle.adjustRadians0To2Pi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in -PI..PI */\r\n  public static adjustRadiansMinusPiPlusPi(radians: number): number {\r\n    if (Math.abs(radians) <= Math.PI)\r\n      return radians;\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      const numPeriods = 1 + Math.floor((radians - Math.PI) / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustRadiansMinusPiPlusPi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Return a (newly allocated) Angle object with value 0 radians */\r\n  public static zero(): Angle {\r\n    return new Angle(0);\r\n  }\r\n  /** Test if the angle is exactly zero. */\r\n  public get isExactZero(): boolean {\r\n    return this.radians === 0;\r\n  }\r\n  /** Test if the angle is almost zero (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostZero(): boolean {\r\n    return Math.abs(this.radians) < Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is almost a north or south pole (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostNorthOrSouthPole(): boolean {\r\n    return Angle.isHalfCircleRadians(this.radians * 2.0);\r\n  }\r\n  /** Create an angle object with degrees adjusted into 0..360. */\r\n  public static createDegreesAdjustPositive(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegrees0To360(degrees));\r\n  }\r\n  /** Create an angle object with degrees adjusted into -180..180. */\r\n  public static createDegreesAdjustSigned180(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegreesSigned180(degrees));\r\n  }\r\n  /**\r\n   * Test if two radians values are equivalent, allowing shift by full circle (i.e. by a multiple of `2*PI`)\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansAllowPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    const delta = Math.abs(radiansA - radiansB);\r\n    if (delta <= radianTol)\r\n      return true;\r\n    const period = Math.PI * 2.0;\r\n    if (Math.abs(delta - period) <= radianTol)\r\n      return true;\r\n    const numPeriod = Math.round(delta / period);\r\n    const delta1 = delta - numPeriod * period;\r\n    return Math.abs(delta1) <= radianTol;\r\n  }\r\n  /**\r\n   * Test if this angle and other are equivalent, allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualAllowPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansAllowPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two angle (in radians)  almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansNoPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Math.abs(radiansA - radiansB) < radianTol;\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualNoPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansNoPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * * This function is same as isAlmostEqualRadiansNoPeriodShift. Please use isAlmostEqualRadiansNoPeriodShift.\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqual(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return this.isAlmostEqualNoPeriodShift(other, radianTol);\r\n  }\r\n  /**\r\n   * Test if dot product values indicate non-zero length perpendicular vectors.\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static isPerpendicularDotSet(dotUU: number, dotVV: number, dotUV: number): boolean {\r\n    return dotUU > Geometry.smallMetricDistanceSquared\r\n      && dotVV > Geometry.smallMetricDistanceSquared\r\n      && dotUV * dotUV <= Geometry.smallAngleRadiansSquared * dotUU * dotVV;\r\n  }\r\n  /**\r\n   * Return cosine, sine, and radians for the half angle of a \"cosine,sine\" pair.\r\n   * * This function assumes the input arguments are related to an angle between -PI and PI\r\n   * * This function returns an angle between -PI and PI\r\n   * @param rCos2A cosine value (scaled by radius) for initial angle.\r\n   * @param rSin2A sine value (scaled by radius) for final angle.\r\n   */\r\n  public static trigValuesToHalfAngleTrigValues(rCos2A: number, rSin2A: number): TrigValues {\r\n    const r = Geometry.hypotenuseXY(rCos2A, rSin2A);\r\n    if (r < Geometry.smallMetricDistance) {\r\n      return { c: 1.0, s: 0.0, radians: 0.0 }; // angle = 0\r\n    } else {\r\n      /* If the caller really gave you sine and cosine values, r should be 1.  However,\r\n       * to allow scaled values -- e.g. the x and y components of any vector -- we normalize\r\n       * right here. This adds an extra sqrt and 2 divides to the whole process, but improves\r\n       * both the usefulness and robustness of the computation.\r\n       */\r\n      let cosA;\r\n      let sinA = 0.0;\r\n      const cos2A = rCos2A / r;\r\n      const sin2A = rSin2A / r;\r\n      // Original angle in NE and SE quadrants. Half angle in same quadrant\r\n      if (cos2A >= 0.0) {\r\n        /*\r\n         * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n         * so 1 + cos2A = 2(cosA)^2 and therefore, cosA = sqrt((1+cos2A)/2)\r\n         * cosine is positive in NE and SE quadrants so we use +sqrt\r\n         */\r\n        cosA = Math.sqrt(0.5 * (1.0 + cos2A));\r\n        // We know sin2A = 2 sinA cosA so sinA = sin2A/(2*cosA)\r\n        sinA = sin2A / (2.0 * cosA);\r\n      } else {\r\n        // Original angle in NW quadrant. Half angle in NE quadrant\r\n        if (sin2A > 0.0) {\r\n          /*\r\n           * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n           * so 1 - cos2A = 2(sinA)^2 and therefore, sinA = sqrt((1-cos2A)/2)\r\n           * sine is positive in NE quadrant so we use +sqrt\r\n           */\r\n          sinA = Math.sqrt(0.5 * (1.0 - cos2A));\r\n          // Original angle in SW quadrant. Half angle in SE quadrant\r\n        } else {\r\n          // sine is negative in SE quadrant so we use -sqrt\r\n          sinA = -Math.sqrt(0.5 * (1.0 - cos2A));\r\n        }\r\n        // We know sin2A = 2 sinA cosA so cosA = sin2A/(2*sinA)\r\n        cosA = sin2A / (2.0 * sinA); // always positive\r\n      }\r\n      return { c: cosA, s: sinA, radians: Math.atan2(sinA, cosA) };\r\n    }\r\n  }\r\n  /** If value is close to -1, -0.5, 0, 0.5, 1, adjust it to the exact value. */\r\n  public static cleanupTrigValue(value: number, tolerance: number = 1.0e-15): number {\r\n    const absValue = Math.abs(value);\r\n    if (absValue <= tolerance)\r\n      return 0;\r\n    let a = Math.abs(absValue - 0.5);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -0.5 : 0.5;\r\n    a = Math.abs(absValue - 1.0);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -1.0 : 1.0;\r\n    return value;\r\n  }\r\n  /**\r\n   * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define\r\n   * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The\r\n   * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.\r\n   * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.\r\n   * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is\r\n   * perpendicular to the ellipse.\r\n   * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)\r\n   * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.\r\n   * math details can be found at docs/learning/geometry/Angle.md\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static dotProductsToHalfAngleTrigValues(\r\n    dotUU: number, dotVV: number, dotUV: number, favorZero: boolean = true\r\n  ): TrigValues {\r\n\r\n    const cos2t0 = dotUU - dotVV;\r\n    const sin2t0 = 2.0 * dotUV;\r\n    if (favorZero && Math.abs(sin2t0) < Geometry.smallAngleRadians * (Math.abs(dotUU) + Math.abs(dotVV)))\r\n      return { c: 1.0, s: 0.0, radians: 0.0 };\r\n    return Angle.trigValuesToHalfAngleTrigValues(cos2t0, sin2t0);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components\r\n   * * The returned angle is between 0 and PI\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   */\r\n  public static radiansBetweenVectorsXYZ(\r\n    ux: number, uy: number, uz: number, vx: number, vy: number, vz: number\r\n  ): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    return Math.atan2(Geometry.crossProductMagnitude(ux, uy, uz, vx, vy, vz), uDotV);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve\r\n   * angle to a full 2PI range.\r\n   * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)\r\n   * * The angle is in the plane of the U and V vectors.\r\n   * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.\r\n   *\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   * @param upVectorX x component of vector to positive side of plane.\r\n   * @param upVectorY y component of vector to positive side of plane.\r\n   * @param upVectorZ z component of vector to positive side of plane.\r\n   * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI). If false, return\r\n   * signed (-PI < radians <= PI)\r\n   */\r\n  public static orientedRadiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number,\r\n    upVectorX: number, upVectorY: number, upVectorZ: number,\r\n    adjustToPositive: boolean = false): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    const wx = uy * vz - uz * vy;\r\n    const wy = uz * vx - ux * vz;\r\n    const wz = ux * vy - uy * vx;\r\n    const upDotW = upVectorX * wx + upVectorY * wy + upVectorZ * wz;\r\n    const crossMagnitude = Geometry.hypotenuseXYZ(wx, wy, wz);\r\n    if (upDotW < 0.0) {\r\n      if (adjustToPositive) {\r\n        // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as\r\n        // addition to PI. This gives the smoothest numerical transition passing PI.\r\n        return Math.PI + Math.atan2(crossMagnitude, -uDotV);\r\n      } else {\r\n        return -Math.atan2(crossMagnitude, uDotV);\r\n      }\r\n    } else {\r\n      return Math.atan2(crossMagnitude, uDotV);\r\n    }\r\n  }\r\n  /**\r\n   * Add a multiple of a full circle angle (360 degrees, 2PI) in place.\r\n   * @param multiple multiplier factor\r\n   */\r\n  public addMultipleOf2PiInPlace(multiple: number) {\r\n    if (this._degrees !== undefined) {\r\n      this._degrees += multiple * 360.0;\r\n      this._radians = Angle.degreesToRadians(this._degrees);\r\n    } else {\r\n      this._radians += multiple * Angle.pi2Radians;\r\n    }\r\n  }\r\n}\r\n"]}