@math.gl/geospatial 5.0.0-alpha.19 → 5.0.0-alpha.20

package/dist/rhumb-path.d.ts

@@ -0,0 +1,28 @@
+import type { NumericArray } from '@math.gl/core/types';
+import { EllipsoidRhumb } from "./ellipsoid-rhumb.js";
+export type RhumbPathOptions = {
+    rhumb?: EllipsoidRhumb;
+    /** Preserve the supplied endpoint interval, including complete turns. */
+    longitudeMode?: 'shortest' | 'unrolled';
+    /** Whole turns added to the shortest endpoint longitude difference. */
+    longitudeTurns?: number;
+};
+/** Constant-bearing endpoint path, in degrees and solver axis units. Height is
+ * linear in surface-distance fraction and does not contribute to length.
+ * Unlike geodesics, rhumbs reject exact poles. Endpoints/options are snapshots.
+ */
+export declare class RhumbPath {
+    readonly length: number;
+    readonly size: 2 | 3;
+    private readonly first;
+    private readonly last;
+    private readonly line;
+    constructor(start: ArrayLike<number>, end: ArrayLike<number>, options?: RhumbPathOptions);
+    get start(): number[];
+    get end(): number[];
+    sample(fraction: number): number[];
+    sample<T extends NumericArray>(fraction: number, result: T): T;
+    sampleByDistance(distance: number): number[];
+    sampleByDistance<T extends NumericArray>(distance: number, result: T): T;
+}
+//# sourceMappingURL=rhumb-path.d.ts.map

package/dist/rhumb-path.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"rhumb-path.d.ts","sourceRoot":"","sources":["../src/rhumb-path.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,qBAAqB,CAAC;AACtD,OAAO,EAAC,cAAc,EAAiB,6BAA0B;AAWjE,MAAM,MAAM,gBAAgB,GAAG;IAC7B,KAAK,CAAC,EAAE,cAAc,CAAC;IACvB,yEAAyE;IACzE,aAAa,CAAC,EAAE,UAAU,GAAG,UAAU,CAAC;IACxC,uEAAuE;IACvE,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB,CAAC;AAEF;;;GAGG;AACH,qBAAa,SAAS;IACpB,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC;IACrB,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAoB;IAC1C,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAoB;IACzC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAY;gBAErB,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,OAAO,GAAE,gBAAqB;IA0C5F,IAAI,KAAK,IAAI,MAAM,EAAE,CAEpB;IACD,IAAI,GAAG,IAAI,MAAM,EAAE,CAElB;IAED,MAAM,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM,EAAE;IAClC,MAAM,CAAC,CAAC,SAAS,YAAY,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC;IAkC9D,gBAAgB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM,EAAE;IAC5C,gBAAgB,CAAC,CAAC,SAAS,YAAY,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,GAAG,CAAC;CAMzE"}

package/dist/rhumb-path.js

@@ -0,0 +1,100 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original distance-parameterized adapter over the attributed math.gl rhumb solver.
+import { EllipsoidRhumb } from "./ellipsoid-rhumb.js";
+import { validateGeographicPosition } from "./geodesic-path.js";
+import { WGS84_RADIUS_X, WGS84_RADIUS_Z } from "./constants.js";
+let wgs84;
+function getDefaultRhumb() {
+    if (!wgs84)
+        wgs84 = new EllipsoidRhumb({ semiMajorAxis: WGS84_RADIUS_X, semiMinorAxis: WGS84_RADIUS_Z });
+    return wgs84;
+}
+/** Constant-bearing endpoint path, in degrees and solver axis units. Height is
+ * linear in surface-distance fraction and does not contribute to length.
+ * Unlike geodesics, rhumbs reject exact poles. Endpoints/options are snapshots.
+ */
+export class RhumbPath {
+    constructor(start, end, options = {}) {
+        validateGeographicPosition(start);
+        validateGeographicPosition(end);
+        if (start.length !== end.length)
+            throw new RangeError('Rhumb endpoint dimensions must match');
+        const solver = options.rhumb || getDefaultRhumb();
+        if ((options.longitudeMode !== undefined &&
+            !['shortest', 'unrolled'].includes(options.longitudeMode)) ||
+            (options.longitudeMode === 'unrolled' && options.longitudeTurns !== undefined))
+            throw new RangeError('Choose explicit turns or unrolled endpoint longitudes');
+        let inverse = solver.inverse(start[1], start[0], end[1], end[0], {
+            longitudeTurns: options.longitudeTurns,
+            unrollLongitude: true
+        });
+        if (options.longitudeMode === 'unrolled') {
+            const turns = Math.round((end[0] - inverse.lon2) / 360);
+            if (turns !== 0)
+                inverse = solver.inverse(start[1], start[0], end[1], end[0], {
+                    longitudeTurns: turns,
+                    unrollLongitude: true
+                });
+        }
+        if (![inverse.s12, inverse.azi12, inverse.lon2].every(Number.isFinite) || inverse.s12 < 0)
+            throw new RangeError('Unresolved rhumb path');
+        this.size = start.length;
+        this.length = inverse.s12;
+        this.first = Object.freeze(Array.from(start));
+        const last = Array.from(end);
+        const turns = Math.round((inverse.lon2 - end[0]) / 360);
+        const endpoint = end[0] + 360 * turns;
+        const retainsPhase = Number.isSafeInteger(turns) &&
+            Number.isFinite(endpoint) &&
+            Math.abs(endpoint - inverse.lon2) <=
+                4 * Number.EPSILON * Math.max(360, Math.abs(inverse.lon2));
+        last[0] = inverse.s12 === 0 ? start[0] : retainsPhase ? endpoint : inverse.lon2;
+        if (inverse.s12 === 0)
+            last[1] = start[1];
+        this.last = Object.freeze(last);
+        this.line = solver.line(start[1], start[0], inverse.azi12, { unrollLongitude: true });
+    }
+    get start() {
+        return this.first.slice();
+    }
+    get end() {
+        return this.last.slice();
+    }
+    sample(fraction, result = new Array(this.size)) {
+        if (!Number.isFinite(fraction) || fraction < 0 || fraction > 1)
+            throw new RangeError('Rhumb fraction must be in [0, 1]');
+        if (result.length !== this.size ||
+            (ArrayBuffer.isView(result) &&
+                !(result instanceof Float32Array) &&
+                !(result instanceof Float64Array)))
+            throw new RangeError('Expected a floating rhumb output with matching dimensions');
+        let longitude = this.first[0], latitude = this.first[1];
+        if (fraction === 1) {
+            longitude = this.last[0];
+            latitude = this.last[1];
+        }
+        else if (fraction !== 0 && this.length !== 0) {
+            const p = this.line.position(this.length * fraction);
+            longitude = p.lon2;
+            latitude = p.lat2;
+        }
+        const height = this.size === 3 ? this.first[2] * (1 - fraction) + this.last[2] * fraction : 0;
+        const values = [longitude, latitude, height];
+        if (!values.every(Number.isFinite) ||
+            (result instanceof Float32Array && !values.every(x => Number.isFinite(Math.fround(x)))))
+            throw new RangeError('Unrepresentable rhumb sample');
+        result[0] = longitude;
+        result[1] = latitude;
+        if (this.size === 3)
+            result[2] = height;
+        return result;
+    }
+    sampleByDistance(distance, result = new Array(this.size)) {
+        if (!Number.isFinite(distance) || distance < 0 || distance > this.length)
+            throw new RangeError('Distance must be within the rhumb path');
+        return this.sample(this.length === 0 ? 0 : distance / this.length, result);
+    }
+}
+//# sourceMappingURL=rhumb-path.js.map

