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

package/dist/geodesic-exact/area.js

@@ -0,0 +1,56 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) Charles Karney (2012-2025)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Area decomposition and I4 integrand adapted from GeographicLib 2.7 GeodesicExact/GeodesicLineExact (MIT). Original derivative-integral divided difference avoids cancellation. See ../../GEOGRAPHICLIB-LICENSE.txt.
+import { integrate } from "./quadrature.js";
+/** Derivative of t(x) = x + sqrt(1+x) asinh(sqrt(x))/sqrt(x) - 1,
+ * analytically continued to negative x. Taylor evaluation avoids subtraction at 0.
+ */
+function derivative(x) {
+    if (Math.abs(x) < 0.125) {
+        let coefficient = 1 / 3, power = 1, sum = 1;
+        for (let n = 1; n <= 24; n++) {
+            sum += n * coefficient * power;
+            coefficient *= (-2 * n) / (2 * n + 3);
+            power *= x;
+        }
+        return sum;
+    }
+    const root = Math.sqrt(Math.abs(x));
+    const ratio = (x > 0 ? Math.asinh(root) : Math.asin(root)) / root;
+    return 1 + (1 - ratio / Math.sqrt(1 + x)) / (2 * x);
+}
+/** Stable t value for a well-separated divided difference. */
+function value(x) {
+    if (Math.abs(x) < 0.125) {
+        let term = x / 3, sum = x;
+        for (let n = 1; n <= 24; n++) {
+            const next = sum + term;
+            if (next === sum)
+                break;
+            sum = next;
+            term *= (-2 * n * x) / (2 * n + 3);
+        }
+        return sum;
+    }
+    const root = Math.sqrt(Math.abs(x));
+    return x + ((Math.sqrt(1 + x) * (x > 0 ? Math.asinh(root) : Math.asin(root))) / root - 1);
+}
+/** Integrate the smooth correction, evaluating t's divided difference as an
+ * integral of its derivative rather than subtracting nearly equal t values.
+ */
+export function areaCorrection(ep2, k2, s1, c1, delta) {
+    // I4 is anti-periodic after pi and has zero integral over a full circuit.
+    const remainder = delta % (2 * Math.PI), tX = value(ep2);
+    return integrate(theta => {
+        const s = s1 * Math.cos(theta) + c1 * Math.sin(theta), y = k2 * s * s;
+        const divided = ep2 === y
+            ? derivative(y)
+            : Math.abs(y) <= 1e-8 * Math.abs(ep2)
+                ? (tX - value(y)) / (ep2 - y)
+                : integrate(u => derivative(y + (ep2 - y) * u), 0, 1);
+        return (-s * divided) / 2;
+    }, 0, remainder);
+}
+//# sourceMappingURL=area.js.map

package/dist/geodesic-exact/area.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"area.js","sourceRoot":"","sources":["../../src/geodesic-exact/area.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,mEAAmE;AACnE,4DAA4D;AAC5D,wOAAwO;AACxO,OAAO,EAAC,SAAS,EAAC,wBAAqB;AACvC;;GAEG;AACH,SAAS,UAAU,CAAC,CAAS;IAC3B,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;QACxB,IAAI,WAAW,GAAG,CAAC,GAAG,CAAC,EACrB,KAAK,GAAG,CAAC,EACT,GAAG,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,GAAG,IAAI,CAAC,GAAG,WAAW,GAAG,KAAK,CAAC;YAC/B,WAAW,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACtC,KAAK,IAAI,CAAC,CAAC;QACb,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC;IAClE,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACtD,CAAC;AACD,8DAA8D;AAC9D,SAAS,KAAK,CAAC,CAAS;IACtB,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;QACxB,IAAI,IAAI,GAAG,CAAC,GAAG,CAAC,EACd,GAAG,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;YACxB,IAAI,IAAI,KAAK,GAAG;gBAAE,MAAM;YACxB,GAAG,GAAG,IAAI,CAAC;YACX,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACrC,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpC,OAAO,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;AAC5F,CAAC;AACD;;GAEG;AACH,MAAM,UAAU,cAAc,CAC5B,GAAW,EACX,EAAU,EACV,EAAU,EACV,EAAU,EACV,KAAa;IAEb,0EAA0E;IAC1E,MAAM,SAAS,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EACrC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;IAClB,OAAO,SAAS,CACd,KAAK,CAAC,EAAE;QACN,MAAM,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EACnD,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,OAAO,GACX,GAAG,KAAK,CAAC;YACP,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC;YACf,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC;gBACnC,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;gBAC7B,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5D,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;IAC5B,CAAC,EACD,CAAC,EACD,SAAS,CACV,CAAC;AACJ,CAAC"}

package/dist/geodesic-exact/inverse-engine.d.ts

@@ -0,0 +1,57 @@
+import type { GeodesicResult } from "../ellipsoid-geodesic.js";
+type State = {
+    sig12: number;
+    eps: number;
+    ssig1: number;
+    csig1: number;
+    ssig2: number;
+    csig2: number;
+    salp2: number;
+    calp2: number;
+    lam12: number;
+    domg12: number;
+    dlam12: number;
+};
+type Lengths = {
+    s12b: number;
+    m12b: number;
+    m0: number;
+    M12: number;
+    M21: number;
+};
+interface Engine {
+    a: number;
+    f: number;
+    _b: number;
+    _f1: number;
+    _ep2: number;
+    _e2: number;
+    _c2: number;
+    Inverse: (a: number, b: number, c: number, d: number, mask: number) => GeodesicResult;
+    Lengths: (eps: number, sig12: number, s1: number, c1: number, d1: number, s2: number, c2: number, d2: number, cb1: number, cb2: number, mask: number, c1a?: number[], c2a?: number[]) => Lengths;
+    Lambda12: (sb1: number, cb1: number, dn1: number, sb2: number, cb2: number, dn2: number, sa1: number, ca1: number, sl: number, cl: number, diff: boolean, c1a: number[], c2a: number[], c3a: number[]) => State;
+    InverseStart: (sb1: number, cb1: number, dn1: number, sb2: number, cb2: number, dn2: number, lam: number, sl: number, cl: number, c1a: number[], c2a: number[]) => {
+        sig12: number;
+        salp1: number;
+        calp1: number;
+    };
+}
+type Namespace = {
+    Geodesic: new (a: number, f: number) => Engine;
+    ALL: number;
+    AREA: number;
+    LONG_UNROLL: number;
+};
+export declare const library: Namespace;
+export declare const angleMath: {
+    sincosd: (x: number) => {
+        s: number;
+        c: number;
+    };
+    AngNormalize: (x: number) => number;
+    atan2d: (y: number, x: number) => number;
+};
+/** Per-instance overrides only: the existing series engine and its prototype remain unchanged. */
+export declare function exactInverseEngine(a: number, b: number): Engine;
+export {};
+//# sourceMappingURL=inverse-engine.d.ts.map

