@math.gl/geospatial 5.0.0-alpha.9 → 5.0.0-beta.2

package/dist/elastic.d.ts

@@ -0,0 +1,5 @@
+/** Optional, published Elastic datasets. Not imported by the main map API. */
+export { elasticI } from "./discontinuous/elastic/elasticI.js";
+export { elasticII } from "./discontinuous/elastic/elasticII.js";
+export { elasticIII } from "./discontinuous/elastic/elasticIII.js";
+//# sourceMappingURL=elastic.d.ts.map

package/dist/elastic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"elastic.d.ts","sourceRoot":"","sources":["../src/elastic.ts"],"names":[],"mappings":"AAGA,8EAA8E;AAC9E,OAAO,EAAC,QAAQ,EAAC,4CAAyC;AAC1D,OAAO,EAAC,SAAS,EAAC,6CAA0C;AAC5D,OAAO,EAAC,UAAU,EAAC,8CAA2C"}

package/dist/elastic.js

@@ -0,0 +1,8 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Optional, published Elastic datasets. Not imported by the main map API. */
+export { elasticI } from "./discontinuous/elastic/elasticI.js";
+export { elasticII } from "./discontinuous/elastic/elasticII.js";
+export { elasticIII } from "./discontinuous/elastic/elasticIII.js";
+//# sourceMappingURL=elastic.js.map

package/dist/elastic.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"elastic.js","sourceRoot":"","sources":["../src/elastic.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,8EAA8E;AAC9E,OAAO,EAAC,QAAQ,EAAC,4CAAyC;AAC1D,OAAO,EAAC,SAAS,EAAC,6CAA0C;AAC5D,OAAO,EAAC,UAAU,EAAC,8CAA2C"}

package/dist/ellipsoid-geodesic.d.ts

@@ -0,0 +1,63 @@
+/** Two-axis revolution geometry. For compatibility, semiMajorAxis means the
+ * equatorial radius a and semiMinorAxis the polar radius b, also when b > a.
+ */
+export type GeodesicParameters = {
+    readonly semiMajorAxis: number;
+    readonly semiMinorAxis: number;
+};
+export type GeodesicOptions = {
+    /** Preserve longitude winding rather than normalizing to [-180, 180]. */
+    unrollLongitude?: boolean;
+    /** Include the signed geodesic edge area S12 in squared axis units. */
+    area?: boolean;
+};
+/** Angles and azimuths are degrees; distances and reduced lengths use the axes' units. */
+export type GeodesicResult = {
+    lat1: number;
+    lon1: number;
+    lat2: number;
+    lon2: number;
+    azi1: number;
+    azi2: number;
+    s12: number;
+    m12: number;
+    M12: number;
+    M21: number;
+    a12: number;
+    S12?: number;
+};
+interface EngineLine {
+    GenPosition(arc: boolean, distance: number, mask: number, output?: GeodesicResult): GeodesicResult;
+    Position(distance: number, mask: number): GeodesicResult;
+    ArcPosition(arc: number, mask: number): GeodesicResult;
+}
+interface Engine {
+    _c2: number;
+    Inverse(lat1: number, lon1: number, lat2: number, lon2: number, mask: number): GeodesicResult;
+    Line(latitude: number, longitude: number, azimuth: number, caps: number): EngineLine;
+}
+/** A geodesic with a fixed starting point and initial azimuth. */
+export interface GeodesicLine {
+    position(distance: number): GeodesicResult;
+    arcPosition(arc: number): GeodesicResult;
+}
+/** Geodesics on an oblate or prolate ellipsoid of revolution, or sphere, using Karney's series solver.
+ * This class snapshots the axes. It does not accept triaxial ellipsoids or infer a CRS.
+ */
+export declare class EllipsoidGeodesic {
+    private readonly engine;
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }: GeodesicParameters);
+    /** Total area in squared axis units, with the same normalization as S12. */
+    get surfaceArea(): number;
+    inverse(lat1: number, lon1: number, lat2: number, lon2: number, options?: GeodesicOptions): GeodesicResult;
+    direct(latitude: number, longitude: number, azimuth: number, distance: number, options?: GeodesicOptions): GeodesicResult;
+    arcDirect(latitude: number, longitude: number, azimuth: number, arc: number, options?: GeodesicOptions): GeodesicResult;
+    line(latitude: number, longitude: number, azimuth: number, options?: GeodesicOptions): GeodesicLine;
+}
+/** @internal Shared engine and non-area mask for projections. Keeping the raw
+ * kernel separate lets projection bundles discard the public optional facade.
+ */
+export declare const nonAreaGeodesicMask: number;
+export declare function createGeodesicEngine({ semiMajorAxis: a, semiMinorAxis: b }: GeodesicParameters): Engine;
+export {};
+//# sourceMappingURL=ellipsoid-geodesic.d.ts.map

package/dist/ellipsoid-geodesic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-geodesic.d.ts","sourceRoot":"","sources":["../src/ellipsoid-geodesic.ts"],"names":[],"mappings":"AAMA;;GAEG;AACH,MAAM,MAAM,kBAAkB,GAAG;IAAC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAA;CAAC,CAAC;AAElG,MAAM,MAAM,eAAe,GAAG;IAC5B,yEAAyE;IACzE,eAAe,CAAC,EAAE,OAAO,CAAC;IAC1B,uEAAuE;IACvE,IAAI,CAAC,EAAE,OAAO,CAAC;CAChB,CAAC;AAEF,0FAA0F;AAC1F,MAAM,MAAM,cAAc,GAAG;IAC3B,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,CAAC,EAAE,MAAM,CAAC;CACd,CAAC;AAEF,UAAU,UAAU;IAClB,WAAW,CACT,GAAG,EAAE,OAAO,EACZ,QAAQ,EAAE,MAAM,EAChB,IAAI,EAAE,MAAM,EACZ,MAAM,CAAC,EAAE,cAAc,GACtB,cAAc,CAAC;IAClB,QAAQ,CAAC,QAAQ,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,cAAc,CAAC;IACzD,WAAW,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,cAAc,CAAC;CACxD;AACD,UAAU,MAAM;IACd,GAAG,EAAE,MAAM,CAAC;IACZ,OAAO,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,cAAc,CAAC;IAC9F,IAAI,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,UAAU,CAAC;CACtF;AAUD,kEAAkE;AAClE,MAAM,WAAW,YAAY;IAC3B,QAAQ,CAAC,QAAQ,EAAE,MAAM,GAAG,cAAc,CAAC;IAC3C,WAAW,CAAC,GAAG,EAAE,MAAM,GAAG,cAAc,CAAC;CAC1C;AAED;;GAEG;AACH,qBAAa,iBAAiB;IAC5B,OAAO,CAAC,QAAQ,CAAC,MAAM,CAAS;gBAEpB,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAC,EAAE,kBAAkB;IAIpE,4EAA4E;IAC5E,IAAI,WAAW,IAAI,MAAM,CAExB;IAED,OAAO,CACL,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,OAAO,GAAE,eAAoB,GAC5B,cAAc;IAMjB,MAAM,CACJ,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,EAAE,MAAM,EACf,QAAQ,EAAE,MAAM,EAChB,OAAO,GAAE,eAAoB,GAC5B,cAAc;IAIjB,SAAS,CACP,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,EAAE,MAAM,EACf,GAAG,EAAE,MAAM,EACX,OAAO,GAAE,eAAoB,GAC5B,cAAc;IAIjB,IAAI,CACF,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,eAAoB,GAC5B,YAAY;CAsBhB;AAmBD;;GAEG;AACH,eAAO,MAAM,mBAAmB,QAA8B,CAAC;AAC/D,wBAAgB,oBAAoB,CAAC,EACnC,aAAa,EAAE,CAAC,EAChB,aAAa,EAAE,CAAC,EACjB,EAAE,kBAAkB,GAAG,MAAM,CAK7B"}