package/dist/rhumb-path.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rhumb-path.js","sourceRoot":"","sources":["../src/rhumb-path.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,sGAAsG;AAGtG,OAAO,EAAC,cAAc,EAAiB,6BAA0B;AACjE,OAAO,EAAC,0BAA0B,EAAC,2BAAwB;AAC3D,OAAO,EAAC,cAAc,EAAE,cAAc,EAAC,uBAAoB;AAE3D,IAAI,KAAiC,CAAC;AACtC,SAAS,eAAe;IACtB,IAAI,CAAC,KAAK;QACR,KAAK,GAAG,IAAI,cAAc,CAAC,EAAC,aAAa,EAAE,cAAc,EAAE,aAAa,EAAE,cAAc,EAAC,CAAC,CAAC;IAC7F,OAAO,KAAK,CAAC;AACf,CAAC;AAUD;;;GAGG;AACH,MAAM,OAAO,SAAS;IAOpB,YAAY,KAAwB,EAAE,GAAsB,EAAE,UAA4B,EAAE;QAC1F,0BAA0B,CAAC,KAAK,CAAC,CAAC;QAClC,0BAA0B,CAAC,GAAG,CAAC,CAAC;QAChC,IAAI,KAAK,CAAC,MAAM,KAAK,GAAG,CAAC,MAAM;YAAE,MAAM,IAAI,UAAU,CAAC,sCAAsC,CAAC,CAAC;QAC9F,MAAM,MAAM,GAAG,OAAO,CAAC,KAAK,IAAI,eAAe,EAAE,CAAC;QAClD,IACE,CAAC,OAAO,CAAC,aAAa,KAAK,SAAS;YAClC,CAAC,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,aAAa,CAAC,CAAC;YAC5D,CAAC,OAAO,CAAC,aAAa,KAAK,UAAU,IAAI,OAAO,CAAC,cAAc,KAAK,SAAS,CAAC;YAE9E,MAAM,IAAI,UAAU,CAAC,uDAAuD,CAAC,CAAC;QAChF,IAAI,OAAO,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE;YAC/D,cAAc,EAAE,OAAO,CAAC,cAAc;YACtC,eAAe,EAAE,IAAI;SACtB,CAAC,CAAC;QACH,IAAI,OAAO,CAAC,aAAa,KAAK,UAAU,EAAE,CAAC;YACzC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;YACxD,IAAI,KAAK,KAAK,CAAC;gBACb,OAAO,GAAG,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE;oBAC3D,cAAc,EAAE,KAAK;oBACrB,eAAe,EAAE,IAAI;iBACtB,CAAC,CAAC;QACP,CAAC;QACD,IAAI,CAAC,CAAC,OAAO,CAAC,GAAG,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,OAAO,CAAC,GAAG,GAAG,CAAC;YACvF,MAAM,IAAI,UAAU,CAAC,uBAAuB,CAAC,CAAC;QAChD,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC,MAAe,CAAC;QAClC,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC;QAC1B,IAAI,CAAC,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC;QACtC,MAAM,YAAY,GAChB,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;YAC3B,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YACzB,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC;gBAC/B,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;QAC/D,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,OAAO,CAAC,IAAI,CAAC;QAChF,IAAI,OAAO,CAAC,GAAG,KAAK,CAAC;YAAE,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QAC1C,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,EAAC,eAAe,EAAE,IAAI,EAAC,CAAC,CAAC;IACtF,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC;IAC5B,CAAC;IACD,IAAI,GAAG;QACL,OAAO,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;IAC3B,CAAC;IAID,MAAM,CAAC,QAAgB,EAAE,SAAuB,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC;QAClE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,GAAG,CAAC,IAAI,QAAQ,GAAG,CAAC;YAC5D,MAAM,IAAI,UAAU,CAAC,kCAAkC,CAAC,CAAC;QAC3D,IACE,MAAM,CAAC,MAAM,KAAK,IAAI,CAAC,IAAI;YAC3B,CAAC,WAAW,CAAC,MAAM,CAAC,MAAM,CAAC;gBACzB,CAAC,CAAC,MAAM,YAAY,YAAY,CAAC;gBACjC,CAAC,CAAC,MAAM,YAAY,YAAY,CAAC,CAAC;YAEpC,MAAM,IAAI,UAAU,CAAC,2DAA2D,CAAC,CAAC;QACpF,IAAI,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAC3B,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;YACnB,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACzB,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC1B,CAAC;aAAM,IAAI,QAAQ,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC;YACrD,SAAS,GAAG,CAAC,CAAC,IAAI,CAAC;YACnB,QAAQ,GAAG,CAAC,CAAC,IAAI,CAAC;QACpB,CAAC;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9F,MAAM,MAAM,GAAG,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;QAC7C,IACE,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAC9B,CAAC,MAAM,YAAY,YAAY,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEvF,MAAM,IAAI,UAAU,CAAC,8BAA8B,CAAC,CAAC;QACvD,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;QACtB,MAAM,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;QACrB,IAAI,IAAI,CAAC,IAAI,KAAK,CAAC;YAAE,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACxC,OAAO,MAAM,CAAC;IAChB,CAAC;IAID,gBAAgB,CAAC,QAAgB,EAAE,SAAuB,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC;QAC5E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,GAAG,CAAC,IAAI,QAAQ,GAAG,IAAI,CAAC,MAAM;YACtE,MAAM,IAAI,UAAU,CAAC,wCAAwC,CAAC,CAAC;QACjE,OAAO,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7E,CAAC;CACF"}