package/dist/geodesic-exact/inverse-engine.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse-engine.d.ts","sourceRoot":"","sources":["../../src/geodesic-exact/inverse-engine.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,cAAc,EAAC,iCAA8B;AAE1D,KAAK,KAAK,GAAG;IACX,KAAK,EAAE,MAAM,CAAC;IACd,GAAG,EAAE,MAAM,CAAC;IACZ,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IACf,MAAM,EAAE,MAAM,CAAC;CAChB,CAAC;AACF,KAAK,OAAO,GAAG;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,MAAM,CAAA;CAAC,CAAC;AAClF,UAAU,MAAM;IACd,CAAC,EAAE,MAAM,CAAC;IACV,CAAC,EAAE,MAAM,CAAC;IACV,EAAE,EAAE,MAAM,CAAC;IACX,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,OAAO,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,KAAK,cAAc,CAAC;IACtF,OAAO,EAAE,CACP,GAAG,EAAE,MAAM,EACX,KAAK,EAAE,MAAM,EACb,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,IAAI,EAAE,MAAM,EACZ,GAAG,CAAC,EAAE,MAAM,EAAE,EACd,GAAG,CAAC,EAAE,MAAM,EAAE,KACX,OAAO,CAAC;IACb,QAAQ,EAAE,CACR,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,IAAI,EAAE,OAAO,EACb,GAAG,EAAE,MAAM,EAAE,EACb,GAAG,EAAE,MAAM,EAAE,EACb,GAAG,EAAE,MAAM,EAAE,KACV,KAAK,CAAC;IACX,YAAY,EAAE,CACZ,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,GAAG,EAAE,MAAM,EAAE,EACb,GAAG,EAAE,MAAM,EAAE,KACV;QAAC,KAAK,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAC,CAAC;CACpD;AACD,KAAK,SAAS,GAAG;IACf,QAAQ,EAAE,KAAK,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;IAC/C,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,WAAW,EAAE,MAAM,CAAC;CACrB,CAAC;AACF,eAAO,MAAM,OAAO,EAAmC,SAAS,CAAC;AACjE,eAAO,MAAM,SAAS;aAGP,CAAC,CAAC,EAAE,MAAM,KAAK;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAA;KAAC;kBAChC,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM;YAC3B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM;CAGxC,CAAC;AACP,kGAAkG;AAClG,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAwE/D"}

package/dist/geodesic-exact/inverse-engine.js

@@ -0,0 +1,56 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) Charles Karney (2012-2025)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Length/Jacobi and Cayley-longitude formulae adapted from GeographicLib 2.7 GeodesicExact.cpp; original quadrature bridge. See ../../GEOGRAPHICLIB-LICENSE.txt.
+import { geodesic } from "../vendor/geographiclib-geodesic/runtime.js";
+import { ellipticIntegrals } from "./quadrature.js";
+export const library = geodesic.Geodesic;
+export const angleMath = geodesic.Math;
+/** Per-instance overrides only: the existing series engine and its prototype remain unchanged. */
+export function exactInverseEngine(a, b) {
+    const engine = new library.Geodesic(a, (a - b) / a);
+    engine.Lengths = function (eps, delta, s1, c1, d1, s2, c2, d2, cb1, cb2) {
+        const k2 = (4 * eps) / (1 - eps) ** 2, I = ellipticIntegrals(k2, this._ep2), start = Math.atan2(s1, c1);
+        const J = I.J(start, delta), t = (this._ep2 * (cb1 - cb2) * (cb1 + cb2)) / (d1 + d2);
+        const cosine = c1 * c2 + s1 * s2;
+        return {
+            s12b: I.E(start, delta),
+            m12b: d2 * (c1 * s2) - d1 * (s1 * c2) - c1 * c2 * J,
+            m0: I.meanJ(),
+            M12: cosine + ((t * s2 - c2 * J) * s1) / d1,
+            M21: cosine - ((t * s1 - c1 * J) * s2) / d2
+        };
+    };
+    // Force the bracketed solve rather than a series short-line approximation.
+    engine.InverseStart = function (sb1, cb1, _d1, sb2, cb2, _d2, _lam, sl, cl) {
+        let s = cb2 * sl, c = cb1 * sb2 - sb1 * cb2 * cl;
+        const length = Math.hypot(s, c);
+        if (length > 0 && s > 0) {
+            s /= length;
+            c /= length;
+        }
+        else {
+            s = 1;
+            c = 0;
+        }
+        return { sig12: -1, salp1: s, calp1: c };
+    };
+    const baseLambda = engine.Lambda12;
+    engine.Lambda12 = function (sb1, cb1, dn1, sb2, cb2, dn2, sa1, ca1, sl, cl, diff, c1a, c2a, c3a) {
+        const result = baseLambda.call(this, sb1, cb1, dn1, sb2, cb2, dn2, sa1, ca1, sl, cl, diff, c1a, c2a, c3a);
+        const sa0 = sa1 * cb1, ca0 = Math.hypot(ca1, sa1 * sb1), k2 = ca0 * ca0 * this._ep2;
+        const om1 = sa0 * sb1, co1 = ca1 * cb1, om2 = sa0 * sb2, co2 = result.calp2 * cb2;
+        const ch1 = this._f1 * dn1 * co1, ch2 = this._f1 * dn2 * co2;
+        const sch = Math.max(0, ch1 * om2 - om1 * ch2), cch = ch1 * ch2 + om1 * om2;
+        const so = Math.max(0, co1 * om2 - om1 * co2), co = co1 * co2 + om1 * om2;
+        const correction = (-this._e2 / this._f1) *
+            sa0 *
+            ellipticIntegrals(k2, this._ep2).H(Math.atan2(result.ssig1, result.csig1), result.sig12);
+        result.lam12 = Math.atan2(sch * cl - cch * sl, cch * cl + sch * sl) + correction;
+        result.domg12 = correction + Math.atan2(sch * co - cch * so, cch * co + sch * so);
+        return result;
+    };
+    return engine;
+}
+//# sourceMappingURL=inverse-engine.js.map

