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

package/dist/ellipsoid-rhumb.js

@@ -0,0 +1,321 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) Charles Karney (2011-2022)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original meridian/isometric divided-integral rhumb solver. Compensated longitude helper adapted from GeographicLib JS 2.2 Math.AngDiff (MIT); rhumbs independently compared with GeographicLib C++ 2.7. See ../GEOGRAPHICLIB-LICENSE.txt.
+import { remainingPoleDistance } from "./ellipsoid-helpers/rhumb-pole-distance.js";
+import { EllipsoidMetrics, integrate } from "./ellipsoid-metrics.js";
+/** A physical rhumb endpoint is undefined here, as distinct from invalid input or failed convergence. */
+export class RhumbDomainError extends RangeError {
+    constructor(reason) {
+        super(reason === 'pole'
+            ? 'Rhumb longitude/bearing is undefined at an exact pole'
+            : 'Rhumb crosses a pole; continuation longitude is undefined');
+        this.reason = reason;
+        this.name = 'RhumbDomainError';
+    }
+}
+const DEG = Math.PI / 180;
+/** Constant-bearing paths on a revolution ellipsoid; explicit degrees and axis units.
+ * Original integral algorithm, qualified for 0.5 <= b/a <= 2. Exact poles are
+ * rejected rather than replaced by an arbitrary epsilon latitude.
+ */
+export class EllipsoidRhumb {
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }) {
+        if (!(Number.isFinite(a) && Number.isFinite(b) && a > 0 && b > 0))
+            throw new RangeError('Rhumbs require finite positive axes');
+        this.a = a;
+        this.b = b;
+        this.r = b / a;
+        this.metrics = new EllipsoidMetrics({ semiMajorAxis: 1, semiMinorAxis: this.r });
+        this.surfaceArea = a * this.metrics.surfaceArea * a;
+        if (!(Number.isFinite(this.surfaceArea) && this.surfaceArea > 0))
+            throw new RangeError('Rhumb surface area is unrepresentable');
+        Object.freeze(this);
+    }
+    inverse(lat1, lon1, lat2, lon2, options = {}) {
+        validateRhumbPoint(lat1, lon1);
+        validateRhumbPoint(lat2, lon2);
+        const turns = options.longitudeTurns ?? 0;
+        if (!Number.isSafeInteger(turns))
+            throw new RangeError('Rhumb longitudeTurns must be a safe integer');
+        const longitude = angleDifference(lon1, lon2) + 360 * turns;
+        if (!Number.isFinite(longitude))
+            throw new RangeError('Unrepresentable rhumb winding');
+        const dphi = (lat2 - lat1) * DEG;
+        const [meanM, meanP] = this.means(lat1, lat2);
+        const dx = (longitude * DEG * meanM) / meanP;
+        const distance = this.a * Math.hypot(dphi * meanM, dx);
+        const azi = Math.atan2(longitude * DEG, dphi * meanP) / DEG;
+        const result = {
+            lat1,
+            lon1,
+            lat2,
+            lon2: options.unrollLongitude ? lon1 + longitude : normalize(lon2),
+            azi12: azi,
+            s12: distance
+        };
+        if (options.area)
+            result.S12 = this.area(lat1, lat2, longitude * DEG);
+        if (!Number.isFinite(distance))
+            throw new RangeError('Unrepresentable rhumb distance');
+        return result;
+    }
+    direct(latitude, longitude, azimuth, distance, options = {}) {
+        return this.line(latitude, longitude, azimuth, options).position(distance);
+    }
+    line(latitude, longitude, azimuth, options = {}) {
+        validateRhumbPoint(latitude, longitude);
+        if (!Number.isFinite(azimuth))
+            throw new RangeError('Rhumb azimuth must be finite');
+        const [sa, ca] = sincosDegrees(azimuth);
+        const area = Boolean(options.area), unroll = Boolean(options.unrollLongitude);
+        const azi = normalize(azimuth);
+        return {
+            position: (distance) => {
+                if (!Number.isFinite(distance))
+                    throw new RangeError('Rhumb distance must be finite');
+                const delta = distance / this.a;
+                if (!Number.isFinite(delta))
+                    throw new RangeError('Unrepresentable rhumb distance');
+                const north = delta * ca;
+                let poleDistance;
+                if (north !== 0) {
+                    // Integrate the pole gap directly: subtracting two absolute meridian
+                    // ordinates discards valid steps close to a pole, even away from it.
+                    const gap = north > 0 ? this.meridianDelta(latitude, 90) : this.meridianDelta(-90, latitude);
+                    poleDistance = gap - Math.abs(north);
+                    if (Math.abs(poleDistance) < 1e-4 * gap)
+                        poleDistance = remainingPoleDistance(this.a, this.b, latitude, azimuth, distance, Math.sign(north));
+                    if (poleDistance < 0)
+                        throw new RhumbDomainError('pole-crossing');
+                    if (poleDistance === 0)
+                        throw new RhumbDomainError('pole');
+                }
+                const polar = north === 0 ? undefined : this.polarEndpoint(latitude, north, poleDistance);
+                const lat2 = north === 0 ? latitude : (polar?.latitude ?? this.inverseMeridian(latitude, north));
+                if (Math.abs(lat2) === 90)
+                    throw new RhumbDomainError('pole');
+                const [meanM, meanP] = polar?.means ?? this.means(latitude, lat2);
+                const dlon = (delta * sa * meanP) / meanM / DEG;
+                const lon2 = unroll ? longitude + dlon : normalize(normalize(longitude) + dlon);
+                if (!Number.isFinite(lon2))
+                    throw new RangeError('Unrepresentable rhumb longitude');
+                const result = {
+                    lat1: latitude,
+                    lon1: longitude,
+                    lat2,
+                    lon2,
+                    azi12: azi,
+                    s12: distance
+                };
+                if (area)
+                    result.S12 = this.area(latitude, lat2, dlon * DEG, polar?.psi);
+                return result;
+            }
+        };
+    }
+    means(lat1, lat2) {
+        const d = (lat2 - lat1) * DEG;
+        const mid = (lat1 + lat2) / 2;
+        const meanM = this.meridianMean(lat1, lat2);
+        const cmin = Math.min(latitudeSinCos(lat1)[1], latitudeSinCos(lat2)[1]);
+        // Divided integral, not subtraction of nearly equal meridian/isometric values.
+        const meanP = Math.abs(d) < 0.25 * cmin
+            ? integrate(t => {
+                const [s, c] = latitudeSinCos(mid + (t * (lat2 - lat1)) / 2);
+                return this.r ** 2 / (c * (c * c + (this.r * s) ** 2));
+            }, -1, 1) / 2
+            : (this.psi(lat2) - this.psi(lat1)) / d;
+        return [meanM, meanP];
+    }
+    psi(latitude) {
+        return this.metrics.latitudeToAuxiliary(latitude, 'isometric') * DEG;
+    }
+    polarEndpoint(latitude, north, poleDistance) {
+        const sign = Math.abs(latitude) > 80 && Math.abs(north) < Math.abs(this.meridianDelta(0, latitude))
+            ? Math.sign(latitude)
+            : Math.sign(north), theta1 = (90 - sign * latitude) * DEG, limit = 10 * DEG;
+        const density = (theta) => this.r ** 2 / Math.hypot(Math.sin(theta), this.r * Math.cos(theta)) ** 3;
+        const gap = (theta) => integrate(density, 0, theta);
+        const target = sign === Math.sign(north)
+            ? poleDistance
+            : this.meridianDelta(sign > 0 ? latitude : -90, sign > 0 ? 90 : latitude) - sign * north;
+        if (!(target > 0 && target < gap(limit)))
+            return undefined;
+        let lo = 0, hi = limit, theta = target * this.r;
+        if (!(theta > lo && theta < hi))
+            theta = theta1;
+        for (let i = 0; i < 64; i++) {
+            const residual = gap(theta) - target;
+            if (residual === 0)
+                break;
+            if (residual < 0)
+                lo = theta;
+            else
+                hi = theta;
+            const next = theta - residual / density(theta);
+            if (next === theta)
+                break;
+            theta = next > lo && next < hi ? next : (lo + hi) / 2;
+            if (theta === lo || theta === hi)
+                break;
+            if (i === 63)
+                throw new Error('Rhumb polar meridian inversion did not converge');
+        }
+        const span = theta - theta1, mid = (theta + theta1) / 2;
+        const meanM = integrate(t => density(mid + (t * span) / 2), -1, 1) / 2;
+        const meanP = Math.abs(span) < 0.25 * Math.min(Math.sin(theta), Math.sin(theta1))
+            ? integrate(t => {
+                const x = mid + (t * span) / 2, c = Math.sin(x), s = Math.cos(x);
+                return this.r ** 2 / (c * (c * c + (this.r * s) ** 2));
+            }, -1, 1) / 2
+            : -(this.polarPsi(theta) - this.polarPsi(theta1)) / span;
+        return {
+            latitude: sign * (90 - theta / DEG),
+            means: [meanM, meanP],
+            psi: sign * this.polarPsi(theta)
+        };
+    }
+    polarPsi(theta) {
+        const s = Math.cos(theta), c = Math.sin(theta), e2 = (1 - this.r) * (1 + this.r), e = Math.sqrt(Math.abs(e2));
+        const I = e === 0 ? s : (e2 > 0 ? Math.atanh(e * s) : Math.atan(e * s)) / e;
+        return Math.asinh(s / c) - e2 * I;
+    }
+    meridianMean(lat1, lat2) {
+        const mid = (lat1 + lat2) / 2;
+        return (integrate(t => {
+            const [s, c] = latitudeSinCos(mid + (t * (lat2 - lat1)) / 2);
+            return this.r ** 2 / Math.hypot(c, this.r * s) ** 3;
+        }, -1, 1) / 2);
+    }
+    meridianDelta(lat1, lat2) {
+        return (lat2 - lat1) * DEG * this.meridianMean(lat1, lat2);
+    }
+    inverseMeridian(start, displacement) {
+        let lo = -90, hi = 90;
+        let phi = start + displacement / this.metrics.meridionalCurvatureRadius(start) / DEG;
+        if (!(phi > lo && phi < hi))
+            phi = (start + Math.sign(displacement) * 90) / 2;
+        for (let i = 0; i < 64; i++) {
+            const residual = this.meridianDelta(start, phi) - displacement;
+            if (residual === 0)
+                return phi;
+            if (residual < 0)
+                lo = phi;
+            else
+                hi = phi;
+            const next = phi - residual / this.metrics.meridionalCurvatureRadius(phi) / DEG;
+            if (next === phi)
+                return phi;
+            phi = next > lo && next < hi ? next : (lo + hi) / 2;
+            if (phi === lo || phi === hi) {
+                // Select a representable non-pole endpoint with the smaller residual.
+                const candidates = [lo, hi].filter(v => Math.abs(v) < 90);
+                return candidates.reduce((best, v) => Math.abs(this.meridianDelta(start, v) - displacement) <
+                    Math.abs(this.meridianDelta(start, best) - displacement)
+                    ? v
+                    : best);
+            }
+        }
+        throw new Error('Rhumb meridian inversion did not converge');
+    }
+    inversePsi(psi) {
+        // Bracketed Newton retains small isometric ordinates instead of bisecting
+        // a fixed 180-degree interval down to a fixed absolute resolution.
+        let lo = -90, hi = 90, phi = Math.atan(Math.sinh(psi)) / DEG;
+        if (Math.abs(phi) === 90)
+            phi = Math.sign(psi) * (90 - 2 ** -46);
+        for (let i = 0; i < 32; i++) {
+            const value = this.psi(phi) - psi;
+            if (value === 0)
+                return phi;
+            if (value < 0)
+                lo = phi;
+            else
+                hi = phi;
+            const [s, c] = latitudeSinCos(phi);
+            const next = phi - (value * c * (c * c + (this.r * s) ** 2)) / this.r ** 2 / DEG;
+            if (next === phi)
+                return phi;
+            phi = next > lo && next < hi ? next : (lo + hi) / 2;
+            if (phi === lo || phi === hi)
+                return phi;
+        }
+        throw new Error('Rhumb isometric inversion did not converge');
+    }
+    potential(latitude) {
+        const [s, c] = latitudeSinCos(latitude), e2 = (1 - this.r) * (1 + this.r), e = Math.sqrt(Math.abs(e2));
+        const I = e === 0 ? s : (e2 > 0 ? Math.atanh(e * s) : Math.atan(e * s)) / e;
+        return (this.r ** 2 * (s / (c * c + (this.r * s) ** 2) + I)) / 2;
+    }
+    area(lat1, lat2, dlon, endpointPsi) {
+        if (dlon === 0 || lat1 === -lat2)
+            return 0;
+        if (lat1 === lat2) {
+            const value = this.a * dlon * (this.a * this.potential(lat1));
+            if (!Number.isFinite(value))
+                throw new RangeError('Unrepresentable rhumb area');
+            return value;
+        }
+        const p1 = this.psi(lat1), p2 = endpointPsi ?? this.psi(lat2), mid = (p1 + p2) / 2, half = (p2 - p1) / 2;
+        // Isometric integration resolves the logarithmic pole approach uniformly.
+        const pieces = Math.max(1, Math.ceil(Math.abs(p2 - p1) / 2));
+        let sum = 0;
+        for (let j = 0; j < pieces; j++)
+            sum += integrate(t => this.potential(this.inversePsi(mid + half * t)), -1 + (2 * j) / pieces, -1 + (2 * (j + 1)) / pieces);
+        const value = this.a * dlon * (this.a * (sum / 2));
+        if (!Number.isFinite(value))
+            throw new RangeError('Unrepresentable rhumb area');
+        return value;
+    }
+}
+/** @internal Shared with RhumbPolygon so even its first vertex rejects exact poles. */
+export function validateRhumbPoint(latitude, longitude) {
+    if (!(Number.isFinite(latitude) && Math.abs(latitude) <= 90 && Number.isFinite(longitude)))
+        throw new RangeError('Rhumb coordinates require finite longitude and latitude in [-90,90]');
+    if (Math.abs(latitude) === 90)
+        throw new RhumbDomainError('pole');
+}
+function latitudeSinCos(phi) {
+    const p = Math.abs(phi);
+    return p <= 45
+        ? [Math.sin(phi * DEG), Math.cos(phi * DEG)]
+        : [Math.sign(phi) * Math.cos((90 - p) * DEG), Math.sin((90 - p) * DEG)];
+}
+function sincosDegrees(angle) {
+    const reduced = angle % 360, q = Math.round(reduced / 90), d = reduced - 90 * q, x = d * DEG;
+    const s = Math.sin(x), c = Math.cos(x);
+    switch (((q % 4) + 4) % 4) {
+        case 0:
+            return [s, c];
+        case 1:
+            return [c, -s];
+        case 2:
+            return [-s, -c];
+        default:
+            return [-c, s];
+    }
+}
+function normalize(angle) {
+    let a = angle % 360;
+    if (a > 180)
+        a -= 360;
+    if (a < -180)
+        a += 360;
+    return a;
+}
+function angleDifference(a, b) {
+    // Retain the error when the first subtraction rounds a tiny dateline gap
+    // to a complete turn. Keep this small helper independent of the runtime.
+    let [d, e] = twoSum(-normalize(a), normalize(b));
+    [d, e] = twoSum(normalize(d), e);
+    if (Math.abs(d) === 180)
+        return (e === 0 ? b < a : e > 0) ? -180 : 180;
+    return d === 0 ? (b < a ? -0 : 0) : d;
+}
+function twoSum(a, b) {
+    const s = a + b, v = s - a;
+    return [s, a - (s - v) + (b - v)];
+}
+//# sourceMappingURL=ellipsoid-rhumb.js.map