package/dist/rhumb-polygon.d.ts

@@ -0,0 +1,11 @@
+import type { GeodesicParameters } from "./ellipsoid-geodesic.js";
+import { SurfacePolygon, type SurfacePointOptions } from "./ellipsoid-helpers/surface-polygon.js";
+export type { SurfacePolygonOptions as RhumbPolygonOptions, SurfacePolygonResult as RhumbPolygonResult, SurfacePointOptions as RhumbPointOptions } from "./ellipsoid-helpers/surface-polygon.js";
+/** Streaming rhumb polygon/polyline. Extra inverse turns are explicit per point. */
+export declare class RhumbPolygon extends SurfacePolygon {
+    constructor(parameters: GeodesicParameters, options?: {
+        polyline?: boolean;
+    });
+    addPoint(latitude: number, longitude: number, options?: SurfacePointOptions): this;
+}
+//# sourceMappingURL=rhumb-polygon.d.ts.map

package/dist/rhumb-polygon.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"rhumb-polygon.d.ts","sourceRoot":"","sources":["../src/rhumb-polygon.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,kBAAkB,EAAC,gCAA6B;AAC7D,OAAO,EAAC,cAAc,EAAE,KAAK,mBAAmB,EAAC,+CAA4C;AAC7F,YAAY,EACV,qBAAqB,IAAI,mBAAmB,EAC5C,oBAAoB,IAAI,kBAAkB,EAC1C,mBAAmB,IAAI,iBAAiB,EACzC,+CAA4C;AAC7C,oFAAoF;AACpF,qBAAa,YAAa,SAAQ,cAAc;gBAClC,UAAU,EAAE,kBAAkB,EAAE,OAAO,GAAE;QAAC,QAAQ,CAAC,EAAE,OAAO,CAAA;KAAM;IAGrE,QAAQ,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,OAAO,GAAE,mBAAwB,GAAG,IAAI;CAMhG"}

package/dist/rhumb-polygon.js

@@ -0,0 +1,19 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original rhumb polygon facade over the original rhumb integrals and separately attributed shared GeographicLib polygon transit/reduction rules.
+import { EllipsoidRhumb, validateRhumbPoint } from "./ellipsoid-rhumb.js";
+import { SurfacePolygon } from "./ellipsoid-helpers/surface-polygon.js";
+/** Streaming rhumb polygon/polyline. Extra inverse turns are explicit per point. */
+export class RhumbPolygon extends SurfacePolygon {
+    constructor(parameters, options = {}) {
+        super(new EllipsoidRhumb(parameters), options);
+    }
+    addPoint(latitude, longitude, options = {}) {
+        validateRhumbPoint(latitude, longitude);
+        if (!Number.isSafeInteger(options.longitudeTurns ?? 0))
+            throw new RangeError('Rhumb longitudeTurns must be a safe integer');
+        return super.addPoint(latitude, longitude, options);
+    }
+}
+//# sourceMappingURL=rhumb-polygon.js.map