package/dist/geodesic-exact/inverse-engine.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse-engine.js","sourceRoot":"","sources":["../../src/geodesic-exact/inverse-engine.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,mEAAmE;AACnE,4DAA4D;AAC5D,mLAAmL;AACnL,OAAO,EAAC,QAAQ,EAAC,oDAAiD;AAElE,OAAO,EAAC,iBAAiB,EAAC,wBAAqB;AA2E/C,MAAM,CAAC,MAAM,OAAO,GAAG,QAAQ,CAAC,QAAgC,CAAC;AACjE,MAAM,CAAC,MAAM,SAAS,GACpB,QAOD,CAAC,IAAI,CAAC;AACP,kGAAkG;AAClG,MAAM,UAAU,kBAAkB,CAAC,CAAS,EAAE,CAAS;IACrD,MAAM,MAAM,GAAG,IAAI,OAAO,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,CAAC,OAAO,GAAG,UAAU,GAAG,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG;QACrE,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,EACnC,CAAC,GAAG,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,EACpC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,EACzB,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QAC1D,MAAM,MAAM,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACjC,OAAO;YACL,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC;YACvB,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;YACnD,EAAE,EAAE,CAAC,CAAC,KAAK,EAAE;YACb,GAAG,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE;YAC3C,GAAG,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE;SAC5C,CAAC;IACJ,CAAC,CAAC;IACF,2EAA2E;IAC3E,MAAM,CAAC,YAAY,GAAG,UAAU,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE;QACxE,IAAI,CAAC,GAAG,GAAG,GAAG,EAAE,EACd,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,IAAI,MAAM,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACxB,CAAC,IAAI,MAAM,CAAC;YACZ,CAAC,IAAI,MAAM,CAAC;QACd,CAAC;aAAM,CAAC;YACN,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC;QACR,CAAC;QACD,OAAO,EAAC,KAAK,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAC,CAAC;IACzC,CAAC,CAAC;IACF,MAAM,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC;IACnC,MAAM,CAAC,QAAQ,GAAG,UAAU,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG;QAC7F,MAAM,MAAM,GAAG,UAAU,CAAC,IAAI,CAC5B,IAAI,EACJ,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,EAAE,EACF,EAAE,EACF,IAAI,EACJ,GAAG,EACH,GAAG,EACH,GAAG,CACJ,CAAC;QACF,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,EACnB,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,EAChC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC;QAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,EACnB,GAAG,GAAG,GAAG,GAAG,GAAG,EACf,GAAG,GAAG,GAAG,GAAG,GAAG,EACf,GAAG,GAAG,MAAM,CAAC,KAAK,GAAG,GAAG,CAAC;QAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,EAC9B,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC7B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,EAC5C,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,EAC3C,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC7B,MAAM,UAAU,GACd,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;YACtB,GAAG;YACH,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QAC3F,MAAM,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,UAAU,CAAC;QACjF,MAAM,CAAC,MAAM,GAAG,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC;QAClF,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC;IACF,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/dist/geodesic-exact/polygon.d.ts

@@ -0,0 +1,12 @@
+import type { GeodesicParameters } from "../ellipsoid-geodesic.js";
+import { SurfacePolygon } from "../ellipsoid-helpers/surface-polygon.js";
+/** Compensated exact-geodesic polygon or open polyline. Angles are degrees. */
+export declare class GeodesicPolygonExact extends SurfacePolygon {
+    constructor(parameters: GeodesicParameters, options?: {
+        polyline?: boolean;
+    });
+    addPoint(latitude: number, longitude: number, options?: {
+        longitudeTurns?: number;
+    }): this;
+}
+//# sourceMappingURL=polygon.d.ts.map

package/dist/geodesic-exact/polygon.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"polygon.d.ts","sourceRoot":"","sources":["../../src/geodesic-exact/polygon.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,kBAAkB,EAAC,iCAA8B;AAC9D,OAAO,EAAC,cAAc,EAAC,gDAA6C;AACpE,+EAA+E;AAC/E,qBAAa,oBAAqB,SAAQ,cAAc;gBAC1C,UAAU,EAAE,kBAAkB,EAAE,OAAO,GAAE;QAAC,QAAQ,CAAC,EAAE,OAAO,CAAA;KAAM;IAarE,QAAQ,CACf,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,GAAE;QAAC,cAAc,CAAC,EAAE,MAAM,CAAA;KAAM,GACtC,IAAI;CASR"}

package/dist/geodesic-exact/polygon.js

@@ -0,0 +1,27 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Exact solver facade over the attributed compensated SurfacePolygon kernel; see ../../GEOGRAPHICLIB-LICENSE.txt.
+import { EllipsoidGeodesicExact } from "../geodesic-exact.js";
+import { SurfacePolygon } from "../ellipsoid-helpers/surface-polygon.js";
+/** Compensated exact-geodesic polygon or open polyline. Angles are degrees. */
+export class GeodesicPolygonExact extends SurfacePolygon {
+    constructor(parameters, options = {}) {
+        const solver = new EllipsoidGeodesicExact(parameters);
+        // Preserve the 360-degree transit parity while bounding internal labels.
+        // External point inputs retain the solver's documented magnitude limit.
+        super({
+            surfaceArea: solver.surfaceArea,
+            inverse: (a, b, c, d, opts) => solver.inverse(a, b % 720, c, d % 720, opts),
+            direct: (a, b, c, d, opts) => solver.direct(a, b % 720, c, d, opts)
+        }, options);
+    }
+    addPoint(latitude, longitude, options = {}) {
+        if (Math.abs(longitude) > 1e6)
+            throw new RangeError('Exact polygon longitude exceeds input profile');
+        if (options.longitudeTurns !== undefined && options.longitudeTurns !== 0)
+            throw new RangeError('Exact geodesic point insertion does not select winding routes; use addEdge');
+        return super.addPoint(latitude, longitude);
+    }
+}
+//# sourceMappingURL=polygon.js.map

package/dist/geodesic-exact/polygon.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"polygon.js","sourceRoot":"","sources":["../../src/geodesic-exact/polygon.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,oIAAoI;AACpI,OAAO,EAAC,sBAAsB,EAAC,6BAA0B;AAEzD,OAAO,EAAC,cAAc,EAAC,gDAA6C;AACpE,+EAA+E;AAC/E,MAAM,OAAO,oBAAqB,SAAQ,cAAc;IACtD,YAAY,UAA8B,EAAE,UAAgC,EAAE;QAC5E,MAAM,MAAM,GAAG,IAAI,sBAAsB,CAAC,UAAU,CAAC,CAAC;QACtD,yEAAyE;QACzE,wEAAwE;QACxE,KAAK,CACH;YACE,WAAW,EAAE,MAAM,CAAC,WAAW;YAC/B,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,IAAI,CAAC;YAC3E,MAAM,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC;SACpE,EACD,OAAO,CACR,CAAC;IACJ,CAAC;IACQ,QAAQ,CACf,QAAgB,EAChB,SAAiB,EACjB,UAAqC,EAAE;QAEvC,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,GAAG,GAAG;YAC3B,MAAM,IAAI,UAAU,CAAC,+CAA+C,CAAC,CAAC;QACxE,IAAI,OAAO,CAAC,cAAc,KAAK,SAAS,IAAI,OAAO,CAAC,cAAc,KAAK,CAAC;YACtE,MAAM,IAAI,UAAU,CAClB,4EAA4E,CAC7E,CAAC;QACJ,OAAO,KAAK,CAAC,QAAQ,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC7C,CAAC;CACF"}

package/dist/geodesic-exact/quadrature.d.ts

@@ -0,0 +1,9 @@
+export declare function integrate(f: (x: number) => number, start: number, length: number): number;
+export declare function ellipticIntegrals(k2: number, ep2: number): {
+    E: (s: number, d: number) => number;
+    J: (s: number, d: number) => number;
+    H: (s: number, d: number) => number;
+    meanE: () => number;
+    meanJ: () => number;
+};
+//# sourceMappingURL=quadrature.d.ts.map

