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

package/dist/triaxial.cjs

@@ -0,0 +1,204 @@
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// dist/triaxial.js
+var triaxial_exports = {};
+__export(triaxial_exports, {
+  TriaxialEllipsoid: () => TriaxialEllipsoid
+});
+module.exports = __toCommonJS(triaxial_exports);
+
+// dist/triaxial/integrate.js
+var C = [0, 1 / 5, 3 / 10, 4 / 5, 8 / 9, 1, 1];
+var A = [
+  [],
+  [1 / 5],
+  [3 / 40, 9 / 40],
+  [44 / 45, -56 / 15, 32 / 9],
+  [19372 / 6561, -25360 / 2187, 64448 / 6561, -212 / 729],
+  [9017 / 3168, -355 / 33, 46732 / 5247, 49 / 176, -5103 / 18656],
+  [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84]
+];
+var B = [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84, 0];
+var ERROR = [
+  35 / 384 - 5179 / 57600,
+  0,
+  500 / 1113 - 7571 / 16695,
+  125 / 192 - 393 / 640,
+  -2187 / 6784 + 92097 / 339200,
+  11 / 84 - 187 / 2100,
+  -1 / 40
+];
+function integrateGeodesic(initial, q, distance) {
+  let state = Array.from(initial), done = 0, step = Math.sign(distance) * Math.min(0.01, Math.abs(distance));
+  const derivative = (y) => {
+    const normal = [q[0] * y[0], q[1] * y[1], q[2] * y[2]];
+    const curvature = (q[0] * y[3] ** 2 + q[1] * y[4] ** 2 + q[2] * y[5] ** 2) / (normal[0] ** 2 + normal[1] ** 2 + normal[2] ** 2);
+    return [
+      y[3],
+      y[4],
+      y[5],
+      -curvature * normal[0],
+      -curvature * normal[1],
+      -curvature * normal[2]
+    ];
+  };
+  for (let attempt = 0; attempt < 5e4; attempt++) {
+    if (done === distance)
+      return state;
+    if (Math.abs(step) > Math.abs(distance - done))
+      step = distance - done;
+    const stages = [];
+    for (let j = 0; j < C.length; j++) {
+      const y = state.map((v, i) => v + step * A[j].reduce((sum, a, k) => sum + a * stages[k][i], 0));
+      stages.push(derivative(y));
+    }
+    const next = state.map((v, i) => v + step * B.reduce((sum, b, j) => sum + b * stages[j][i], 0));
+    let error = 0;
+    for (let i = 0; i < 6; i++)
+      error = Math.max(error, Math.abs(step * ERROR.reduce((sum, e, j) => sum + e * stages[j][i], 0)));
+    const tolerance = 2e-15;
+    if (error <= tolerance) {
+      state = next;
+      done += step;
+      const surface = Math.sqrt(q[0] * state[0] ** 2 + q[1] * state[1] ** 2 + q[2] * state[2] ** 2);
+      for (let i = 0; i < 3; i++)
+        state[i] /= surface;
+      const normal = state.slice(0, 3).map((x, i) => q[i] * x), nn = normal.reduce((s, x) => s + x * x, 0);
+      const dot = normal.reduce((s, x, i) => s + x * state[i + 3], 0) / nn;
+      for (let i = 0; i < 3; i++)
+        state[i + 3] -= normal[i] * dot;
+      const speed = Math.hypot(state[3], state[4], state[5]);
+      for (let i = 3; i < 6; i++)
+        state[i] /= speed;
+    }
+    step *= Math.max(0.2, Math.min(4, 0.9 * (tolerance / Math.max(error, 1e-300)) ** 0.2));
+    if (done + step === done)
+      throw new RangeError("Triaxial integration step underflow");
+  }
+  throw new RangeError("Triaxial integration work limit exceeded");
+}
+
+// dist/triaxial.js
+function unit(v) {
+  const scale = Math.max(Math.abs(v[0]), Math.abs(v[1]), Math.abs(v[2]));
+  if (!(scale > 0 && Number.isFinite(scale)))
+    throw new RangeError("Expected nonzero finite direction");
+  const x = v[0] / scale, y = v[1] / scale, z = v[2] / scale, n = Math.hypot(x, y, z);
+  return [x / n, y / n, z / n];
+}
+function vector(v) {
+  if (v.length !== 3 || !v.every(Number.isFinite))
+    throw new RangeError("Expected three finite coordinates");
+}
+function system(s) {
+  if (!["geodetic", "parametric", "geocentric"].includes(s))
+    throw new RangeError("Unknown triaxial coordinate system");
+}
+function trig(degrees) {
+  const x = (degrees % 360 + 540) % 360 - 180;
+  return {
+    s: x === 0 || Math.abs(x) === 180 ? 0 : Math.sin(x * Math.PI / 180),
+    c: Math.abs(x) === 90 ? 0 : Math.cos(x * Math.PI / 180)
+  };
+}
+var TriaxialEllipsoid = class {
+  constructor({ radii }) {
+    if (radii.length !== 3 || !radii.every((x) => Number.isFinite(x) && x >= 1e-60 && x <= 1e60) || Math.min(...radii) / Math.max(...radii) < 0.25)
+      throw new RangeError("Triaxial axes require [1e-60,1e60] and min/max >= 0.25");
+    this.snapshot = Object.freeze([...radii]);
+    this.scale = Math.max(...radii);
+    this.axes = Object.freeze(radii.map((r) => r / this.scale));
+  }
+  get radii() {
+    return this.snapshot;
+  }
+  toCartesian(latitude, longitude, coordinates = "geodetic") {
+    system(coordinates);
+    if (!Number.isFinite(latitude) || Math.abs(latitude) > 90 || !Number.isFinite(longitude))
+      throw new RangeError("Invalid triaxial latitude/longitude");
+    const p = trig(latitude), l = trig(longitude), n = [p.c * l.c, p.c * l.s, p.s], r = this.axes;
+    if (coordinates === "parametric")
+      return n.map((x, i) => x * this.radii[i]);
+    const denominator = coordinates === "geodetic" ? Math.hypot(...n.map((x, i) => x * r[i])) : Math.hypot(...n.map((x, i) => x / r[i]));
+    return n.map((x, i) => this.scale * x * (coordinates === "geodetic" ? r[i] * r[i] : 1) / denominator);
+  }
+  fromCartesian(position, coordinates = "geodetic") {
+    system(coordinates);
+    const p = this.surface(position);
+    const v = p.map((x, i) => coordinates === "geodetic" ? x / this.axes[i] ** 2 : coordinates === "parametric" ? x / this.axes[i] : x);
+    return {
+      latitude: Math.atan2(v[2], Math.hypot(v[0], v[1])) * 180 / Math.PI,
+      longitude: v[0] === 0 && v[1] === 0 ? 0 : Math.atan2(v[1], v[0]) * 180 / Math.PI
+    };
+  }
+  directionToCartesian(position, azimuth) {
+    if (!Number.isFinite(azimuth))
+      throw new RangeError("Invalid triaxial bearing");
+    const { north, east } = this.basis(position), a = trig(azimuth);
+    return unit(north.map((x, i) => x * a.c + east[i] * a.s));
+  }
+  directionFromCartesian(position, direction) {
+    vector(direction);
+    const { north, east, normal } = this.basis(position), v = unit(direction);
+    if (Math.abs(v.reduce((s, x, i) => s + x * normal[i], 0)) > 2e-12)
+      throw new RangeError("Direction must be tangent");
+    return Math.atan2(v.reduce((s, x, i) => s + x * east[i], 0), v.reduce((s, x, i) => s + x * north[i], 0)) * 180 / Math.PI;
+  }
+  surface(position) {
+    vector(position);
+    const p = Array.from(position, (x) => x / this.scale);
+    const residual = p.reduce((s, x, i) => s + (x / this.axes[i]) ** 2, 0);
+    if (!Number.isFinite(residual) || Math.abs(residual - 1) > 2e-12)
+      throw new RangeError("Point must lie on the ellipsoid surface");
+    return p;
+  }
+  basis(position) {
+    const p = this.surface(position), normal = unit(p.map((x, i) => x / this.axes[i] ** 2));
+    const horizontal = Math.hypot(p[0], p[1]), meridian = horizontal ? [-p[1] / horizontal, p[0] / horizontal, 0] : [0, 1, 0];
+    const cross = (a, b) => [
+      a[1] * b[2] - a[2] * b[1],
+      a[2] * b[0] - a[0] * b[2],
+      a[0] * b[1] - a[1] * b[0]
+    ];
+    const north = unit(cross(normal, meridian)), east = unit(cross(north, normal));
+    return { north, east, normal };
+  }
+  /** Cartesian initial-value geodesic. No global triaxial inverse is implied. */
+  line(position, direction) {
+    const p = this.surface(position);
+    this.directionFromCartesian(position, direction);
+    const v = unit(direction), initial = Object.freeze([...p, ...v]), q = this.axes.map((r) => 1 / (r * r)), scale = this.scale;
+    return {
+      position: (distance) => {
+        if (!Number.isFinite(distance) || Math.abs(distance) > 2 * Math.PI * scale)
+          throw new RangeError("Triaxial distance exceeds one max-axis circuit");
+        const y = distance === 0 ? Array.from(initial) : integrateGeodesic(initial, q, distance / scale);
+        return {
+          position: y.slice(0, 3).map((x) => x * scale),
+          direction: y.slice(3),
+          distance
+        };
+      }
+    };
+  }
+  direct(position, direction, distance) {
+    return this.line(position, direction).position(distance);
+  }
+};
+//# sourceMappingURL=triaxial.cjs.map