package/dist/ellipsoid-geodesic.js

@@ -0,0 +1,70 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original typed facade over the separately attributed GeographicLib geodesic solver. See ../GEOGRAPHICLIB-LICENSE.txt.
+import { geodesic } from "./vendor/geographiclib-geodesic/runtime.js";
+const library = geodesic.Geodesic;
+/** Geodesics on an oblate or prolate ellipsoid of revolution, or sphere, using Karney's series solver.
+ * This class snapshots the axes. It does not accept triaxial ellipsoids or infer a CRS.
+ */
+export class EllipsoidGeodesic {
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }) {
+        this.engine = createGeodesicEngine({ semiMajorAxis: a, semiMinorAxis: b });
+    }
+    /** Total area in squared axis units, with the same normalization as S12. */
+    get surfaceArea() {
+        return 4 * Math.PI * this.engine._c2;
+    }
+    inverse(lat1, lon1, lat2, lon2, options = {}) {
+        validatePoint(lat1, lon1);
+        validatePoint(lat2, lon2);
+        return this.engine.Inverse(lat1, lon1, lat2, lon2, outputMask(options));
+    }
+    direct(latitude, longitude, azimuth, distance, options = {}) {
+        return this.line(latitude, longitude, azimuth, options).position(distance);
+    }
+    arcDirect(latitude, longitude, azimuth, arc, options = {}) {
+        return this.line(latitude, longitude, azimuth, options).arcPosition(arc);
+    }
+    line(latitude, longitude, azimuth, options = {}) {
+        validatePoint(latitude, longitude);
+        validateFinite(azimuth);
+        const selectedMask = outputMask(options);
+        const line = this.engine.Line(latitude, longitude, azimuth, selectedMask | library.DISTANCE_IN);
+        return {
+            // Fourth argument is an internal result request; ordinary Array.map arguments remain ignored.
+            position(distance, _index, _array, output) {
+                validateFinite(distance);
+                return line.GenPosition(false, distance, selectedMask, output);
+            },
+            arcPosition(arc) {
+                validateFinite(arc);
+                return line.ArcPosition(arc, selectedMask);
+            }
+        };
+    }
+}
+function validatePoint(latitude, longitude) {
+    if (!(Number.isFinite(latitude) && Math.abs(latitude) <= 90 && Number.isFinite(longitude))) {
+        throw new Error('Invalid geodesic coordinates');
+    }
+}
+function outputMask(options) {
+    return ((library.ALL & ~(options.area ? 0 : library.AREA)) |
+        (options.unrollLongitude ? library.LONG_UNROLL : 0));
+}
+function validateFinite(value) {
+    if (!Number.isFinite(value))
+        throw new Error('Nonfinite axis, azimuth, distance or arc');
+}
+/** @internal Shared engine and non-area mask for projections. Keeping the raw
+ * kernel separate lets projection bundles discard the public optional facade.
+ */
+export const nonAreaGeodesicMask = library.ALL & ~library.AREA;
+export function createGeodesicEngine({ semiMajorAxis: a, semiMinorAxis: b }) {
+    validateFinite(a);
+    validateFinite(b);
+    // The engine validates positivity of both derived axes.
+    return new library.Geodesic(a, (a - b) / a);
+}
+//# sourceMappingURL=ellipsoid-geodesic.js.map

package/dist/ellipsoid-geodesic.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-geodesic.js","sourceRoot":"","sources":["../src/ellipsoid-geodesic.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,0IAA0I;AAC1I,OAAO,EAAC,QAAQ,EAAC,mDAAgD;AAoDjE,MAAM,OAAO,GAAG,QAAQ,CAAC,QAAgC,CAAC;AAQ1D;;GAEG;AACH,MAAM,OAAO,iBAAiB;IAG5B,YAAY,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAqB;QAClE,IAAI,CAAC,MAAM,GAAG,oBAAoB,CAAC,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAC,CAAC,CAAC;IAC3E,CAAC;IAED,4EAA4E;IAC5E,IAAI,WAAW;QACb,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC;IACvC,CAAC;IAED,OAAO,CACL,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,UAA2B,EAAE;QAE7B,aAAa,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC1B,aAAa,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC1B,OAAO,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,UAAU,CAAC,OAAO,CAAC,CAAC,CAAC;IAC1E,CAAC;IAED,MAAM,CACJ,QAAgB,EAChB,SAAiB,EACjB,OAAe,EACf,QAAgB,EAChB,UAA2B,EAAE;QAE7B,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC7E,CAAC;IAED,SAAS,CACP,QAAgB,EAChB,SAAiB,EACjB,OAAe,EACf,GAAW,EACX,UAA2B,EAAE;QAE7B,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC;IAC3E,CAAC;IAED,IAAI,CACF,QAAgB,EAChB,SAAiB,EACjB,OAAe,EACf,UAA2B,EAAE;QAE7B,aAAa,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACnC,cAAc,CAAC,OAAO,CAAC,CAAC;QACxB,MAAM,YAAY,GAAG,UAAU,CAAC,OAAO,CAAC,CAAC;QACzC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,YAAY,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC;QAChG,OAAO;YACL,8FAA8F;YAC9F,QAAQ,CACN,QAAgB,EAChB,MAAgB,EAChB,MAAgB,EAChB,MAAuB;gBAEvB,cAAc,CAAC,QAAQ,CAAC,CAAC;gBACzB,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,QAAQ,EAAE,YAAY,EAAE,MAAM,CAAC,CAAC;YACjE,CAAC;YACD,WAAW,CAAC,GAAW;gBACrB,cAAc,CAAC,GAAG,CAAC,CAAC;gBACpB,OAAO,IAAI,CAAC,WAAW,CAAC,GAAG,EAAE,YAAY,CAAC,CAAC;YAC7C,CAAC;SACF,CAAC;IACJ,CAAC;CACF;AAED,SAAS,aAAa,CAAC,QAAgB,EAAE,SAAiB;IACxD,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,EAAE,IAAI,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC;QAC3F,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,CAAC;AACH,CAAC;AAED,SAAS,UAAU,CAAC,OAAwB;IAC1C,OAAO,CACL,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;QAClD,CAAC,OAAO,CAAC,eAAe,CAAC,CAAC,CAAC,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CACpD,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CAAC,KAAa;IACnC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;AAC3F,CAAC;AAED;;GAEG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAG,OAAO,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC;AAC/D,MAAM,UAAU,oBAAoB,CAAC,EACnC,aAAa,EAAE,CAAC,EAChB,aAAa,EAAE,CAAC,EACG;IACnB,cAAc,CAAC,CAAC,CAAC,CAAC;IAClB,cAAc,CAAC,CAAC,CAAC,CAAC;IAClB,wDAAwD;IACxD,OAAO,IAAI,OAAO,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AAC9C,CAAC"}