package/dist/rhumb-polygon.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rhumb-polygon.js","sourceRoot":"","sources":["../src/rhumb-polygon.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,oKAAoK;AACpK,OAAO,EAAC,cAAc,EAAE,kBAAkB,EAAC,6BAA0B;AAErE,OAAO,EAAC,cAAc,EAA2B,+CAA4C;AAM7F,oFAAoF;AACpF,MAAM,OAAO,YAAa,SAAQ,cAAc;IAC9C,YAAY,UAA8B,EAAE,UAAgC,EAAE;QAC5E,KAAK,CAAC,IAAI,cAAc,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC,CAAC;IACjD,CAAC;IACQ,QAAQ,CAAC,QAAgB,EAAE,SAAiB,EAAE,UAA+B,EAAE;QACtF,kBAAkB,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACxC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,CAAC,cAAc,IAAI,CAAC,CAAC;YACpD,MAAM,IAAI,UAAU,CAAC,6CAA6C,CAAC,CAAC;QACtE,OAAO,KAAK,CAAC,QAAQ,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,CAAC,CAAC;IACtD,CAAC;CACF"}

package/dist/subdivide-geodesic-polyline.d.ts

@@ -0,0 +1,15 @@
+import type { SubdividePolylineOptions, SubdividedPolyline } from '@math.gl/geometry/polygon';
+import { type GeodesicPathOptions } from "./geodesic-path.js";
+export type SubdivideGeodesicPolylineOptions = Omit<SubdividePolylineOptions, 'createEdgeSampler' | 'maxSegmentLength'> & GeodesicPathOptions & {
+    /** Positive maximum surface distance per output edge, in solver axis units.
+     * Height is excluded. Infinity by default; sampled target error is still checked.
+     */
+    maxSurfaceDistance?: number;
+};
+/** Adaptively sample shortest endpoint geodesics through a supplied continuous transform.
+ * Geographic source angles are degrees; heights interpolate by surface-distance fraction.
+ * Output longitudes are unrolled across successive edges. No seam splitting, clipping,
+ * rhumb semantics or universal error bound is implied. Original-edge provenance is retained.
+ */
+export declare function subdivideGeodesicPolyline(input: ArrayLike<number>, options: SubdivideGeodesicPolylineOptions): SubdividedPolyline;
+//# sourceMappingURL=subdivide-geodesic-polyline.d.ts.map

package/dist/subdivide-geodesic-polyline.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"subdivide-geodesic-polyline.d.ts","sourceRoot":"","sources":["../src/subdivide-geodesic-polyline.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,wBAAwB,EAAE,kBAAkB,EAAC,MAAM,2BAA2B,CAAC;AAC5F,OAAO,EAA2C,KAAK,mBAAmB,EAAC,2BAAwB;AAEnG,MAAM,MAAM,gCAAgC,GAAG,IAAI,CACjD,wBAAwB,EACxB,mBAAmB,GAAG,kBAAkB,CACzC,GACC,mBAAmB,GAAG;IACpB;;OAEG;IACH,kBAAkB,CAAC,EAAE,MAAM,CAAC;CAC7B,CAAC;AAEJ;;;;GAIG;AACH,wBAAgB,yBAAyB,CACvC,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,EACxB,OAAO,EAAE,gCAAgC,GACxC,kBAAkB,CA6CpB"}

package/dist/subdivide-geodesic-polyline.js

@@ -0,0 +1,52 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original geodesic edge adapter over the shared adaptive polyline subdivision engine.
+import { subdividePolyline } from '@math.gl/geometry/polygon';
+import { GeodesicPath, validateGeographicPosition } from "./geodesic-path.js";
+/** Adaptively sample shortest endpoint geodesics through a supplied continuous transform.
+ * Geographic source angles are degrees; heights interpolate by surface-distance fraction.
+ * Output longitudes are unrolled across successive edges. No seam splitting, clipping,
+ * rhumb semantics or universal error bound is implied. Original-edge provenance is retained.
+ */
+export function subdivideGeodesicPolyline(input, options) {
+    const { geodesic, size = 2, maxSurfaceDistance = Infinity, ...subdivision } = options;
+    if ((size !== 2 && size !== 3) ||
+        !Number.isSafeInteger(input.length) ||
+        input.length < 0 ||
+        input.length % size !== 0) {
+        throw new RangeError('Expected packed longitude/latitude[/height] positions');
+    }
+    // Reject size/vertex-limit violations before constructing potentially many line objects.
+    const maxVertices = subdivision.maxVertices ?? 65536;
+    if (!Number.isSafeInteger(maxVertices) ||
+        maxVertices < 1 ||
+        maxVertices > 0xffffffff ||
+        input.length / size > maxVertices) {
+        throw new RangeError('Invalid or exceeded maxVertices');
+    }
+    if (!(maxSurfaceDistance > 0) ||
+        (maxSurfaceDistance !== Infinity && !Number.isFinite(maxSurfaceDistance))) {
+        throw new RangeError('maxSurfaceDistance must be positive');
+    }
+    const positions = Array.from(input);
+    for (let i = 0; i < positions.length; i += size)
+        validateGeographicPosition(positions.slice(i, i + size));
+    const paths = [];
+    for (let i = size; i < positions.length; i += size) {
+        const path = new GeodesicPath(positions.slice(i - size, i), positions.slice(i, i + size), {
+            geodesic
+        });
+        paths.push(path);
+        const end = path.end;
+        for (let component = 0; component < size; component++)
+            positions[i + component] = end[component];
+    }
+    return subdividePolyline(positions, {
+        ...subdivision,
+        size,
+        maxSegmentLength: maxSurfaceDistance,
+        createEdgeSampler: (_start, _end, index) => paths[index]
+    });
+}
+//# sourceMappingURL=subdivide-geodesic-polyline.js.map