package/dist/triaxial.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../src/triaxial.ts", "../src/triaxial/integrate.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAKA;;;;;;;ACDA,IAAM,IAAI,CAAC,GAAG,IAAI,GAAG,IAAI,IAAI,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC;AAC/C,IAAM,IAAI;EACR,CAAA;EACA,CAAC,IAAI,CAAC;EACN,CAAC,IAAI,IAAI,IAAI,EAAE;EACf,CAAC,KAAK,IAAI,MAAM,IAAI,KAAK,CAAC;EAC1B,CAAC,QAAQ,MAAM,SAAS,MAAM,QAAQ,MAAM,OAAO,GAAG;EACtD,CAAC,OAAO,MAAM,OAAO,IAAI,QAAQ,MAAM,KAAK,KAAK,QAAQ,KAAK;EAC9D,CAAC,KAAK,KAAK,GAAG,MAAM,MAAM,MAAM,KAAK,QAAQ,MAAM,KAAK,EAAE;;AAE5D,IAAM,IAAI,CAAC,KAAK,KAAK,GAAG,MAAM,MAAM,MAAM,KAAK,QAAQ,MAAM,KAAK,IAAI,CAAC;AACvE,IAAM,QAAQ;EACZ,KAAK,MAAM,OAAO;EAClB;EACA,MAAM,OAAO,OAAO;EACpB,MAAM,MAAM,MAAM;EAClB,QAAQ,OAAO,QAAQ;EACvB,KAAK,KAAK,MAAM;EAChB,KAAK;;AAED,SAAU,kBACd,SACA,GACA,UAAgB;AAEhB,MAAI,QAAQ,MAAM,KAAK,OAAO,GAC5B,OAAO,GACP,OAAO,KAAK,KAAK,QAAQ,IAAI,KAAK,IAAI,MAAM,KAAK,IAAI,QAAQ,CAAC;AAChE,QAAM,aAAa,CAAC,MAAwB;AAC1C,UAAM,SAAS,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC;AACrD,UAAM,aACH,EAAE,CAAC,IAAI,EAAE,CAAC,KAAK,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,KAAK,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,KAAK,MACrD,OAAO,CAAC,KAAK,IAAI,OAAO,CAAC,KAAK,IAAI,OAAO,CAAC,KAAK;AAClD,WAAO;MACL,EAAE,CAAC;MACH,EAAE,CAAC;MACH,EAAE,CAAC;MACH,CAAC,YAAY,OAAO,CAAC;MACrB,CAAC,YAAY,OAAO,CAAC;MACrB,CAAC,YAAY,OAAO,CAAC;;EAEzB;AACA,WAAS,UAAU,GAAG,UAAU,KAAO,WAAW;AAChD,QAAI,SAAS;AAAU,aAAO;AAC9B,QAAI,KAAK,IAAI,IAAI,IAAI,KAAK,IAAI,WAAW,IAAI;AAAG,aAAO,WAAW;AAClE,UAAM,SAAqB,CAAA;AAC3B,aAAS,IAAI,GAAG,IAAI,EAAE,QAAQ,KAAK;AACjC,YAAM,IAAI,MAAM,IACd,CAAC,GAAG,MAAM,IAAI,OAAO,EAAE,CAAC,EAAE,OAAO,CAAC,KAAK,GAAG,MAAM,MAAM,IAAI,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;AAE5E,aAAO,KAAK,WAAW,CAAC,CAAC;IAC3B;AACA,UAAM,OAAO,MAAM,IAAI,CAAC,GAAG,MAAM,IAAI,OAAO,EAAE,OAAO,CAAC,KAAK,GAAG,MAAM,MAAM,IAAI,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;AAC9F,QAAI,QAAQ;AACZ,aAAS,IAAI,GAAG,IAAI,GAAG;AACrB,cAAQ,KAAK,IACX,OACA,KAAK,IAAI,OAAO,MAAM,OAAO,CAAC,KAAK,GAAG,MAAM,MAAM,IAAI,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AAE3E,UAAM,YAAY;AAClB,QAAI,SAAS,WAAW;AACtB,cAAQ;AACR,cAAQ;AAER,YAAM,UAAU,KAAK,KAAK,EAAE,CAAC,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,CAAC,IAAI,MAAM,CAAC,KAAK,IAAI,EAAE,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC;AAC5F,eAAS,IAAI,GAAG,IAAI,GAAG;AAAK,cAAM,CAAC,KAAK;AACxC,YAAM,SAAS,MAAM,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,IAAI,CAAC,GACrD,KAAK,OAAO,OAAO,CAAC,GAAG,MAAM,IAAI,IAAI,GAAG,CAAC;AAC3C,YAAM,MAAM,OAAO,OAAO,CAAC,GAAG,GAAG,MAAM,IAAI,IAAI,MAAM,IAAI,CAAC,GAAG,CAAC,IAAI;AAClE,eAAS,IAAI,GAAG,IAAI,GAAG;AAAK,cAAM,IAAI,CAAC,KAAK,OAAO,CAAC,IAAI;AACxD,YAAM,QAAQ,KAAK,MAAM,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC;AACrD,eAAS,IAAI,GAAG,IAAI,GAAG;AAAK,cAAM,CAAC,KAAK;IAC1C;AACA,YAAQ,KAAK,IAAI,KAAK,KAAK,IAAI,GAAG,OAAO,YAAY,KAAK,IAAI,OAAO,MAAM,MAAM,GAAG,CAAC;AACrF,QAAI,OAAO,SAAS;AAAM,YAAM,IAAI,WAAW,qCAAqC;EACtF;AACA,QAAM,IAAI,WAAW,0CAA0C;AACjE;;;ADhEA,SAAS,KAAK,GAAM;AAClB,QAAM,QAAQ,KAAK,IAAI,KAAK,IAAI,EAAE,CAAC,CAAC,GAAG,KAAK,IAAI,EAAE,CAAC,CAAC,GAAG,KAAK,IAAI,EAAE,CAAC,CAAC,CAAC;AACrE,MAAI,EAAE,QAAQ,KAAK,OAAO,SAAS,KAAK;AACtC,UAAM,IAAI,WAAW,mCAAmC;AAC1D,QAAM,IAAI,EAAE,CAAC,IAAI,OACf,IAAI,EAAE,CAAC,IAAI,OACX,IAAI,EAAE,CAAC,IAAI,OACX,IAAI,KAAK,MAAM,GAAG,GAAG,CAAC;AACxB,SAAO,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;AAC7B;AACA,SAAS,OAAO,GAAM;AACpB,MAAI,EAAE,WAAW,KAAK,CAAC,EAAE,MAAM,OAAO,QAAQ;AAC5C,UAAM,IAAI,WAAW,mCAAmC;AAC5D;AACA,SAAS,OAAO,GAAsB;AACpC,MAAI,CAAC,CAAC,YAAY,cAAc,YAAY,EAAE,SAAS,CAAC;AACtD,UAAM,IAAI,WAAW,oCAAoC;AAC7D;AACA,SAAS,KAAK,SAAe;AAC3B,QAAM,KAAO,UAAU,MAAO,OAAO,MAAO;AAC5C,SAAO;IACL,GAAG,MAAM,KAAK,KAAK,IAAI,CAAC,MAAM,MAAM,IAAI,KAAK,IAAK,IAAI,KAAK,KAAM,GAAG;IACpE,GAAG,KAAK,IAAI,CAAC,MAAM,KAAK,IAAI,KAAK,IAAK,IAAI,KAAK,KAAM,GAAG;;AAE5D;AAIM,IAAO,oBAAP,MAAwB;EAI5B,YAAY,EAAC,MAAK,GAAqB;AACrC,QACE,MAAM,WAAW,KACjB,CAAC,MAAM,MAAM,OAAK,OAAO,SAAS,CAAC,KAAK,KAAK,SAAS,KAAK,IAAI,KAC/D,KAAK,IAAI,GAAG,KAAK,IAAI,KAAK,IAAI,GAAG,KAAK,IAAI;AAE1C,YAAM,IAAI,WAAW,wDAAwD;AAC/E,SAAK,WAAW,OAAO,OAAO,CAAC,GAAG,KAAK,CAAC;AACxC,SAAK,QAAQ,KAAK,IAAI,GAAG,KAAK;AAC9B,SAAK,OAAO,OAAO,OAAO,MAAM,IAAI,OAAK,IAAI,KAAK,KAAK,CAAC;EAC1D;EACA,IAAI,QAAK;AACP,WAAO,KAAK;EACd;EACA,YACE,UACA,WACA,cAAmC,YAAU;AAE7C,WAAO,WAAW;AAClB,QAAI,CAAC,OAAO,SAAS,QAAQ,KAAK,KAAK,IAAI,QAAQ,IAAI,MAAM,CAAC,OAAO,SAAS,SAAS;AACrF,YAAM,IAAI,WAAW,qCAAqC;AAC5D,UAAM,IAAI,KAAK,QAAQ,GACrB,IAAI,KAAK,SAAS,GAClB,IAAI,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,GAC9B,IAAI,KAAK;AACX,QAAI,gBAAgB;AAClB,aAAO,EAAE,IAAI,CAAC,GAAG,MAAM,IAAI,KAAK,MAAM,CAAC,CAAC;AAC1C,UAAM,cACJ,gBAAgB,aACZ,KAAK,MAAM,GAAG,EAAE,IAAI,CAAC,GAAG,MAAM,IAAI,EAAE,CAAC,CAAC,CAAC,IACvC,KAAK,MAAM,GAAG,EAAE,IAAI,CAAC,GAAG,MAAM,IAAI,EAAE,CAAC,CAAC,CAAC;AAC7C,WAAO,EAAE,IACP,CAAC,GAAG,MAAO,KAAK,QAAQ,KAAK,gBAAgB,aAAa,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,KAAM,WAAW;EAE7F;EACA,cACE,UACA,cAAmC,YAAU;AAE7C,WAAO,WAAW;AAClB,UAAM,IAAI,KAAK,QAAQ,QAAQ;AAC/B,UAAM,IAAI,EAAE,IAAI,CAAC,GAAG,MAClB,gBAAgB,aACZ,IAAI,KAAK,KAAK,CAAC,KAAK,IACpB,gBAAgB,eACd,IAAI,KAAK,KAAK,CAAC,IACf,CAAC;AAET,WAAO;MACL,UAAW,KAAK,MAAM,EAAE,CAAC,GAAG,KAAK,MAAM,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,MAAO,KAAK;MAClE,WAAW,EAAE,CAAC,MAAM,KAAK,EAAE,CAAC,MAAM,IAAI,IAAK,KAAK,MAAM,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,MAAO,KAAK;;EAEpF;EACA,qBAAqB,UAAe,SAAe;AACjD,QAAI,CAAC,OAAO,SAAS,OAAO;AAAG,YAAM,IAAI,WAAW,0BAA0B;AAC9E,UAAM,EAAC,OAAO,KAAI,IAAI,KAAK,MAAM,QAAQ,GACvC,IAAI,KAAK,OAAO;AAClB,WAAO,KAAK,MAAM,IAAI,CAAC,GAAG,MAAM,IAAI,EAAE,IAAI,KAAK,CAAC,IAAI,EAAE,CAAC,CAAC;EAC1D;EACA,uBAAuB,UAAe,WAAc;AAClD,WAAO,SAAS;AAChB,UAAM,EAAC,OAAO,MAAM,OAAM,IAAI,KAAK,MAAM,QAAQ,GAC/C,IAAI,KAAK,SAAS;AACpB,QAAI,KAAK,IAAI,EAAE,OAAO,CAAC,GAAG,GAAG,MAAM,IAAI,IAAI,OAAO,CAAC,GAAG,CAAC,CAAC,IAAI;AAC1D,YAAM,IAAI,WAAW,2BAA2B;AAClD,WACG,KAAK,MACJ,EAAE,OAAO,CAAC,GAAG,GAAG,MAAM,IAAI,IAAI,KAAK,CAAC,GAAG,CAAC,GACxC,EAAE,OAAO,CAAC,GAAG,GAAG,MAAM,IAAI,IAAI,MAAM,CAAC,GAAG,CAAC,CAAC,IAE1C,MACF,KAAK;EAET;EACQ,QAAQ,UAAa;AAC3B,WAAO,QAAQ;AACf,UAAM,IAAI,MAAM,KAAK,UAAU,OAAK,IAAI,KAAK,KAAK;AAClD,UAAM,WAAW,EAAE,OAAO,CAAC,GAAG,GAAG,MAAM,KAAK,IAAI,KAAK,KAAK,CAAC,MAAM,GAAG,CAAC;AACrE,QAAI,CAAC,OAAO,SAAS,QAAQ,KAAK,KAAK,IAAI,WAAW,CAAC,IAAI;AACzD,YAAM,IAAI,WAAW,yCAAyC;AAChE,WAAO;EACT;EACQ,MAAM,UAAa;AACzB,UAAM,IAAI,KAAK,QAAQ,QAAQ,GAC7B,SAAS,KAAK,EAAE,IAAI,CAAC,GAAG,MAAM,IAAI,KAAK,KAAK,CAAC,KAAK,CAAC,CAAC;AACtD,UAAM,aAAa,KAAK,MAAM,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,GACtC,WAAW,aAAa,CAAC,CAAC,EAAE,CAAC,IAAI,YAAY,EAAE,CAAC,IAAI,YAAY,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;AAC/E,UAAM,QAAQ,CAAC,GAAQ,MAAW;MAChC,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC;MACxB,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC;MACxB,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC,IAAI,EAAE,CAAC;;AAE1B,UAAM,QAAQ,KAAK,MAAM,QAAQ,QAAQ,CAAC,GACxC,OAAO,KAAK,MAAM,OAAO,MAAM,CAAC;AAClC,WAAO,EAAC,OAAO,MAAM,OAAM;EAC7B;;EAEA,KAAK,UAAe,WAAc;AAChC,UAAM,IAAI,KAAK,QAAQ,QAAQ;AAC/B,SAAK,uBAAuB,UAAU,SAAS;AAC/C,UAAM,IAAI,KAAK,SAAS,GACtB,UAAU,OAAO,OAAO,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GACpC,IAAI,KAAK,KAAK,IAAI,OAAK,KAAK,IAAI,EAAE,GAClC,QAAQ,KAAK;AACf,WAAO;MACL,UAAU,CAAC,aAAoB;AAC7B,YAAI,CAAC,OAAO,SAAS,QAAQ,KAAK,KAAK,IAAI,QAAQ,IAAI,IAAI,KAAK,KAAK;AACnE,gBAAM,IAAI,WAAW,gDAAgD;AACvE,cAAM,IACJ,aAAa,IAAI,MAAM,KAAK,OAAO,IAAI,kBAAkB,SAAS,GAAG,WAAW,KAAK;AACvF,eAAO;UACL,UAAU,EAAE,MAAM,GAAG,CAAC,EAAE,IAAI,OAAK,IAAI,KAAK;UAC1C,WAAW,EAAE,MAAM,CAAC;UACpB;;MAEJ;;EAEJ;EACA,OAAO,UAAe,WAAgB,UAAgB;AACpD,WAAO,KAAK,KAAK,UAAU,SAAS,EAAE,SAAS,QAAQ;EACzD;;",
+  "names": []
+}