package/dist/ellipsoid-helpers/ellipsoid-transform.d.ts

@@ -1,4 +1,4 @@
-import { NumericArray } from '@math.gl/types';
+import { NumericArray } from '@math.gl/core';
 import type { Ellipsoid } from "../ellipsoid.js";
 export type AxisDirection = 'up' | 'down' | 'north' | 'east' | 'south' | 'west';
 export declare function localFrameToFixedFrame(ellipsoid: Ellipsoid, firstAxis: AxisDirection, secondAxis: AxisDirection, thirdAxis: AxisDirection, cartesianOrigin: Readonly<NumericArray>, result: number[]): number[];

package/dist/ellipsoid-helpers/ellipsoid-transform.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ellipsoid-transform.d.ts","sourceRoot":"","sources":["../../src/ellipsoid-helpers/ellipsoid-transform.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,YAAY,EAAC,MAAM,gBAAgB,CAAC;AAG5C,OAAO,KAAK,EAAC,SAAS,EAAC,wBAAqB;AAM5C,MAAM,MAAM,aAAa,GAAG,IAAI,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,CAAC;AAsEhF,wBAAgB,sBAAsB,CACpC,SAAS,EAAE,SAAS,EACpB,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,eAAe,EAAE,QAAQ,CAAC,YAAY,CAAC,EACvC,MAAM,EAAE,MAAM,EAAE,GACf,MAAM,EAAE,CAuEV"}
+{"version":3,"file":"ellipsoid-transform.d.ts","sourceRoot":"","sources":["../../src/ellipsoid-helpers/ellipsoid-transform.ts"],"names":[],"mappings":"AAMA,OAAO,EAAC,YAAY,EAAC,MAAM,eAAe,CAAC;AAQ3C,OAAO,KAAK,EAAC,SAAS,EAAC,wBAAqB;AAO5C,MAAM,MAAM,aAAa,GAAG,IAAI,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,GAAG,OAAO,GAAG,MAAM,CAAC;AAsEhF,wBAAgB,sBAAsB,CACpC,SAAS,EAAE,SAAS,EACpB,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,eAAe,EAAE,QAAQ,CAAC,YAAY,CAAC,EACvC,MAAM,EAAE,MAAM,EAAE,GACf,MAAM,EAAE,CA0HV"}

package/dist/ellipsoid-helpers/ellipsoid-transform.js

@@ -1,9 +1,13 @@
 // math.gl
-// SPDX-License-Identifier: MIT and Apache-2.0
-// Copyright (c) vis.gl contributors
-import { Vector3, assert, equals as equalsEpsilon } from '@math.gl/core';
+// SPDX-License-Identifier: MIT AND Apache-2.0
+// SPDX-FileCopyrightText: Copyright 2011-2018 CesiumJS Contributors
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Derived from Cesium. See the repository LICENSE for upstream attribution and Apache-2.0 terms. Scaled east normalization and the small-unit pole neighborhood are original math.gl fixes.
+import { Vector3, assert } from '@math.gl/core';
+import { createLocalFrameBasis, eastNorthUpBasisFromDirections, localFrameToMatrix } from '@math.gl/core';
 const EPSILON14 = 1e-14;
 const scratchOrigin = new Vector3();