package/dist/geodesic-exact/quadrature.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"quadrature.d.ts","sourceRoot":"","sources":["../../src/geodesic-exact/quadrature.ts"],"names":[],"mappings":"AAcA,wBAAgB,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM,CAmBzF;AACD,wBAAgB,iBAAiB,CAAC,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;WAmB9C,MAAM,KAAK,MAAM;WACjB,MAAM,KAAK,MAAM;WACjB,MAAM,KAAK,MAAM;;;EAI3B"}

package/dist/geodesic-exact/quadrature.js

@@ -0,0 +1,62 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original bounded adaptive Gauss-Legendre evaluation of elliptic geodesic integrals.
+const X16 = [
+    0.09501250983763744, 0.2816035507792589, 0.4580167776572274, 0.6178762444026438,
+    0.755404408355003, 0.8656312023878318, 0.9445750230732326, 0.9894009349916499
+];
+const W16 = [
+    0.1894506104550685, 0.1826034150449236, 0.16915651939500254, 0.14959598881657673,
+    0.12462897125553387, 0.09515851168249278, 0.06225352393864789, 0.027152459411754096
+];
+const X8 = [0.1834346424956498, 0.525532409916329, 0.7966664774136267, 0.9602898564975363];
+const W8 = [0.362683783378362, 0.3137066458778873, 0.2223810344533745, 0.1012285362903763];
+export function integrate(f, start, length) {
+    if (length === 0)
+        return 0;
+    function rule(x, w, a, h) {
+        const mid = a + h / 2, half = h / 2;
+        let sum = 0;
+        for (let i = 0; i < x.length; i++)
+            sum += w[i] * (f(mid - half * x[i]) + f(mid + half * x[i]));
+        return half * sum;
+    }
+    function visit(a, h, depth) {
+        const high = rule(X16, W16, a, h), low = rule(X8, W8, a, h);
+        if (!Number.isFinite(high))
+            throw new RangeError('Nonfinite elliptic integral');
+        if (Math.abs(high - low) <= 2e-15 * Math.abs(h) + 8 * Number.EPSILON * Math.abs(high))
+            return high;
+        if (depth === 16)
+            throw new RangeError('Elliptic quadrature did not converge');
+        return visit(a, h / 2, depth + 1) + visit(a + h / 2, h / 2, depth + 1);
+    }
+    return visit(start, length, 0);
+}
+export function ellipticIntegrals(k2, ep2) {
+    const dn = (x) => Math.sqrt(1 + k2 * Math.sin(x) ** 2);
+    const densityD = (x) => Math.sin(x) ** 2 / dn(x);
+    const densityH = (x) => Math.cos(x) ** 2 / ((1 + ep2 * Math.sin(x) ** 2) * dn(x));
+    const complete = new Map();
+    const total = (f) => {
+        let v = complete.get(f);
+        if (v === undefined) {
+            v = integrate(f, 0, Math.PI);
+            complete.set(f, v);
+        }
+        return v;
+    };
+    const periodic = (f, start, delta) => {
+        const turns = Math.trunc(delta / Math.PI), remainder = delta - turns * Math.PI;
+        return (turns ? turns * total(f) : 0) + integrate(f, start % Math.PI, remainder);
+    };
+    return {
+        E: (s, d) => periodic(dn, s, d),
+        J: (s, d) => k2 * periodic(densityD, s, d),
+        H: (s, d) => periodic(densityH, s, d),
+        meanE: () => total(dn) / Math.PI,
+        meanJ: () => (k2 * total(densityD)) / Math.PI
+    };
+}
+//# sourceMappingURL=quadrature.js.map

package/dist/geodesic-exact/quadrature.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"quadrature.js","sourceRoot":"","sources":["../../src/geodesic-exact/quadrature.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,wGAAwG;AACxG,MAAM,GAAG,GAAG;IACV,mBAAmB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;IAC/E,iBAAiB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC9E,CAAC;AACF,MAAM,GAAG,GAAG;IACV,kBAAkB,EAAE,kBAAkB,EAAE,mBAAmB,EAAE,mBAAmB;IAChF,mBAAmB,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,oBAAoB;CACpF,CAAC;AACF,MAAM,EAAE,GAAG,CAAC,kBAAkB,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAC3F,MAAM,EAAE,GAAG,CAAC,iBAAiB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAC3F,MAAM,UAAU,SAAS,CAAC,CAAwB,EAAE,KAAa,EAAE,MAAc;IAC/E,IAAI,MAAM,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAC3B,SAAS,IAAI,CAAC,CAAoB,EAAE,CAAoB,EAAE,CAAS,EAAE,CAAS;QAC5E,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EACnB,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;QACf,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/F,OAAO,IAAI,GAAG,GAAG,CAAC;IACpB,CAAC;IACD,SAAS,KAAK,CAAC,CAAS,EAAE,CAAS,EAAE,KAAa;QAChD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAC/B,GAAG,GAAG,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,CAAC,CAAC;QAChF,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;YACnF,OAAO,IAAI,CAAC;QACd,IAAI,KAAK,KAAK,EAAE;YAAE,MAAM,IAAI,UAAU,CAAC,sCAAsC,CAAC,CAAC;QAC/E,OAAO,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;IACzE,CAAC;IACD,OAAO,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;AACjC,CAAC;AACD,MAAM,UAAU,iBAAiB,CAAC,EAAU,EAAE,GAAW;IACvD,MAAM,EAAE,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IAC/D,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACzD,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1F,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAoB,CAAC;IAC7C,MAAM,KAAK,GAAG,CAAC,CAAwB,EAAE,EAAE;QACzC,IAAI,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;YACpB,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;YAC7B,QAAQ,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACrB,CAAC;QACD,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;IACF,MAAM,QAAQ,GAAG,CAAC,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,EAAE;QAC1E,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,EACvC,SAAS,GAAG,KAAK,GAAG,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QACtC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC,EAAE,EAAE,SAAS,CAAC,CAAC;IACnF,CAAC,CAAC;IACF,OAAO;QACL,CAAC,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,QAAQ,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;QAC/C,CAAC,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,EAAE,GAAG,QAAQ,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC;QAC1D,CAAC,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,QAAQ,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC;QACrD,KAAK,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;QAChC,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE;KAC9C,CAAC;AACJ,CAAC"}

package/dist/geodesic-exact.d.ts