package/dist/triaxial.d.ts

@@ -0,0 +1,38 @@
+import type { NumericArray } from '@math.gl/core/types';
+export type TriaxialCoordinates = 'geodetic' | 'parametric' | 'geocentric';
+export type TriaxialParameters = {
+    radii: readonly [number, number, number];
+};
+export type CartesianGeodesicState = {
+    position: [number, number, number];
+    direction: [number, number, number];
+    distance: number;
+};
+export interface TriaxialGeodesicLine {
+    position(distance: number): CartesianGeodesicState;
+}
+type Vec = Readonly<NumericArray>;
+/** Surface-only coordinate/direction conversions on three-radius ellipsoids.
+ * Angles are degrees, pole longitude is zero. Shapes need not be axis-sorted.
+ */
+export declare class TriaxialEllipsoid {
+    private readonly snapshot;
+    private readonly scale;
+    private readonly axes;
+    constructor({ radii }: TriaxialParameters);
+    get radii(): readonly [number, number, number];
+    toCartesian(latitude: number, longitude: number, coordinates?: TriaxialCoordinates): [number, number, number];
+    fromCartesian(position: Vec, coordinates?: TriaxialCoordinates): {
+        latitude: number;
+        longitude: number;
+    };
+    directionToCartesian(position: Vec, azimuth: number): [number, number, number];
+    directionFromCartesian(position: Vec, direction: Vec): number;
+    private surface;
+    private basis;
+    /** Cartesian initial-value geodesic. No global triaxial inverse is implied. */
+    line(position: Vec, direction: Vec): TriaxialGeodesicLine;
+    direct(position: Vec, direction: Vec, distance: number): CartesianGeodesicState;
+}
+export {};
+//# sourceMappingURL=triaxial.d.ts.map