package/dist/ellipsoid-rhumb.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-rhumb.js","sourceRoot":"","sources":["../src/ellipsoid-rhumb.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,mEAAmE;AACnE,4DAA4D;AAC5D,8PAA8P;AAC9P,OAAO,EAAC,qBAAqB,EAAC,mDAAgD;AAC9E,OAAO,EAAC,gBAAgB,EAAE,SAAS,EAAC,+BAA4B;AAsBhE,yGAAyG;AACzG,MAAM,OAAO,gBAAiB,SAAQ,UAAU;IAC9C,YAAqB,MAAgC;QACnD,KAAK,CACH,MAAM,KAAK,MAAM;YACf,CAAC,CAAC,uDAAuD;YACzD,CAAC,CAAC,2DAA2D,CAChE,CAAC;QALiB,WAAM,GAAN,MAAM,CAA0B;QAMnD,IAAI,CAAC,IAAI,GAAG,kBAAkB,CAAC;IACjC,CAAC;CACF;AACD,MAAM,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;AAE1B;;;GAGG;AACH,MAAM,OAAO,cAAc;IAMzB,YAAY,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAqB;QAClE,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YAC/D,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;QAC9D,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACf,IAAI,CAAC,OAAO,GAAG,IAAI,gBAAgB,CAAC,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,IAAI,CAAC,CAAC,EAAC,CAAC,CAAC;QAC/E,IAAI,CAAC,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,WAAW,GAAG,CAAC,CAAC;QACpD,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,WAAW,CAAC,IAAI,IAAI,CAAC,WAAW,GAAG,CAAC,CAAC;YAC9D,MAAM,IAAI,UAAU,CAAC,uCAAuC,CAAC,CAAC;QAChE,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IACtB,CAAC;IAED,OAAO,CACL,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,IAAY,EACZ,UAAwB,EAAE;QAE1B,kBAAkB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC/B,kBAAkB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC/B,MAAM,KAAK,GAAG,OAAO,CAAC,cAAc,IAAI,CAAC,CAAC;QAC1C,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC;YAC9B,MAAM,IAAI,UAAU,CAAC,6CAA6C,CAAC,CAAC;QACtE,MAAM,SAAS,GAAG,eAAe,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC;QAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,+BAA+B,CAAC,CAAC;QACvF,MAAM,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;QACjC,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,CAAC,SAAS,GAAG,GAAG,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;QAC7C,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,KAAK,EAAE,EAAE,CAAC,CAAC;QACvD,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,GAAG,EAAE,IAAI,GAAG,KAAK,CAAC,GAAG,GAAG,CAAC;QAC5D,MAAM,MAAM,GAAgB;YAC1B,IAAI;YACJ,IAAI;YACJ,IAAI;YACJ,IAAI,EAAE,OAAO,CAAC,eAAe,CAAC,CAAC,CAAC,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC;YAClE,KAAK,EAAE,GAAG;YACV,GAAG,EAAE,QAAQ;SACd,CAAC;QACF,IAAI,OAAO,CAAC,IAAI;YAAE,MAAM,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,EAAE,SAAS,GAAG,GAAG,CAAC,CAAC;QACtE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,gCAAgC,CAAC,CAAC;QACvF,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,MAAM,CACJ,QAAgB,EAChB,SAAiB,EACjB,OAAe,EACf,QAAgB,EAChB,UAAwB,EAAE;QAE1B,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC7E,CAAC;IAED,IAAI,CACF,QAAgB,EAChB,SAAiB,EACjB,OAAe,EACf,UAAwB,EAAE;QAE1B,kBAAkB,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACxC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,8BAA8B,CAAC,CAAC;QACpF,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,aAAa,CAAC,OAAO,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,OAAO,CAAC,OAAO,CAAC,IAAI,CAAC,EAChC,MAAM,GAAG,OAAO,CAAC,OAAO,CAAC,eAAe,CAAC,CAAC;QAC5C,MAAM,GAAG,GAAG,SAAS,CAAC,OAAO,CAAC,CAAC;QAC/B,OAAO;YACL,QAAQ,EAAE,CAAC,QAAgB,EAAe,EAAE;gBAC1C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;oBAAE,MAAM,IAAI,UAAU,CAAC,+BAA+B,CAAC,CAAC;gBACtF,MAAM,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC;gBAChC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;oBAAE,MAAM,IAAI,UAAU,CAAC,gCAAgC,CAAC,CAAC;gBACpF,MAAM,KAAK,GAAG,KAAK,GAAG,EAAE,CAAC;gBACzB,IAAI,YAAgC,CAAC;gBACrC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;oBAChB,qEAAqE;oBACrE,qEAAqE;oBACrE,MAAM,GAAG,GACP,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,EAAE,QAAQ,CAAC,CAAC;oBACnF,YAAY,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;oBACrC,IAAI,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,GAAG,IAAI,GAAG,GAAG;wBACrC,YAAY,GAAG,qBAAqB,CAClC,IAAI,CAAC,CAAC,EACN,IAAI,CAAC,CAAC,EACN,QAAQ,EACR,OAAO,EACP,QAAQ,EACR,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CACjB,CAAC;oBACJ,IAAI,YAAY,GAAG,CAAC;wBAAE,MAAM,IAAI,gBAAgB,CAAC,eAAe,CAAC,CAAC;oBAClE,IAAI,YAAY,KAAK,CAAC;wBAAE,MAAM,IAAI,gBAAgB,CAAC,MAAM,CAAC,CAAC;gBAC7D,CAAC;gBACD,MAAM,KAAK,GAAG,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,QAAQ,EAAE,KAAK,EAAE,YAAa,CAAC,CAAC;gBAC3F,MAAM,IAAI,GACR,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,QAAQ,IAAI,IAAI,CAAC,eAAe,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC;gBACtF,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,EAAE;oBAAE,MAAM,IAAI,gBAAgB,CAAC,MAAM,CAAC,CAAC;gBAC9D,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,GAAG,KAAK,EAAE,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;gBAClE,MAAM,IAAI,GAAG,CAAC,KAAK,GAAG,EAAE,GAAG,KAAK,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC;gBAChD,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,SAAS,CAAC,SAAS,CAAC,GAAG,IAAI,CAAC,CAAC;gBAChF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;oBAAE,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;gBACpF,MAAM,MAAM,GAAgB;oBAC1B,IAAI,EAAE,QAAQ;oBACd,IAAI,EAAE,SAAS;oBACf,IAAI;oBACJ,IAAI;oBACJ,KAAK,EAAE,GAAG;oBACV,GAAG,EAAE,QAAQ;iBACd,CAAC;gBACF,IAAI,IAAI;oBAAE,MAAM,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,EAAE,IAAI,GAAG,GAAG,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;gBACzE,OAAO,MAAM,CAAC;YAChB,CAAC;SACF,CAAC;IACJ,CAAC;IAEO,KAAK,CAAC,IAAY,EAAE,IAAY;QACtC,MAAM,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;QAC9B,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,YAAY,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxE,+EAA+E;QAC/E,MAAM,KAAK,GACT,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI;YACvB,CAAC,CAAC,SAAS,CACP,CAAC,CAAC,EAAE;gBACF,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC7D,OAAO,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACzD,CAAC,EACD,CAAC,CAAC,EACF,CAAC,CACF,GAAG,CAAC;YACP,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QAC5C,OAAO,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC;IAEO,GAAG,CAAC,QAAgB;QAC1B,OAAO,IAAI,CAAC,OAAO,CAAC,mBAAmB,CAAC,QAAQ,EAAE,WAAW,CAAC,GAAG,GAAG,CAAC;IACvE,CAAC;IAEO,aAAa,CACnB,QAAgB,EAChB,KAAa,EACb,YAAoB;QAEpB,MAAM,IAAI,GACN,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,EAAE,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;YACpF,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC;YACrB,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EACtB,MAAM,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,QAAQ,CAAC,GAAG,GAAG,EACrC,KAAK,GAAG,EAAE,GAAG,GAAG,CAAC;QACnB,MAAM,OAAO,GAAG,CAAC,KAAa,EAAE,EAAE,CAChC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,CAAC;QAC3E,MAAM,GAAG,GAAG,CAAC,KAAa,EAAE,EAAE,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QAC5D,MAAM,MAAM,GACV,IAAI,KAAK,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC;YACvB,CAAC,CAAC,YAAY;YACd,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,IAAI,GAAG,KAAK,CAAC;QAC7F,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC;YAAE,OAAO,SAAS,CAAC;QAC3D,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,KAAK,EACV,KAAK,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;QAC1B,IAAI,CAAC,CAAC,KAAK,GAAG,EAAE,IAAI,KAAK,GAAG,EAAE,CAAC;YAAE,KAAK,GAAG,MAAM,CAAC;QAChD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,QAAQ,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC;YACrC,IAAI,QAAQ,KAAK,CAAC;gBAAE,MAAM;YAC1B,IAAI,QAAQ,GAAG,CAAC;gBAAE,EAAE,GAAG,KAAK,CAAC;;gBACxB,EAAE,GAAG,KAAK,CAAC;YAChB,MAAM,IAAI,GAAG,KAAK,GAAG,QAAQ,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC;YAC/C,IAAI,IAAI,KAAK,KAAK;gBAAE,MAAM;YAC1B,KAAK,GAAG,IAAI,GAAG,EAAE,IAAI,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACtD,IAAI,KAAK,KAAK,EAAE,IAAI,KAAK,KAAK,EAAE;gBAAE,MAAM;YACxC,IAAI,CAAC,KAAK,EAAE;gBAAE,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;QACnF,CAAC;QACD,MAAM,IAAI,GAAG,KAAK,GAAG,MAAM,EACzB,GAAG,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QACvE,MAAM,KAAK,GACT,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;YACjE,CAAC,CAAC,SAAS,CACP,CAAC,CAAC,EAAE;gBACF,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAC5B,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EACf,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAClB,OAAO,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACzD,CAAC,EACD,CAAC,CAAC,EACF,CAAC,CACF,GAAG,CAAC;YACP,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC;QAC7D,OAAO;YACL,QAAQ,EAAE,IAAI,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,GAAG,CAAC;YACnC,KAAK,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC;YACrB,GAAG,EAAE,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC;SACjC,CAAC;IACJ,CAAC;IAEO,QAAQ,CAAC,KAAa;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EACvB,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EACnB,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAChC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC5E,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACpC,CAAC;IAEO,YAAY,CAAC,IAAY,EAAE,IAAY;QAC7C,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC9B,OAAO,CACL,SAAS,CACP,CAAC,CAAC,EAAE;YACF,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC7D,OAAO,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;QACtD,CAAC,EACD,CAAC,CAAC,EACF,CAAC,CACF,GAAG,CAAC,CACN,CAAC;IACJ,CAAC;IAEO,aAAa,CAAC,IAAY,EAAE,IAAY;QAC9C,OAAO,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,YAAY,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IAC7D,CAAC;IAEO,eAAe,CAAC,KAAa,EAAE,YAAoB;QACzD,IAAI,EAAE,GAAG,CAAC,EAAE,EACV,EAAE,GAAG,EAAE,CAAC;QACV,IAAI,GAAG,GAAG,KAAK,GAAG,YAAY,GAAG,IAAI,CAAC,OAAO,CAAC,yBAAyB,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QACrF,IAAI,CAAC,CAAC,GAAG,GAAG,EAAE,IAAI,GAAG,GAAG,EAAE,CAAC;YAAE,GAAG,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC9E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,QAAQ,GAAG,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,YAAY,CAAC;YAC/D,IAAI,QAAQ,KAAK,CAAC;gBAAE,OAAO,GAAG,CAAC;YAC/B,IAAI,QAAQ,GAAG,CAAC;gBAAE,EAAE,GAAG,GAAG,CAAC;;gBACtB,EAAE,GAAG,GAAG,CAAC;YACd,MAAM,IAAI,GAAG,GAAG,GAAG,QAAQ,GAAG,IAAI,CAAC,OAAO,CAAC,yBAAyB,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;YAChF,IAAI,IAAI,KAAK,GAAG;gBAAE,OAAO,GAAG,CAAC;YAC7B,GAAG,GAAG,IAAI,GAAG,EAAE,IAAI,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACpD,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE,EAAE,CAAC;gBAC7B,sEAAsE;gBACtE,MAAM,UAAU,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;gBAC1D,OAAO,UAAU,CAAC,MAAM,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CACnC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,YAAY,CAAC;oBACrD,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,YAAY,CAAC;oBACtD,CAAC,CAAC,CAAC;oBACH,CAAC,CAAC,IAAI,CACT,CAAC;YACJ,CAAC;QACH,CAAC;QACD,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,CAAC;IAEO,UAAU,CAAC,GAAW;QAC5B,0EAA0E;QAC1E,mEAAmE;QACnE,IAAI,EAAE,GAAG,CAAC,EAAE,EACV,EAAE,GAAG,EAAE,EACP,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACxC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,EAAE;YAAE,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;YAClC,IAAI,KAAK,KAAK,CAAC;gBAAE,OAAO,GAAG,CAAC;YAC5B,IAAI,KAAK,GAAG,CAAC;gBAAE,EAAE,GAAG,GAAG,CAAC;;gBACnB,EAAE,GAAG,GAAG,CAAC;YACd,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,GAAG,CAAC,CAAC;YACnC,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;YACjF,IAAI,IAAI,KAAK,GAAG;gBAAE,OAAO,GAAG,CAAC;YAC7B,GAAG,GAAG,IAAI,GAAG,EAAE,IAAI,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACpD,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE;gBAAE,OAAO,GAAG,CAAC;QAC3C,CAAC;QACD,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;IAChE,CAAC;IAEO,SAAS,CAAC,QAAgB;QAChC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,QAAQ,CAAC,EACrC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAChC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC5E,OAAO,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACnE,CAAC;IAEO,IAAI,CAAC,IAAY,EAAE,IAAY,EAAE,IAAY,EAAE,WAAoB;QACzE,IAAI,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,IAAI;YAAE,OAAO,CAAC,CAAC;QAC3C,IAAI,IAAI,KAAK,IAAI,EAAE,CAAC;YAClB,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC;YAC9D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;YAChF,OAAO,KAAK,CAAC;QACf,CAAC;QACD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EACvB,EAAE,GAAG,WAAW,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAClC,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EACnB,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QACvB,0EAA0E;QAC1E,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7D,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;YAC7B,GAAG,IAAI,SAAS,CACd,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,UAAU,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,EACpD,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,EACrB,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAC5B,CAAC;QACJ,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACnD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;QAChF,OAAO,KAAK,CAAC;IACf,CAAC;CACF;AAED,uFAAuF;AACvF,MAAM,UAAU,kBAAkB,CAAC,QAAgB,EAAE,SAAiB;IACpE,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,EAAE,IAAI,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC;QACxF,MAAM,IAAI,UAAU,CAAC,qEAAqE,CAAC,CAAC;IAC9F,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE;QAAE,MAAM,IAAI,gBAAgB,CAAC,MAAM,CAAC,CAAC;AACpE,CAAC;AACD,SAAS,cAAc,CAAC,GAAW;IACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxB,OAAO,CAAC,IAAI,EAAE;QACZ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAC5C,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AAC5E,CAAC;AACD,SAAS,aAAa,CAAC,KAAa;IAClC,MAAM,OAAO,GAAG,KAAK,GAAG,GAAG,EACzB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,EAAE,CAAC,EAC5B,CAAC,GAAG,OAAO,GAAG,EAAE,GAAG,CAAC,EACpB,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;IACd,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAClB,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QAC1B,KAAK,CAAC;YACJ,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChB,KAAK,CAAC;YACJ,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjB,KAAK,CAAC;YACJ,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAClB;YACE,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACnB,CAAC;AACH,CAAC;AACD,SAAS,SAAS,CAAC,KAAa;IAC9B,IAAI,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC;IACpB,IAAI,CAAC,GAAG,GAAG;QAAE,CAAC,IAAI,GAAG,CAAC;IACtB,IAAI,CAAC,GAAG,CAAC,GAAG;QAAE,CAAC,IAAI,GAAG,CAAC;IACvB,OAAO,CAAC,CAAC;AACX,CAAC;AACD,SAAS,eAAe,CAAC,CAAS,EAAE,CAAS;IAC3C,yEAAyE;IACzE,yEAAyE;IACzE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG;QAAE,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IACvE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACxC,CAAC;AACD,SAAS,MAAM,CAAC,CAAS,EAAE,CAAS;IAClC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EACb,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACZ,OAAO,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACpC,CAAC"}

