@math.gl/geospatial 5.0.0-alpha.8 → 5.0.0-beta.1

package/dist/projection/experimental/kernels/poly.d.ts

@@ -0,0 +1,10 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    ml0: number;
+    meridianPole: number;
+    poleRadius: number;
+};
+export declare function forward(state: State, p: Point): Point | null | undefined | number;
+export declare function inverse(state: State, p: Point): Point | null | undefined | number;
+export declare function createState(base: KernelParameters, options?: Partial<State>): State;
+//# sourceMappingURL=poly.d.ts.map

package/dist/projection/experimental/kernels/poly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"poly.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/poly.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,GAAG,EAAE,MAAM,CAAC;IACZ,YAAY,EAAE,MAAM,CAAC;IACrB,UAAU,EAAE,MAAM,CAAC;CACpB,CAAC;AAmBF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CA4BjF;AAID,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAyFjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAUvF"}

package/dist/projection/experimental/kernels/poly.js

@@ -0,0 +1,193 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original math.gl polar inverse tangent-plane limit and shared meridional quadrature remove cancellation, truncation and inverse residuals.  Adapted from proj4js 2.22.0 with corrected spherical meridian distance across pole-relative latitudes. See ../../../../projection/PROJ4-LICENSE.md for upstream attribution. Independent corrections acknowledge PROJ 9.5.1 src/projections/poly.cpp; see ../../../../projection/PROJ-LICENSE.txt.
+import adjust_lon from "../common/adjust_lon.js";
+import { EPSLN } from "../common/constants.js";
+import gN from "../common/gN.js";
+import { meridianDistance } from "../common/meridian-distance.js";
+var MAX_ITER = 20;
+function initialize(state) {
+    /* Place parameters in static storage for common use
+        ------------------------------------------------- */
+    state.meridianPole = meridianDistance(state.a, state.es, 0, Math.PI / 2);
+    state.poleRadius = state.a / Math.sqrt(1 - state.es);
+    state.ml0 = meridianDistance(state.a, state.es, 0, state.lat0);
+}
+/* Polyconic forward equations--mapping lat,long to x,y
+    --------------------------------------------------- */
+export function forward(state, p) {
+    var lon = p.x;
+    var lat = p.y;
+    var x, y, el;
+    var dlon = adjust_lon(lon - state.long0, state.over);
+    el = dlon * Math.sin(lat);
+    if (state.sphere) {
+        if (lat === 0) {
+            x = state.a * dlon;
+            y = -1 * state.a * state.lat0;
+        }
+        else {
+            x = (state.a * Math.sin(el)) / Math.tan(lat);
+            // This is a meridian distance, not a latitude to wrap at either pole.
+            y = state.a * (lat - state.lat0 + (2 * Math.sin(el / 2) ** 2) / Math.tan(lat));
+        }
+    }
+    else {
+        if (lat === 0) {
+            x = state.a * dlon;
+            y = -1 * state.ml0;
+        }
+        else {
+            var nl = gN(state.a, state.e, Math.sin(lat)) / Math.tan(lat);
+            x = nl * Math.sin(el);
+            y = meridianDistance(state.a, state.es, state.lat0, lat) + nl * (2 * Math.sin(el / 2) ** 2);
+        }
+    }
+    p.x = x + state.x0;
+    p.y = y + state.y0;
+    return p;
+}
+/* Inverse equations
+  ----------------- */
+export function inverse(state, p) {
+    var lon, lat, x, y, i;
+    var al, bl;
+    var phi, dphi;
+    x = p.x - state.x0;
+    y = p.y - state.y0;
+    const meridian = y + (state.sphere ? state.a * state.lat0 : state.ml0);
+    const sign = meridian < 0 ? -1 : 1;
+    const pole = state.sphere ? (state.a * Math.PI) / 2 : state.meridianPole;
+    const poleOffset = pole - sign * meridian;
+    const radius = Math.hypot(x, poleOffset);
+    // The polar limit is a tangent plane: x=N*u*sin(lambda),
+    // y_pole-y=sign*N*u*cos(lambda). Neglected terms are O(a*u^3),
+    // below a nanometre for this branch, before binary64 rounding.
+    if (radius / state.a < 1e-6) {
+        p.x = adjust_lon(state.long0 + Math.atan2(x, poleOffset), state.over);
+        p.y = sign * (Math.PI / 2 - radius / state.poleRadius);
+        return p;
+    }
+    if (state.sphere) {
+        if (Math.abs(y + state.a * state.lat0) <= EPSLN) {
+            lon = adjust_lon(x / state.a + state.long0, state.over);
+            lat = 0;
+        }
+        else {
+            al = state.lat0 + y / state.a;
+            bl = (x * x) / state.a / state.a + al * al;
+            phi = al;
+            var tanphi;
+            for (i = MAX_ITER; i; --i) {
+                tanphi = Math.tan(phi);
+                dphi =
+                    (-1 * (al * (phi * tanphi + 1) - phi - 0.5 * (phi * phi + bl) * tanphi)) /
+                        ((phi - al) / tanphi - 1);
+                phi += dphi;
+                if (Math.abs(dphi) <= 1e-14) {
+                    lat = phi;
+                    break;
+                }
+            }
+            lat = refineInverseLatitude(state, lat === undefined ? phi : lat, x, y);
+            if (lat === undefined)
+                return null;
+            phi = lat;
+            lon = adjust_lon(state.long0 + Math.asin((x * Math.tan(phi)) / state.a) / Math.sin(lat), state.over);
+        }
+    }
+    else {
+        if (Math.abs(y + state.ml0) <= EPSLN) {
+            lat = 0;
+            lon = adjust_lon(state.long0 + x / state.a, state.over);
+        }
+        else {
+            al = (state.ml0 + y) / state.a;
+            bl = (x * x) / state.a / state.a + al * al;
+            phi = al;
+            var cl, mln, mlnp, ma;
+            var con;
+            for (i = MAX_ITER; i; --i) {
+                con = state.e * Math.sin(phi);
+                cl = Math.sqrt(1 - con * con) * Math.tan(phi);
+                mln = meridianDistance(state.a, state.es, 0, phi);
+                mlnp = (1 - state.es) / Math.pow(1 - con * con, 1.5);
+                ma = mln / state.a;
+                dphi =
+                    (al * (cl * ma + 1) - ma - 0.5 * cl * (ma * ma + bl)) /
+                        ((state.es * Math.sin(2 * phi) * (ma * ma + bl - 2 * al * ma)) / (4 * cl) +
+                            (al - ma) * (cl * mlnp - 2 / Math.sin(2 * phi)) -
+                            mlnp);
+                phi -= dphi;
+                if (Math.abs(dphi) <= 1e-14) {
+                    lat = phi;
+                    break;
+                }
+            }
+            lat = refineInverseLatitude(state, lat === undefined ? phi : lat, x, y);
+            if (lat === undefined)
+                return null;
+            phi = lat;
+            // lat=phi4z(state.e,state.e0,state.e1,state.e2,state.e3,al,bl,0,0);
+            cl = Math.sqrt(1 - state.es * Math.pow(Math.sin(lat), 2)) * Math.tan(lat);
+            lon = adjust_lon(state.long0 + Math.asin((x * cl) / state.a) / Math.sin(lat), state.over);
+        }
+    }
+    p.x = lon;
+    p.y = lat;
+    return p;
+}
+export function createState(base, options = {}) {
+    const state = {
+        ...base,
+        ml0: 0,
+        meridianPole: 0,
+        poleRadius: 0,
+        ...options
+    };
+    initialize(state);
+    return state;
+}
+/** Refine a stalled iterate using the original northing at fixed easting.
+ * Acceptance is measured in planar axis units, not an angular step threshold.
+ */
+function refineInverseLatitude(state, latitude, x, y) {
+    for (let i = 0; i < 8; i++) {
+        if (!Number.isFinite(latitude) || Math.abs(latitude) >= Math.PI / 2)
+            return undefined;
+        const sin = Math.sin(latitude), cos = Math.cos(latitude);
+        const denominator = 1 - state.es * sin * sin;
+        const normal = state.sphere ? state.a : gN(state.a, state.e, sin);
+        const radius = normal / Math.tan(latitude);
+        const arc = Math.asin(x / radius);
+        const easting = radius * Math.sin(arc);
+        const oneMinusCos = 2 * Math.sin(arc / 2) ** 2;
+        const meridian = state.sphere
+            ? state.a * (latitude - state.lat0)
+            : meridianDistance(state.a, state.es, state.lat0, latitude);
+        const correction = radius * oneMinusCos;
+        const northing = meridian + correction;
+        const error = northing - y;
+        // Scale each ordinate by the terms forming it. A full-axis tolerance on
+        // tiny equatorial northings would accept a large relative latitude error.
+        if (Math.abs(easting - x) <= 64 * Number.EPSILON * Math.max(state.a, Math.abs(x)) &&
+            Math.abs(error) <=
+                64 *
+                    Number.EPSILON *
+                    Math.max(Math.abs(y), Math.abs(meridian), Math.abs(correction), state.a * Math.max(Math.abs(latitude), Math.abs(state.lat0))))
+            return latitude;
+        const meridianDerivative = state.sphere
+            ? state.a
+            : (state.a * (1 - state.es)) / Math.pow(denominator, 1.5);
+        const radiusDerivative = normal * ((state.sphere ? 0 : (state.es * cos * cos) / denominator) - 1 / (sin * sin));
+        const derivative = meridianDerivative + radiusDerivative * (-oneMinusCos / Math.cos(arc));
+        const step = error / derivative;
+        if (!Number.isFinite(step) || Math.abs(step) > 1e-6)
+            return undefined;
+        latitude -= step;
+    }
+    return undefined;
+}
+//# sourceMappingURL=poly.js.map