@@ -0,0 +1,23 @@
+import type { GeodesicParameters, GeodesicResult, GeodesicLine } from "./ellipsoid-geodesic.js";
+export type ExactGeodesicOptions = {
+    unrollLongitude?: boolean;
+    area?: boolean;
+};
+/** Elliptic-integral revolution geodesics, independently selected from the series solver.
+ * Angles are degrees. Qualified axis ratio b/a is [0.25,4], axis scales [1e-60,1e60].
+ * Signed edge area is opt-in and uses squared axis units.
+ */
+export declare class EllipsoidGeodesicExact {
+    private readonly engine;
+    private readonly axes;
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }: GeodesicParameters);
+    get semiMajorAxis(): number;
+    get semiMinorAxis(): number;
+    get surfaceArea(): number;
+    inverse(lat1: number, lon1: number, lat2: number, lon2: number, options?: ExactGeodesicOptions): GeodesicResult;
+    direct(lat: number, lon: number, azimuth: number, distance: number, options?: ExactGeodesicOptions): GeodesicResult;
+    line(lat: number, lon: number, azimuth: number, options?: ExactGeodesicOptions): GeodesicLine;
+}
+export { GeodesicPolygonExact } from "./geodesic-exact/polygon.js";
+export type { SurfacePolygonOptions as ExactPolygonOptions, SurfacePolygonResult as ExactPolygonResult } from "./ellipsoid-helpers/surface-polygon.js";
+//# sourceMappingURL=geodesic-exact.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"geodesic-exact.d.ts","sourceRoot":"","sources":["../src/geodesic-exact.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,kBAAkB,EAAE,cAAc,EAAE,YAAY,EAAC,gCAA6B;AAI3F,MAAM,MAAM,oBAAoB,GAAG;IAAC,eAAe,CAAC,EAAE,OAAO,CAAC;IAAC,IAAI,CAAC,EAAE,OAAO,CAAA;CAAC,CAAC;AAS/E;;;GAGG;AACH,qBAAa,sBAAsB;IACjC,OAAO,CAAC,QAAQ,CAAC,MAAM,CAAwC;IAC/D,OAAO,CAAC,QAAQ,CAAC,IAAI,CAA4B;gBACrC,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAC,EAAE,kBAAkB;IAapE,IAAI,aAAa,IAAI,MAAM,CAE1B;IACD,IAAI,aAAa,IAAI,MAAM,CAE1B;IACD,IAAI,WAAW,IAAI,MAAM,CAExB;IACD,OAAO,CACL,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,OAAO,GAAE,oBAAyB,GACjC,cAAc;IA8EjB,MAAM,CACJ,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,OAAO,EAAE,MAAM,EACf,QAAQ,EAAE,MAAM,EAChB,OAAO,GAAE,oBAAyB,GACjC,cAAc;IAGjB,IAAI,CACF,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,EACX,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,oBAAyB,GACjC,YAAY;CA0IhB;AAED,OAAO,EAAC,oBAAoB,EAAC,oCAAiC;AAC9D,YAAY,EACV,qBAAqB,IAAI,mBAAmB,EAC5C,oBAAoB,IAAI,kBAAkB,EAC3C,+CAA4C"}

package/dist/geodesic-exact.js