package/dist/ellipsoid-tangent-plane.d.ts

@@ -1,28 +1,50 @@
 import { Vector2, Vector3 } from '@math.gl/core';
+import type { NumericArray } from '@math.gl/core/types';
 import { Plane } from '@math.gl/geometry/culling';
 import { Ellipsoid } from "./ellipsoid.js";
+export type TangentPlaneProjectionMode = 'normal' | 'radial';
+export type TangentPlaneProjection = {
+    positions: Float64Array;
+    /** One byte per input point: 1 for success, 0 for no radial intersection. */
+    valid: Uint8Array;
+};
 /** A two-dimensional east-north plane tangent to an ellipsoid. */
 export declare class EllipsoidTangentPlane {
-    private _origin;
-    private _xAxis;
-    private _yAxis;
+    private anchor;
+    private east;
+    private north;
     private _plane;
+    private readonly radii;
+    private readonly radius?;
     /**
-     * Creates a new plane tangent to the WGS84 ellipsoid at the provided origin.
+     * Creates a new plane tangent to an ellipsoid at the provided origin.
      * If origin is not on the surface of the ellipsoid, its surface projection will be used.
      *
      * @param origin A nonzero Cartesian point on or near the ellipsoid.
      * @param ellipsoid The ellipsoid to which the origin belongs. Defaults to WGS84.
      */
-    constructor(origin: number[], ellipsoid?: Ellipsoid);
+    constructor(origin: Readonly<NumericArray>, ellipsoid?: Ellipsoid);
     /**
      * Computes the projection of the provided 3D position onto the 2D plane, along the plane normal.
      *
-     * @param cartesian The WGS84 Cartesian point to project.
+     * @param cartesian The fixed Cartesian point to project.
      * @param result The optional object in which to store the local east-north coordinates.
      * @returns The supplied result or a new Vector2.
      */
-    projectPointToNearestOnPlane(cartesian: Vector3, result?: Vector2): Vector2;
+    projectPointToNearestOnPlane(cartesian: Readonly<NumericArray>, result?: Vector2): Vector2;
+    /** Project along the outward ray from the ellipsoid center.
+     * Back-facing, horizon and nonfinite intersections return undefined without writing result.
+     */
+    projectPointRadiallyOntoPlane(cartesian: Readonly<NumericArray>, result?: Vector2): Vector2 | undefined;
+    /** Lift a local east/north point to the ellipsoid. Radial is the inverse of radial
+     * projection. Normal selects the nearest-facing line intersection and can miss.
+     * Neither operation restores a source point's height. Failure leaves result untouched.
+     */
+    projectPointOntoEllipsoid(point: Readonly<NumericArray>, mode?: TangentPlaneProjectionMode, result?: Vector3): Vector3 | undefined;
+    /** Packed XYZ to packed east/north. Invalid radial points produce NaN pairs and a zero mask.
+     * Optional output must have exactly two components per point and must not overlap input.
+     */
+    projectPointsFlat(input: Readonly<NumericArray>, mode?: TangentPlaneProjectionMode, output?: Float64Array): TangentPlaneProjection;
     get plane(): Plane;
     get origin(): Vector3;
     get xAxis(): Vector3;

package/dist/ellipsoid-tangent-plane.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ellipsoid-tangent-plane.d.ts","sourceRoot":"","sources":["../src/ellipsoid-tangent-plane.ts"],"names":[],"mappings":"AASA,OAAO,EAAC,OAAO,EAAE,OAAO,EAAU,MAAM,eAAe,CAAC;AACxD,OAAO,EAAC,KAAK,EAAC,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAC,SAAS,EAAC,uBAAoB;AAOtC,kEAAkE;AAClE,qBAAa,qBAAqB;IAChC,OAAO,CAAC,OAAO,CAAU;IACzB,OAAO,CAAC,MAAM,CAAU;IACxB,OAAO,CAAC,MAAM,CAAU;IACxB,OAAO,CAAC,MAAM,CAAQ;IAEtB;;;;;;OAMG;gBACS,MAAM,EAAE,MAAM,EAAE,EAAE,SAAS,GAAE,SAA2B;IAapE;;;;;;OAMG;IACH,4BAA4B,CAAC,SAAS,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAiB3E,IAAI,KAAK,IAAI,KAAK,CAEjB;IAED,IAAI,MAAM,IAAI,OAAO,CAEpB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;CACF"}
+{"version":3,"file":"ellipsoid-tangent-plane.d.ts","sourceRoot":"","sources":["../src/ellipsoid-tangent-plane.ts"],"names":[],"mappings":"AASA,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAC/C,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,qBAAqB,CAAC;AACtD,OAAO,EAAC,KAAK,EAAC,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAC,SAAS,EAAC,uBAAoB;AAGtC,MAAM,MAAM,0BAA0B,GAAG,QAAQ,GAAG,QAAQ,CAAC;AAC7D,MAAM,MAAM,sBAAsB,GAAG;IACnC,SAAS,EAAE,YAAY,CAAC;IACxB,6EAA6E;IAC7E,KAAK,EAAE,UAAU,CAAC;CACnB,CAAC;AAEF,kEAAkE;AAClE,qBAAa,qBAAqB;IAChC,OAAO,CAAC,MAAM,CAAU;IACxB,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,MAAM,CAAQ;IACtB,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAoB;IAC1C,OAAO,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAS;IAEjC;;;;;;OAMG;gBACS,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,SAAS,GAAE,SAA2B;IA4BlF;;;;;;OAMG;IACH,4BAA4B,CAAC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,UAAgB,GAAG,OAAO;IAWhG;;OAEG;IACH,6BAA6B,CAC3B,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EACjC,MAAM,UAAgB,GACrB,OAAO,GAAG,SAAS;IAyBtB;;;OAGG;IACH,yBAAyB,CACvB,KAAK,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC7B,IAAI,GAAE,0BAAqC,EAC3C,MAAM,UAAgB,GACrB,OAAO,GAAG,SAAS;IAgEtB;;OAEG;IACH,iBAAiB,CACf,KAAK,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC7B,IAAI,GAAE,0BAAqC,EAC3C,MAAM,CAAC,EAAE,YAAY,GACpB,sBAAsB;IAkCzB,IAAI,KAAK,IAAI,KAAK,CAEjB;IAED,IAAI,MAAM,IAAI,OAAO,CAEpB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;IAED,IAAI,KAAK,IAAI,OAAO,CAEnB;CACF"}