+const scratchBasis = createLocalFrameBasis();
 // Caclulate third axis from given two axii
 const VECTOR_PRODUCT_LOCAL_FRAME = {
     up: {
@@ -74,7 +78,38 @@ export function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAx
     let thirdAxisVector;
     const origin = scratchOrigin.copy(cartesianOrigin);
     // If x and y are zero, assume origin is at a pole, which is a special case.
-    const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);
+    if (!Number.isFinite(origin.x) ||
+        !Number.isFinite(origin.y) ||
+        !Number.isFinite(origin.z) ||
+        (!Array.isArray(result) && result.length < 16))
+        throw new Error('Unsupported local frame');
+    // Retain the legacy orientation neighborhood for ordinary units, but never
+    // classify an entire tiny ellipsoid as a pole merely because of its unit.
+    const poleTolerance = EPSILON14 * Math.min(1, ellipsoid.maximumRadius);
+    const atPole = Math.abs(origin.x) <= poleTolerance && Math.abs(origin.y) <= poleTolerance;
+    // Only the qualified sphere/oblate basis is shared. Three-radius and center
+    // conventions keep the retained Cesium path, without changing its normal choice.
+    const radii = ellipsoid.radii;
+    if (radii.x === radii.y &&
+        radii.x > 0 &&
+        radii.z > 0 &&
+        radii.z <= radii.x &&
+        (!atPole || origin.z !== 0)) {
+        const { east, up } = scratchAxisVectors;
+        if (atPole) {
+            east.set(0, 1, 0);
+            up.set(0, 0, Math.sign(origin.z));
+        }
+        else {
+            const scale = Math.max(Math.abs(origin.x), Math.abs(origin.y));
+            east.set(-origin.y / scale, origin.x / scale, 0).normalize();
+            ellipsoid.geodeticSurfaceNormal(origin, up);
+        }
+        if (!eastNorthUpBasisFromDirections(east, up, scratchBasis) ||
+            !localFrameToMatrix(scratchBasis, origin, firstAxis, secondAxis, thirdAxis, result))
+            throw new Error('Unsupported local frame');
+        return result;
+    }
     if (atPole) {
         // Look up axis value and adjust
         const sign = Math.sign(origin.z);
@@ -94,7 +129,8 @@ export function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAx
     else {
         // Calculate all axis
         const { up, east, north } = scratchAxisVectors;
-        east.set(-origin.y, origin.x, 0.0).normalize();
+        const scale = Math.max(Math.abs(origin.x), Math.abs(origin.y));
+        east.set(-origin.y / scale, origin.x / scale, 0).normalize();
         ellipsoid.geodeticSurfaceNormal(origin, up);
         north.copy(up).cross(east);
         const { down, west, south } = scratchAxisVectors;
@@ -106,22 +142,27 @@ export function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAx
         secondAxisVector = scratchAxisVectors[secondAxis];
         thirdAxisVector = scratchAxisVectors[thirdAxis];
     }
-    // TODO - assuming the result is column-major
-    result[0] = firstAxisVector.x;
-    result[1] = firstAxisVector.y;
-    result[2] = firstAxisVector.z;
+    // Snapshot the retained three-radius/center path too: application result setters
+    // can recursively invoke either path and overwrite all module scratch vectors.
+    const ax = firstAxisVector.x, ay = firstAxisVector.y, az = firstAxisVector.z;
+    const bx = secondAxisVector.x, by = secondAxisVector.y, bz = secondAxisVector.z;
+    const cx = thirdAxisVector.x, cy = thirdAxisVector.y, cz = thirdAxisVector.z;
+    const x = origin.x, y = origin.y, z = origin.z;
+    result[0] = ax;
+    result[1] = ay;
+    result[2] = az;
     result[3] = 0.0;
-    result[4] = secondAxisVector.x;
-    result[5] = secondAxisVector.y;
-    result[6] = secondAxisVector.z;
+    result[4] = bx;
+    result[5] = by;
+    result[6] = bz;
     result[7] = 0.0;
-    result[8] = thirdAxisVector.x;
-    result[9] = thirdAxisVector.y;
-    result[10] = thirdAxisVector.z;
+    result[8] = cx;
+    result[9] = cy;
+    result[10] = cz;
     result[11] = 0.0;
-    result[12] = origin.x;
-    result[13] = origin.y;
-    result[14] = origin.z;
+    result[12] = x;
+    result[13] = y;
+    result[14] = z;
     result[15] = 1.0;
     return result;
 }

package/dist/ellipsoid-helpers/ellipsoid-transform.js.map

@@ -1 +1 @@
-{"version":3,"file":"ellipsoid-transform.js","sourceRoot":"","sources":["../../src/ellipsoid-helpers/ellipsoid-transform.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oCAAoC;AAGpC,OAAO,EAAC,OAAO,EAAE,MAAM,EAAE,MAAM,IAAI,aAAa,EAAC,MAAM,eAAe,CAAC;AAIvE,MAAM,SAAS,GAAG,KAAK,CAAC;AAExB,MAAM,aAAa,GAAG,IAAI,OAAO,EAAE,CAAC;AAIpC,2CAA2C;AAC3C,MAAM,0BAA0B,GAG5B;IACF,EAAE,EAAE;QACF,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,MAAM;QACb,IAAI,EAAE,OAAO;QACb,IAAI,EAAE,OAAO;KACd;IACD,IAAI,EAAE;QACJ,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,MAAM;QACb,IAAI,EAAE,OAAO;QACb,IAAI,EAAE,OAAO;KACd;IACD,KAAK,EAAE;QACL,EAAE,EAAE,MAAM;QACV,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,IAAI;KACX;IACD,KAAK,EAAE;QACL,EAAE,EAAE,MAAM;QACV,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,IAAI;QACV,IAAI,EAAE,MAAM;KACb;IACD,IAAI,EAAE;QACJ,EAAE,EAAE,OAAO;QACX,IAAI,EAAE,OAAO;QACb,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,IAAI;KACZ;IACD,IAAI,EAAE;QACJ,EAAE,EAAE,OAAO;QACX,IAAI,EAAE,OAAO;QACb,KAAK,EAAE,IAAI;QACX,KAAK,EAAE,MAAM;KACd;CACO,CAAC;AAEX,MAAM,4BAA4B,GAAG;IACnC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACb,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;CACR,CAAC;AAEX,MAAM,kBAAkB,GAAG;IACzB,IAAI,EAAE,IAAI,OAAO,EAAE;IACnB,KAAK,EAAE,IAAI,OAAO,EAAE;IACpB,EAAE,EAAE,IAAI,OAAO,EAAE;IACjB,IAAI,EAAE,IAAI,OAAO,EAAE;IACnB,KAAK,EAAE,IAAI,OAAO,EAAE;IACpB,IAAI,EAAE,IAAI,OAAO,EAAE;CACpB,CAAC;AAEF,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AACrC,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AACrC,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AAErC,8DAA8D;AAC9D,qFAAqF;AACrF,kEAAkE;AAClE,MAAM,UAAU,sBAAsB,CACpC,SAAoB,EACpB,SAAwB,EACxB,UAAyB,EACzB,SAAwB,EACxB,eAAuC,EACvC,MAAgB;IAEhB,MAAM,iBAAiB,GACrB,0BAA0B,CAAC,SAAS,CAAC,IAAI,0BAA0B,CAAC,SAAS,CAAC,CAAC,UAAU,CAAC,CAAC;IAC7F,4EAA4E;IAC5E,MAAM,CAAC,iBAAiB,IAAI,CAAC,CAAC,SAAS,IAAI,SAAS,KAAK,iBAAiB,CAAC,CAAC,CAAC;IAE7E,IAAI,eAAwB,CAAC;IAC7B,IAAI,gBAAyB,CAAC;IAC9B,IAAI,eAAwB,CAAC;IAE7B,MAAM,MAAM,GAAG,aAAa,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC;IAEnD,4EAA4E;IAC5E,MAAM,MAAM,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,SAAS,CAAC,IAAI,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,SAAS,CAAC,CAAC;IAElG,IAAI,MAAM,EAAE,CAAC;QACX,gCAAgC;QAChC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAEjC,eAAe,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,SAAS,CAAC,CAAC,CAAC;QACpF,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,MAAM,EAAE,CAAC;YACjD,eAAe,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9B,CAAC;QAED,gBAAgB,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,UAAU,CAAC,CAAC,CAAC;QACtF,IAAI,UAAU,KAAK,MAAM,IAAI,UAAU,KAAK,MAAM,EAAE,CAAC;YACnD,gBAAgB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC/B,CAAC;QAED,eAAe,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,SAAS,CAAC,CAAC,CAAC;QACpF,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,MAAM,EAAE,CAAC;YACjD,eAAe,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9B,CAAC;IACH,CAAC;SAAM,CAAC;QACN,qBAAqB;QACrB,MAAM,EAAC,EAAE,EAAE,IAAI,EAAE,KAAK,EAAC,GAAG,kBAAkB,CAAC;QAE7C,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QAC/C,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;QAC5C,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAE3B,MAAM,EAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAC,GAAG,kBAAkB,CAAC;QAE/C,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1B,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAE5B,+CAA+C;QAC/C,eAAe,GAAG,kBAAkB,CAAC,SAAS,CAAC,CAAC;QAChD,gBAAgB,GAAG,kBAAkB,CAAC,UAAU,CAAC,CAAC;QAClD,eAAe,GAAG,kBAAkB,CAAC,SAAS,CAAC,CAAC;IAClD,CAAC;IAED,6CAA6C;IAC7C,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAChB,MAAM,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAChB,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjB,MAAM,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACtB,MAAM,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACtB,MAAM,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACtB,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjB,OAAO,MAAM,CAAC;AAChB,CAAC"}
+{"version":3,"file":"ellipsoid-transform.js","sourceRoot":"","sources":["../../src/ellipsoid-helpers/ellipsoid-transform.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oEAAoE;AACpE,4DAA4D;AAC5D,8MAA8M;AAG9M,OAAO,EAAC,OAAO,EAAE,MAAM,EAAC,MAAM,eAAe,CAAC;AAC9C,OAAO,EACL,qBAAqB,EACrB,8BAA8B,EAC9B,kBAAkB,EACnB,MAAM,eAAe,CAAC;AAIvB,MAAM,SAAS,GAAG,KAAK,CAAC;AAExB,MAAM,aAAa,GAAG,IAAI,OAAO,EAAE,CAAC;AACpC,MAAM,YAAY,GAAG,qBAAqB,EAAE,CAAC;AAI7C,2CAA2C;AAC3C,MAAM,0BAA0B,GAG5B;IACF,EAAE,EAAE;QACF,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,MAAM;QACb,IAAI,EAAE,OAAO;QACb,IAAI,EAAE,OAAO;KACd;IACD,IAAI,EAAE;QACJ,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,MAAM;QACb,IAAI,EAAE,OAAO;QACb,IAAI,EAAE,OAAO;KACd;IACD,KAAK,EAAE;QACL,EAAE,EAAE,MAAM;QACV,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,IAAI;KACX;IACD,KAAK,EAAE;QACL,EAAE,EAAE,MAAM;QACV,IAAI,EAAE,MAAM;QACZ,IAAI,EAAE,IAAI;QACV,IAAI,EAAE,MAAM;KACb;IACD,IAAI,EAAE;QACJ,EAAE,EAAE,OAAO;QACX,IAAI,EAAE,OAAO;QACb,KAAK,EAAE,MAAM;QACb,KAAK,EAAE,IAAI;KACZ;IACD,IAAI,EAAE;QACJ,EAAE,EAAE,OAAO;QACX,IAAI,EAAE,OAAO;QACb,KAAK,EAAE,IAAI;QACX,KAAK,EAAE,MAAM;KACd;CACO,CAAC;AAEX,MAAM,4BAA4B,GAAG;IACnC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACf,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACb,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;CACR,CAAC;AAEX,MAAM,kBAAkB,GAAG;IACzB,IAAI,EAAE,IAAI,OAAO,EAAE;IACnB,KAAK,EAAE,IAAI,OAAO,EAAE;IACpB,EAAE,EAAE,IAAI,OAAO,EAAE;IACjB,IAAI,EAAE,IAAI,OAAO,EAAE;IACnB,KAAK,EAAE,IAAI,OAAO,EAAE;IACpB,IAAI,EAAE,IAAI,OAAO,EAAE;CACpB,CAAC;AAEF,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AACrC,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AACrC,MAAM,cAAc,GAAG,IAAI,OAAO,EAAE,CAAC;AAErC,8DAA8D;AAC9D,qFAAqF;AACrF,kEAAkE;AAClE,MAAM,UAAU,sBAAsB,CACpC,SAAoB,EACpB,SAAwB,EACxB,UAAyB,EACzB,SAAwB,EACxB,eAAuC,EACvC,MAAgB;IAEhB,MAAM,iBAAiB,GACrB,0BAA0B,CAAC,SAAS,CAAC,IAAI,0BAA0B,CAAC,SAAS,CAAC,CAAC,UAAU,CAAC,CAAC;IAC7F,4EAA4E;IAC5E,MAAM,CAAC,iBAAiB,IAAI,CAAC,CAAC,SAAS,IAAI,SAAS,KAAK,iBAAiB,CAAC,CAAC,CAAC;IAE7E,IAAI,eAAwB,CAAC;IAC7B,IAAI,gBAAyB,CAAC;IAC9B,IAAI,eAAwB,CAAC;IAE7B,MAAM,MAAM,GAAG,aAAa,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC;IAEnD,4EAA4E;IAC5E,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;QAC1B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;QAC1B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;QAC1B,CAAC,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,IAAK,MAAuB,CAAC,MAAM,GAAG,EAAE,CAAC;QAEhE,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;IAC7C,2EAA2E;IAC3E,0EAA0E;IAC1E,MAAM,aAAa,GAAG,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,CAAC,aAAa,CAAC,CAAC;IACvE,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,aAAa,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,aAAa,CAAC;IAE1F,4EAA4E;IAC5E,iFAAiF;IACjF,MAAM,KAAK,GAAG,SAAS,CAAC,KAAK,CAAC;IAC9B,IACE,KAAK,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC;QACnB,KAAK,CAAC,CAAC,GAAG,CAAC;QACX,KAAK,CAAC,CAAC,GAAG,CAAC;QACX,KAAK,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC;QAClB,CAAC,CAAC,MAAM,IAAI,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,EAC3B,CAAC;QACD,MAAM,EAAC,IAAI,EAAE,EAAE,EAAC,GAAG,kBAAkB,CAAC;QACtC,IAAI,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAClB,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACpC,CAAC;aAAM,CAAC;YACN,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/D,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,EAAE,MAAM,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;YAC7D,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;QAC9C,CAAC;QACD,IACE,CAAC,8BAA8B,CAAC,IAAI,EAAE,EAAE,EAAE,YAAY,CAAC;YACvD,CAAC,kBAAkB,CAAC,YAAY,EAAE,MAAM,EAAE,SAAS,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,CAAC;YAEnF,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;QAC7C,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,IAAI,MAAM,EAAE,CAAC;QACX,gCAAgC;QAChC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAEjC,eAAe,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,SAAS,CAAC,CAAC,CAAC;QACpF,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,MAAM,EAAE,CAAC;YACjD,eAAe,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9B,CAAC;QAED,gBAAgB,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,UAAU,CAAC,CAAC,CAAC;QACtF,IAAI,UAAU,KAAK,MAAM,IAAI,UAAU,KAAK,MAAM,EAAE,CAAC;YACnD,gBAAgB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC/B,CAAC;QAED,eAAe,GAAG,cAAc,CAAC,SAAS,CAAC,4BAA4B,CAAC,SAAS,CAAC,CAAC,CAAC;QACpF,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,MAAM,EAAE,CAAC;YACjD,eAAe,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9B,CAAC;IACH,CAAC;SAAM,CAAC;QACN,qBAAqB;QACrB,MAAM,EAAC,EAAE,EAAE,IAAI,EAAE,KAAK,EAAC,GAAG,kBAAkB,CAAC;QAE7C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/D,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,EAAE,MAAM,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QAC7D,SAAS,CAAC,qBAAqB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;QAC5C,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAE3B,MAAM,EAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAC,GAAG,kBAAkB,CAAC;QAE/C,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1B,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAE5B,+CAA+C;QAC/C,eAAe,GAAG,kBAAkB,CAAC,SAAS,CAAC,CAAC;QAChD,gBAAgB,GAAG,kBAAkB,CAAC,UAAU,CAAC,CAAC;QAClD,eAAe,GAAG,kBAAkB,CAAC,SAAS,CAAC,CAAC;IAClD,CAAC;IAED,iFAAiF;IACjF,+EAA+E;IAC/E,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,EAC1B,EAAE,GAAG,eAAe,CAAC,CAAC,EACtB,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAC3B,EAAE,GAAG,gBAAgB,CAAC,CAAC,EACvB,EAAE,GAAG,gBAAgB,CAAC,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,EAC1B,EAAE,GAAG,eAAe,CAAC,CAAC,EACtB,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC;IACzB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,EAChB,CAAC,GAAG,MAAM,CAAC,CAAC,EACZ,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAChB,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAChB,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACf,MAAM,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IAChB,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjB,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACf,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACf,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACf,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;IACjB,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/dist/ellipsoid-helpers/exact-nearest-surface.d.ts

@@ -0,0 +1,16 @@
+/** Rare fallback: exact signs of the multiplier equation, not arbitrary precision
+ * transcendental arithmetic. BigInt/string allocations are confined to this path.
+ * Input dyadics share a unit; t uses at least 64 fraction bits. For a nonzero
+ * shortest-axis coordinate the root is at least one squared input unit, so this
+ * resolves its relative bracket to substantially better than binary64 precision.
+ */
+export declare function exactNearestSurface(x: number, y: number, z: number, a: number, b: number, c: number, result: {
+    x: number;
+    y: number;
+    z: number;
+}): boolean;
+/** Round a positive rational dyadic to a binary64 significand without converting
+ * unbounded integers to Number first. Split the exponent for subnormal outputs.
+ */
+export declare function ratio(n: bigint, d: bigint, exponent: number): number;
+//# sourceMappingURL=exact-nearest-surface.d.ts.map

package/dist/ellipsoid-helpers/exact-nearest-surface.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"exact-nearest-surface.d.ts","sourceRoot":"","sources":["../../src/ellipsoid-helpers/exact-nearest-surface.ts"],"names":[],"mappings":"AAKA;;;;;GAKG;AACH,wBAAgB,mBAAmB,CACjC,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,MAAM,EAAE;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAC,GACxC,OAAO,CAkFT;AAED;;GAEG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM,CAWpE"}

package/dist/ellipsoid-helpers/exact-nearest-surface.js

@@ -0,0 +1,108 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original exact-dyadic fallback for ill-conditioned nearest-point geometry; no upstream code.
+/** Rare fallback: exact signs of the multiplier equation, not arbitrary precision
+ * transcendental arithmetic. BigInt/string allocations are confined to this path.
+ * Input dyadics share a unit; t uses at least 64 fraction bits. For a nonzero
+ * shortest-axis coordinate the root is at least one squared input unit, so this
+ * resolves its relative bracket to substantially better than binary64 precision.
+ */
+export function exactNearestSurface(x, y, z, a, b, c, result) {
+    const values = [a, b, c, Math.abs(x), Math.abs(y), Math.abs(z)];
+    const parts = values.map(value => {
+        const bits = value.toString(2);
+        const dot = bits.indexOf('.');
+        return {
+            integer: BigInt('0b' + bits.replace('.', '')),
+            exponent: dot < 0 ? 0 : dot + 1 - bits.length
+        };
+    });
+    const exponent = Math.min(...parts.filter(p => p.integer !== 0n).map(p => p.exponent));
+    const integers = parts.map(p => p.integer << BigInt(p.exponent - exponent));
+    const squares = integers.slice(0, 3).map(v => v * v);
+    const minimum = squares.reduce((p, v) => (p < v ? p : v));
+    const offsets = squares.map(v => v - minimum);
+    const products = squares.map((_, i) => integers[i] * integers[i + 3]);
+    let precision = 64;
+    if (!products.some((u, i) => offsets[i] === 0n && u !== 0n)) {
+        // Exact boundary classification. Negative means multiple minima; equality
+        // means the unique cusp point has zero shortest-axis coordinates.
+        let numerator = -1n, denominator = 1n;
+        for (let i = 0; i < 3; i++) {
+            if (offsets[i] === 0n)
+                continue;
+            const d2 = offsets[i] * offsets[i];
+            numerator = numerator * d2 + products[i] * products[i] * denominator;
+            denominator *= d2;
+        }
+        if (numerator < 0n)
+            return false;
+        if (numerator === 0n) {
+            return commit(0n, 1n);
+        }
+        // Convexity bounds the positive root below by F(0)/abs(F'(0)).
+        // abs(F'(0)) <= 2*(1+F(0))/minPositiveOffset. Use this exact
+        // rational bound when a zero shortest-axis coordinate admits a tiny root.
+        const positive = offsets.filter(d => d > 0n);
+        const least = positive.reduce((p, v) => (p < v ? p : v));
+        precision = Math.max(64, (2n * (denominator + numerator)).toString(2).length -
+            (numerator * least).toString(2).length +
+            128);
+    }
+    const unit = 1n << BigInt(precision);
+    const d = offsets.map(v => v * unit);
+    const u = products.map(v => v * unit);
+    const u2 = u.map(v => v * v);
+    let lower = 0n, upper = u.reduce((sum, v) => sum + v, 0n);
+    if (upper === 0n)
+        return false;
+    for (let iteration = 0; iteration < 4096; iteration++) {
+        const t = (lower + upper) >> 1n;
+        if (t === lower)
+            return false;
+        const dx = d[0] + t, dy = d[1] + t, dz = d[2] + t;
+        const xx = dx * dx, yy = dy * dy, zz = dz * dz;
+        const sign = u2[0] * yy * zz + u2[1] * xx * zz + u2[2] * xx * yy - xx * yy * zz;
+        if (sign === 0n)
+            return commit(t, unit);
+        if (sign > 0n)
+            lower = t;
+        else
+            upper = t;
+        if (lower > 0n && upper - lower <= lower >> 60n)
+            return commit((lower + upper) >> 1n, unit);
+    }
+    return false;
+    function commit(t, unit) {
+        const output = squares.map((r2, i) => {
+            const denominator = offsets[i] * unit + t;
+            return denominator === 0n ? 0 : ratio(r2 * integers[i + 3] * unit, denominator, exponent);
+        });
+        if (!output.every(Number.isFinite))
+            return false;
+        const qx = Math.sign(x) * output[0], qy = Math.sign(y) * output[1], qz = Math.sign(z) * output[2];
+        result.x = qx;
+        result.y = qy;
+        result.z = qz;
+        return true;
+    }
+}
+/** Round a positive rational dyadic to a binary64 significand without converting
+ * unbounded integers to Number first. Split the exponent for subnormal outputs.
+ */
+export function ratio(n, d, exponent) {
+    if (n === 0n)
+        return 0;
+    const shift = 52 - (n.toString(2).length - d.toString(2).length);
+    const numerator = shift >= 0 ? n << BigInt(shift) : n;
+    const denominator = shift >= 0 ? d : d << BigInt(-shift);
+    let quotient = numerator / denominator;
+    const remainder = numerator % denominator;
+    if (2n * remainder > denominator || (2n * remainder === denominator && quotient & 1n))
+        quotient++;
+    const power = exponent - shift;
+    const first = Math.max(-1022, Math.min(1023, power));
+    return Number(quotient) * 2 ** first * 2 ** (power - first);
+}
+//# sourceMappingURL=exact-nearest-surface.js.map

package/dist/ellipsoid-helpers/exact-nearest-surface.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"exact-nearest-surface.js","sourceRoot":"","sources":["../../src/ellipsoid-helpers/exact-nearest-surface.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,iHAAiH;AAEjH;;;;;GAKG;AACH,MAAM,UAAU,mBAAmB,CACjC,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAyC;IAEzC,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE;QAC/B,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;QAC9B,OAAO;YACL,OAAO,EAAE,MAAM,CAAC,IAAI,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;YAC7C,QAAQ,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM;SAC9C,CAAC;IACJ,CAAC,CAAC,CAAC;IACH,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC;IACvF,MAAM,QAAQ,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,IAAI,MAAM,CAAC,CAAC,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC;IAC5E,MAAM,OAAO,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrD,MAAM,OAAO,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1D,MAAM,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;IAC9C,MAAM,QAAQ,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACtE,IAAI,SAAS,GAAG,EAAE,CAAC;IACnB,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC;QAC5D,0EAA0E;QAC1E,kEAAkE;QAClE,IAAI,SAAS,GAAG,CAAC,EAAE,EACjB,WAAW,GAAG,EAAE,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,EAAE;gBAAE,SAAS;YAChC,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;YACnC,SAAS,GAAG,SAAS,GAAG,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC;YACrE,WAAW,IAAI,EAAE,CAAC;QACpB,CAAC;QACD,IAAI,SAAS,GAAG,EAAE;YAAE,OAAO,KAAK,CAAC;QACjC,IAAI,SAAS,KAAK,EAAE,EAAE,CAAC;YACrB,OAAO,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACxB,CAAC;QACD,+DAA+D;QAC/D,6DAA6D;QAC7D,0EAA0E;QAC1E,MAAM,QAAQ,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QAC7C,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzD,SAAS,GAAG,IAAI,CAAC,GAAG,CAClB,EAAE,EACF,CAAC,EAAE,GAAG,CAAC,WAAW,GAAG,SAAS,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM;YACjD,CAAC,SAAS,GAAG,KAAK,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM;YACtC,GAAG,CACN,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,GAAG,EAAE,IAAI,MAAM,CAAC,SAAS,CAAC,CAAC;IACrC,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IACrC,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7B,IAAI,KAAK,GAAG,EAAE,EACZ,KAAK,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC;IAC5C,IAAI,KAAK,KAAK,EAAE;QAAE,OAAO,KAAK,CAAC;IAC/B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI,EAAE,SAAS,EAAE,EAAE,CAAC;QACtD,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC;QAChC,IAAI,CAAC,KAAK,KAAK;YAAE,OAAO,KAAK,CAAC;QAC9B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACjB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EACb,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,EAChB,EAAE,GAAG,EAAE,GAAG,EAAE,EACZ,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACf,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAChF,IAAI,IAAI,KAAK,EAAE;YAAE,OAAO,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;QACxC,IAAI,IAAI,GAAG,EAAE;YAAE,KAAK,GAAG,CAAC,CAAC;;YACpB,KAAK,GAAG,CAAC,CAAC;QACf,IAAI,KAAK,GAAG,EAAE,IAAI,KAAK,GAAG,KAAK,IAAI,KAAK,IAAI,GAAG;YAAE,OAAO,MAAM,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,EAAE,EAAE,IAAI,CAAC,CAAC;IAC9F,CAAC;IACD,OAAO,KAAK,CAAC;IAEb,SAAS,MAAM,CAAC,CAAS,EAAE,IAAY;QACrC,MAAM,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;YACnC,MAAM,WAAW,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YAC1C,OAAO,WAAW,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,EAAE,WAAW,EAAE,QAAQ,CAAC,CAAC;QAC5F,CAAC,CAAC,CAAC;QACH,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QACjD,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EACjC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAC7B,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;QACd,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;QACd,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;QACd,OAAO,IAAI,CAAC;IACd,CAAC;AACH,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,KAAK,CAAC,CAAS,EAAE,CAAS,EAAE,QAAgB;IAC1D,IAAI,CAAC,KAAK,EAAE;QAAE,OAAO,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACjE,MAAM,SAAS,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtD,MAAM,WAAW,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC;IACzD,IAAI,QAAQ,GAAG,SAAS,GAAG,WAAW,CAAC;IACvC,MAAM,SAAS,GAAG,SAAS,GAAG,WAAW,CAAC;IAC1C,IAAI,EAAE,GAAG,SAAS,GAAG,WAAW,IAAI,CAAC,EAAE,GAAG,SAAS,KAAK,WAAW,IAAI,QAAQ,GAAG,EAAE,CAAC;QAAE,QAAQ,EAAE,CAAC;IAClG,MAAM,KAAK,GAAG,QAAQ,GAAG,KAAK,CAAC;IAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC;IACrD,OAAO,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,KAAK,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;AAC9D,CAAC"}

package/dist/ellipsoid-helpers/intersection-bounds.d.ts

@@ -0,0 +1,18 @@
+import type { EllipsoidGeodesic, GeodesicLine } from "../ellipsoid-geodesic.js";
+/** @internal Native characteristic distances, in normalized axis units. */
+export declare class IntersectionBounds {
+    private readonly solver;
+    readonly margin: number;
+    readonly capture: number;
+    readonly spacing: number;
+    readonly halfCircumference: number;
+    constructor(solver: EllipsoidGeodesic, ratio: number, radius: number);
+    conjugate(line: GeodesicLine, guess: number, semi?: boolean, reference?: {
+        m12: number;
+        M12: number;
+        M21: number;
+    }): number;
+    private polarDistance;
+    private polarMinimum;
+}
+//# sourceMappingURL=intersection-bounds.d.ts.map

package/dist/ellipsoid-helpers/intersection-bounds.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"intersection-bounds.d.ts","sourceRoot":"","sources":["../../src/ellipsoid-helpers/intersection-bounds.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,iBAAiB,EAAE,YAAY,EAAC,iCAA8B;AAE3E,2EAA2E;AAC3E,qBAAa,kBAAkB;IAM3B,OAAO,CAAC,QAAQ,CAAC,MAAM;IALzB,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB,QAAQ,CAAC,iBAAiB,EAAE,MAAM,CAAC;gBAEhB,MAAM,EAAE,iBAAiB,EAC1C,KAAK,EAAE,MAAM,EACb,MAAM,EAAE,MAAM;IAUhB,SAAS,CACP,IAAI,EAAE,YAAY,EAClB,KAAK,EAAE,MAAM,EACb,IAAI,UAAQ,EACZ,SAAS;;;;KAA2B,GACnC,MAAM;IAiBT,OAAO,CAAC,aAAa;IAOrB,OAAO,CAAC,YAAY;CA6BrB"}

package/dist/ellipsoid-helpers/intersection-bounds.js

@@ -0,0 +1,60 @@
+/** @internal Native characteristic distances, in normalized axis units. */
+export class IntersectionBounds {
+    constructor(solver, ratio, radius) {
+        this.solver = solver;
+        this.halfCircumference = Math.PI * radius;
+        this.margin = this.halfCircumference * Number.EPSILON ** (1 / 5);
+        this.capture = ratio <= 1 ? Math.PI * ratio : this.polarMinimum(ratio);
+        this.spacing = this.capture - this.margin;
+        if (!(Number.isFinite(this.spacing) && this.spacing > 0))
+            throw new Error('Unresolved geodesic capture bound');
+    }
+    conjugate(line, guess, semi = false, reference = { m12: 0, M12: 1, M21: 1 }) {
+        const { m12, M12, M21 } = reference;
+        let s = guess;
+        for (let i = 0; i < 100; i++) {
+            const p = line.position(s);
+            const m = p.m12 * M12 - m12 * p.M12;
+            const forward = p.M12 * M21 + (m12 === 0 ? 0 : ((1 - M12 * M21) * p.m12) / m12);
+            const reverse = p.M21 * M12 + (p.m12 === 0 ? 0 : ((1 - p.M12 * p.M21) * m12) / p.m12);
+            const step = semi ? (m * forward) / (1 - forward * reverse) : -m / reverse;
+            if (!Number.isFinite(step))
+                throw new Error('Unresolved conjugate distance');
+            const next = s + step;
+            if (next === s || Math.abs(step) <= 5e-14)
+                return next;
+            s = next;
+        }
+        throw new Error('Conjugate distance did not converge');
+    }
+    polarDistance(latitude, ratio) {
+        return this.conjugate(this.solver.line(latitude, 0, 0), ((1 + (1 - ratio) / 2) * Math.PI) / 2, true);
+    }
+    polarMinimum(ratio) {
+        let lat0 = 63, lat1 = 65, lat2 = 64;
+        let s0 = this.polarDistance(lat0, ratio), s1 = this.polarDistance(lat1, ratio), s2 = this.polarDistance(lat2, ratio);
+        let best = s2;
+        for (let i = 0; i < 10; i++) {
+            const denominator = (lat1 - lat0) * s2 + (lat0 - lat2) * s1 + (lat2 - lat1) * s0;
+            if (denominator === 0)
+                break;
+            const next = ((lat1 - lat0) * (lat1 + lat0) * s2 +
+                (lat0 - lat2) * (lat0 + lat2) * s1 +
+                (lat2 - lat1) * (lat2 + lat1) * s0) /
+                (2 * denominator);
+            // Near a sphere the quadratic becomes numerically flat. The capture
+            // safety margin is much larger than the roundoff in these distances.
+            if (!(Number.isFinite(next) && next >= 0 && next <= 90))
+                break;
+            lat0 = lat1;
+            s0 = s1;
+            lat1 = lat2;
+            s1 = s2;
+            lat2 = next;
+            s2 = this.polarDistance(lat2, ratio);
+            best = Math.min(best, s2);
+        }
+        return 2 * best;
+    }
+}
+//# sourceMappingURL=intersection-bounds.js.map

package/dist/ellipsoid-helpers/intersection-bounds.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"intersection-bounds.js","sourceRoot":"","sources":["../../src/ellipsoid-helpers/intersection-bounds.ts"],"names":[],"mappings":"AAOA,2EAA2E;AAC3E,MAAM,OAAO,kBAAkB;IAK7B,YACmB,MAAyB,EAC1C,KAAa,EACb,MAAc;QAFG,WAAM,GAAN,MAAM,CAAmB;QAI1C,IAAI,CAAC,iBAAiB,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC;QAC1C,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,iBAAiB,GAAG,MAAM,CAAC,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACjE,IAAI,CAAC,OAAO,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC;QACvE,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,MAAM,CAAC;QAC1C,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,OAAO,GAAG,CAAC,CAAC;YACtD,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;IACzD,CAAC;IAED,SAAS,CACP,IAAkB,EAClB,KAAa,EACb,IAAI,GAAG,KAAK,EACZ,SAAS,GAAG,EAAC,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAC;QAEpC,MAAM,EAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAC,GAAG,SAAS,CAAC;QAClC,IAAI,CAAC,GAAG,KAAK,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC;YACpC,MAAM,OAAO,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;YAChF,MAAM,OAAO,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACtF,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YAC3E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;gBAAE,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;YAC7E,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;YACtB,IAAI,IAAI,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK;gBAAE,OAAO,IAAI,CAAC;YACvD,CAAC,GAAG,IAAI,CAAC;QACX,CAAC;QACD,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IACzD,CAAC;IAEO,aAAa,CAAC,QAAgB,EAAE,KAAa;QACnD,OAAO,IAAI,CAAC,SAAS,CACnB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC,EAChC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EACrC,IAAI,CACL,CAAC;IACJ,CAAC;IACO,YAAY,CAAC,KAAa;QAChC,IAAI,IAAI,GAAG,EAAE,EACX,IAAI,GAAG,EAAE,EACT,IAAI,GAAG,EAAE,CAAC;QACZ,IAAI,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,IAAI,EAAE,KAAK,CAAC,EACtC,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,IAAI,EAAE,KAAK,CAAC,EACpC,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;QACvC,IAAI,IAAI,GAAG,EAAE,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,WAAW,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;YACjF,IAAI,WAAW,KAAK,CAAC;gBAAE,MAAM;YAC7B,MAAM,IAAI,GACR,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE;gBACjC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE;gBAClC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;gBACrC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;YACpB,oEAAoE;YACpE,qEAAqE;YACrE,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,IAAI,IAAI,IAAI,EAAE,CAAC;gBAAE,MAAM;YAC/D,IAAI,GAAG,IAAI,CAAC;YACZ,EAAE,GAAG,EAAE,CAAC;YACR,IAAI,GAAG,IAAI,CAAC;YACZ,EAAE,GAAG,EAAE,CAAC;YACR,IAAI,GAAG,IAAI,CAAC;YACZ,EAAE,GAAG,IAAI,CAAC,aAAa,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;YACrC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC;QAC5B,CAAC;QACD,OAAO,CAAC,GAAG,IAAI,CAAC;IAClB,CAAC;CACF"}

package/dist/ellipsoid-helpers/nearest-surface.d.ts

@@ -0,0 +1,12 @@
+/** Unique nearest point on a positive three-axis ellipsoid. No output on failure.
+ * With t = lambda + min(radii)^2 > 0 the constrained Hessian is positive
+ * definite, certifying a global minimum in exact arithmetic. Unique boundary
+ * roots and ill-conditioned residuals use exact dyadic classification/solving. Numerical qualification
+ * is sampled; unrepresentable products and exhausted convergence reject.
+ */
+export declare function nearestSurface(x: number, y: number, z: number, a: number, b: number, c: number, result: {
+    x: number;
+    y: number;
+    z: number;
+}): boolean;
+//# sourceMappingURL=nearest-surface.d.ts.map

package/dist/ellipsoid-helpers/nearest-surface.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"nearest-surface.d.ts","sourceRoot":"","sources":["../../src/ellipsoid-helpers/nearest-surface.ts"],"names":[],"mappings":"AAOA;;;;;GAKG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,MAAM,EAAE;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAC,GACxC,OAAO,CAyGT"}