package/dist/subdivide-geodesic-polyline.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"subdivide-geodesic-polyline.js","sourceRoot":"","sources":["../src/subdivide-geodesic-polyline.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,yGAAyG;AAEzG,OAAO,EAAC,iBAAiB,EAAC,MAAM,2BAA2B,CAAC;AAE5D,OAAO,EAAC,YAAY,EAAE,0BAA0B,EAA2B,2BAAwB;AAanG;;;;GAIG;AACH,MAAM,UAAU,yBAAyB,CACvC,KAAwB,EACxB,OAAyC;IAEzC,MAAM,EAAC,QAAQ,EAAE,IAAI,GAAG,CAAC,EAAE,kBAAkB,GAAG,QAAQ,EAAE,GAAG,WAAW,EAAC,GAAG,OAAO,CAAC;IACpF,IACE,CAAC,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC;QAC1B,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC;QACnC,KAAK,CAAC,MAAM,GAAG,CAAC;QAChB,KAAK,CAAC,MAAM,GAAG,IAAI,KAAK,CAAC,EACzB,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,uDAAuD,CAAC,CAAC;IAChF,CAAC;IACD,yFAAyF;IACzF,MAAM,WAAW,GAAG,WAAW,CAAC,WAAW,IAAI,KAAK,CAAC;IACrD,IACE,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC;QAClC,WAAW,GAAG,CAAC;QACf,WAAW,GAAG,UAAU;QACxB,KAAK,CAAC,MAAM,GAAG,IAAI,GAAG,WAAW,EACjC,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;IAC1D,CAAC;IACD,IACE,CAAC,CAAC,kBAAkB,GAAG,CAAC,CAAC;QACzB,CAAC,kBAAkB,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,kBAAkB,CAAC,CAAC,EACzE,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;IAC9D,CAAC;IACD,MAAM,SAAS,GAAG,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,IAAI,IAAI;QAC7C,0BAA0B,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IAC3D,MAAM,KAAK,GAAmB,EAAE,CAAC;IACjC,KAAK,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC;QACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;YACxF,QAAQ;SACT,CAAC,CAAC;QACH,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QACrB,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI,EAAE,SAAS,EAAE;YACnD,SAAS,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC;IAC9C,CAAC;IACD,OAAO,iBAAiB,CAAC,SAAS,EAAE;QAClC,GAAG,WAAW;QACd,IAAI;QACJ,gBAAgB,EAAE,kBAAkB;QACpC,iBAAiB,EAAE,CAAC,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC;KACzD,CAAC,CAAC;AACL,CAAC"}

package/dist/subdivide-rhumb-polyline.d.ts

@@ -0,0 +1,12 @@
+import type { SubdividePolylineOptions, SubdividedPolyline } from '@math.gl/geometry/polygon';
+import { type RhumbPathOptions } from "./rhumb-path.js";
+export type SubdivideRhumbPolylineOptions = Omit<SubdividePolylineOptions, 'createEdgeSampler' | 'maxSegmentLength'> & Omit<RhumbPathOptions, 'longitudeTurns'> & {
+    /** One safe integer turn count per original edge; copied before sampling. */
+    longitudeTurns?: readonly number[];
+    maxSurfaceDistance?: number;
+};
+/** Rhumb sampling through a continuous transform. This is a sampled error test,
+ * not seam clipping or a universal error certificate. Source provenance is kept.
+ */
+export declare function subdivideRhumbPolyline(input: ArrayLike<number>, options: SubdivideRhumbPolylineOptions): SubdividedPolyline;
+//# sourceMappingURL=subdivide-rhumb-polyline.d.ts.map

package/dist/subdivide-rhumb-polyline.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"subdivide-rhumb-polyline.d.ts","sourceRoot":"","sources":["../src/subdivide-rhumb-polyline.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EACV,wBAAwB,EACxB,kBAAkB,EAEnB,MAAM,2BAA2B,CAAC;AACnC,OAAO,EAAY,KAAK,gBAAgB,EAAC,wBAAqB;AAG9D,MAAM,MAAM,6BAA6B,GAAG,IAAI,CAC9C,wBAAwB,EACxB,mBAAmB,GAAG,kBAAkB,CACzC,GACC,IAAI,CAAC,gBAAgB,EAAE,gBAAgB,CAAC,GAAG;IACzC,6EAA6E;IAC7E,cAAc,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IACnC,kBAAkB,CAAC,EAAE,MAAM,CAAC;CAC7B,CAAC;AAEJ;;GAEG;AACH,wBAAgB,sBAAsB,CACpC,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,EACxB,OAAO,EAAE,6BAA6B,GACrC,kBAAkB,CA6EpB"}

package/dist/subdivide-rhumb-polyline.js