package/dist/projection/experimental/kernels/poly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"poly.js","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/poly.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,mcAAmc;AAQnc,OAAO,UAAU,gCAA6B;AAC9C,OAAO,EAAC,KAAK,EAAC,+BAA4B;AAE1C,OAAO,EAAE,wBAAqB;AAC9B,OAAO,EAAC,gBAAgB,EAAC,uCAAoC;AAE7D,IAAI,QAAQ,GAAG,EAAE,CAAC;AAElB,SAAS,UAAU,CAAC,KAAY;IAC9B;4DACwD;IACxD,KAAK,CAAC,YAAY,GAAG,gBAAgB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACzE,KAAK,CAAC,UAAU,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;IACrD,KAAK,CAAC,GAAG,GAAG,gBAAgB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;AACjE,CAAC;AAED;0DAC0D;AAC1D,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;IACb,IAAI,IAAI,GAAG,UAAU,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACrD,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC1B,IAAI,KAAK,CAAC,MAAM,EAAE,CAAC;QACjB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;YACd,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC;YACnB,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC;QAChC,CAAC;aAAM,CAAC;YACN,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC7C,sEAAsE;YACtE,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACjF,CAAC;IACH,CAAC;SAAM,CAAC;QACN,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;YACd,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC;YACnB,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC;QACrB,CAAC;aAAM,CAAC;YACN,IAAI,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC7D,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;YACtB,CAAC,GAAG,gBAAgB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;QAC9F,CAAC;IACH,CAAC;IACD,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IACnB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IACnB,OAAO,CAAC,CAAC;AACX,CAAC;AAED;sBACsB;AACtB,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACtB,IAAI,EAAE,EAAE,EAAE,CAAC;IACX,IAAI,GAAG,EAAE,IAAI,CAAC;IACd,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IACnB,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAEnB,MAAM,QAAQ,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACvE,MAAM,IAAI,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,IAAI,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,YAAY,CAAC;IACzE,MAAM,UAAU,GAAG,IAAI,GAAG,IAAI,GAAG,QAAQ,CAAC;IAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IACzC,yDAAyD;IACzD,+DAA+D;IAC/D,+DAA+D;IAC/D,IAAI,MAAM,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;QAC5B,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,UAAU,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACtE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;QACvD,OAAO,CAAC,CAAC;IACX,CAAC;IAED,IAAI,KAAK,CAAC,MAAM,EAAE,CAAC;QACjB,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;YAChD,GAAG,GAAG,UAAU,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;YACxD,GAAG,GAAG,CAAC,CAAC;QACV,CAAC;aAAM,CAAC;YACN,EAAE,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YAC9B,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;YAC3C,GAAG,GAAG,EAAE,CAAC;YACT,IAAI,MAAM,CAAC;YACX,KAAK,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;gBAC1B,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACvB,IAAI;oBACF,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;wBACxE,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC;gBAC5B,GAAG,IAAI,IAAI,CAAC;gBACZ,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;oBAC5B,GAAG,GAAG,GAAG,CAAC;oBACV,MAAM;gBACR,CAAC;YACH,CAAC;YACD,GAAG,GAAG,qBAAqB,CAAC,KAAK,EAAE,GAAG,KAAK,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACxE,IAAI,GAAG,KAAK,SAAS;gBAAE,OAAO,IAAI,CAAC;YACnC,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,GAAG,UAAU,CACd,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EACtE,KAAK,CAAC,IAAI,CACX,CAAC;QACJ,CAAC;IACH,CAAC;SAAM,CAAC;QACN,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,KAAK,EAAE,CAAC;YACrC,GAAG,GAAG,CAAC,CAAC;YACR,GAAG,GAAG,UAAU,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAC1D,CAAC;aAAM,CAAC;YACN,EAAE,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YAC/B,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;YAC3C,GAAG,GAAG,EAAE,CAAC;YACT,IAAI,EAAE,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;YACtB,IAAI,GAAG,CAAC;YACR,KAAK,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;gBAC1B,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC9B,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC9C,GAAG,GAAG,gBAAgB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;gBAClD,IAAI,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC;gBACrD,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;gBACnB,IAAI;oBACF,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;wBACrD,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;4BACvE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;4BAC/C,IAAI,CAAC,CAAC;gBACV,GAAG,IAAI,IAAI,CAAC;gBACZ,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;oBAC5B,GAAG,GAAG,GAAG,CAAC;oBACV,MAAM;gBACR,CAAC;YACH,CAAC;YAED,GAAG,GAAG,qBAAqB,CAAC,KAAK,EAAE,GAAG,KAAK,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACxE,IAAI,GAAG,KAAK,SAAS;gBAAE,OAAO,IAAI,CAAC;YACnC,GAAG,GAAG,GAAG,CAAC;YACV,oEAAoE;YACpE,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC1E,GAAG,GAAG,UAAU,CAAC,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAC5F,CAAC;IACH,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU;QACnB,GAAG,IAAI;QACP,GAAG,EAAE,CAAC;QACN,YAAY,EAAE,CAAC;QACf,UAAU,EAAE,CAAC;QACb,GAAG,OAAO;KACX,CAAC;IACF,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;GAEG;AACH,SAAS,qBAAqB,CAC5B,KAAY,EACZ,QAAgB,EAChB,CAAS,EACT,CAAS;IAET,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,EAAE,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QACtF,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;QAC3B,MAAM,WAAW,GAAG,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;QAC7C,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAClE,MAAM,MAAM,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;QAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QAClC,MAAM,OAAO,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvC,MAAM,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;QAC/C,MAAM,QAAQ,GAAG,KAAK,CAAC,MAAM;YAC3B,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC;YACnC,CAAC,CAAC,gBAAgB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAC9D,MAAM,UAAU,GAAG,MAAM,GAAG,WAAW,CAAC;QACxC,MAAM,QAAQ,GAAG,QAAQ,GAAG,UAAU,CAAC;QACvC,MAAM,KAAK,GAAG,QAAQ,GAAG,CAAC,CAAC;QAC3B,wEAAwE;QACxE,0EAA0E;QAC1E,IACE,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,MAAM,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7E,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;gBACb,EAAE;oBACA,MAAM,CAAC,OAAO;oBACd,IAAI,CAAC,GAAG,CACN,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EACX,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAClB,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,EACpB,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAC7D;YAEL,OAAO,QAAQ,CAAC;QAClB,MAAM,kBAAkB,GAAG,KAAK,CAAC,MAAM;YACrC,CAAC,CAAC,KAAK,CAAC,CAAC;YACT,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,EAAE,GAAG,CAAC,CAAC;QAC5D,MAAM,gBAAgB,GACpB,MAAM,GAAG,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QACzF,MAAM,UAAU,GAAG,kBAAkB,GAAG,gBAAgB,GAAG,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1F,MAAM,IAAI,GAAG,KAAK,GAAG,UAAU,CAAC;QAChC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI;YAAE,OAAO,SAAS,CAAC;QACtE,QAAQ,IAAI,IAAI,CAAC;IACnB,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/dist/projection/experimental/kernels/qsc.d.ts