package/dist/triaxial.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"triaxial.d.ts","sourceRoot":"","sources":["../src/triaxial.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,qBAAqB,CAAC;AAEtD,MAAM,MAAM,mBAAmB,GAAG,UAAU,GAAG,YAAY,GAAG,YAAY,CAAC;AAC3E,MAAM,MAAM,kBAAkB,GAAG;IAAC,KAAK,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;CAAC,CAAC;AAC5E,MAAM,MAAM,sBAAsB,GAAG;IACnC,QAAQ,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACnC,SAAS,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACpC,QAAQ,EAAE,MAAM,CAAC;CAClB,CAAC;AACF,MAAM,WAAW,oBAAoB;IACnC,QAAQ,CAAC,QAAQ,EAAE,MAAM,GAAG,sBAAsB,CAAC;CACpD;AACD,KAAK,GAAG,GAAG,QAAQ,CAAC,YAAY,CAAC,CAAC;AA0BlC;;GAEG;AACH,qBAAa,iBAAiB;IAC5B,OAAO,CAAC,QAAQ,CAAC,QAAQ,CAAoC;IAC7D,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAS;IAC/B,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAoB;gBAC7B,EAAC,KAAK,EAAC,EAAE,kBAAkB;IAWvC,IAAI,KAAK,IAAI,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAE7C;IACD,WAAW,CACT,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,WAAW,GAAE,mBAAgC,GAC5C,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;IAkB3B,aAAa,CACX,QAAQ,EAAE,GAAG,EACb,WAAW,GAAE,mBAAgC,GAC5C;QAAC,QAAQ,EAAE,MAAM,CAAC;QAAC,SAAS,EAAE,MAAM,CAAA;KAAC;IAexC,oBAAoB,CAAC,QAAQ,EAAE,GAAG,EAAE,OAAO,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;IAM9E,sBAAsB,CAAC,QAAQ,EAAE,GAAG,EAAE,SAAS,EAAE,GAAG,GAAG,MAAM;IAe7D,OAAO,CAAC,OAAO;IAQf,OAAO,CAAC,KAAK;IAcb,+EAA+E;IAC/E,IAAI,CAAC,QAAQ,EAAE,GAAG,EAAE,SAAS,EAAE,GAAG,GAAG,oBAAoB;IAqBzD,MAAM,CAAC,QAAQ,EAAE,GAAG,EAAE,SAAS,EAAE,GAAG,EAAE,QAAQ,EAAE,MAAM,GAAG,sBAAsB;CAGhF"}