@@ -0,0 +1,82 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original rhumb adapter for the existing adaptive subdivision engine.
+import { subdividePolyline } from '@math.gl/geometry/polygon';
+import { RhumbPath } from "./rhumb-path.js";
+import { validateGeographicPosition } from "./geodesic-path.js";
+/** Rhumb sampling through a continuous transform. This is a sampled error test,
+ * not seam clipping or a universal error certificate. Source provenance is kept.
+ */
+export function subdivideRhumbPolyline(input, options) {
+    const { rhumb, longitudeTurns, longitudeMode, size = 2, maxSurfaceDistance = Infinity, ...subdivision } = options;
+    const maxVertices = subdivision.maxVertices ?? 65536;
+    if ((longitudeMode !== undefined && !['shortest', 'unrolled'].includes(longitudeMode)) ||
+        (longitudeMode === 'unrolled' && longitudeTurns !== undefined) ||
+        (size !== 2 && size !== 3) ||
+        !Number.isSafeInteger(input.length) ||
+        input.length < 0 ||
+        input.length % size !== 0 ||
+        !Number.isSafeInteger(maxVertices) ||
+        maxVertices < 1 ||
+        maxVertices > 0xffffffff ||
+        input.length / size > maxVertices ||
+        !(maxSurfaceDistance > 0) ||
+        (maxSurfaceDistance !== Infinity && !Number.isFinite(maxSurfaceDistance)))
+        throw new RangeError('Invalid rhumb subdivision dimensions or limits');
+    const count = Math.max(0, input.length / size - 1);
+    const turns = longitudeTurns ? Array.from(longitudeTurns) : undefined;
+    if (turns && (turns.length !== count || !turns.every(Number.isSafeInteger)))
+        throw new RangeError('Expected one integer longitudeTurns value per rhumb edge');
+    const positions = Array.from(input);
+    for (let i = 0; i < positions.length; i += size)
+        validateGeographicPosition(positions.slice(i, i + size));
+    const original = positions.slice();
+    const paths = [];
+    for (let i = size; i < positions.length; i += size) {
+        const path = new RhumbPath(original.slice(i - size, i), original.slice(i, i + size), {
+            rhumb,
+            longitudeMode,
+            longitudeTurns: turns?.[paths.length]
+        });
+        // Select each signed interval before translating it into the accumulated
+        // world. Otherwise a preceding wrap can reverse an exact 180-degree tie.
+        const first = positions.slice(i - size, i);
+        const worldTurns = Math.round((first[0] - original[i - size]) / 360);
+        const offset = 360 * worldTurns;
+        if (!Number.isSafeInteger(worldTurns) || original[i - size] + offset !== first[0])
+            throw new RangeError('Unrepresentable rhumb longitude chart');
+        const last = path.end;
+        const endTurns = Math.round((last[0] - original[i]) / 360);
+        const translated = last[0] + offset;
+        const endpoint = original[i] + (offset + 360 * endTurns);
+        last[0] =
+            Number.isSafeInteger(endTurns) &&
+                Number.isFinite(endpoint) &&
+                Math.abs(endpoint - translated) <= 4 * Number.EPSILON * Math.max(360, Math.abs(translated))
+                ? endpoint
+                : translated;
+        if (path.length === 0) {
+            last[0] = first[0];
+            last[1] = first[1];
+        }
+        const sample = (fraction) => {
+            if (fraction === 0)
+                return first.slice();
+            if (fraction === 1)
+                return last.slice();
+            const p = path.sample(fraction);
+            p[0] += offset;
+            return p;
+        };
+        paths.push({ length: path.length, sample });
+        positions.splice(i, size, ...sample(1));
+    }
+    return subdividePolyline(positions, {
+        ...subdivision,
+        size,
+        maxSegmentLength: maxSurfaceDistance,
+        createEdgeSampler: (_start, _end, index) => paths[index]
+    });
+}
+//# sourceMappingURL=subdivide-rhumb-polyline.js.map