@@ -0,0 +1,210 @@
+import { exactInverseEngine, angleMath, library } from "./geodesic-exact/inverse-engine.js";
+import { ellipticIntegrals } from "./geodesic-exact/quadrature.js";
+import { areaCorrection } from "./geodesic-exact/area.js";
+const DEG = Math.PI / 180;
+function point(lat, lon) {
+    if (!Number.isFinite(lat) || Math.abs(lat) > 90 || !Number.isFinite(lon) || Math.abs(lon) > 1e6)
+        throw new RangeError('Expected latitude in [-90,90] and longitude in [-1e6,1e6]');
+}
+function finite(x) {
+    if (!Number.isFinite(x))
+        throw new RangeError('Expected finite geodesic argument');
+}
+/** Elliptic-integral revolution geodesics, independently selected from the series solver.
+ * Angles are degrees. Qualified axis ratio b/a is [0.25,4], axis scales [1e-60,1e60].
+ * Signed edge area is opt-in and uses squared axis units.
+ */
+export class EllipsoidGeodesicExact {
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }) {
+        if (!Number.isFinite(a) ||
+            !Number.isFinite(b) ||
+            Math.min(a, b) < 1e-60 ||
+            Math.max(a, b) > 1e60 ||
+            b / a < 0.25 ||
+            b / a > 4)
+            throw new RangeError('Exact geodesics require axes in [1e-60,1e60] and b/a in [0.25,4]');
+        this.axes = Object.freeze([a, b]);
+        this.engine = exactInverseEngine(a, b);
+    }
+    get semiMajorAxis() {
+        return this.axes[0];
+    }
+    get semiMinorAxis() {
+        return this.axes[1];
+    }
+    get surfaceArea() {
+        return 4 * Math.PI * this.engine._c2;
+    }
+    inverse(lat1, lon1, lat2, lon2, options = {}) {
+        point(lat1, lon1);
+        point(lat2, lon2);
+        const result = this.engine.Inverse(lat1, lon1, lat2, lon2, (library.ALL & ~library.AREA) | (options.unrollLongitude ? library.LONG_UNROLL : 0));
+        if (![result.s12, result.azi1, result.azi2, result.m12, result.M12, result.M21, result.a12].every(Number.isFinite))
+            throw new RangeError('Exact inverse did not converge');
+        let azimuth = result.azi1, distance = result.s12;
+        const target = angleMath.sincosd(lat2), f1 = this.semiMinorAxis / this.semiMajorAxis;
+        const denominator = Math.sqrt(1 - (1 - f1 * f1) * target.s * target.s);
+        const meridian = (this.semiMajorAxis * f1 * f1) / denominator ** 3;
+        const parallel = (this.semiMajorAxis * target.c) / denominator;
+        for (let iteration = 0; iteration < 8; iteration++) {
+            const endpoint = this.direct(lat1, lon1, azimuth, distance, { area: options.area });
+            const north = meridian * (lat2 - endpoint.lat2) * DEG;
+            const east = parallel * angleMath.AngNormalize(angleMath.AngNormalize(lon2) - endpoint.lon2) * DEG;
+            // Propagate angular subtraction roundoff into physical endpoint units.
+            // Longitude labels near 180 degrees can lose several ulps even on short routes.
+            const rounding = 8 *
+                Number.EPSILON *
+                (options.area &&
+                    distance !== 0 &&
+                    Math.abs(distance) < 1e-5 * Math.max(this.semiMajorAxis, this.semiMinorAxis)
+                    ? Math.abs(distance)
+                    : Math.max(this.semiMajorAxis, this.semiMinorAxis)) +
+                4 *
+                    Number.EPSILON *
+                    DEG *
+                    (Math.abs(meridian) * (Math.abs(lat2) + Math.abs(endpoint.lat2)) +
+                        Math.abs(parallel) *
+                            (Math.abs(angleMath.AngNormalize(lon2)) + Math.abs(endpoint.lon2)));
+            if (Math.hypot(north, east) <= rounding) {
+                let edgeArea = endpoint.S12;
+                if (options.area && Math.abs(lat1) < Math.abs(lat2) && Math.abs(lat2) < 90) {
+                    // Evaluate area from the more polar endpoint, so terminal bearing
+                    // roundoff is not amplified at a nearly polar endpoint. Reversal is
+                    // algebraically exact; the reversed call cannot take this branch.
+                    edgeArea = -this.inverse(lat2, lon2, lat1, lon1, { area: true }).S12;
+                }
+                return {
+                    ...endpoint,
+                    lat2,
+                    lon1: result.lon1,
+                    lon2: result.lon2,
+                    ...(options.area ? { S12: edgeArea } : {}),
+                    ...(options.area && Math.abs(lat2) === 90
+                        ? { S12: this.engine._c2 * angleMath.AngNormalize(result.azi2 - result.azi1) * DEG }
+                        : {}),
+                    // At an exact pole, preserve the inverse framework's requested-meridian
+                    // bearing; a rounded direct endpoint cannot determine that convention.
+                    azi2: Math.abs(lat2) === 90 ? result.azi2 : endpoint.azi2
+                };
+            }
+            if (distance === 0 ||
+                (options.area ? endpoint.m12 === 0 : Math.abs(endpoint.m12) < 1e-15 * this.semiMajorAxis))
+                break;
+            const heading = angleMath.sincosd(endpoint.azi2);
+            distance += heading.c * north + heading.s * east;
+            azimuth += (heading.c * east - heading.s * north) / endpoint.m12 / DEG;
+        }
+        throw new RangeError('Exact inverse endpoint residual');
+    }
+    direct(lat, lon, azimuth, distance, options = {}) {
+        return this.line(lat, lon, azimuth, options).position(distance);
+    }
+    line(lat, lon, azimuth, options = {}) {
+        point(lat, lon);
+        finite(azimuth);
+        const a = this.semiMajorAxis, b = this.semiMinorAxis, f1 = b / a, ep2 = (a / b) ** 2 - 1, e2 = 1 - f1 * f1;
+        const phi = angleMath.sincosd(lat), alpha = angleMath.sincosd(angleMath.AngNormalize(azimuth));
+        const norm = Math.hypot(phi.s * f1, phi.c), sb1 = (phi.s * f1) / norm, cb1 = Math.max(1e-150, phi.c / norm);
+        const sa0 = alpha.s * cb1, ca0 = Math.hypot(alpha.c, alpha.s * sb1);
+        let s1 = sb1, c1 = sb1 !== 0 || alpha.c !== 0 ? cb1 * alpha.c : 1;
+        const sigmaNorm = Math.hypot(s1, c1);
+        s1 /= sigmaNorm;
+        c1 /= sigmaNorm;
+        const sig1 = Math.atan2(s1, c1);
+        const k2 = ca0 * ca0 * ep2, I = ellipticIntegrals(k2, ep2), dn1 = Math.sqrt(1 + k2 * s1 * s1);
+        const sign = sa0 < 0 || Object.is(sa0, -0) ? -1 : 1;
+        const includeArea = Boolean(options.area);
+        const unroll = Boolean(options.unrollLongitude), azi1 = angleMath.AngNormalize(azimuth);
+        const distanceLimit = b * I.E(sig1, 4 * Math.PI);
+        const evaluate = (delta, distance, arcDegrees) => {
+            const phase = arcDegrees === undefined
+                ? { s: Math.sin(delta), c: Math.cos(delta) }
+                : angleMath.sincosd(arcDegrees);
+            const s2 = s1 * phase.c + c1 * phase.s, c2 = c1 * phase.c - s1 * phase.s;
+            const dn2 = Math.sqrt(1 + k2 * s2 * s2), J = I.J(sig1, delta);
+            const chi2 = Math.atan2(sign * sa0 * s2, f1 * dn2 * c2), signedChi1 = Math.atan2(sign * sa0 * s1, f1 * dn1 * c1);
+            const t = (k2 * (s2 - s1) * (s2 + s1)) / (dn1 + dn2);
+            // Form short angular differences from cross/dot products, avoiding subtraction
+            // of near-pi phases on tiny edges. The dn difference is rationalized above.
+            const turns = Math.trunc(arcDegrees === undefined ? delta / Math.PI : arcDegrees / 180), remainder = arcDegrees === undefined ? delta - turns * Math.PI : arcDegrees - turns * 180;
+            const parity = Math.abs(turns % 2) === 1 ? -1 : 1;
+            const sr = phase.s * parity;
+            const endS = s2 * parity, endC = c2 * parity;
+            const chi = sa0 === 0
+                ? Math.round((sign * (delta - (Math.atan2(s2, c2) - Math.atan2(s1, c1)) + (chi2 - signedChi1))) /
+                    Math.PI) * Math.PI
+                : sign *
+                    (turns * Math.PI +
+                        (remainder === 0
+                            ? 0
+                            : Math.atan2(sign * sa0 * f1 * (dn1 * sr - t * s1 * endC), f1 * f1 * dn1 * dn2 * c1 * endC + sa0 * sa0 * s1 * endS)));
+            const longitude = (unroll ? lon : angleMath.AngNormalize(lon)) +
+                (chi - (e2 / f1) * sa0 * I.H(sig1, delta)) / DEG;
+            const cosine = phase.c;
+            return {
+                lat1: lat,
+                lon1: unroll ? lon : angleMath.AngNormalize(lon),
+                azi1,
+                lat2: angleMath.atan2d(ca0 * s2, f1 * Math.hypot(sa0, ca0 * c2)),
+                lon2: unroll ? longitude : angleMath.AngNormalize(longitude),
+                azi2: angleMath.atan2d(sa0, ca0 * c2),
+                s12: distance === undefined ? b * I.E(sig1, delta) : distance,
+                a12: delta / DEG,
+                m12: b * (dn2 * (c1 * s2) - dn1 * (s1 * c2) - c1 * c2 * J),
+                M12: cosine + ((t * s2 - c2 * J) * s1) / dn1,
+                M21: cosine - ((t * s1 - c1 * J) * s2) / dn2,
+                ...(includeArea
+                    ? {
+                        S12: (() => {
+                            const sd = phase.s, cd = phase.c;
+                            const sine = ca0 *
+                                sa0 *
+                                (cd <= 0 ? c1 * (1 - cd) + sd * s1 : sd * ((c1 * sd) / (1 + cd) + s1));
+                            const cosine = sa0 * sa0 + ca0 * ca0 * c1 * c2;
+                            const angle = ca0 === 0 || sa0 === 0
+                                ? Math.atan2(sa0 * alpha.c - ca0 * c2 * alpha.s, ca0 * c2 * alpha.c + sa0 * alpha.s)
+                                : Math.atan2(sine, cosine);
+                            return (this.engine._c2 * angle +
+                                (a * a * e2 * ca0 * sa0 === 0
+                                    ? 0
+                                    : a * a * e2 * ca0 * sa0 * areaCorrection(ep2, k2, s1, c1, delta)));
+                        })()
+                    }
+                    : {})
+            };
+        };
+        return {
+            arcPosition: (arc) => {
+                finite(arc);
+                if (Math.abs(arc) > 720)
+                    throw new RangeError('Exact line arc exceeds 720 degrees');
+                return evaluate(arc * DEG, undefined, arc);
+            },
+            position: (distance) => {
+                finite(distance);
+                if (Math.abs(distance) > distanceLimit)
+                    throw new RangeError('Exact line distance exceeds support profile');
+                const target = distance / b, lowDn = Math.min(1, a / b), highDn = Math.max(1, a / b);
+                let lower = Math.max(-4 * Math.PI, Math.min(target / lowDn, target / highDn)), upper = Math.min(4 * Math.PI, Math.max(target / lowDn, target / highDn));
+                let delta = Math.max(lower, Math.min(upper, target / I.meanE()));
+                for (let i = 0; i < 80; i++) {
+                    const residual = I.E(sig1, delta) - target;
+                    if (Math.abs(residual) <= Number.EPSILON * Math.max(Math.abs(target), 1e-300))
+                        return evaluate(delta, distance);
+                    if (residual > 0)
+                        upper = delta;
+                    else
+                        lower = delta;
+                    const next = delta - residual / Math.sqrt(1 + k2 * Math.sin(sig1 + delta) ** 2);
+                    const value = next > lower && next < upper ? next : (lower + upper) / 2;
+                    if (value === delta)
+                        return evaluate(delta, distance);
+                    delta = value;
+                }
+                throw new RangeError('Exact distance inversion did not converge');
+            }
+        };
+    }
+}
+export { GeodesicPolygonExact } from "./geodesic-exact/polygon.js";
+//# sourceMappingURL=geodesic-exact.js.map