package/dist/triaxial.js

@@ -0,0 +1,121 @@
+import { integrateGeodesic } from "./triaxial/integrate.js";
+function unit(v) {
+    const scale = Math.max(Math.abs(v[0]), Math.abs(v[1]), Math.abs(v[2]));
+    if (!(scale > 0 && Number.isFinite(scale)))
+        throw new RangeError('Expected nonzero finite direction');
+    const x = v[0] / scale, y = v[1] / scale, z = v[2] / scale, n = Math.hypot(x, y, z);
+    return [x / n, y / n, z / n];
+}
+function vector(v) {
+    if (v.length !== 3 || !v.every(Number.isFinite))
+        throw new RangeError('Expected three finite coordinates');
+}
+function system(s) {
+    if (!['geodetic', 'parametric', 'geocentric'].includes(s))
+        throw new RangeError('Unknown triaxial coordinate system');
+}
+function trig(degrees) {
+    const x = (((degrees % 360) + 540) % 360) - 180;
+    return {
+        s: x === 0 || Math.abs(x) === 180 ? 0 : Math.sin((x * Math.PI) / 180),
+        c: Math.abs(x) === 90 ? 0 : Math.cos((x * Math.PI) / 180)
+    };
+}
+/** Surface-only coordinate/direction conversions on three-radius ellipsoids.
+ * Angles are degrees, pole longitude is zero. Shapes need not be axis-sorted.
+ */
+export class TriaxialEllipsoid {
+    constructor({ radii }) {
+        if (radii.length !== 3 ||
+            !radii.every(x => Number.isFinite(x) && x >= 1e-60 && x <= 1e60) ||
+            Math.min(...radii) / Math.max(...radii) < 0.25)
+            throw new RangeError('Triaxial axes require [1e-60,1e60] and min/max >= 0.25');
+        this.snapshot = Object.freeze([...radii]);
+        this.scale = Math.max(...radii);
+        this.axes = Object.freeze(radii.map(r => r / this.scale));
+    }
+    get radii() {
+        return this.snapshot;
+    }
+    toCartesian(latitude, longitude, coordinates = 'geodetic') {
+        system(coordinates);
+        if (!Number.isFinite(latitude) || Math.abs(latitude) > 90 || !Number.isFinite(longitude))
+            throw new RangeError('Invalid triaxial latitude/longitude');
+        const p = trig(latitude), l = trig(longitude), n = [p.c * l.c, p.c * l.s, p.s], r = this.axes;
+        if (coordinates === 'parametric')
+            return n.map((x, i) => x * this.radii[i]);
+        const denominator = coordinates === 'geodetic'
+            ? Math.hypot(...n.map((x, i) => x * r[i]))
+            : Math.hypot(...n.map((x, i) => x / r[i]));
+        return n.map((x, i) => (this.scale * x * (coordinates === 'geodetic' ? r[i] * r[i] : 1)) / denominator);
+    }
+    fromCartesian(position, coordinates = 'geodetic') {
+        system(coordinates);
+        const p = this.surface(position);
+        const v = p.map((x, i) => coordinates === 'geodetic'
+            ? x / this.axes[i] ** 2
+            : coordinates === 'parametric'
+                ? x / this.axes[i]
+                : x);
+        return {
+            latitude: (Math.atan2(v[2], Math.hypot(v[0], v[1])) * 180) / Math.PI,
+            longitude: v[0] === 0 && v[1] === 0 ? 0 : (Math.atan2(v[1], v[0]) * 180) / Math.PI
+        };
+    }
+    directionToCartesian(position, azimuth) {
+        if (!Number.isFinite(azimuth))
+            throw new RangeError('Invalid triaxial bearing');
+        const { north, east } = this.basis(position), a = trig(azimuth);
+        return unit(north.map((x, i) => x * a.c + east[i] * a.s));
+    }
+    directionFromCartesian(position, direction) {
+        vector(direction);
+        const { north, east, normal } = this.basis(position), v = unit(direction);
+        if (Math.abs(v.reduce((s, x, i) => s + x * normal[i], 0)) > 2e-12)
+            throw new RangeError('Direction must be tangent');
+        return ((Math.atan2(v.reduce((s, x, i) => s + x * east[i], 0), v.reduce((s, x, i) => s + x * north[i], 0)) *
+            180) /
+            Math.PI);
+    }
+    surface(position) {
+        vector(position);
+        const p = Array.from(position, x => x / this.scale);
+        const residual = p.reduce((s, x, i) => s + (x / this.axes[i]) ** 2, 0);
+        if (!Number.isFinite(residual) || Math.abs(residual - 1) > 2e-12)
+            throw new RangeError('Point must lie on the ellipsoid surface');
+        return p;
+    }
+    basis(position) {
+        const p = this.surface(position), normal = unit(p.map((x, i) => x / this.axes[i] ** 2));
+        const horizontal = Math.hypot(p[0], p[1]), meridian = horizontal ? [-p[1] / horizontal, p[0] / horizontal, 0] : [0, 1, 0];
+        const cross = (a, b) => [
+            a[1] * b[2] - a[2] * b[1],
+            a[2] * b[0] - a[0] * b[2],
+            a[0] * b[1] - a[1] * b[0]
+        ];
+        const north = unit(cross(normal, meridian)), east = unit(cross(north, normal));
+        return { north, east, normal };
+    }
+    /** Cartesian initial-value geodesic. No global triaxial inverse is implied. */
+    line(position, direction) {
+        const p = this.surface(position);
+        this.directionFromCartesian(position, direction);
+        const v = unit(direction), initial = Object.freeze([...p, ...v]), q = this.axes.map(r => 1 / (r * r)), scale = this.scale;
+        return {
+            position: (distance) => {
+                if (!Number.isFinite(distance) || Math.abs(distance) > 2 * Math.PI * scale)
+                    throw new RangeError('Triaxial distance exceeds one max-axis circuit');
+                const y = distance === 0 ? Array.from(initial) : integrateGeodesic(initial, q, distance / scale);
+                return {
+                    position: y.slice(0, 3).map(x => x * scale),
+                    direction: y.slice(3),
+                    distance
+                };
+            }
+        };
+    }
+    direct(position, direction, distance) {
+        return this.line(position, direction).position(distance);
+    }
+}
+//# sourceMappingURL=triaxial.js.map