package/dist/subdivide-rhumb-polyline.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"subdivide-rhumb-polyline.js","sourceRoot":"","sources":["../src/subdivide-rhumb-polyline.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,yFAAyF;AAEzF,OAAO,EAAC,iBAAiB,EAAC,MAAM,2BAA2B,CAAC;AAM5D,OAAO,EAAC,SAAS,EAAwB,wBAAqB;AAC9D,OAAO,EAAC,0BAA0B,EAAC,2BAAwB;AAY3D;;GAEG;AACH,MAAM,UAAU,sBAAsB,CACpC,KAAwB,EACxB,OAAsC;IAEtC,MAAM,EACJ,KAAK,EACL,cAAc,EACd,aAAa,EACb,IAAI,GAAG,CAAC,EACR,kBAAkB,GAAG,QAAQ,EAC7B,GAAG,WAAW,EACf,GAAG,OAAO,CAAC;IACZ,MAAM,WAAW,GAAG,WAAW,CAAC,WAAW,IAAI,KAAK,CAAC;IACrD,IACE,CAAC,aAAa,KAAK,SAAS,IAAI,CAAC,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC,QAAQ,CAAC,aAAa,CAAC,CAAC;QAClF,CAAC,aAAa,KAAK,UAAU,IAAI,cAAc,KAAK,SAAS,CAAC;QAC9D,CAAC,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC;QAC1B,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC;QACnC,KAAK,CAAC,MAAM,GAAG,CAAC;QAChB,KAAK,CAAC,MAAM,GAAG,IAAI,KAAK,CAAC;QACzB,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC;QAClC,WAAW,GAAG,CAAC;QACf,WAAW,GAAG,UAAU;QACxB,KAAK,CAAC,MAAM,GAAG,IAAI,GAAG,WAAW;QACjC,CAAC,CAAC,kBAAkB,GAAG,CAAC,CAAC;QACzB,CAAC,kBAAkB,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,kBAAkB,CAAC,CAAC;QAEzE,MAAM,IAAI,UAAU,CAAC,gDAAgD,CAAC,CAAC;IACzE,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,MAAM,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;IACnD,MAAM,KAAK,GAAG,cAAc,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;IACtE,IAAI,KAAK,IAAI,CAAC,KAAK,CAAC,MAAM,KAAK,KAAK,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC;QACzE,MAAM,IAAI,UAAU,CAAC,0DAA0D,CAAC,CAAC;IACnF,MAAM,SAAS,GAAG,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,IAAI,IAAI;QAC7C,0BAA0B,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IAC3D,MAAM,QAAQ,GAAG,SAAS,CAAC,KAAK,EAAE,CAAC;IACnC,MAAM,KAAK,GAA0B,EAAE,CAAC;IACxC,KAAK,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC;QACnD,MAAM,IAAI,GAAG,IAAI,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;YACnF,KAAK;YACL,aAAa;YACb,cAAc,EAAE,KAAK,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC;SACtC,CAAC,CAAC;QACH,yEAAyE;QACzE,yEAAyE;QACzE,MAAM,KAAK,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QACrE,MAAM,MAAM,GAAG,GAAG,GAAG,UAAU,CAAC;QAChC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,UAAU,CAAC,IAAI,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,MAAM,KAAK,KAAK,CAAC,CAAC,CAAC;YAC/E,MAAM,IAAI,UAAU,CAAC,uCAAuC,CAAC,CAAC;QAChE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC;QACtB,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC3D,MAAM,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACpC,MAAM,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,QAAQ,CAAC,CAAC;QACzD,IAAI,CAAC,CAAC,CAAC;YACL,MAAM,CAAC,aAAa,CAAC,QAAQ,CAAC;gBAC9B,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;gBACzB,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,UAAU,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;gBACzF,CAAC,CAAC,QAAQ;gBACV,CAAC,CAAC,UAAU,CAAC;QACjB,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YACtB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACnB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACrB,CAAC;QACD,MAAM,MAAM,GAAG,CAAC,QAAgB,EAAY,EAAE;YAC5C,IAAI,QAAQ,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC,KAAK,EAAE,CAAC;YACzC,IAAI,QAAQ,KAAK,CAAC;gBAAE,OAAO,IAAI,CAAC,KAAK,EAAE,CAAC;YACxC,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;YAChC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;YACf,OAAO,CAAC,CAAC;QACX,CAAC,CAAC;QACF,KAAK,CAAC,IAAI,CAAC,EAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAC,CAAC,CAAC;QAC1C,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD,OAAO,iBAAiB,CAAC,SAAS,EAAE;QAClC,GAAG,WAAW;QACd,IAAI;QACJ,gBAAgB,EAAE,kBAAkB;QACpC,iBAAiB,EAAE,CAAC,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC;KACzD,CAAC,CAAC;AACL,CAAC"}

package/dist/triaxial/integrate.d.ts

@@ -0,0 +1,2 @@
+export declare function integrateGeodesic(initial: readonly number[], q: readonly number[], distance: number): number[];
+//# sourceMappingURL=integrate.d.ts.map

package/dist/triaxial/integrate.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"integrate.d.ts","sourceRoot":"","sources":["../../src/triaxial/integrate.ts"],"names":[],"mappings":"AAwBA,wBAAgB,iBAAiB,CAC/B,OAAO,EAAE,SAAS,MAAM,EAAE,EAC1B,CAAC,EAAE,SAAS,MAAM,EAAE,EACpB,QAAQ,EAAE,MAAM,GACf,MAAM,EAAE,CAqDV"}

package/dist/triaxial/integrate.js