package/dist/geodesic-exact.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geodesic-exact.js","sourceRoot":"","sources":["../src/geodesic-exact.ts"],"names":[],"mappings":"AAMA,OAAO,EAAC,kBAAkB,EAAE,SAAS,EAAE,OAAO,EAAC,2CAAwC;AACvF,OAAO,EAAC,iBAAiB,EAAC,uCAAoC;AAC9D,OAAO,EAAC,cAAc,EAAC,iCAA8B;AAErD,MAAM,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;AAC1B,SAAS,KAAK,CAAC,GAAW,EAAE,GAAW;IACrC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG;QAC7F,MAAM,IAAI,UAAU,CAAC,2DAA2D,CAAC,CAAC;AACtF,CAAC;AACD,SAAS,MAAM,CAAC,CAAS;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;AACrF,CAAC;AACD;;;GAGG;AACH,MAAM,OAAO,sBAAsB;IAGjC,YAAY,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAqB;QAClE,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YACnB,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YACnB,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK;YACtB,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI;YACrB,CAAC,GAAG,CAAC,GAAG,IAAI;YACZ,CAAC,GAAG,CAAC,GAAG,CAAC;YAET,MAAM,IAAI,UAAU,CAAC,kEAAkE,CAAC,CAAC;QAC3F,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAU,CAAC,CAAC;QAC3C,IAAI,CAAC,MAAM,GAAG,kBAAkB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,IAAI,aAAa;QACf,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,aAAa;QACf,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,IAAI,WAAW;QACb,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC;IACvC,CAAC;IACD,OAAO,CACL,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,UAAgC,EAAE;QAElC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAClB,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAClB,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAChC,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,eAAe,CAAC,CAAC,CAAC,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CACpF,CAAC;QACF,IACE,CAAC,CAAC,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,IAAI,EAAE,MAAM,CAAC,IAAI,EAAE,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,CAC3F,MAAM,CAAC,QAAQ,CAChB;YAED,MAAM,IAAI,UAAU,CAAC,gCAAgC,CAAC,CAAC;QACzD,IAAI,OAAO,GAAG,MAAM,CAAC,IAAI,EACvB,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC;QACxB,MAAM,MAAM,GAAG,SAAS,CAAC,OAAO,CAAC,IAAI,CAAC,EACpC,EAAE,GAAG,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,aAAa,CAAC;QAC/C,MAAM,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACvE,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,aAAa,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,WAAW,IAAI,CAAC,CAAC;QACnE,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,aAAa,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC;QAC/D,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE,CAAC;YACnD,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,IAAI,EAAE,OAAO,EAAE,QAAQ,EAAE,EAAC,IAAI,EAAE,OAAO,CAAC,IAAI,EAAC,CAAC,CAAC;YAClF,MAAM,KAAK,GAAG,QAAQ,GAAG,CAAC,IAAI,GAAG,QAAQ,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;YACtD,MAAM,IAAI,GACR,QAAQ,GAAG,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC,YAAY,CAAC,IAAI,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;YACxF,uEAAuE;YACvE,gFAAgF;YAChF,MAAM,QAAQ,GACZ,CAAC;gBACC,MAAM,CAAC,OAAO;gBACd,CAAC,OAAO,CAAC,IAAI;oBACb,QAAQ,KAAK,CAAC;oBACd,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,EAAE,IAAI,CAAC,aAAa,CAAC;oBAC1E,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;oBACpB,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;gBACvD,CAAC;oBACC,MAAM,CAAC,OAAO;oBACd,GAAG;oBACH,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;wBAC9D,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;4BAChB,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC5E,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,QAAQ,EAAE,CAAC;gBACxC,IAAI,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC;gBAC5B,IAAI,OAAO,CAAC,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;oBAC3E,kEAAkE;oBAClE,oEAAoE;oBACpE,kEAAkE;oBAClE,QAAQ,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EAAC,IAAI,EAAE,IAAI,EAAC,CAAC,CAAC,GAAI,CAAC;gBACtE,CAAC;gBACD,OAAO;oBACL,GAAG,QAAQ;oBACX,IAAI;oBACJ,IAAI,EAAE,MAAM,CAAC,IAAI;oBACjB,IAAI,EAAE,MAAM,CAAC,IAAI;oBACjB,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAC,GAAG,EAAE,QAAQ,EAAC,CAAC,CAAC,CAAC,EAAE,CAAC;oBACxC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,EAAE;wBACvC,CAAC,CAAC,EAAC,GAAG,EAAE,IAAI,CAAC,MAAM,CAAC,GAAG,GAAG,SAAS,CAAC,YAAY,CAAC,MAAM,CAAC,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,GAAG,EAAC;wBAClF,CAAC,CAAC,EAAE,CAAC;oBACP,wEAAwE;oBACxE,uEAAuE;oBACvE,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI;iBAC1D,CAAC;YACJ,CAAC;YACD,IACE,QAAQ,KAAK,CAAC;gBACd,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,aAAa,CAAC;gBAEzF,MAAM;YACR,MAAM,OAAO,GAAG,SAAS,CAAC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YACjD,QAAQ,IAAI,OAAO,CAAC,CAAC,GAAG,KAAK,GAAG,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC;YACjD,OAAO,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,IAAI,GAAG,OAAO,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,QAAQ,CAAC,GAAG,GAAG,GAAG,CAAC;QACzE,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;IAC1D,CAAC;IAED,MAAM,CACJ,GAAW,EACX,GAAW,EACX,OAAe,EACf,QAAgB,EAChB,UAAgC,EAAE;QAElC,OAAO,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAClE,CAAC;IACD,IAAI,CACF,GAAW,EACX,GAAW,EACX,OAAe,EACf,UAAgC,EAAE;QAElC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAChB,MAAM,CAAC,OAAO,CAAC,CAAC;QAChB,MAAM,CAAC,GAAG,IAAI,CAAC,aAAa,EAC1B,CAAC,GAAG,IAAI,CAAC,aAAa,EACtB,EAAE,GAAG,CAAC,GAAG,CAAC,EACV,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EACtB,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;QACnB,MAAM,GAAG,GAAG,SAAS,CAAC,OAAO,CAAC,GAAG,CAAC,EAChC,KAAK,GAAG,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC;QAC7D,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,EACxC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,EACzB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QACvC,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,EACvB,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC3C,IAAI,EAAE,GAAG,GAAG,EACV,EAAE,GAAG,GAAG,KAAK,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtD,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACrC,EAAE,IAAI,SAAS,CAAC;QAChB,EAAE,IAAI,SAAS,CAAC;QAChB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChC,MAAM,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EACxB,CAAC,GAAG,iBAAiB,CAAC,EAAE,EAAE,GAAG,CAAC,EAC9B,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACpC,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,IAAI,MAAM,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpD,MAAM,WAAW,GAAG,OAAO,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;QAC1C,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,CAAC,eAAe,CAAC,EAC7C,IAAI,GAAG,SAAS,CAAC,YAAY,CAAC,OAAO,CAAC,CAAC;QACzC,MAAM,aAAa,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACjD,MAAM,QAAQ,GAAG,CAAC,KAAa,EAAE,QAAiB,EAAE,UAAmB,EAAkB,EAAE;YACzF,MAAM,KAAK,GACT,UAAU,KAAK,SAAS;gBACtB,CAAC,CAAC,EAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAC;gBAC1C,CAAC,CAAC,SAAS,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC;YACpC,MAAM,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,EACpC,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;YACnC,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EACrC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;YACvB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,GAAG,GAAG,EAAE,EAAE,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,EACrD,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,GAAG,GAAG,EAAE,EAAE,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC;YAC1D,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YACrD,+EAA+E;YAC/E,4EAA4E;YAC5E,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,UAAU,GAAG,GAAG,CAAC,EACrF,SAAS,GAAG,UAAU,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,UAAU,GAAG,KAAK,GAAG,GAAG,CAAC;YAC5F,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClD,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC;YAC5B,MAAM,IAAI,GAAG,EAAE,GAAG,MAAM,EACtB,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC;YACrB,MAAM,GAAG,GACP,GAAG,KAAK,CAAC;gBACP,CAAC,CAAC,IAAI,CAAC,KAAK,CACR,CAAC,IAAI,GAAG,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,UAAU,CAAC,CAAC,CAAC;oBAChF,IAAI,CAAC,EAAE,CACV,GAAG,IAAI,CAAC,EAAE;gBACb,CAAC,CAAC,IAAI;oBACJ,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE;wBACd,CAAC,SAAS,KAAK,CAAC;4BACd,CAAC,CAAC,CAAC;4BACH,CAAC,CAAC,IAAI,CAAC,KAAK,CACR,IAAI,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,EAC5C,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,CACxD,CAAC,CAAC,CAAC;YAChB,MAAM,SAAS,GACb,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;gBAC5C,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC;YACnD,MAAM,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC;YACvB,OAAO;gBACL,IAAI,EAAE,GAAG;gBACT,IAAI,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,YAAY,CAAC,GAAG,CAAC;gBAChD,IAAI;gBACJ,IAAI,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,GAAG,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,GAAG,EAAE,CAAC,CAAC;gBAChE,IAAI,EAAE,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC;gBAC5D,IAAI,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,EAAE,GAAG,GAAG,EAAE,CAAC;gBACrC,GAAG,EAAE,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ;gBAC7D,GAAG,EAAE,KAAK,GAAG,GAAG;gBAChB,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;gBAC1D,GAAG,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG;gBAC5C,GAAG,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG;gBAC5C,GAAG,CAAC,WAAW;oBACb,CAAC,CAAC;wBACE,GAAG,EAAE,CAAC,GAAG,EAAE;4BACT,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAChB,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;4BACf,MAAM,IAAI,GACR,GAAG;gCACH,GAAG;gCACH,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;4BACzE,MAAM,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;4BAC/C,MAAM,KAAK,GACT,GAAG,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC;gCACpB,CAAC,CAAC,IAAI,CAAC,KAAK,CACR,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,EAClC,GAAG,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CACnC;gCACH,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,MAAM,CAAC,CAAC;4BAC/B,OAAO,CACL,IAAI,CAAC,MAAM,CAAC,GAAG,GAAG,KAAK;gCACvB,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,KAAK,CAAC;oCAC3B,CAAC,CAAC,CAAC;oCACH,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC,CACrE,CAAC;wBACJ,CAAC,CAAC,EAAE;qBACL;oBACH,CAAC,CAAC,EAAE,CAAC;aACR,CAAC;QACJ,CAAC,CAAC;QACF,OAAO;YACL,WAAW,EAAE,CAAC,GAAW,EAAE,EAAE;gBAC3B,MAAM,CAAC,GAAG,CAAC,CAAC;gBACZ,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG;oBAAE,MAAM,IAAI,UAAU,CAAC,oCAAoC,CAAC,CAAC;gBACpF,OAAO,QAAQ,CAAC,GAAG,GAAG,GAAG,EAAE,SAAS,EAAE,GAAG,CAAC,CAAC;YAC7C,CAAC;YACD,QAAQ,EAAE,CAAC,QAAgB,EAAE,EAAE;gBAC7B,MAAM,CAAC,QAAQ,CAAC,CAAC;gBACjB,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,aAAa;oBACpC,MAAM,IAAI,UAAU,CAAC,6CAA6C,CAAC,CAAC;gBACtE,MAAM,MAAM,GAAG,QAAQ,GAAG,CAAC,EACzB,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAC1B,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9B,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,KAAK,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC,EAC3E,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,KAAK,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;gBAC3E,IAAI,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,MAAM,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,MAAM,CAAC;oBAC3C,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,MAAM,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,MAAM,CAAC;wBAC3E,OAAO,QAAQ,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;oBACnC,IAAI,QAAQ,GAAG,CAAC;wBAAE,KAAK,GAAG,KAAK,CAAC;;wBAC3B,KAAK,GAAG,KAAK,CAAC;oBACnB,MAAM,IAAI,GAAG,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;oBAChF,MAAM,KAAK,GAAG,IAAI,GAAG,KAAK,IAAI,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;oBACxE,IAAI,KAAK,KAAK,KAAK;wBAAE,OAAO,QAAQ,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;oBACtD,KAAK,GAAG,KAAK,CAAC;gBAChB,CAAC;gBACD,MAAM,IAAI,UAAU,CAAC,2CAA2C,CAAC,CAAC;YACpE,CAAC;SACF,CAAC;IACJ,CAAC;CACF;AAED,OAAO,EAAC,oBAAoB,EAAC,oCAAiC"}