package/dist/triaxial.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"triaxial.js","sourceRoot":"","sources":["../src/triaxial.ts"],"names":[],"mappings":"AAKA,OAAO,EAAC,iBAAiB,EAAC,gCAA6B;AAYvD,SAAS,IAAI,CAAC,CAAM;IAClB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvE,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QACxC,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;IAC5D,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EACpB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAChB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAChB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1B,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;AAC/B,CAAC;AACD,SAAS,MAAM,CAAC,CAAM;IACpB,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAC7C,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;AAC9D,CAAC;AACD,SAAS,MAAM,CAAC,CAAsB;IACpC,IAAI,CAAC,CAAC,UAAU,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC;QACvD,MAAM,IAAI,UAAU,CAAC,oCAAoC,CAAC,CAAC;AAC/D,CAAC;AACD,SAAS,IAAI,CAAC,OAAe;IAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IAChD,OAAO;QACL,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;QACrE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;KAC1D,CAAC;AACJ,CAAC;AACD;;GAEG;AACH,MAAM,OAAO,iBAAiB;IAI5B,YAAY,EAAC,KAAK,EAAqB;QACrC,IACE,KAAK,CAAC,MAAM,KAAK,CAAC;YAClB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC;YAChE,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI;YAE9C,MAAM,IAAI,UAAU,CAAC,wDAAwD,CAAC,CAAC;QACjF,IAAI,CAAC,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,KAAK,CAAC,CAAsC,CAAC;QAC/E,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;QAChC,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;IAC5D,CAAC;IACD,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IACD,WAAW,CACT,QAAgB,EAChB,SAAiB,EACjB,cAAmC,UAAU;QAE7C,MAAM,CAAC,WAAW,CAAC,CAAC;QACpB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,EAAE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;YACtF,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;QAC9D,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,EACtB,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,EACnB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAC/B,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QAChB,IAAI,WAAW,KAAK,YAAY;YAC9B,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAA6B,CAAC;QACxE,MAAM,WAAW,GACf,WAAW,KAAK,UAAU;YACxB,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/C,OAAO,CAAC,CAAC,GAAG,CACV,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,WAAW,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAC9D,CAAC;IAChC,CAAC;IACD,aAAa,CACX,QAAa,EACb,cAAmC,UAAU;QAE7C,MAAM,CAAC,WAAW,CAAC,CAAC;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CACvB,WAAW,KAAK,UAAU;YACxB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC;YACvB,CAAC,CAAC,WAAW,KAAK,YAAY;gBAC5B,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;gBAClB,CAAC,CAAC,CAAC,CACR,CAAC;QACF,OAAO;YACL,QAAQ,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE;YACpE,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE;SACnF,CAAC;IACJ,CAAC;IACD,oBAAoB,CAAC,QAAa,EAAE,OAAe;QACjD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,0BAA0B,CAAC,CAAC;QAChF,MAAM,EAAC,KAAK,EAAE,IAAI,EAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,EACxC,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;QACpB,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,CAAC;IACD,sBAAsB,CAAC,QAAa,EAAE,SAAc;QAClD,MAAM,CAAC,SAAS,CAAC,CAAC;QAClB,MAAM,EAAC,KAAK,EAAE,IAAI,EAAE,MAAM,EAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,EAChD,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC;QACtB,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK;YAC/D,MAAM,IAAI,UAAU,CAAC,2BAA2B,CAAC,CAAC;QACpD,OAAO,CACL,CAAC,IAAI,CAAC,KAAK,CACT,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EACzC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAC3C;YACC,GAAG,CAAC;YACN,IAAI,CAAC,EAAE,CACR,CAAC;IACJ,CAAC;IACO,OAAO,CAAC,QAAa;QAC3B,MAAM,CAAC,QAAQ,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACpD,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACvE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK;YAC9D,MAAM,IAAI,UAAU,CAAC,yCAAyC,CAAC,CAAC;QAClE,OAAO,CAAC,CAAC;IACX,CAAC;IACO,KAAK,CAAC,QAAa;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,EAC9B,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EACvC,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACjF,MAAM,KAAK,GAAG,CAAC,CAAM,EAAE,CAAM,EAAE,EAAE,CAAC;YAChC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACzB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;SAC1B,CAAC;QACF,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,EACzC,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC;QACpC,OAAO,EAAC,KAAK,EAAE,IAAI,EAAE,MAAM,EAAC,CAAC;IAC/B,CAAC;IACD,+EAA+E;IAC/E,IAAI,CAAC,QAAa,EAAE,SAAc;QAChC,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QACjC,IAAI,CAAC,sBAAsB,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;QACjD,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,EACvB,OAAO,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EACrC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EACnC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,OAAO;YACL,QAAQ,EAAE,CAAC,QAAgB,EAAE,EAAE;gBAC7B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK;oBACxE,MAAM,IAAI,UAAU,CAAC,gDAAgD,CAAC,CAAC;gBACzE,MAAM,CAAC,GACL,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,iBAAiB,CAAC,OAAO,EAAE,CAAC,EAAE,QAAQ,GAAG,KAAK,CAAC,CAAC;gBACzF,OAAO;oBACL,QAAQ,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAA6B;oBACvE,SAAS,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAA6B;oBACjD,QAAQ;iBACT,CAAC;YACJ,CAAC;SACF,CAAC;IACJ,CAAC;IACD,MAAM,CAAC,QAAa,EAAE,SAAc,EAAE,QAAgB;QACpD,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC3D,CAAC;CACF"}

package/dist/vendor/geographiclib-geodesic/runtime.d.ts

@@ -0,0 +1,9 @@
+export namespace geodesic {
+    let Constants: {};
+    let Math: {};
+    let Accumulator: {};
+    let Geodesic: {};
+    let GeodesicLine: {};
+    let PolygonArea: {};
+}
+//# sourceMappingURL=runtime.d.ts.map

package/dist/vendor/geographiclib-geodesic/runtime.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"runtime.d.ts","sourceRoot":"","sources":["../../../src/vendor/geographiclib-geodesic/runtime.js"],"names":[],"mappings":""}