@@ -0,0 +1,76 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adaptive Dormand-Prince integration of the constrained Cartesian geodesic equation; mathematical coefficients, no third-party implementation copied.
+const C = [0, 1 / 5, 3 / 10, 4 / 5, 8 / 9, 1, 1];
+const A = [
+    [],
+    [1 / 5],
+    [3 / 40, 9 / 40],
+    [44 / 45, -56 / 15, 32 / 9],
+    [19372 / 6561, -25360 / 2187, 64448 / 6561, -212 / 729],
+    [9017 / 3168, -355 / 33, 46732 / 5247, 49 / 176, -5103 / 18656],
+    [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84]
+];
+const B = [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84, 0];
+const ERROR = [
+    35 / 384 - 5179 / 57600,
+    0,
+    500 / 1113 - 7571 / 16695,
+    125 / 192 - 393 / 640,
+    -2187 / 6784 + 92097 / 339200,
+    11 / 84 - 187 / 2100,
+    -1 / 40
+];
+export function integrateGeodesic(initial, q, distance) {
+    let state = Array.from(initial), done = 0, step = Math.sign(distance) * Math.min(0.01, Math.abs(distance));
+    const derivative = (y) => {
+        const normal = [q[0] * y[0], q[1] * y[1], q[2] * y[2]];
+        const curvature = (q[0] * y[3] ** 2 + q[1] * y[4] ** 2 + q[2] * y[5] ** 2) /
+            (normal[0] ** 2 + normal[1] ** 2 + normal[2] ** 2);
+        return [
+            y[3],
+            y[4],
+            y[5],
+            -curvature * normal[0],
+            -curvature * normal[1],
+            -curvature * normal[2]
+        ];
+    };
+    for (let attempt = 0; attempt < 50000; attempt++) {
+        if (done === distance)
+            return state;
+        if (Math.abs(step) > Math.abs(distance - done))
+            step = distance - done;
+        const stages = [];
+        for (let j = 0; j < C.length; j++) {
+            const y = state.map((v, i) => v + step * A[j].reduce((sum, a, k) => sum + a * stages[k][i], 0));
+            stages.push(derivative(y));
+        }
+        const next = state.map((v, i) => v + step * B.reduce((sum, b, j) => sum + b * stages[j][i], 0));
+        let error = 0;
+        for (let i = 0; i < 6; i++)
+            error = Math.max(error, Math.abs(step * ERROR.reduce((sum, e, j) => sum + e * stages[j][i], 0)));
+        const tolerance = 2e-15;
+        if (error <= tolerance) {
+            state = next;
+            done += step;
+            // Restore the two constraints and unit speed only at accepted steps.
+            const surface = Math.sqrt(q[0] * state[0] ** 2 + q[1] * state[1] ** 2 + q[2] * state[2] ** 2);
+            for (let i = 0; i < 3; i++)
+                state[i] /= surface;
+            const normal = state.slice(0, 3).map((x, i) => q[i] * x), nn = normal.reduce((s, x) => s + x * x, 0);
+            const dot = normal.reduce((s, x, i) => s + x * state[i + 3], 0) / nn;
+            for (let i = 0; i < 3; i++)
+                state[i + 3] -= normal[i] * dot;
+            const speed = Math.hypot(state[3], state[4], state[5]);
+            for (let i = 3; i < 6; i++)
+                state[i] /= speed;
+        }
+        step *= Math.max(0.2, Math.min(4, 0.9 * (tolerance / Math.max(error, 1e-300)) ** 0.2));
+        if (done + step === done)
+            throw new RangeError('Triaxial integration step underflow');
+    }
+    throw new RangeError('Triaxial integration work limit exceeded');
+}
+//# sourceMappingURL=integrate.js.map

package/dist/triaxial/integrate.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"integrate.js","sourceRoot":"","sources":["../../src/triaxial/integrate.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,kLAAkL;AAClL,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACjD,MAAM,CAAC,GAAG;IACR,EAAE;IACF,CAAC,CAAC,GAAG,CAAC,CAAC;IACP,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC;IAChB,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAC3B,CAAC,KAAK,GAAG,IAAI,EAAE,CAAC,KAAK,GAAG,IAAI,EAAE,KAAK,GAAG,IAAI,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC;IACvD,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,KAAK,GAAG,IAAI,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,KAAK,CAAC;IAC/D,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,CAAC;CAC5D,CAAC;AACF,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;AACzE,MAAM,KAAK,GAAG;IACZ,EAAE,GAAG,GAAG,GAAG,IAAI,GAAG,KAAK;IACvB,CAAC;IACD,GAAG,GAAG,IAAI,GAAG,IAAI,GAAG,KAAK;IACzB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG;IACrB,CAAC,IAAI,GAAG,IAAI,GAAG,KAAK,GAAG,MAAM;IAC7B,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,IAAI;IACpB,CAAC,CAAC,GAAG,EAAE;CACR,CAAC;AACF,MAAM,UAAU,iBAAiB,CAC/B,OAA0B,EAC1B,CAAoB,EACpB,QAAgB;IAEhB,IAAI,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,EAC7B,IAAI,GAAG,CAAC,EACR,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;IAClE,MAAM,UAAU,GAAG,CAAC,CAAoB,EAAE,EAAE;QAC1C,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvD,MAAM,SAAS,GACb,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YACxD,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;QACrD,OAAO;YACL,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;YACtB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;YACtB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;SACvB,CAAC;IACJ,CAAC,CAAC;IACF,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,EAAE,OAAO,EAAE,EAAE,CAAC;QACjD,IAAI,IAAI,KAAK,QAAQ;YAAE,OAAO,KAAK,CAAC;QACpC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC;YAAE,IAAI,GAAG,QAAQ,GAAG,IAAI,CAAC;QACvE,MAAM,MAAM,GAAe,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CACjB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAC3E,CAAC;YACF,MAAM,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAChG,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,KAAK,GAAG,IAAI,CAAC,GAAG,CACd,KAAK,EACL,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CACxE,CAAC;QACJ,MAAM,SAAS,GAAG,KAAK,CAAC;QACxB,IAAI,KAAK,IAAI,SAAS,EAAE,CAAC;YACvB,KAAK,GAAG,IAAI,CAAC;YACb,IAAI,IAAI,IAAI,CAAC;YACb,qEAAqE;YACrE,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YAC9F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC;YAChD,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EACtD,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7C,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;YACrE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAC5D,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;QAChD,CAAC;QACD,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QACvF,IAAI,IAAI,GAAG,IAAI,KAAK,IAAI;YAAE,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;IACxF,CAAC;IACD,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;AACnE,CAAC"}