@@ -0,0 +1,10 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    face: number;
+    one_minus_f: number;
+    one_minus_f_squared: number;
+};
+export declare function forward(state: State, p: Point): Point | null | undefined | number;
+export declare function inverse(state: State, p: Point): Point | null | undefined | number;
+export declare function createState(base: KernelParameters, options?: Partial<State>): State;
+//# sourceMappingURL=qsc.d.ts.map

package/dist/projection/experimental/kernels/qsc.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"qsc.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/qsc.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,IAAI,EAAE,MAAM,CAAC;IACb,WAAW,EAAE,MAAM,CAAC;IACpB,mBAAmB,EAAE,MAAM,CAAC;CAC7B,CAAC;AAqDF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAuIjF;AAID,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAqIjF;AAkBD,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAIvF"}

package/dist/projection/experimental/kernels/qsc.js

@@ -0,0 +1,364 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original math.gl half-angle and atan2 stability corrections. Direct TypeScript port of proj4js 2.22.0, modified to reuse the caller-owned point and represent area/temporary coordinates as local numbers. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
+// QSC projection rewritten from the original PROJ4
+// https://github.com/OSGeo/proj.4/blob/master/src/PJ_qsc.c
+import { TWO_PI, HALF_PI, FORTPI } from "../common/constants.js";
+/* constants */
+var FACE_ENUM = {
+    FRONT: 1,
+    RIGHT: 2,
+    BACK: 3,
+    LEFT: 4,
+    TOP: 5,
+    BOTTOM: 6
+};
+var AREA_ENUM = {
+    AREA_0: 1,
+    AREA_1: 2,
+    AREA_2: 3,
+    AREA_3: 4
+};
+function initialize(state) {
+    state.x0 = state.x0 || 0;
+    state.y0 = state.y0 || 0;
+    state.lat0 = state.lat0 || 0;
+    state.long0 = state.long0 || 0;
+    state.lat_ts = state.lat_ts || 0;
+    /* Determine the cube face from the center of projection. */
+    if (state.lat0 >= HALF_PI - FORTPI / 2.0) {
+        state.face = FACE_ENUM.TOP;
+    }
+    else if (state.lat0 <= -(HALF_PI - FORTPI / 2.0)) {
+        state.face = FACE_ENUM.BOTTOM;
+    }
+    else if (Math.abs(state.long0) <= FORTPI) {
+        state.face = FACE_ENUM.FRONT;
+    }
+    else if (Math.abs(state.long0) <= HALF_PI + FORTPI) {
+        state.face = state.long0 > 0.0 ? FACE_ENUM.RIGHT : FACE_ENUM.LEFT;
+    }
+    else {
+        state.face = FACE_ENUM.BACK;
+    }
+    /* Fill in useful values for the ellipsoid <-> sphere shift
+     * described in [LK12]. */
+    if (state.es !== 0) {
+        state.one_minus_f = 1 - (state.a - state.b) / state.a;
+        state.one_minus_f_squared = state.one_minus_f * state.one_minus_f;
+    }
+}
+// QSC forward equations--mapping lat,long to x,y
+// -----------------------------------------------------------------
+export function forward(state, p) {
+    var xy = p;
+    var lat, lon;
+    var theta, phi;
+    var t, mu;
+    /* nu; */
+    var area = 0;
+    // move lon according to projection's lon
+    p.x -= state.long0;
+    /* Convert the geodetic latitude to a geocentric latitude.
+     * This corresponds to the shift from the ellipsoid to the sphere
+     * described in [LK12]. */
+    if (state.es !== 0) {
+        // if (P->es != 0) {
+        lat = Math.atan(state.one_minus_f_squared * Math.tan(p.y));
+    }
+    else {
+        lat = p.y;
+    }
+    /* Convert the input lat, lon into theta, phi as used by QSC.
+     * This depends on the cube face and the area on it.
+     * For the top and bottom face, we can compute theta and phi
+     * directly from phi, lam. For the other faces, we must use
+     * unit sphere cartesian coordinates as an intermediate step. */
+    lon = p.x; // lon = lam;
+    if (state.face === FACE_ENUM.TOP) {
+        phi = HALF_PI - lat;
+        if (lon >= FORTPI && lon <= HALF_PI + FORTPI) {
+            area = AREA_ENUM.AREA_0;
+            theta = lon - HALF_PI;
+        }
+        else if (lon > HALF_PI + FORTPI || lon <= -(HALF_PI + FORTPI)) {
+            area = AREA_ENUM.AREA_1;
+            theta = lon > 0.0 ? lon - Math.PI : lon + Math.PI;
+        }
+        else if (lon > -(HALF_PI + FORTPI) && lon <= -FORTPI) {
+            area = AREA_ENUM.AREA_2;
+            theta = lon + HALF_PI;
+        }
+        else {
+            area = AREA_ENUM.AREA_3;
+            theta = lon;
+        }
+    }
+    else if (state.face === FACE_ENUM.BOTTOM) {
+        phi = HALF_PI + lat;
+        if (lon >= FORTPI && lon <= HALF_PI + FORTPI) {
+            area = AREA_ENUM.AREA_0;
+            theta = -lon + HALF_PI;
+        }
+        else if (lon < FORTPI && lon >= -FORTPI) {
+            area = AREA_ENUM.AREA_1;
+            theta = -lon;
+        }
+        else if (lon < -FORTPI && lon >= -(HALF_PI + FORTPI)) {
+            area = AREA_ENUM.AREA_2;
+            theta = -lon - HALF_PI;
+        }
+        else {
+            area = AREA_ENUM.AREA_3;
+            theta = lon > 0.0 ? -lon + Math.PI : -lon - Math.PI;
+        }
+    }
+    else {
+        var q, r, s;
+        var sinlat, coslat;
+        var sinlon, coslon;
+        if (state.face === FACE_ENUM.RIGHT) {
+            lon = qsc_shift_lon_origin(lon, +HALF_PI);
+        }
+        else if (state.face === FACE_ENUM.BACK) {
+            lon = qsc_shift_lon_origin(lon, +Math.PI);
+        }
+        else if (state.face === FACE_ENUM.LEFT) {
+            lon = qsc_shift_lon_origin(lon, -HALF_PI);
+        }
+        sinlat = Math.sin(lat);
+        coslat = Math.cos(lat);
+        sinlon = Math.sin(lon);
+        coslon = Math.cos(lon);
+        q = coslat * coslon;
+        r = coslat * sinlon;
+        s = sinlat;
+        if (state.face === FACE_ENUM.FRONT) {
+            phi = Math.atan2(Math.hypot(r, s), q);
+            theta = qsc_fwd_equat_face_theta(phi, s, r);
+        }
+        else if (state.face === FACE_ENUM.RIGHT) {
+            phi = Math.atan2(Math.hypot(q, s), r);
+            theta = qsc_fwd_equat_face_theta(phi, s, -q);
+        }
+        else if (state.face === FACE_ENUM.BACK) {
+            phi = Math.atan2(Math.hypot(r, s), -q);
+            theta = qsc_fwd_equat_face_theta(phi, s, -r);
+        }
+        else if (state.face === FACE_ENUM.LEFT) {
+            phi = Math.atan2(Math.hypot(q, s), -r);
+            theta = qsc_fwd_equat_face_theta(phi, s, q);
+        }
+        else {
+            /* Impossible */
+            phi = theta = 0;
+        }
+        if (Math.abs(theta) <= FORTPI) {
+            area = AREA_ENUM.AREA_0;
+        }
+        else if (theta > FORTPI && theta <= HALF_PI + FORTPI) {
+            area = AREA_ENUM.AREA_1;
+            theta -= HALF_PI;
+        }
+        else if (theta > HALF_PI + FORTPI || theta <= -(HALF_PI + FORTPI)) {
+            area = AREA_ENUM.AREA_2;
+            theta = theta >= 0.0 ? theta - Math.PI : theta + Math.PI;
+        }
+        else {
+            area = AREA_ENUM.AREA_3;
+            theta += HALF_PI;
+        }
+    }
+    /* Compute mu and nu for the area of definition.
+     * For mu, see Eq. (3-21) in [OL76], but note the typos:
+     * compare with Eq. (3-14). For nu, see Eq. (3-38). */
+    mu = Math.atan((12 / Math.PI) * (theta + Math.acos(Math.sin(theta) * Math.cos(FORTPI)) - HALF_PI));
+    t = Math.sqrt((2 * Math.sin(phi / 2) ** 2) /
+        (Math.cos(mu) * Math.cos(mu)) /
+        (1 - Math.cos(Math.atan(1 / Math.cos(theta)))));
+    /* Apply the result to the real area. */
+    if (area === AREA_ENUM.AREA_1) {
+        mu += HALF_PI;
+    }
+    else if (area === AREA_ENUM.AREA_2) {
+        mu += Math.PI;
+    }
+    else if (area === AREA_ENUM.AREA_3) {
+        mu += 1.5 * Math.PI;
+    }
+    /* Now compute x, y from mu and nu */
+    xy.x = t * Math.cos(mu);
+    xy.y = t * Math.sin(mu);
+    xy.x = xy.x * state.a + state.x0;
+    xy.y = xy.y * state.a + state.y0;
+    return p;
+}
+// QSC inverse equations--mapping x,y to lat/long
+// -----------------------------------------------------------------
+export function inverse(state, p) {
+    var lam = 0, latitude = 0;
+    var mu, cosmu, tannu;
+    var tantheta, theta, cosphi, phi;
+    var t;
+    var area = 0;
+    /* de-offset */
+    p.x = (p.x - state.x0) / state.a;
+    p.y = (p.y - state.y0) / state.a;
+    /* Convert the input x, y to the mu and nu angles as used by QSC.
+     * This depends on the area of the cube face. */
+    mu = Math.atan2(p.y, p.x);
+    if (p.x >= 0.0 && p.x >= Math.abs(p.y)) {
+        area = AREA_ENUM.AREA_0;
+    }
+    else if (p.y >= 0.0 && p.y >= Math.abs(p.x)) {
+        area = AREA_ENUM.AREA_1;
+        mu -= HALF_PI;
+    }
+    else if (p.x < 0.0 && -p.x >= Math.abs(p.y)) {
+        area = AREA_ENUM.AREA_2;
+        mu = mu < 0.0 ? mu + Math.PI : mu - Math.PI;
+    }
+    else {
+        area = AREA_ENUM.AREA_3;
+        mu += HALF_PI;
+    }
+    /* Compute phi and theta for the area of definition.
+     * The inverse projection is not described in the original paper, but some
+     * good hints can be found here (as of 2011-12-14):
+     * http://fits.gsfc.nasa.gov/fitsbits/saf.93/saf.9302
+     * (search for "Message-Id: <9302181759.AA25477 at fits.cv.nrao.edu>") */
+    t = (Math.PI / 12) * Math.tan(mu);
+    tantheta = Math.sin(t) / (Math.cos(t) - 1 / Math.sqrt(2));
+    theta = Math.atan(tantheta);
+    cosmu = Math.cos(mu);
+    tannu = Math.hypot(p.x, p.y);
+    const oneMinusCosPhi = cosmu * cosmu * tannu * tannu * (1 - Math.cos(Math.atan(1 / Math.cos(theta))));
+    cosphi = 1 - oneMinusCosPhi;
+    if (cosphi < -1) {
+        cosphi = -1;
+    }
+    else if (cosphi > +1) {
+        cosphi = +1;
+    }
+    /* Apply the result to the real area on the cube face.
+     * For the top and bottom face, we can compute phi and lam directly.
+     * For the other faces, we must use unit sphere cartesian coordinates
+     * as an intermediate step. */
+    if (state.face === FACE_ENUM.TOP) {
+        phi = 2 * Math.asin(Math.sqrt(Math.min(1, oneMinusCosPhi / 2)));
+        latitude = HALF_PI - phi;
+        if (area === AREA_ENUM.AREA_0) {
+            lam = theta + HALF_PI;
+        }
+        else if (area === AREA_ENUM.AREA_1) {
+            lam = theta < 0.0 ? theta + Math.PI : theta - Math.PI;
+        }
+        else if (area === AREA_ENUM.AREA_2) {
+            lam = theta - HALF_PI;
+        } /* area == AREA_ENUM.AREA_3 */
+        else {
+            lam = theta;
+        }
+    }
+    else if (state.face === FACE_ENUM.BOTTOM) {
+        phi = 2 * Math.asin(Math.sqrt(Math.min(1, oneMinusCosPhi / 2)));
+        latitude = phi - HALF_PI;
+        if (area === AREA_ENUM.AREA_0) {
+            lam = -theta + HALF_PI;
+        }
+        else if (area === AREA_ENUM.AREA_1) {
+            lam = -theta;
+        }
+        else if (area === AREA_ENUM.AREA_2) {
+            lam = -theta - HALF_PI;
+        } /* area == AREA_ENUM.AREA_3 */
+        else {
+            lam = theta < 0.0 ? -theta - Math.PI : -theta + Math.PI;
+        }
+    }
+    else {
+        /* Compute phi and lam via cartesian unit sphere coordinates. */
+        var q, r, s;
+        q = cosphi;
+        // Retain sin(phi) from the half-angle complement. Computing 1-q*q
+        // discards face-centre offsets after q has rounded to one.
+        const sinphi = Math.sqrt(Math.max(0, oneMinusCosPhi * (2 - oneMinusCosPhi)));
+        s = sinphi * Math.sin(theta);
+        r = sinphi * Math.cos(theta);
+        /* Rotate q,r,s into the correct area. */
+        if (area === AREA_ENUM.AREA_1) {
+            t = r;
+            r = -s;
+            s = t;
+        }
+        else if (area === AREA_ENUM.AREA_2) {
+            r = -r;
+            s = -s;
+        }
+        else if (area === AREA_ENUM.AREA_3) {
+            t = r;
+            r = s;
+            s = -t;
+        }
+        /* Rotate q,r,s into the correct cube face. */
+        if (state.face === FACE_ENUM.RIGHT) {
+            t = q;
+            q = -r;
+            r = t;
+        }
+        else if (state.face === FACE_ENUM.BACK) {
+            q = -q;
+            r = -r;
+        }
+        else if (state.face === FACE_ENUM.LEFT) {
+            t = q;
+            q = r;
+            r = -t;
+        }
+        /* Now compute phi and lam from the unit sphere coordinates. */
+        latitude = Math.atan2(s, Math.hypot(q, r));
+        lam = Math.atan2(r, q);
+        if (state.face === FACE_ENUM.RIGHT) {
+            lam = qsc_shift_lon_origin(lam, -HALF_PI);
+        }
+        else if (state.face === FACE_ENUM.BACK) {
+            lam = qsc_shift_lon_origin(lam, -Math.PI);
+        }
+        else if (state.face === FACE_ENUM.LEFT) {
+            lam = qsc_shift_lon_origin(lam, +HALF_PI);
+        }
+    }
+    /* Apply the shift from the sphere to the ellipsoid as described
+     * in [LK12]. */
+    if (state.es !== 0) {
+        // Geocentric to geodetic latitude without a cancelling difference of squares.
+        latitude = Math.atan2(Math.sin(latitude), state.one_minus_f_squared * Math.cos(latitude));
+    }
+    lam += state.long0;
+    p.x = lam;
+    p.y = latitude;
+    return p;
+}
+/* Equatorial-face theta before scalar area selection in forward. */
+function qsc_fwd_equat_face_theta(phi, y, x) {
+    return phi === 0 ? 0.0 : Math.atan2(y, x);
+}
+/* Helper function: shift the longitude. */
+function qsc_shift_lon_origin(lon, offset) {
+    var slon = lon + offset;
+    if (slon < -Math.PI) {
+        slon += TWO_PI;
+    }
+    else if (slon > +Math.PI) {
+        slon -= TWO_PI;
+    }
+    return slon;
+}
+export function createState(base, options = {}) {
+    const state = { ...base, face: 0, one_minus_f: 0, one_minus_f_squared: 0, ...options };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=qsc.js.map

package/dist/projection/experimental/kernels/qsc.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"qsc.js","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/qsc.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,qTAAqT;AAQrT,mDAAmD;AACnD,2DAA2D;AAE3D,OAAO,EAAC,MAAM,EAAE,OAAO,EAAE,MAAM,EAAC,+BAA4B;AAE5D,eAAe;AACf,IAAI,SAAS,GAAG;IACd,KAAK,EAAE,CAAC;IACR,KAAK,EAAE,CAAC;IACR,IAAI,EAAE,CAAC;IACP,IAAI,EAAE,CAAC;IACP,GAAG,EAAE,CAAC;IACN,MAAM,EAAE,CAAC;CACV,CAAC;AAEF,IAAI,SAAS,GAAG;IACd,MAAM,EAAE,CAAC;IACT,MAAM,EAAE,CAAC;IACT,MAAM,EAAE,CAAC;IACT,MAAM,EAAE,CAAC;CACV,CAAC;AAEF,SAAS,UAAU,CAAC,KAAY;IAC9B,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC;IACzB,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC;IACzB,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,IAAI,CAAC,CAAC;IAC7B,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,IAAI,CAAC,CAAC;IAC/B,KAAK,CAAC,MAAM,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,CAAC;IAEjC,4DAA4D;IAC5D,IAAI,KAAK,CAAC,IAAI,IAAI,OAAO,GAAG,MAAM,GAAG,GAAG,EAAE,CAAC;QACzC,KAAK,CAAC,IAAI,GAAG,SAAS,CAAC,GAAG,CAAC;IAC7B,CAAC;SAAM,IAAI,KAAK,CAAC,IAAI,IAAI,CAAC,CAAC,OAAO,GAAG,MAAM,GAAG,GAAG,CAAC,EAAE,CAAC;QACnD,KAAK,CAAC,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;IAChC,CAAC;SAAM,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,MAAM,EAAE,CAAC;QAC3C,KAAK,CAAC,IAAI,GAAG,SAAS,CAAC,KAAK,CAAC;IAC/B,CAAC;SAAM,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,OAAO,GAAG,MAAM,EAAE,CAAC;QACrD,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC;IACpE,CAAC;SAAM,CAAC;QACN,KAAK,CAAC,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC;IAC9B,CAAC;IAED;8BAC0B;IAC1B,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,KAAK,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QACtD,KAAK,CAAC,mBAAmB,GAAG,KAAK,CAAC,WAAW,GAAG,KAAK,CAAC,WAAW,CAAC;IACpE,CAAC;AACH,CAAC;AAED,iDAAiD;AACjD,oEAAoE;AACpE,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,GAAG,EAAE,GAAG,CAAC;IACb,IAAI,KAAK,EAAE,GAAG,CAAC;IACf,IAAI,CAAC,EAAE,EAAE,CAAC;IACV,SAAS;IACT,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,yCAAyC;IACzC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,KAAK,CAAC;IAEnB;;8BAE0B;IAC1B,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,oBAAoB;QACpB,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,CAAC;SAAM,CAAC;QACN,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACZ,CAAC;IAED;;;;oEAIgE;IAChE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,aAAa;IACxB,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,GAAG,EAAE,CAAC;QACjC,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;QACpB,IAAI,GAAG,IAAI,MAAM,IAAI,GAAG,IAAI,OAAO,GAAG,MAAM,EAAE,CAAC;YAC7C,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,GAAG,GAAG,OAAO,CAAC;QACxB,CAAC;aAAM,IAAI,GAAG,GAAG,OAAO,GAAG,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,OAAO,GAAG,MAAM,CAAC,EAAE,CAAC;YAChE,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;QACpD,CAAC;aAAM,IAAI,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,MAAM,CAAC,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;YACvD,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,GAAG,GAAG,OAAO,CAAC;QACxB,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,GAAG,CAAC;QACd,CAAC;IACH,CAAC;SAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;QAC3C,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;QACpB,IAAI,GAAG,IAAI,MAAM,IAAI,GAAG,IAAI,OAAO,GAAG,MAAM,EAAE,CAAC;YAC7C,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC;QACzB,CAAC;aAAM,IAAI,GAAG,GAAG,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;YAC1C,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,CAAC,GAAG,CAAC;QACf,CAAC;aAAM,IAAI,GAAG,GAAG,CAAC,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,OAAO,GAAG,MAAM,CAAC,EAAE,CAAC;YACvD,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC;QACzB,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;QACtD,CAAC;IACH,CAAC;SAAM,CAAC;QACN,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,IAAI,MAAM,EAAE,MAAM,CAAC;QACnB,IAAI,MAAM,EAAE,MAAM,CAAC;QAEnB,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;YACnC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC;QAC5C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAC5C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC;QAC5C,CAAC;QACD,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC;QACpB,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC;QACpB,CAAC,GAAG,MAAM,CAAC;QAEX,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;YACnC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtC,KAAK,GAAG,wBAAwB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;YAC1C,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACtC,KAAK,GAAG,wBAAwB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YACvC,KAAK,GAAG,wBAAwB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YACvC,KAAK,GAAG,wBAAwB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,CAAC;aAAM,CAAC;YACN,gBAAgB;YAChB,GAAG,GAAG,KAAK,GAAG,CAAC,CAAC;QAClB,CAAC;QACD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,MAAM,EAAE,CAAC;YAC9B,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;QAC1B,CAAC;aAAM,IAAI,KAAK,GAAG,MAAM,IAAI,KAAK,IAAI,OAAO,GAAG,MAAM,EAAE,CAAC;YACvD,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,IAAI,OAAO,CAAC;QACnB,CAAC;aAAM,IAAI,KAAK,GAAG,OAAO,GAAG,MAAM,IAAI,KAAK,IAAI,CAAC,CAAC,OAAO,GAAG,MAAM,CAAC,EAAE,CAAC;YACpE,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,GAAG,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QAC3D,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;YACxB,KAAK,IAAI,OAAO,CAAC;QACnB,CAAC;IACH,CAAC;IAED;;0DAEsD;IACtD,EAAE,GAAG,IAAI,CAAC,IAAI,CACZ,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,OAAO,CAAC,CACnF,CAAC;IACF,CAAC,GAAG,IAAI,CAAC,IAAI,CACX,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;QAC1B,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAC7B,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CACjD,CAAC;IAEF,wCAAwC;IACxC,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;QAC9B,EAAE,IAAI,OAAO,CAAC;IAChB,CAAC;SAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;QACrC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;IAChB,CAAC;SAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;QACrC,EAAE,IAAI,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;IACtB,CAAC;IAED,qCAAqC;IACrC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACxB,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IACjC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAEjC,OAAO,CAAC,CAAC;AACX,CAAC;AAED,iDAAiD;AACjD,oEAAoE;AACpE,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,CAAC,EACT,QAAQ,GAAG,CAAC,CAAC;IACf,IAAI,EAAE,EAAE,KAAK,EAAE,KAAK,CAAC;IACrB,IAAI,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,CAAC;IACjC,IAAI,CAAC,CAAC;IACN,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,eAAe;IACf,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAEjC;oDACgD;IAChD,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1B,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QACvC,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;IAC1B,CAAC;SAAM,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAC9C,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;QACxB,EAAE,IAAI,OAAO,CAAC;IAChB,CAAC;SAAM,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAC9C,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;QACxB,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;IAC9C,CAAC;SAAM,CAAC;QACN,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC;QACxB,EAAE,IAAI,OAAO,CAAC;IAChB,CAAC;IAED;;;;6EAIyE;IACzE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAClC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1D,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC5B,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACrB,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,cAAc,GAClB,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACjF,MAAM,GAAG,CAAC,GAAG,cAAc,CAAC;IAC5B,IAAI,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC;QAChB,MAAM,GAAG,CAAC,CAAC,CAAC;IACd,CAAC;SAAM,IAAI,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC;QACvB,MAAM,GAAG,CAAC,CAAC,CAAC;IACd,CAAC;IAED;;;kCAG8B;IAC9B,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,GAAG,EAAE,CAAC;QACjC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,cAAc,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChE,QAAQ,GAAG,OAAO,GAAG,GAAG,CAAC;QACzB,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YAC9B,GAAG,GAAG,KAAK,GAAG,OAAO,CAAC;QACxB,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,GAAG,GAAG,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QACxD,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,GAAG,GAAG,KAAK,GAAG,OAAO,CAAC;QACxB,CAAC,CAAC,8BAA8B;aAAM,CAAC;YACrC,GAAG,GAAG,KAAK,CAAC;QACd,CAAC;IACH,CAAC;SAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;QAC3C,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,cAAc,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChE,QAAQ,GAAG,GAAG,GAAG,OAAO,CAAC;QACzB,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YAC9B,GAAG,GAAG,CAAC,KAAK,GAAG,OAAO,CAAC;QACzB,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,GAAG,GAAG,CAAC,KAAK,CAAC;QACf,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,GAAG,GAAG,CAAC,KAAK,GAAG,OAAO,CAAC;QACzB,CAAC,CAAC,8BAA8B;aAAM,CAAC;YACrC,GAAG,GAAG,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QAC1D,CAAC;IACH,CAAC;SAAM,CAAC;QACN,gEAAgE;QAChE,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACZ,CAAC,GAAG,MAAM,CAAC;QACX,kEAAkE;QAClE,2DAA2D;QAC3D,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,cAAc,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;QAC7E,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7B,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7B,yCAAyC;QACzC,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YAC9B,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC,CAAC;YACP,CAAC,GAAG,CAAC,CAAC;QACR,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,CAAC,GAAG,CAAC,CAAC,CAAC;YACP,CAAC,GAAG,CAAC,CAAC,CAAC;QACT,CAAC;aAAM,IAAI,IAAI,KAAK,SAAS,CAAC,MAAM,EAAE,CAAC;YACrC,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC,CAAC;QACT,CAAC;QACD,8CAA8C;QAC9C,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;YACnC,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC,CAAC;YACP,CAAC,GAAG,CAAC,CAAC;QACR,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,CAAC,GAAG,CAAC,CAAC,CAAC;YACP,CAAC,GAAG,CAAC,CAAC,CAAC;QACT,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC,CAAC;QACT,CAAC;QACD,+DAA+D;QAC/D,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC3C,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvB,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,KAAK,EAAE,CAAC;YACnC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC;QAC5C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAC5C,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,EAAE,CAAC;YACzC,GAAG,GAAG,oBAAoB,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC;QAC5C,CAAC;IACH,CAAC;IAED;oBACgB;IAChB,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,8EAA8E;QAC9E,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,KAAK,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;IAC5F,CAAC;IAED,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC;IACnB,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;IACf,OAAO,CAAC,CAAC;AACX,CAAC;AAED,oEAAoE;AACpE,SAAS,wBAAwB,CAAC,GAAW,EAAE,CAAS,EAAE,CAAS;IACjE,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAC5C,CAAC;AAED,2CAA2C;AAC3C,SAAS,oBAAoB,CAAC,GAAW,EAAE,MAAc;IACvD,IAAI,IAAI,GAAG,GAAG,GAAG,MAAM,CAAC;IACxB,IAAI,IAAI,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC;QACpB,IAAI,IAAI,MAAM,CAAC;IACjB,CAAC;SAAM,IAAI,IAAI,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,IAAI,MAAM,CAAC;IACjB,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU,EAAC,GAAG,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,mBAAmB,EAAE,CAAC,EAAE,GAAG,OAAO,EAAC,CAAC;IAC5F,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/projection/experimental/kernels/robin.d.ts

@@ -0,0 +1,6 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {};
+export declare function forward(state: State, ll: Point): Point | null | undefined | number;
+export declare function inverse(state: State, xy: Point): Point | null | undefined | number;
+export declare function createState(base: KernelParameters, options?: Partial<State>): State;
+//# sourceMappingURL=robin.d.ts.map

package/dist/projection/experimental/kernels/robin.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"robin.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/robin.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG,EAAE,CAAC;AAgG1C,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAsBlF;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAoDlF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAQvF"}