@math.gl/projection 0.0.0-stage → 5.0.0-alpha.12

package/dist/experimental/kernels/eqearth.js

@@ -0,0 +1,103 @@
+// math.gl
+// SPDX-License-Identifier: Apache-2.0
+// SPDX-FileCopyrightText: 2018 Bernie Jenny, Monash University, Melbourne, Australia (Equal Earth, ported through proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0 lib/projections/eqearth.js. Modified to use explicit TypeScript state. Original Apache-2.0 notice retained below. See ../../../APACHE-2.0-LICENSE.txt and ../../../THIRD-PARTY-NOTICES.md.
+/**
+ * Copyright 2018 Bernie Jenny, Monash University, Melbourne, Australia.
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ * http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ *
+ * Equal Earth is a projection inspired by the Robinson projection, but unlike
+ * the Robinson projection retains the relative size of areas. The projection
+ * was designed in 2018 by Bojan Savric, Tom Patterson and Bernhard Jenny.
+ *
+ * Publication:
+ * Bojan Savric, Tom Patterson & Bernhard Jenny (2018). The Equal Earth map
+ * projection, International Journal of Geographical Information Science,
+ * DOI: 10.1080/13658816.2018.1504949
+ *
+ * Code released August 2018
+ * Ported to JavaScript and adapted for mapshaper-proj by Matthew Bloch August 2018
+ * Modified for proj4js by Andreas Hocevar by Andreas Hocevar March 2024
+ */
+import adjust_lon from "../common/adjust_lon.js";
+import qsfnz from "../common/qsfnz.js";
+import authset from "../common/authset.js";
+import authlat from "../common/authlat.js";
+var A1 = 1.340264, A2 = -0.081106, A3 = 0.000893, A4 = 0.003796, M = Math.sqrt(3) / 2.0;
+function initialize(state) {
+    state.long0 = state.long0 !== undefined ? state.long0 : 0;
+    state.x0 = state.x0 !== undefined ? state.x0 : 0;
+    state.y0 = state.y0 !== undefined ? state.y0 : 0;
+    if (state.es !== 0) {
+        state.apa = authset(state.es);
+        state.qp = qsfnz(state.e, 1);
+        state.rqda = Math.sqrt(0.5 * state.qp);
+    }
+}
+export function forward(state, p) {
+    var lam = adjust_lon(p.x - state.long0, state.over);
+    var phi = p.y;
+    var sinphi = Math.sin(phi);
+    if (state.es !== 0) {
+        sinphi = qsfnz(state.e, sinphi) / state.qp;
+    }
+    var paramLat = Math.asin(M * sinphi), paramLatSq = paramLat * paramLat, paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
+    p.x =
+        (lam * Math.cos(paramLat)) /
+            (M * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq)));
+    p.y = paramLat * (A1 + A2 * paramLatSq + paramLatPow6 * (A3 + A4 * paramLatSq));
+    if (state.es !== 0) {
+        p.x *= state.rqda;
+        p.y *= state.rqda;
+    }
+    p.x = state.a * p.x + state.x0;
+    p.y = state.a * p.y + state.y0;
+    return p;
+}
+export function inverse(state, p) {
+    p.x = (p.x - state.x0) / state.a;
+    p.y = (p.y - state.y0) / state.a;
+    if (state.es !== 0) {
+        p.x /= state.rqda;
+        p.y /= state.rqda;
+    }
+    var EPS = 1e-9, NITER = 12, paramLat = p.y, paramLatSq, paramLatPow6, fy, fpy, dlat, i;
+    for (i = 0; i < NITER; ++i) {
+        paramLatSq = paramLat * paramLat;
+        paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
+        fy = paramLat * (A1 + A2 * paramLatSq + paramLatPow6 * (A3 + A4 * paramLatSq)) - p.y;
+        fpy = A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq);
+        paramLat -= dlat = fy / fpy;
+        if (Math.abs(dlat) < EPS) {
+            break;
+        }
+    }
+    paramLatSq = paramLat * paramLat;
+    paramLatPow6 = paramLatSq * paramLatSq * paramLatSq;
+    p.x =
+        (M * p.x * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq))) /
+            Math.cos(paramLat);
+    p.y = Math.asin(Math.sin(paramLat) / M);
+    if (state.es !== 0) {
+        p.y = authlat(p.y, state.apa);
+    }
+    p.x = adjust_lon(p.x + state.long0, state.over);
+    return p;
+}
+export function createState(base, options = {}) {
+    const state = { ...base, apa: [], qp: 0, rqda: 0, ...options };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=eqearth.js.map

package/dist/experimental/kernels/eqearth.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"eqearth.js","sourceRoot":"","sources":["../../../src/experimental/kernels/eqearth.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,sCAAsC;AACtC,2HAA2H;AAC3H,wPAAwP;AAQxP;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AAEH,OAAO,UAAU,gCAA6B;AAC9C,OAAO,KAAK,2BAAwB;AACpC,OAAO,OAAO,6BAA0B;AACxC,OAAO,OAAO,6BAA0B;AAExC,IAAI,EAAE,GAAG,QAAQ,EACf,EAAE,GAAG,CAAC,QAAQ,EACd,EAAE,GAAG,QAAQ,EACb,EAAE,GAAG,QAAQ,EACb,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;AAEzB,SAAS,UAAU,CAAC,KAAY;IAC9B,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1D,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,KAAK,CAAC,GAAG,GAAG,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QAC9B,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7B,KAAK,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;IACzC,CAAC;AACH,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACpD,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3B,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAC7C,CAAC;IACD,IAAI,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,EAClC,UAAU,GAAG,QAAQ,GAAG,QAAQ,EAChC,YAAY,GAAG,UAAU,GAAG,UAAU,GAAG,UAAU,CAAC;IACtD,CAAC,CAAC,CAAC;QACD,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YAC1B,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,GAAG,YAAY,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACnF,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,UAAU,GAAG,YAAY,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC;IAEhF,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC;QAClB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC;IACpB,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAC/B,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAC/B,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,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,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC;QAClB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC;IACpB,CAAC;IAED,IAAI,GAAG,GAAG,IAAI,EACZ,KAAK,GAAG,EAAE,EACV,QAAQ,GAAG,CAAC,CAAC,CAAC,EACd,UAAU,EACV,YAAY,EACZ,EAAE,EACF,GAAG,EACH,IAAI,EACJ,CAAC,CAAC;IAEJ,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC;QAC3B,UAAU,GAAG,QAAQ,GAAG,QAAQ,CAAC;QACjC,YAAY,GAAG,UAAU,GAAG,UAAU,GAAG,UAAU,CAAC;QACpD,EAAE,GAAG,QAAQ,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,UAAU,GAAG,YAAY,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACrF,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,GAAG,YAAY,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,CAAC,CAAC;QAC/E,QAAQ,IAAI,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;QAC5B,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;YACzB,MAAM;QACR,CAAC;IACH,CAAC;IACD,UAAU,GAAG,QAAQ,GAAG,QAAQ,CAAC;IACjC,YAAY,GAAG,UAAU,GAAG,UAAU,GAAG,UAAU,CAAC;IACpD,CAAC,CAAC,CAAC;QACD,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,GAAG,YAAY,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC;YACtF,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IACrB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC;IAExC,IAAI,KAAK,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACnB,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC;IAChC,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IAChD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU,EAAC,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,GAAG,OAAO,EAAC,CAAC;IACpE,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/experimental/kernels/equi.d.ts

@@ -0,0 +1,9 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    t2: number;
+    t1: 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=equi.d.ts.map

package/dist/experimental/kernels/equi.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"equi.d.ts","sourceRoot":"","sources":["../../../src/experimental/kernels/equi.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ,CAAC;AAWF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAajF;AAID,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CASjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAIvF"}

package/dist/experimental/kernels/equi.js

@@ -0,0 +1,41 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+import adjust_lon from "../common/adjust_lon.js";
+function initialize(state) {
+    state.x0 = state.x0 || 0;
+    state.y0 = state.y0 || 0;
+    state.lat0 = state.lat0 || 0;
+    state.long0 = state.long0 || 0;
+    /// state.t2;
+}
+export function forward(state, p) {
+    var lon = p.x;
+    var lat = p.y;
+    var dlon = adjust_lon(lon - state.long0, state.over);
+    var x = state.x0 + state.a * dlon * Math.cos(state.lat0);
+    var y = state.y0 + state.a * lat;
+    state.t1 = x;
+    state.t2 = Math.cos(state.lat0);
+    p.x = x;
+    p.y = y;
+    return p;
+}
+/* Equirectangular inverse equations--mapping x,y to lat/long
+  --------------------------------------------------------- */
+export function inverse(state, p) {
+    p.x -= state.x0;
+    p.y -= state.y0;
+    var lat = p.y / state.a;
+    var lon = adjust_lon(state.long0 + p.x / (state.a * Math.cos(state.lat0)), state.over);
+    p.x = lon;
+    p.y = lat;
+    return p;
+}
+export function createState(base, options = {}) {
+    const state = { ...base, t2: 0, t1: 0, ...options };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=equi.js.map

package/dist/experimental/kernels/equi.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"equi.js","sourceRoot":"","sources":["../../../src/experimental/kernels/equi.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,sIAAsI;AAOtI,OAAO,UAAU,gCAA6B;AAE9C,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,aAAa;AACf,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAEd,IAAI,IAAI,GAAG,UAAU,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACrD,IAAI,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACzD,IAAI,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC;IAEjC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;IACb,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IAChC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACR,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACR,OAAO,CAAC,CAAC;AACX,CAAC;AAED;8DAC8D;AAC9D,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;IAChB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;IAChB,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAExB,IAAI,GAAG,GAAG,UAAU,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACvF,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,EAAC,GAAG,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,OAAO,EAAC,CAAC;IACzD,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/experimental/kernels/etmerc.d.ts

@@ -0,0 +1,13 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    cgb: number[];
+    cbg: number[];
+    utg: number[];
+    gtu: number[];
+    Qn: number;
+    Zb: 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): State;
+//# sourceMappingURL=etmerc.d.ts.map

package/dist/experimental/kernels/etmerc.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"etmerc.d.ts","sourceRoot":"","sources":["../../../src/experimental/kernels/etmerc.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,GAAG,EAAE,MAAM,EAAE,CAAC;IACd,GAAG,EAAE,MAAM,EAAE,CAAC;IACd,GAAG,EAAE,MAAM,EAAE,CAAC;IACd,GAAG,EAAE,MAAM,EAAE,CAAC;IACd,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ,CAAC;AA0FF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAkCjF;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAoCjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,GAAG,KAAK,CAYzD"}

package/dist/experimental/kernels/etmerc.js

@@ -0,0 +1,145 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// Heavily based on this etmerc projection implementation
+// https://github.com/mbloch/mapshaper-proj/blob/master/src/projections/etmerc.js
+import sinh from "../common/sinh.js";
+import hypot from "../common/hypot.js";
+import asinhy from "../common/asinhy.js";
+import gatg from "../common/gatg.js";
+import clens from "../common/clens.js";
+import clens_cmplx from "../common/clens_cmplx.js";
+import adjust_lon from "../common/adjust_lon.js";
+function initialize(state) {
+    state.x0 = state.x0 !== undefined ? state.x0 : 0;
+    state.y0 = state.y0 !== undefined ? state.y0 : 0;
+    state.long0 = state.long0 !== undefined ? state.long0 : 0;
+    state.lat0 = state.lat0 !== undefined ? state.lat0 : 0;
+    state.k0 = state.k0 !== undefined ? state.k0 : 1;
+    state.cgb = [];
+    state.cbg = [];
+    state.utg = [];
+    state.gtu = [];
+    var f = state.es / (1 + Math.sqrt(1 - state.es));
+    var n = f / (2 - f);
+    var np = n;
+    state.cgb[0] =
+        n * (2 + n * (-2 / 3 + n * (-2 + n * (116 / 45 + n * (26 / 45 + n * (-2854 / 675))))));
+    state.cbg[0] =
+        n * (-2 + n * (2 / 3 + n * (4 / 3 + n * (-82 / 45 + n * (32 / 45 + n * (4642 / 4725))))));
+    np = np * n;
+    state.cgb[1] =
+        np * (7 / 3 + n * (-8 / 5 + n * (-227 / 45 + n * (2704 / 315 + n * (2323 / 945)))));
+    state.cbg[1] =
+        np * (5 / 3 + n * (-16 / 15 + n * (-13 / 9 + n * (904 / 315 + n * (-1522 / 945)))));
+    np = np * n;
+    state.cgb[2] = np * (56 / 15 + n * (-136 / 35 + n * (-1262 / 105 + n * (73814 / 2835))));
+    state.cbg[2] = np * (-26 / 15 + n * (34 / 21 + n * (8 / 5 + n * (-12686 / 2835))));
+    np = np * n;
+    state.cgb[3] = np * (4279 / 630 + n * (-332 / 35 + n * (-399572 / 14175)));
+    state.cbg[3] = np * (1237 / 630 + n * (-12 / 5 + n * (-24832 / 14175)));
+    np = np * n;
+    state.cgb[4] = np * (4174 / 315 + n * (-144838 / 6237));
+    state.cbg[4] = np * (-734 / 315 + n * (109598 / 31185));
+    np = np * n;
+    state.cgb[5] = np * (601676 / 22275);
+    state.cbg[5] = np * (444337 / 155925);
+    np = Math.pow(n, 2);
+    state.Qn = (state.k0 / (1 + n)) * (1 + np * (1 / 4 + np * (1 / 64 + np / 256)));
+    state.utg[0] =
+        n *
+            (-0.5 + n * (2 / 3 + n * (-37 / 96 + n * (1 / 360 + n * (81 / 512 + n * (-96199 / 604800))))));
+    state.gtu[0] =
+        n *
+            (0.5 + n * (-2 / 3 + n * (5 / 16 + n * (41 / 180 + n * (-127 / 288 + n * (7891 / 37800))))));
+    state.utg[1] =
+        np * (-1 / 48 + n * (-1 / 15 + n * (437 / 1440 + n * (-46 / 105 + n * (1118711 / 3870720)))));
+    state.gtu[1] =
+        np * (13 / 48 + n * (-3 / 5 + n * (557 / 1440 + n * (281 / 630 + n * (-1983433 / 1935360)))));
+    np = np * n;
+    state.utg[2] = np * (-17 / 480 + n * (37 / 840 + n * (209 / 4480 + n * (-5569 / 90720))));
+    state.gtu[2] = np * (61 / 240 + n * (-103 / 140 + n * (15061 / 26880 + n * (167603 / 181440))));
+    np = np * n;
+    state.utg[3] = np * (-4397 / 161280 + n * (11 / 504 + n * (830251 / 7257600)));
+    state.gtu[3] = np * (49561 / 161280 + n * (-179 / 168 + n * (6601661 / 7257600)));
+    np = np * n;
+    state.utg[4] = np * (-4583 / 161280 + n * (108847 / 3991680));
+    state.gtu[4] = np * (34729 / 80640 + n * (-3418889 / 1995840));
+    np = np * n;
+    state.utg[5] = np * (-20648693 / 638668800);
+    state.gtu[5] = np * (212378941 / 319334400);
+    var Z = gatg(state.cbg, state.lat0);
+    state.Zb = -state.Qn * (Z + clens(state.gtu, 2 * Z));
+}
+export function forward(state, p) {
+    var Ce = adjust_lon(p.x - state.long0, state.over);
+    var Cn = p.y;
+    Cn = gatg(state.cbg, Cn);
+    var sin_Cn = Math.sin(Cn);
+    var cos_Cn = Math.cos(Cn);
+    var sin_Ce = Math.sin(Ce);
+    var cos_Ce = Math.cos(Ce);
+    Cn = Math.atan2(sin_Cn, cos_Ce * cos_Cn);
+    Ce = Math.atan2(sin_Ce * cos_Cn, hypot(sin_Cn, cos_Cn * cos_Ce));
+    Ce = asinhy(Math.tan(Ce));
+    clens_cmplx(state.gtu, 2 * Cn, 2 * Ce, p);
+    Cn = Cn + p.x;
+    Ce = Ce + p.y;
+    var x;
+    var y;
+    if (Math.abs(Ce) <= 2.623395162778) {
+        x = state.a * (state.Qn * Ce) + state.x0;
+        y = state.a * (state.Qn * Cn + state.Zb) + state.y0;
+    }
+    else {
+        x = Infinity;
+        y = Infinity;
+    }
+    p.x = x;
+    p.y = y;
+    return p;
+}
+export function inverse(state, p) {
+    var Ce = (p.x - state.x0) * (1 / state.a);
+    var Cn = (p.y - state.y0) * (1 / state.a);
+    Cn = (Cn - state.Zb) / state.Qn;
+    Ce = Ce / state.Qn;
+    var lon;
+    var lat;
+    if (Math.abs(Ce) <= 2.623395162778) {
+        clens_cmplx(state.utg, 2 * Cn, 2 * Ce, p);
+        Cn = Cn + p.x;
+        Ce = Ce + p.y;
+        Ce = Math.atan(sinh(Ce));
+        var sin_Cn = Math.sin(Cn);
+        var cos_Cn = Math.cos(Cn);
+        var sin_Ce = Math.sin(Ce);
+        var cos_Ce = Math.cos(Ce);
+        Cn = Math.atan2(sin_Cn * cos_Ce, hypot(sin_Ce, cos_Ce * cos_Cn));
+        Ce = Math.atan2(sin_Ce, cos_Ce * cos_Cn);
+        lon = adjust_lon(Ce + state.long0, state.over);
+        lat = gatg(state.cgb, Cn);
+    }
+    else {
+        lon = Infinity;
+        lat = Infinity;
+    }
+    p.x = lon;
+    p.y = lat;
+    return p;
+}
+export function createState(base) {
+    const state = {
+        ...base,
+        cgb: [],
+        cbg: [],
+        utg: [],
+        gtu: [],
+        Qn: 0,
+        Zb: 0
+    };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=etmerc.js.map

package/dist/experimental/kernels/etmerc.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"etmerc.js","sourceRoot":"","sources":["../../../src/experimental/kernels/etmerc.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,sIAAsI;AAWtI,yDAAyD;AACzD,iFAAiF;AAEjF,OAAO,IAAI,0BAAuB;AAClC,OAAO,KAAK,2BAAwB;AACpC,OAAO,MAAM,4BAAyB;AACtC,OAAO,IAAI,0BAAuB;AAClC,OAAO,KAAK,2BAAwB;AACpC,OAAO,WAAW,iCAA8B;AAChD,OAAO,UAAU,gCAA6B;AAE9C,SAAS,UAAU,CAAC,KAAY;IAC9B,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1D,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEjD,KAAK,CAAC,GAAG,GAAG,EAAE,CAAC;IACf,KAAK,CAAC,GAAG,GAAG,EAAE,CAAC;IACf,KAAK,CAAC,GAAG,GAAG,EAAE,CAAC;IACf,KAAK,CAAC,GAAG,GAAG,EAAE,CAAC;IAEf,IAAI,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IACjD,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpB,IAAI,EAAE,GAAG,CAAC,CAAC;IAEX,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE5F,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtF,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEtF,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAEnF,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAC3E,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAExE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;IACxD,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC;IAExD,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC;IACrC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;IAEtC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,KAAK,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAEhF,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,CAAC;YACD,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACjG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,CAAC;YACD,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE/F,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACV,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAEhG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1F,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAEhG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAC/E,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAElF,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC;IAC9D,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC;IAE/D,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACZ,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,QAAQ,GAAG,SAAS,CAAC,CAAC;IAC5C,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC,CAAC;IAE5C,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACpC,KAAK,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACvD,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IACnD,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IAEb,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;IACzB,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAE1B,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC;IACzC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,EAAE,KAAK,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;IACjE,EAAE,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAE1B,WAAW,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IAE1C,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IAEd,IAAI,CAAC,CAAC;IACN,IAAI,CAAC,CAAC;IAEN,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC;QACnC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;QACzC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IACtD,CAAC;SAAM,CAAC;QACN,CAAC,GAAG,QAAQ,CAAC;QACb,CAAC,GAAG,QAAQ,CAAC;IACf,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACR,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAER,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAC1C,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAE1C,EAAE,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;IAChC,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAEnB,IAAI,GAAG,CAAC;IACR,IAAI,GAAG,CAAC;IAER,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC;QACnC,WAAW,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;QAE1C,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzB,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAC1B,IAAI,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAE1B,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,EAAE,KAAK,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;QACjE,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC,CAAC;QAEzC,GAAG,GAAG,UAAU,CAAC,EAAE,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAC/C,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;IAC5B,CAAC;SAAM,CAAC;QACN,GAAG,GAAG,QAAQ,CAAC;QACf,GAAG,GAAG,QAAQ,CAAC;IACjB,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAEV,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB;IAChD,MAAM,KAAK,GAAU;QACnB,GAAG,IAAI;QACP,GAAG,EAAE,EAAE;QACP,GAAG,EAAE,EAAE;QACP,GAAG,EAAE,EAAE;QACP,GAAG,EAAE,EAAE;QACP,EAAE,EAAE,CAAC;QACL,EAAE,EAAE,CAAC;KACN,CAAC;IACF,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/experimental/kernels/gauss.d.ts

@@ -0,0 +1,12 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    rc: number;
+    C: number;
+    phic0: number;
+    ratexp: number;
+    K: 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): State;
+//# sourceMappingURL=gauss.d.ts.map

package/dist/experimental/kernels/gauss.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gauss.d.ts","sourceRoot":"","sources":["../../../src/experimental/kernels/gauss.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,EAAE,EAAE,MAAM,CAAC;IACX,CAAC,EAAE,MAAM,CAAC;IACV,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IACf,CAAC,EAAE,MAAM,CAAC;CACX,CAAC;AAkBF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAcjF;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAmBjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,GAAG,KAAK,CAWzD"}

package/dist/experimental/kernels/gauss.js

@@ -0,0 +1,64 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+import srat from "../common/srat.js";
+var MAX_ITER = 20;
+import { HALF_PI, FORTPI } from "../common/constants.js";
+function initialize(state) {
+    var sphi = Math.sin(state.lat0);
+    var cphi = Math.cos(state.lat0);
+    cphi *= cphi;
+    state.rc = Math.sqrt(1 - state.es) / (1 - state.es * sphi * sphi);
+    state.C = Math.sqrt(1 + (state.es * cphi * cphi) / (1 - state.es));
+    state.phic0 = Math.asin(sphi / state.C);
+    state.ratexp = 0.5 * state.C * state.e;
+    state.K =
+        Math.tan(0.5 * state.phic0 + FORTPI) /
+            (Math.pow(Math.tan(0.5 * state.lat0 + FORTPI), state.C) * srat(state.e * sphi, state.ratexp));
+}
+export function forward(state, p) {
+    var lon = p.x;
+    var lat = p.y;
+    p.y =
+        2 *
+            Math.atan(state.K *
+                Math.pow(Math.tan(0.5 * lat + FORTPI), state.C) *
+                srat(state.e * Math.sin(lat), state.ratexp)) -
+            HALF_PI;
+    p.x = state.C * lon;
+    return p;
+}
+export function inverse(state, p) {
+    var DEL_TOL = 1e-14;
+    var lon = p.x / state.C;
+    var lat = p.y;
+    var num = Math.pow(Math.tan(0.5 * lat + FORTPI) / state.K, 1 / state.C);
+    for (var i = MAX_ITER; i > 0; --i) {
+        lat = 2 * Math.atan(num * srat(state.e * Math.sin(p.y), -0.5 * state.e)) - HALF_PI;
+        if (Math.abs(lat - p.y) < DEL_TOL) {
+            break;
+        }
+        p.y = lat;
+    }
+    /* convergence failed */
+    if (!i) {
+        return null;
+    }
+    p.x = lon;
+    p.y = lat;
+    return p;
+}
+export function createState(base) {
+    const state = {
+        ...base,
+        rc: 0,
+        C: 0,
+        phic0: 0,
+        ratexp: 0,
+        K: 0
+    };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=gauss.js.map

package/dist/experimental/kernels/gauss.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gauss.js","sourceRoot":"","sources":["../../../src/experimental/kernels/gauss.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,sIAAsI;AAUtI,OAAO,IAAI,0BAAuB;AAClC,IAAI,QAAQ,GAAG,EAAE,CAAC;AAClB,OAAO,EAAC,OAAO,EAAE,MAAM,EAAC,+BAA4B;AAEpD,SAAS,UAAU,CAAC,KAAY;IAC9B,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IAChC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IAChC,IAAI,IAAI,IAAI,CAAC;IACb,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;IAClE,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IACnE,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACxC,KAAK,CAAC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACvC,KAAK,CAAC,CAAC;QACL,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC;YACpC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,IAAI,GAAG,MAAM,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;AAClG,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAEd,CAAC,CAAC,CAAC;QACD,CAAC;YACC,IAAI,CAAC,IAAI,CACP,KAAK,CAAC,CAAC;gBACL,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,MAAM,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;gBAC/C,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,CAC9C;YACH,OAAO,CAAC;IACV,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC;IACpB,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,OAAO,GAAG,KAAK,CAAC;IACpB,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACxB,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,KAAK,IAAI,CAAC,GAAG,QAAQ,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;QAClC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QACnF,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,CAAC;YAClC,MAAM;QACR,CAAC;QACD,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACZ,CAAC;IACD,wBAAwB;IACxB,IAAI,CAAC,CAAC,EAAE,CAAC;QACP,OAAO,IAAI,CAAC;IACd,CAAC;IACD,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;IAChD,MAAM,KAAK,GAAU;QACnB,GAAG,IAAI;QACP,EAAE,EAAE,CAAC;QACL,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,CAAC;QACR,MAAM,EAAE,CAAC;QACT,CAAC,EAAE,CAAC;KACL,CAAC;IACF,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/experimental/kernels/geos.d.ts

@@ -0,0 +1,17 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    flip_axis: number;
+    sweep: string;
+    h: number;
+    radius_g_1: number;
+    radius_g: number;
+    C: number;
+    radius_p: number;
+    radius_p2: number;
+    radius_p_inv2: number;
+    shape: string;
+};
+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=geos.d.ts.map

package/dist/experimental/kernels/geos.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geos.d.ts","sourceRoot":"","sources":["../../../src/experimental/kernels/geos.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,SAAS,EAAE,MAAM,CAAC;IAClB,KAAK,EAAE,MAAM,CAAC;IACd,CAAC,EAAE,MAAM,CAAC;IACV,UAAU,EAAE,MAAM,CAAC;IACnB,QAAQ,EAAE,MAAM,CAAC;IACjB,CAAC,EAAE,MAAM,CAAC;IACV,QAAQ,EAAE,MAAM,CAAC;IACjB,SAAS,EAAE,MAAM,CAAC;IAClB,aAAa,EAAE,MAAM,CAAC;IACtB,KAAK,EAAE,MAAM,CAAC;CACf,CAAC;AAiCF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CA8CjF;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAiEjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAiBvF"}

package/dist/experimental/kernels/geos.js

@@ -0,0 +1,153 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+import hypot from "../common/hypot.js";
+function initialize(state) {
+    state.flip_axis = state.sweep === 'x' ? 1 : 0;
+    state.h = Number(state.h);
+    state.radius_g_1 = state.h / state.a;
+    if (state.radius_g_1 <= 0 || state.radius_g_1 > 1e10) {
+        throw new Error();
+    }
+    state.radius_g = 1.0 + state.radius_g_1;
+    state.C = state.radius_g * state.radius_g - 1.0;
+    if (state.es !== 0.0) {
+        var one_es = 1.0 - state.es;
+        var rone_es = 1 / one_es;
+        state.radius_p = Math.sqrt(one_es);
+        state.radius_p2 = one_es;
+        state.radius_p_inv2 = rone_es;
+        state.shape = 'ellipse'; // Use as a condition in the forward and inverse functions.
+    }
+    else {
+        state.radius_p = 1.0;
+        state.radius_p2 = 1.0;
+        state.radius_p_inv2 = 1.0;
+        state.shape = 'sphere'; // Use as a condition in the forward and inverse functions.
+    }
+}
+export function forward(state, p) {
+    var lon = p.x;
+    var lat = p.y;
+    var tmp, v_x, v_y, v_z;
+    lon = lon - state.long0;
+    if (state.shape === 'ellipse') {
+        lat = Math.atan(state.radius_p2 * Math.tan(lat));
+        var r = state.radius_p / hypot(state.radius_p * Math.cos(lat), Math.sin(lat));
+        v_x = r * Math.cos(lon) * Math.cos(lat);
+        v_y = r * Math.sin(lon) * Math.cos(lat);
+        v_z = r * Math.sin(lat);
+        if ((state.radius_g - v_x) * v_x - v_y * v_y - v_z * v_z * state.radius_p_inv2 < 0.0) {
+            p.x = Number.NaN;
+            p.y = Number.NaN;
+            return p;
+        }
+        tmp = state.radius_g - v_x;
+        if (state.flip_axis) {
+            p.x = state.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
+            p.y = state.radius_g_1 * Math.atan(v_z / tmp);
+        }
+        else {
+            p.x = state.radius_g_1 * Math.atan(v_y / tmp);
+            p.y = state.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
+        }
+    }
+    else if (state.shape === 'sphere') {
+        tmp = Math.cos(lat);
+        v_x = Math.cos(lon) * tmp;
+        v_y = Math.sin(lon) * tmp;
+        v_z = Math.sin(lat);
+        tmp = state.radius_g - v_x;
+        if (state.flip_axis) {
+            p.x = state.radius_g_1 * Math.atan(v_y / hypot(v_z, tmp));
+            p.y = state.radius_g_1 * Math.atan(v_z / tmp);
+        }
+        else {
+            p.x = state.radius_g_1 * Math.atan(v_y / tmp);
+            p.y = state.radius_g_1 * Math.atan(v_z / hypot(v_y, tmp));
+        }
+    }
+    p.x = p.x * state.a;
+    p.y = p.y * state.a;
+    return p;
+}
+export function inverse(state, p) {
+    var v_x = -1.0;
+    var v_y = 0.0;
+    var v_z = 0.0;
+    var a, b, det, k;
+    p.x = p.x / state.a;
+    p.y = p.y / state.a;
+    if (state.shape === 'ellipse') {
+        if (state.flip_axis) {
+            v_z = Math.tan(p.y / state.radius_g_1);
+            v_y = Math.tan(p.x / state.radius_g_1) * hypot(1.0, v_z);
+        }
+        else {
+            v_y = Math.tan(p.x / state.radius_g_1);
+            v_z = Math.tan(p.y / state.radius_g_1) * hypot(1.0, v_y);
+        }
+        var v_zp = v_z / state.radius_p;
+        a = v_y * v_y + v_zp * v_zp + v_x * v_x;
+        b = 2 * state.radius_g * v_x;
+        det = b * b - 4 * a * state.C;
+        if (det < 0.0) {
+            p.x = Number.NaN;
+            p.y = Number.NaN;
+            return p;
+        }
+        k = (-b - Math.sqrt(det)) / (2.0 * a);
+        v_x = state.radius_g + k * v_x;
+        v_y *= k;
+        v_z *= k;
+        p.x = Math.atan2(v_y, v_x);
+        p.y = Math.atan((v_z * Math.cos(p.x)) / v_x);
+        p.y = Math.atan(state.radius_p_inv2 * Math.tan(p.y));
+    }
+    else if (state.shape === 'sphere') {
+        if (state.flip_axis) {
+            v_z = Math.tan(p.y / state.radius_g_1);
+            v_y = Math.tan(p.x / state.radius_g_1) * Math.sqrt(1.0 + v_z * v_z);
+        }
+        else {
+            v_y = Math.tan(p.x / state.radius_g_1);
+            v_z = Math.tan(p.y / state.radius_g_1) * Math.sqrt(1.0 + v_y * v_y);
+        }
+        a = v_y * v_y + v_z * v_z + v_x * v_x;
+        b = 2 * state.radius_g * v_x;
+        det = b * b - 4 * a * state.C;
+        if (det < 0.0) {
+            p.x = Number.NaN;
+            p.y = Number.NaN;
+            return p;
+        }
+        k = (-b - Math.sqrt(det)) / (2.0 * a);
+        v_x = state.radius_g + k * v_x;
+        v_y *= k;
+        v_z *= k;
+        p.x = Math.atan2(v_y, v_x);
+        p.y = Math.atan((v_z * Math.cos(p.x)) / v_x);
+    }
+    p.x = p.x + state.long0;
+    return p;
+}
+export function createState(base, options = {}) {
+    const state = {
+        ...base,
+        flip_axis: 0,
+        sweep: '',
+        h: 0,
+        radius_g_1: 0,
+        radius_g: 0,
+        C: 0,
+        radius_p: 0,
+        radius_p2: 0,
+        radius_p_inv2: 0,
+        shape: '',
+        ...options
+    };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=geos.js.map

package/dist/experimental/kernels/geos.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geos.js","sourceRoot":"","sources":["../../../src/experimental/kernels/geos.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,sIAAsI;AAetI,OAAO,KAAK,2BAAwB;AAEpC,SAAS,UAAU,CAAC,KAAY;IAC9B,KAAK,CAAC,SAAS,GAAG,KAAK,CAAC,KAAK,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9C,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC1B,KAAK,CAAC,UAAU,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAErC,IAAI,KAAK,CAAC,UAAU,IAAI,CAAC,IAAI,KAAK,CAAC,UAAU,GAAG,IAAI,EAAE,CAAC;QACrD,MAAM,IAAI,KAAK,EAAE,CAAC;IACpB,CAAC;IAED,KAAK,CAAC,QAAQ,GAAG,GAAG,GAAG,KAAK,CAAC,UAAU,CAAC;IACxC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;IAEhD,IAAI,KAAK,CAAC,EAAE,KAAK,GAAG,EAAE,CAAC;QACrB,IAAI,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC;QAC5B,IAAI,OAAO,GAAG,CAAC,GAAG,MAAM,CAAC;QAEzB,KAAK,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACnC,KAAK,CAAC,SAAS,GAAG,MAAM,CAAC;QACzB,KAAK,CAAC,aAAa,GAAG,OAAO,CAAC;QAE9B,KAAK,CAAC,KAAK,GAAG,SAAS,CAAC,CAAC,2DAA2D;IACtF,CAAC;SAAM,CAAC;QACN,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;QACrB,KAAK,CAAC,SAAS,GAAG,GAAG,CAAC;QACtB,KAAK,CAAC,aAAa,GAAG,GAAG,CAAC;QAE1B,KAAK,CAAC,KAAK,GAAG,QAAQ,CAAC,CAAC,2DAA2D;IACrF,CAAC;AACH,CAAC;AAED,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,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;IACvB,GAAG,GAAG,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC;IAExB,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;QAC9B,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACjD,IAAI,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,KAAK,CAAC,KAAK,CAAC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAE9E,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACxC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACxC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,KAAK,CAAC,aAAa,GAAG,GAAG,EAAE,CAAC;YACrF,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,OAAO,CAAC,CAAC;QACX,CAAC;QAED,GAAG,GAAG,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;QAC3B,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;YACpB,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;YAC1D,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAChD,CAAC;aAAM,CAAC;YACN,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YAC9C,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;SAAM,IAAI,KAAK,CAAC,KAAK,KAAK,QAAQ,EAAE,CAAC;QACpC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACpB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QAC1B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QAC1B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACpB,GAAG,GAAG,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;QAE3B,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;YACpB,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;YAC1D,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAChD,CAAC;aAAM,CAAC;YACN,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YAC9C,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IACD,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACpB,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC;IACf,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAEjB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAEpB,IAAI,KAAK,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;QAC9B,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;YACpB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;YACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3D,CAAC;aAAM,CAAC;YACN,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;YACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3D,CAAC;QAED,IAAI,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC,QAAQ,CAAC;QAChC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;QACxC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;QAC7B,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAE9B,IAAI,GAAG,GAAG,GAAG,EAAE,CAAC;YACd,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,OAAO,CAAC,CAAC;QACX,CAAC;QAED,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACtC,GAAG,GAAG,KAAK,CAAC,QAAQ,GAAG,CAAC,GAAG,GAAG,CAAC;QAC/B,GAAG,IAAI,CAAC,CAAC;QACT,GAAG,IAAI,CAAC,CAAC;QAET,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC7C,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvD,CAAC;SAAM,IAAI,KAAK,CAAC,KAAK,KAAK,QAAQ,EAAE,CAAC;QACpC,IAAI,KAAK,CAAC,SAAS,EAAE,CAAC;YACpB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;YACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACtE,CAAC;aAAM,CAAC;YACN,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;YACvC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACtE,CAAC;QAED,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACtC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,QAAQ,GAAG,GAAG,CAAC;QAC7B,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC9B,IAAI,GAAG,GAAG,GAAG,EAAE,CAAC;YACd,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACjB,OAAO,CAAC,CAAC;QACX,CAAC;QAED,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACtC,GAAG,GAAG,KAAK,CAAC,QAAQ,GAAG,CAAC,GAAG,GAAG,CAAC;QAC/B,GAAG,IAAI,CAAC,CAAC;QACT,GAAG,IAAI,CAAC,CAAC;QAET,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAC/C,CAAC;IACD,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC;IACxB,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU;QACnB,GAAG,IAAI;QACP,SAAS,EAAE,CAAC;QACZ,KAAK,EAAE,EAAE;QACT,CAAC,EAAE,CAAC;QACJ,UAAU,EAAE,CAAC;QACb,QAAQ,EAAE,CAAC;QACX,CAAC,EAAE,CAAC;QACJ,QAAQ,EAAE,CAAC;QACX,SAAS,EAAE,CAAC;QACZ,aAAa,EAAE,CAAC;QAChB,KAAK,EAAE,EAAE;QACT,GAAG,OAAO;KACX,CAAC;IACF,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}