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

package/dist/spheroid.cjs

@@ -23,10 +23,11 @@ __export(spheroid_exports, {
   spheroidToCartesian: () => spheroidToCartesian
 });
 module.exports = __toCommonJS(spheroid_exports);
+var EPSILON = Number.EPSILON;
 function spheroidToCartesian(point, geometry, result = point) {
   const longitude = point.x, latitude = point.y, height = point.z;
   const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
-  const sin = Math.sin(latitude), cos = Math.cos(latitude);
+  const sin = Math.sin(latitude), cos = Math.abs(latitude) === Math.PI / 2 ? 0 : Math.cos(latitude);
   const ratio = b / a;
   const radius = es < 0.9 ? a / Math.sqrt(1 - es * sin * sin) : a / Math.hypot(cos, ratio * sin);
   const x = (radius + height) * cos * Math.cos(longitude);
@@ -37,18 +38,19 @@ function spheroidToCartesian(point, geometry, result = point) {
 function cartesianToSpheroid(point, geometry, result = point) {
   const x = point.x, y = point.y, z = point.z;
   const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
+  const longitude = Math.atan2(y, x);
   const p = Math.hypot(x, y), rr = Math.hypot(x, y, z);
   if (!(rr > 0) || !Number.isFinite(rr))
     return false;
   if (es === 0)
-    return commit(result, p === 0 ? 0 : Math.atan2(y, x), Math.atan2(z, p), rr - a);
+    return commit(result, p === 0 ? 0 : longitude, Math.atan2(z, p), rr - a);
   if (p === 0)
     return commit(result, 0, Math.sign(z) * Math.PI / 2, Math.abs(z) - b);
   const horizontal = p / a, vertical = z / b;
-  if (horizontal * horizontal + vertical * vertical < 1 - 16 * Number.EPSILON)
+  if (horizontal * horizontal + vertical * vertical < 1 - 16 * EPSILON)
     return interiorInverse(x, y, z, p, a, b, result);
   if (es >= 0.9)
-    return exteriorInverse(x, y, z, p, a, b, result);
+    return interiorInverse(x, y, z, p, a, b, result);
   const ct = z / rr, st = p / rr;
   let rx = 1 / Math.sqrt(1 - es * (2 - es) * st * st);
   let cos = st * (1 - es) * rx, sin = ct * rx;
@@ -59,64 +61,65 @@ function cartesianToSpheroid(point, geometry, result = point) {
     rx = 1 / Math.sqrt(1 - rk * (2 - rk) * st * st);
     const nextCos = st * (1 - rk) * rx, nextSin = ct * rx;
     if (Math.abs(nextSin * cos - nextCos * sin) <= 1e-12) {
-      if (!Number.isFinite(height) || !Number.isFinite(nextSin) || !Number.isFinite(nextCos))
+      if (!Number.isFinite(height))
         break;
-      return commit(result, Math.atan2(y, x), Math.atan2(nextSin, Math.abs(nextCos)), height);
+      return commit(result, longitude, Math.atan2(nextSin, Math.abs(nextCos)), height);
     }
     cos = nextCos;
     sin = nextSin;
   }
-  return exteriorInverse(x, y, z, p, a, b, result);
+  return interiorInverse(x, y, z, p, a, b, result);
 }
 function interiorInverse(x, y, z, p, a, b, result) {
+  const longitude = Math.atan2(y, x);
   const ratio = b / a, ratio2 = ratio * ratio;
   const horizontal = p / a, vertical = z / a, delta = 1 - ratio2;
   if (!(ratio2 > 0) || ratio > 1)
     return false;
-  if (z === 0) {
-    if (horizontal < delta)
-      return false;
-    return commit(result, Math.atan2(y, x), 0, p - a);
-  }
   const aa = a * a, bb = b * b;
   const ratioError = Number.isFinite(aa) && aa > 0 ? (bb - ratio2 * aa + productError(b, b, bb) - ratio2 * productError(a, a, aa) - productError(ratio2, aa, ratio2 * aa)) / aa : 0;
   const pp = p * p, xx = x * x, yy = y * y, sum = xx + yy;
   const radiusError = Number.isFinite(pp) && pp > 0 ? (sum - pp + sumError(xx, yy, sum) + productError(x, x, xx) + productError(y, y, yy) - productError(p, p, pp)) / (2 * p) : 0;
-  const shift = a * ratio2, base = p - a + shift;
-  const correction = radiusError + sumError(p - a, shift, base) + productError(a, ratio2, shift) + a * ratioError;
+  const shift = a * ratio2, distance = p - a, base = distance + shift;
+  const correction = radiusError + sumError(p, -a, distance) + sumError(distance, shift, base) + productError(a, ratio2, shift) + a * ratioError;
+  if (z === 0) {
+    if (base + correction < 0)
+      return false;
+    return commit(result, longitude, 0, distance);
+  }
   const weighted = Math.abs(vertical) * ratio;
-  if (!(weighted > 0))
+  let lower = 0, upper = ratio2 + Math.hypot(horizontal, vertical * ratio);
+  if (!Number.isFinite(upper))
     return false;
-  let lower = 0, upper = ratio2;
   let t = Math.max(weighted, horizontal - delta);
   if (!(t > lower && t < upper))
     t = upper / 2;
-  for (let i = 0; i < 128; i++) {
+  const gap = -(base + correction) / a;
+  const limiting = Math.sqrt(gap * (horizontal + delta)) / delta;
+  let latitude = gap > 0 && Math.abs(z) / limiting / a * ratio <= EPSILON * gap ? Math.atan2(Math.sign(z) * limiting / ratio, horizontal / delta) : void 0;
+  for (let i = 0; i < 1024 && latitude === void 0; i++) {
     const equatorial = delta + t;
-    const u = horizontal / equatorial, v = weighted / t;
+    const u = horizontal / equatorial, verticalAt = z / t / a, v = verticalAt * ratio;
     const at = a * t, difference = base - at;
     const offset = (difference + (sumError(base, -at, difference) + correction - productError(a, t, at))) / (a * equatorial);
     const term = offset * (u + 1);
     const residual = term + v * v;
-    if (Math.abs(residual) <= 4 * Number.EPSILON * (Math.abs(term) + v * v) || upper - lower <= 4 * Number.EPSILON * t) {
-      const latitude = Math.atan2(vertical / t, horizontal / equatorial);
-      const sin = Math.sin(latitude), cos = Math.cos(latitude);
-      return commit(result, Math.atan2(y, x), latitude, p * cos + z * sin - a * Math.hypot(cos, ratio * sin));
-    }
+    if (!Number.isFinite(residual))
+      return false;
     if (residual > 0)
       lower = t;
     else
       upper = t;
     const derivative = -2 * (u * u / equatorial + v * v / t);
     const next = t - residual / derivative;
-    const candidate = Number.isFinite(next) && next > lower && next < upper ? next : lower + (upper - lower) / 2;
-    if (candidate === t) {
-      const latitude = Math.atan2(vertical / t, horizontal / equatorial);
-      return commit(result, Math.atan2(y, x), latitude, p * Math.cos(latitude) + z * Math.sin(latitude) - a * Math.hypot(Math.cos(latitude), ratio * Math.sin(latitude)));
+    const candidate = next > lower && next < upper ? next : lower + (upper - lower) / 2;
+    if (Math.abs(residual) <= 4 * EPSILON * (Math.abs(term) + v * v) || upper - lower <= 4 * EPSILON * t) {
+      latitude = Math.atan2(verticalAt, u);
+      break;
     }
     t = candidate;
   }
-  return false;
+  return latitude !== void 0 && commitNormal(result, longitude, latitude, p, z, a, ratio);
 }
 function sumError(x, y, sum) {
   const part = sum - x;
@@ -129,35 +132,9 @@ function productError(x, y, product) {
   const hx = sx - (sx - x), hy = sy - (sy - y), lx = x - hx, ly = y - hy;
   return hx * hy - product + hx * ly + lx * hy + lx * ly;
 }
-function exteriorInverse(x, y, z, p, a, b, result) {
-  const ratio = b / a, ratio2 = ratio * ratio, horizontal = p / a, vertical = z / a;
-  const norm = Math.hypot(horizontal, z / b);
-  if (!Number.isFinite(norm) || norm < 1 - 8 * Number.EPSILON || !(ratio2 > 0))
-    return false;
-  let lower = 0, upper = Math.hypot(horizontal, vertical * ratio);
-  if (!Number.isFinite(upper))
-    return false;
-  let multiplier = norm <= 1 + 8 * Number.EPSILON ? 0 : upper / 2;
-  for (let i = 0; i < 96; i++) {
-    const equatorial = 1 + multiplier, polar = ratio2 + multiplier;
-    const u = horizontal / equatorial, v = vertical * ratio / polar;
-    const length = Math.hypot(u, v);
-    const residual = length - 1;
-    if (Math.abs(residual) <= 4 * Number.EPSILON) {
-      const latitude = Math.atan2(vertical / polar, horizontal / equatorial);
-      const sin = Math.sin(latitude), cos = Math.cos(latitude);
-      const height = p * cos + z * sin - a * Math.hypot(cos, ratio * sin);
-      return commit(result, Math.atan2(y, x), latitude, height);
-    }
-    if (residual > 0)
-      lower = multiplier;
-    else
-      upper = multiplier;
-    const derivative = -(u * u / equatorial + v * v / polar) / length;
-    const next = multiplier - residual / derivative;
-    multiplier = Number.isFinite(next) && next > lower && next < upper ? next : (lower + upper) / 2;
-  }
-  return false;
+function commitNormal(result, longitude, latitude, p, z, a, ratio) {
+  const sin = Math.sin(latitude), cos = Math.cos(latitude);
+  return commit(result, longitude, latitude, p * cos + z * sin - a * Math.hypot(cos, ratio * sin));
 }
 function commit(result, x, y, z) {
   if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))

package/dist/spheroid.cjs.map

@@ -1,6 +1,6 @@
 {
   "version": 3,
   "sources": ["../src/spheroid.ts"],
-  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;;;;;;AAYM,SAAU,oBACd,OACA,UACA,SAAwB,OAAK;AAE7B,QAAM,YAAY,MAAM,GACtB,WAAW,MAAM,GACjB,SAAS,MAAM;AACjB,QAAM,EAAC,eAAe,GAAG,eAAe,GAAG,qBAAqB,GAAE,IAAI;AACtE,QAAM,MAAM,KAAK,IAAI,QAAQ,GAC3B,MAAM,KAAK,IAAI,QAAQ;AAEzB,QAAM,QAAQ,IAAI;AAClB,QAAM,SAAS,KAAK,MAAM,IAAI,KAAK,KAAK,IAAI,KAAK,MAAM,GAAG,IAAI,IAAI,KAAK,MAAM,KAAK,QAAQ,GAAG;AAC7F,QAAM,KAAK,SAAS,UAAU,MAAM,KAAK,IAAI,SAAS;AACtD,QAAM,KAAK,SAAS,UAAU,MAAM,KAAK,IAAI,SAAS;AACtD,QAAM,KAAK,UAAU,KAAK,MAAM,IAAI,KAAK,QAAQ,SAAS,UAAU;AACpE,SAAO,OAAO,QAAQ,GAAG,GAAG,CAAC;AAC/B;AASM,SAAU,oBACd,OACA,UACA,SAAwB,OAAK;AAE7B,QAAM,IAAI,MAAM,GACd,IAAI,MAAM,GACV,IAAI,MAAM;AACZ,QAAM,EAAC,eAAe,GAAG,eAAe,GAAG,qBAAqB,GAAE,IAAI;AACtE,QAAM,IAAI,KAAK,MAAM,GAAG,CAAC,GACvB,KAAK,KAAK,MAAM,GAAG,GAAG,CAAC;AACzB,MAAI,EAAE,KAAK,MAAM,CAAC,OAAO,SAAS,EAAE;AAAG,WAAO;AAC9C,MAAI,OAAO;AAAG,WAAO,OAAO,QAAQ,MAAM,IAAI,IAAI,KAAK,MAAM,GAAG,CAAC,GAAG,KAAK,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;AAC5F,MAAI,MAAM;AAAG,WAAO,OAAO,QAAQ,GAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAM,GAAG,KAAK,IAAI,CAAC,IAAI,CAAC;AACnF,QAAM,aAAa,IAAI,GACrB,WAAW,IAAI;AACjB,MAAI,aAAa,aAAa,WAAW,WAAW,IAAI,KAAK,OAAO;AAClE,WAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AACjD,MAAI,MAAM;AAAK,WAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AAC9D,QAAM,KAAK,IAAI,IACb,KAAK,IAAI;AACX,MAAI,KAAK,IAAI,KAAK,KAAK,IAAI,MAAM,IAAI,MAAM,KAAK,EAAE;AAClD,MAAI,MAAM,MAAM,IAAI,MAAM,IACxB,MAAM,KAAK;AACb,WAAS,IAAI,GAAG,IAAI,IAAI,KAAK;AAC3B,UAAM,KAAK,IAAI,KAAK,KAAK,IAAI,KAAK,MAAM,GAAG;AAC3C,UAAM,SAAS,IAAI,MAAM,IAAI,MAAM,MAAM,IAAI,KAAK,MAAM;AACxD,UAAM,KAAM,KAAK,MAAO,KAAK;AAC7B,SAAK,IAAI,KAAK,KAAK,IAAI,MAAM,IAAI,MAAM,KAAK,EAAE;AAC9C,UAAM,UAAU,MAAM,IAAI,MAAM,IAC9B,UAAU,KAAK;AACjB,QAAI,KAAK,IAAI,UAAU,MAAM,UAAU,GAAG,KAAK,OAAO;AACpD,UAAI,CAAC,OAAO,SAAS,MAAM,KAAK,CAAC,OAAO,SAAS,OAAO,KAAK,CAAC,OAAO,SAAS,OAAO;AAAG;AACxF,aAAO,OAAO,QAAQ,KAAK,MAAM,GAAG,CAAC,GAAG,KAAK,MAAM,SAAS,KAAK,IAAI,OAAO,CAAC,GAAG,MAAM;IACxF;AACA,UAAM;AACN,UAAM;EACR;AACA,SAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AACjD;AAKA,SAAS,gBACP,GACA,GACA,GACA,GACA,GACA,GACA,QAAqB;AAErB,QAAM,QAAQ,IAAI,GAChB,SAAS,QAAQ;AACnB,QAAM,aAAa,IAAI,GACrB,WAAW,IAAI,GACf,QAAQ,IAAI;AACd,MAAI,EAAE,SAAS,MAAM,QAAQ;AAAG,WAAO;AACvC,MAAI,MAAM,GAAG;AACX,QAAI,aAAa;AAAO,aAAO;AAC/B,WAAO,OAAO,QAAQ,KAAK,MAAM,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC;EAClD;AACA,QAAM,KAAK,IAAI,GACb,KAAK,IAAI;AACX,QAAM,aACJ,OAAO,SAAS,EAAE,KAAK,KAAK,KACvB,KACC,SAAS,KACT,aAAa,GAAG,GAAG,EAAE,IACrB,SAAS,aAAa,GAAG,GAAG,EAAE,IAC9B,aAAa,QAAQ,IAAI,SAAS,EAAE,KACtC,KACA;AACN,QAAM,KAAK,IAAI,GACb,KAAK,IAAI,GACT,KAAK,IAAI,GACT,MAAM,KAAK;AACb,QAAM,cACJ,OAAO,SAAS,EAAE,KAAK,KAAK,KACvB,MACC,KACA,SAAS,IAAI,IAAI,GAAG,IACpB,aAAa,GAAG,GAAG,EAAE,IACrB,aAAa,GAAG,GAAG,EAAE,IACrB,aAAa,GAAG,GAAG,EAAE,MACtB,IAAI,KACL;AACN,QAAM,QAAQ,IAAI,QAChB,OAAO,IAAI,IAAI;AACjB,QAAM,aACJ,cAAc,SAAS,IAAI,GAAG,OAAO,IAAI,IAAI,aAAa,GAAG,QAAQ,KAAK,IAAI,IAAI;AACpF,QAAM,WAAW,KAAK,IAAI,QAAQ,IAAI;AACtC,MAAI,EAAE,WAAW;AAAI,WAAO;AAC5B,MAAI,QAAQ,GACV,QAAQ;AAGV,MAAI,IAAI,KAAK,IAAI,UAAU,aAAa,KAAK;AAC7C,MAAI,EAAE,IAAI,SAAS,IAAI;AAAQ,QAAI,QAAQ;AAC3C,WAAS,IAAI,GAAG,IAAI,KAAK,KAAK;AAC5B,UAAM,aAAa,QAAQ;AAC3B,UAAM,IAAI,aAAa,YACrB,IAAI,WAAW;AACjB,UAAM,KAAK,IAAI,GACb,aAAa,OAAO;AACtB,UAAM,UACH,cAAc,SAAS,MAAM,CAAC,IAAI,UAAU,IAAI,aAAa,aAAa,GAAG,GAAG,EAAE,OAClF,IAAI;AACP,UAAM,OAAO,UAAU,IAAI;AAC3B,UAAM,WAAW,OAAO,IAAI;AAC5B,QACE,KAAK,IAAI,QAAQ,KAAK,IAAI,OAAO,WAAW,KAAK,IAAI,IAAI,IAAI,IAAI,MACjE,QAAQ,SAAS,IAAI,OAAO,UAAU,GACtC;AACA,YAAM,WAAW,KAAK,MAAM,WAAW,GAAG,aAAa,UAAU;AACjE,YAAM,MAAM,KAAK,IAAI,QAAQ,GAC3B,MAAM,KAAK,IAAI,QAAQ;AACzB,aAAO,OACL,QACA,KAAK,MAAM,GAAG,CAAC,GACf,UACA,IAAI,MAAM,IAAI,MAAM,IAAI,KAAK,MAAM,KAAK,QAAQ,GAAG,CAAC;IAExD;AACA,QAAI,WAAW;AAAG,cAAQ;;AACrB,cAAQ;AACb,UAAM,aAAa,MAAO,IAAI,IAAK,aAAc,IAAI,IAAK;AAC1D,UAAM,OAAO,IAAI,WAAW;AAC5B,UAAM,YACJ,OAAO,SAAS,IAAI,KAAK,OAAO,SAAS,OAAO,QAAQ,OAAO,SAAS,QAAQ,SAAS;AAC3F,QAAI,cAAc,GAAG;AACnB,YAAM,WAAW,KAAK,MAAM,WAAW,GAAG,aAAa,UAAU;AACjE,aAAO,OACL,QACA,KAAK,MAAM,GAAG,CAAC,GACf,UACA,IAAI,KAAK,IAAI,QAAQ,IACnB,IAAI,KAAK,IAAI,QAAQ,IACrB,IAAI,KAAK,MAAM,KAAK,IAAI,QAAQ,GAAG,QAAQ,KAAK,IAAI,QAAQ,CAAC,CAAC;IAEpE;AACA,QAAI;EACN;AACA,SAAO;AACT;AAGA,SAAS,SAAS,GAAW,GAAW,KAAW;AACjD,QAAM,OAAO,MAAM;AACnB,SAAO,KAAK,MAAM,SAAS,IAAI;AACjC;AACA,SAAS,aAAa,GAAW,GAAW,SAAe;AACzD,QAAM,KAAK,YAAY,GACrB,KAAK,YAAY;AACnB,MAAI,CAAC,OAAO,SAAS,EAAE,KAAK,CAAC,OAAO,SAAS,EAAE;AAAG,WAAO;AACzD,QAAM,KAAK,MAAM,KAAK,IACpB,KAAK,MAAM,KAAK,IAChB,KAAK,IAAI,IACT,KAAK,IAAI;AACX,SAAO,KAAK,KAAK,UAAU,KAAK,KAAK,KAAK,KAAK,KAAK;AACtD;AAGA,SAAS,gBACP,GACA,GACA,GACA,GACA,GACA,GACA,QAAqB;AAErB,QAAM,QAAQ,IAAI,GAChB,SAAS,QAAQ,OACjB,aAAa,IAAI,GACjB,WAAW,IAAI;AACjB,QAAM,OAAO,KAAK,MAAM,YAAY,IAAI,CAAC;AACzC,MAAI,CAAC,OAAO,SAAS,IAAI,KAAK,OAAO,IAAI,IAAI,OAAO,WAAW,EAAE,SAAS;AAAI,WAAO;AACrF,MAAI,QAAQ,GACV,QAAQ,KAAK,MAAM,YAAY,WAAW,KAAK;AACjD,MAAI,CAAC,OAAO,SAAS,KAAK;AAAG,WAAO;AACpC,MAAI,aAAa,QAAQ,IAAI,IAAI,OAAO,UAAU,IAAI,QAAQ;AAC9D,WAAS,IAAI,GAAG,IAAI,IAAI,KAAK;AAC3B,UAAM,aAAa,IAAI,YACrB,QAAQ,SAAS;AACnB,UAAM,IAAI,aAAa,YACrB,IAAK,WAAW,QAAS;AAC3B,UAAM,SAAS,KAAK,MAAM,GAAG,CAAC;AAC9B,UAAM,WAAW,SAAS;AAC1B,QAAI,KAAK,IAAI,QAAQ,KAAK,IAAI,OAAO,SAAS;AAC5C,YAAM,WAAW,KAAK,MAAM,WAAW,OAAO,aAAa,UAAU;AACrE,YAAM,MAAM,KAAK,IAAI,QAAQ,GAC3B,MAAM,KAAK,IAAI,QAAQ;AACzB,YAAM,SAAS,IAAI,MAAM,IAAI,MAAM,IAAI,KAAK,MAAM,KAAK,QAAQ,GAAG;AAClE,aAAO,OAAO,QAAQ,KAAK,MAAM,GAAG,CAAC,GAAG,UAAU,MAAM;IAC1D;AACA,QAAI,WAAW;AAAG,cAAQ;;AACrB,cAAQ;AACb,UAAM,aAAa,EAAG,IAAI,IAAK,aAAc,IAAI,IAAK,SAAS;AAC/D,UAAM,OAAO,aAAa,WAAW;AACrC,iBAAa,OAAO,SAAS,IAAI,KAAK,OAAO,SAAS,OAAO,QAAQ,QAAQ,QAAQ,SAAS;EAChG;AACA,SAAO;AACT;AAGA,SAAS,OAAO,QAAuB,GAAW,GAAW,GAAS;AACpE,MAAI,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC;AAAG,WAAO;AAC9E,SAAO,IAAI;AACX,SAAO,IAAI;AACX,SAAO,IAAI;AACX,SAAO;AACT;",
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;;;;;;AAWA,IAAM,UAAU,OAAO;AAGjB,SAAU,oBACd,OACA,UACA,SAAwB,OAAK;AAE7B,QAAM,YAAY,MAAM,GACtB,WAAW,MAAM,GACjB,SAAS,MAAM;AACjB,QAAM,EAAC,eAAe,GAAG,eAAe,GAAG,qBAAqB,GAAE,IAAI;AACtE,QAAM,MAAM,KAAK,IAAI,QAAQ,GAC3B,MAAM,KAAK,IAAI,QAAQ,MAAM,KAAK,KAAK,IAAI,IAAI,KAAK,IAAI,QAAQ;AAElE,QAAM,QAAQ,IAAI;AAClB,QAAM,SAAS,KAAK,MAAM,IAAI,KAAK,KAAK,IAAI,KAAK,MAAM,GAAG,IAAI,IAAI,KAAK,MAAM,KAAK,QAAQ,GAAG;AAC7F,QAAM,KAAK,SAAS,UAAU,MAAM,KAAK,IAAI,SAAS;AACtD,QAAM,KAAK,SAAS,UAAU,MAAM,KAAK,IAAI,SAAS;AACtD,QAAM,KAAK,UAAU,KAAK,MAAM,IAAI,KAAK,QAAQ,SAAS,UAAU;AACpE,SAAO,OAAO,QAAQ,GAAG,GAAG,CAAC;AAC/B;AASM,SAAU,oBACd,OACA,UACA,SAAwB,OAAK;AAE7B,QAAM,IAAI,MAAM,GACd,IAAI,MAAM,GACV,IAAI,MAAM;AACZ,QAAM,EAAC,eAAe,GAAG,eAAe,GAAG,qBAAqB,GAAE,IAAI;AACtE,QAAM,YAAY,KAAK,MAAM,GAAG,CAAC;AACjC,QAAM,IAAI,KAAK,MAAM,GAAG,CAAC,GACvB,KAAK,KAAK,MAAM,GAAG,GAAG,CAAC;AACzB,MAAI,EAAE,KAAK,MAAM,CAAC,OAAO,SAAS,EAAE;AAAG,WAAO;AAC9C,MAAI,OAAO;AAAG,WAAO,OAAO,QAAQ,MAAM,IAAI,IAAI,WAAW,KAAK,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;AACrF,MAAI,MAAM;AAAG,WAAO,OAAO,QAAQ,GAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAM,GAAG,KAAK,IAAI,CAAC,IAAI,CAAC;AACnF,QAAM,aAAa,IAAI,GACrB,WAAW,IAAI;AACjB,MAAI,aAAa,aAAa,WAAW,WAAW,IAAI,KAAK;AAC3D,WAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AACjD,MAAI,MAAM;AAAK,WAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AAC9D,QAAM,KAAK,IAAI,IACb,KAAK,IAAI;AACX,MAAI,KAAK,IAAI,KAAK,KAAK,IAAI,MAAM,IAAI,MAAM,KAAK,EAAE;AAClD,MAAI,MAAM,MAAM,IAAI,MAAM,IACxB,MAAM,KAAK;AACb,WAAS,IAAI,GAAG,IAAI,IAAI,KAAK;AAC3B,UAAM,KAAK,IAAI,KAAK,KAAK,IAAI,KAAK,MAAM,GAAG;AAC3C,UAAM,SAAS,IAAI,MAAM,IAAI,MAAM,MAAM,IAAI,KAAK,MAAM;AACxD,UAAM,KAAM,KAAK,MAAO,KAAK;AAC7B,SAAK,IAAI,KAAK,KAAK,IAAI,MAAM,IAAI,MAAM,KAAK,EAAE;AAC9C,UAAM,UAAU,MAAM,IAAI,MAAM,IAC9B,UAAU,KAAK;AACjB,QAAI,KAAK,IAAI,UAAU,MAAM,UAAU,GAAG,KAAK,OAAO;AACpD,UAAI,CAAC,OAAO,SAAS,MAAM;AAAG;AAC9B,aAAO,OAAO,QAAQ,WAAW,KAAK,MAAM,SAAS,KAAK,IAAI,OAAO,CAAC,GAAG,MAAM;IACjF;AACA,UAAM;AACN,UAAM;EACR;AACA,SAAO,gBAAgB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM;AACjD;AAKA,SAAS,gBACP,GACA,GACA,GACA,GACA,GACA,GACA,QAAqB;AAErB,QAAM,YAAY,KAAK,MAAM,GAAG,CAAC;AACjC,QAAM,QAAQ,IAAI,GAChB,SAAS,QAAQ;AACnB,QAAM,aAAa,IAAI,GACrB,WAAW,IAAI,GACf,QAAQ,IAAI;AACd,MAAI,EAAE,SAAS,MAAM,QAAQ;AAAG,WAAO;AACvC,QAAM,KAAK,IAAI,GACb,KAAK,IAAI;AACX,QAAM,aACJ,OAAO,SAAS,EAAE,KAAK,KAAK,KACvB,KACC,SAAS,KACT,aAAa,GAAG,GAAG,EAAE,IACrB,SAAS,aAAa,GAAG,GAAG,EAAE,IAC9B,aAAa,QAAQ,IAAI,SAAS,EAAE,KACtC,KACA;AACN,QAAM,KAAK,IAAI,GACb,KAAK,IAAI,GACT,KAAK,IAAI,GACT,MAAM,KAAK;AACb,QAAM,cACJ,OAAO,SAAS,EAAE,KAAK,KAAK,KACvB,MACC,KACA,SAAS,IAAI,IAAI,GAAG,IACpB,aAAa,GAAG,GAAG,EAAE,IACrB,aAAa,GAAG,GAAG,EAAE,IACrB,aAAa,GAAG,GAAG,EAAE,MACtB,IAAI,KACL;AACN,QAAM,QAAQ,IAAI,QAChB,WAAW,IAAI,GACf,OAAO,WAAW;AACpB,QAAM,aACJ,cACA,SAAS,GAAG,CAAC,GAAG,QAAQ,IACxB,SAAS,UAAU,OAAO,IAAI,IAC9B,aAAa,GAAG,QAAQ,KAAK,IAC7B,IAAI;AACN,MAAI,MAAM,GAAG;AAGX,QAAI,OAAO,aAAa;AAAG,aAAO;AAClC,WAAO,OAAO,QAAQ,WAAW,GAAG,QAAQ;EAC9C;AACA,QAAM,WAAW,KAAK,IAAI,QAAQ,IAAI;AACtC,MAAI,QAAQ,GACV,QAAQ,SAAS,KAAK,MAAM,YAAY,WAAW,KAAK;AAC1D,MAAI,CAAC,OAAO,SAAS,KAAK;AAAG,WAAO;AAGpC,MAAI,IAAI,KAAK,IAAI,UAAU,aAAa,KAAK;AAC7C,MAAI,EAAE,IAAI,SAAS,IAAI;AAAQ,QAAI,QAAQ;AAI3C,QAAM,MAAM,EAAE,OAAO,cAAc;AACnC,QAAM,WAAW,KAAK,KAAK,OAAO,aAAa,MAAM,IAAI;AACzD,MAAI,WACF,MAAM,KAAM,KAAK,IAAI,CAAC,IAAI,WAAW,IAAK,SAAS,UAAU,MACzD,KAAK,MAAO,KAAK,KAAK,CAAC,IAAI,WAAY,OAAO,aAAa,KAAK,IAChE;AACN,WAAS,IAAI,GAAG,IAAI,QAAQ,aAAa,QAAW,KAAK;AACvD,UAAM,aAAa,QAAQ;AAC3B,UAAM,IAAI,aAAa,YAGrB,aAAa,IAAI,IAAI,GACrB,IAAI,aAAa;AACnB,UAAM,KAAK,IAAI,GACb,aAAa,OAAO;AACtB,UAAM,UACH,cAAc,SAAS,MAAM,CAAC,IAAI,UAAU,IAAI,aAAa,aAAa,GAAG,GAAG,EAAE,OAClF,IAAI;AACP,UAAM,OAAO,UAAU,IAAI;AAC3B,UAAM,WAAW,OAAO,IAAI;AAC5B,QAAI,CAAC,OAAO,SAAS,QAAQ;AAAG,aAAO;AACvC,QAAI,WAAW;AAAG,cAAQ;;AACrB,cAAQ;AACb,UAAM,aAAa,MAAO,IAAI,IAAK,aAAc,IAAI,IAAK;AAC1D,UAAM,OAAO,IAAI,WAAW;AAC5B,UAAM,YAAY,OAAO,SAAS,OAAO,QAAQ,OAAO,SAAS,QAAQ,SAAS;AAClF,QACE,KAAK,IAAI,QAAQ,KAAK,IAAI,WAAW,KAAK,IAAI,IAAI,IAAI,IAAI,MAC1D,QAAQ,SAAS,IAAI,UAAU,GAC/B;AACA,iBAAW,KAAK,MAAM,YAAY,CAAC;AACnC;IACF;AACA,QAAI;EACN;AACA,SAAO,aAAa,UAAa,aAAa,QAAQ,WAAW,UAAU,GAAG,GAAG,GAAG,KAAK;AAC3F;AAGA,SAAS,SAAS,GAAW,GAAW,KAAW;AACjD,QAAM,OAAO,MAAM;AACnB,SAAO,KAAK,MAAM,SAAS,IAAI;AACjC;AACA,SAAS,aAAa,GAAW,GAAW,SAAe;AACzD,QAAM,KAAK,YAAY,GACrB,KAAK,YAAY;AACnB,MAAI,CAAC,OAAO,SAAS,EAAE,KAAK,CAAC,OAAO,SAAS,EAAE;AAAG,WAAO;AACzD,QAAM,KAAK,MAAM,KAAK,IACpB,KAAK,MAAM,KAAK,IAChB,KAAK,IAAI,IACT,KAAK,IAAI;AACX,SAAO,KAAK,KAAK,UAAU,KAAK,KAAK,KAAK,KAAK,KAAK;AACtD;AAGA,SAAS,aACP,QACA,WACA,UACA,GACA,GACA,GACA,OAAa;AAEb,QAAM,MAAM,KAAK,IAAI,QAAQ,GAC3B,MAAM,KAAK,IAAI,QAAQ;AACzB,SAAO,OAAO,QAAQ,WAAW,UAAU,IAAI,MAAM,IAAI,MAAM,IAAI,KAAK,MAAM,KAAK,QAAQ,GAAG,CAAC;AACjG;AAGA,SAAS,OAAO,QAAuB,GAAW,GAAW,GAAS;AACpE,MAAI,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC;AAAG,WAAO;AAC9E,SAAO,IAAI;AACX,SAAO,IAAI;AACX,SAAO,IAAI;AACX,SAAO;AACT;",
   "names": []
 }

package/dist/spheroid.d.ts

@@ -12,7 +12,7 @@ export type SpheroidPoint = {
 export declare function spheroidToCartesian(point: SpheroidPoint, geometry: SpheroidGeometry, result?: SpheroidPoint): boolean;
 /**
  * Bounded inverse. Sphere is analytic; ordinary spheroids use the retained Hannover
- * iteration on exterior points. Interiors use a safeguarded nearest-footpoint solve.
+ * iteration on exterior points. Interiors and thin-shape exteriors share a safeguarded nearest-footpoint solve.
  * Equatorial interiors below the evolute cusp have two nearest normals: reject.
  * false leaves output untouched, including the center and unresolved binary64 cases.
  * Only exact poles canonicalize longitude; adapters may retain other conventions.

package/dist/spheroid.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"spheroid.d.ts","sourceRoot":"","sources":["../src/spheroid.ts"],"names":[],"mappings":"AAMA,kGAAkG;AAClG,OAAO,KAAK,EAAC,kBAAkB,EAAC,mBAAgB;AAChD,MAAM,MAAM,gBAAgB,GAAG,kBAAkB,GAAG;IAAC,QAAQ,CAAC,mBAAmB,EAAE,MAAM,CAAA;CAAC,CAAC;AAC3F,MAAM,MAAM,aAAa,GAAG;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAC,CAAC;AAE9D,6FAA6F;AAC7F,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,aAAa,EACpB,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,GAAE,aAAqB,GAC5B,OAAO,CAcT;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,aAAa,EACpB,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,GAAE,aAAqB,GAC5B,OAAO,CAmCT"}
+{"version":3,"file":"spheroid.d.ts","sourceRoot":"","sources":["../src/spheroid.ts"],"names":[],"mappings":"AAMA,kGAAkG;AAClG,OAAO,KAAK,EAAC,kBAAkB,EAAC,mBAAgB;AAChD,MAAM,MAAM,gBAAgB,GAAG,kBAAkB,GAAG;IAAC,QAAQ,CAAC,mBAAmB,EAAE,MAAM,CAAA;CAAC,CAAC;AAC3F,MAAM,MAAM,aAAa,GAAG;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAC,CAAC;AAI9D,6FAA6F;AAC7F,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,aAAa,EACpB,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,GAAE,aAAqB,GAC5B,OAAO,CAcT;AAED;;;;;;GAMG;AACH,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,aAAa,EACpB,QAAQ,EAAE,gBAAgB,EAC1B,MAAM,GAAE,aAAqB,GAC5B,OAAO,CAoCT"}

package/dist/spheroid.js

@@ -2,12 +2,13 @@
 // 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: Forward geocentric equations and Hannover inverse moved/adapted from math.gl projection's proj4js 2.22.0 datumUtils.js port. See ../PROJ4-LICENSE.md. Safeguarded nearest-footpoint inverse and numeric output commit are original math.gl code.
+// SPDX-FileComment: Forward geocentric equations and Hannover inverse moved/adapted from math.gl projection's proj4js 2.22.0 datumUtils.js port. See ../PROJ4-LICENSE.md. Safeguarded compensated nearest-footpoint inverse, thin-exterior/pole fixes and numeric output commit are original math.gl code, qualified with independent mathematics and offline GeographicLib 2.7 (Charles Karney, MIT).
+const EPSILON = Number.EPSILON;
 /** Caller supplies validated positive finite axes and eccentricity. No latitude wrapping. */
 export function spheroidToCartesian(point, geometry, result = point) {
     const longitude = point.x, latitude = point.y, height = point.z;
     const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
-    const sin = Math.sin(latitude), cos = Math.cos(latitude);
+    const sin = Math.sin(latitude), cos = Math.abs(latitude) === Math.PI / 2 ? 0 : Math.cos(latitude);
     // Use the axis ratio for very flat shapes: 1 - es*sin² loses polar precision.
     const ratio = b / a;
     const radius = es < 0.9 ? a / Math.sqrt(1 - es * sin * sin) : a / Math.hypot(cos, ratio * sin);
@@ -18,7 +19,7 @@ export function spheroidToCartesian(point, geometry, result = point) {
 }
 /**
  * Bounded inverse. Sphere is analytic; ordinary spheroids use the retained Hannover
- * iteration on exterior points. Interiors use a safeguarded nearest-footpoint solve.
+ * iteration on exterior points. Interiors and thin-shape exteriors share a safeguarded nearest-footpoint solve.
  * Equatorial interiors below the evolute cusp have two nearest normals: reject.
  * false leaves output untouched, including the center and unresolved binary64 cases.
  * Only exact poles canonicalize longitude; adapters may retain other conventions.
@@ -26,18 +27,19 @@ export function spheroidToCartesian(point, geometry, result = point) {
 export function cartesianToSpheroid(point, geometry, result = point) {
     const x = point.x, y = point.y, z = point.z;
     const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
+    const longitude = Math.atan2(y, x);
     const p = Math.hypot(x, y), rr = Math.hypot(x, y, z);
     if (!(rr > 0) || !Number.isFinite(rr))
         return false;
     if (es === 0)
-        return commit(result, p === 0 ? 0 : Math.atan2(y, x), Math.atan2(z, p), rr - a);
+        return commit(result, p === 0 ? 0 : longitude, Math.atan2(z, p), rr - a);
     if (p === 0)
         return commit(result, 0, (Math.sign(z) * Math.PI) / 2, Math.abs(z) - b);
     const horizontal = p / a, vertical = z / b;
-    if (horizontal * horizontal + vertical * vertical < 1 - 16 * Number.EPSILON)
+    if (horizontal * horizontal + vertical * vertical < 1 - 16 * EPSILON)
         return interiorInverse(x, y, z, p, a, b, result);
     if (es >= 0.9)
-        return exteriorInverse(x, y, z, p, a, b, result);
+        return interiorInverse(x, y, z, p, a, b, result);
     const ct = z / rr, st = p / rr;
     let rx = 1 / Math.sqrt(1 - es * (2 - es) * st * st);
     let cos = st * (1 - es) * rx, sin = ct * rx;
@@ -48,28 +50,24 @@ export function cartesianToSpheroid(point, geometry, result = point) {
         rx = 1 / Math.sqrt(1 - rk * (2 - rk) * st * st);
         const nextCos = st * (1 - rk) * rx, nextSin = ct * rx;
         if (Math.abs(nextSin * cos - nextCos * sin) <= 1e-12) {
-            if (!Number.isFinite(height) || !Number.isFinite(nextSin) || !Number.isFinite(nextCos))
+            if (!Number.isFinite(height))
                 break;
-            return commit(result, Math.atan2(y, x), Math.atan2(nextSin, Math.abs(nextCos)), height);
+            return commit(result, longitude, Math.atan2(nextSin, Math.abs(nextCos)), height);
         }
         cos = nextCos;
         sin = nextSin;
     }
-    return exteriorInverse(x, y, z, p, a, b, result);
+    return interiorInverse(x, y, z, p, a, b, result);
 }
-/** Original nearest-normal solve in t = (lambda + b²)/a², avoiding interior cancellation.
+/** Original nearest-normal solve in t = (lambda + b²)/a², avoiding interior cancellation and thin-exterior residual loss.
  * The positive interval selects the globally nearest branch, rather than a stationary normal.
  */
 function interiorInverse(x, y, z, p, a, b, result) {
+    const longitude = Math.atan2(y, x);
     const ratio = b / a, ratio2 = ratio * ratio;
     const horizontal = p / a, vertical = z / a, delta = 1 - ratio2;
     if (!(ratio2 > 0) || ratio > 1)
         return false;
-    if (z === 0) {
-        if (horizontal < delta)
-            return false;
-        return commit(result, Math.atan2(y, x), 0, p - a);
-    }
     const aa = a * a, bb = b * b;
     const ratioError = Number.isFinite(aa) && aa > 0
         ? (bb -
@@ -89,47 +87,64 @@ function interiorInverse(x, y, z, p, a, b, result) {
             productError(p, p, pp)) /
             (2 * p)
         : 0;
-    const shift = a * ratio2, base = p - a + shift;
-    const correction = radiusError + sumError(p - a, shift, base) + productError(a, ratio2, shift) + a * ratioError;
+    const shift = a * ratio2, distance = p - a, base = distance + shift;
+    const correction = radiusError +
+        sumError(p, -a, distance) +
+        sumError(distance, shift, base) +
+        productError(a, ratio2, shift) +
+        a * ratioError;
+    if (z === 0) {
+        // A rounded hypot can equal the cusp for points just inside it. Include
+        // the radius/axis corrections before selecting a unique equatorial normal.
+        if (base + correction < 0)
+            return false;
+        return commit(result, longitude, 0, distance);
+    }
     const weighted = Math.abs(vertical) * ratio;
-    if (!(weighted > 0))
+    let lower = 0, upper = ratio2 + Math.hypot(horizontal, vertical * ratio);
+    if (!Number.isFinite(upper))
         return false;
-    let lower = 0, upper = ratio2;
     // A lower bound near the short-axis singularity prevents thousands of bisections
     // for deep, nearly equatorial interiors. No subtraction of lambda and -b².
     let t = Math.max(weighted, horizontal - delta);
     if (!(t > lower && t < upper))
         t = upper / 2;
-    for (let i = 0; i < 128; i++) {
+    // Below the cusp, a vanishing polar input can make t subnormal. The
+    // limiting vertical footpoint remains well conditioned. Use it only when
+    // its bound t <= abs(z)*ratio/(a*v) is negligible relative to the cusp gap.
+    const gap = -(base + correction) / a;
+    const limiting = Math.sqrt(gap * (horizontal + delta)) / delta;
+    let latitude = gap > 0 && (Math.abs(z) / limiting / a) * ratio <= EPSILON * gap
+        ? Math.atan2((Math.sign(z) * limiting) / ratio, horizontal / delta)
+        : undefined;
+    for (let i = 0; i < 1024 && latitude === undefined; i++) {
         const equatorial = delta + t;
-        const u = horizontal / equatorial, v = weighted / t;
+        const u = horizontal / equatorial, 
+        // Divide by t first: both z/a and a*t can lose precision in the
+        // subnormal range although the near-cusp latitude is representable.
+        verticalAt = z / t / a, v = verticalAt * ratio;
         const at = a * t, difference = base - at;
         const offset = (difference + (sumError(base, -at, difference) + correction - productError(a, t, at))) /
             (a * equatorial);
         const term = offset * (u + 1);
         const residual = term + v * v;
-        if (Math.abs(residual) <= 4 * Number.EPSILON * (Math.abs(term) + v * v) ||
-            upper - lower <= 4 * Number.EPSILON * t) {
-            const latitude = Math.atan2(vertical / t, horizontal / equatorial);
-            const sin = Math.sin(latitude), cos = Math.cos(latitude);
-            return commit(result, Math.atan2(y, x), latitude, p * cos + z * sin - a * Math.hypot(cos, ratio * sin));
-        }
+        if (!Number.isFinite(residual))
+            return false;
         if (residual > 0)
             lower = t;
         else
             upper = t;
         const derivative = -2 * ((u * u) / equatorial + (v * v) / t);
         const next = t - residual / derivative;
-        const candidate = Number.isFinite(next) && next > lower && next < upper ? next : lower + (upper - lower) / 2;
-        if (candidate === t) {
-            const latitude = Math.atan2(vertical / t, horizontal / equatorial);
-            return commit(result, Math.atan2(y, x), latitude, p * Math.cos(latitude) +
-                z * Math.sin(latitude) -
-                a * Math.hypot(Math.cos(latitude), ratio * Math.sin(latitude)));
+        const candidate = next > lower && next < upper ? next : lower + (upper - lower) / 2;
+        if (Math.abs(residual) <= 4 * EPSILON * (Math.abs(term) + v * v) ||
+            upper - lower <= 4 * EPSILON * t) {
+            latitude = Math.atan2(verticalAt, u);
+            break;
         }
         t = candidate;
     }
-    return false;
+    return latitude !== undefined && commitNormal(result, longitude, latitude, p, z, a, ratio);
 }
 /** Error-free transforms on the finite squared-axis domain used by geometry adapters. */
 function sumError(x, y, sum) {
@@ -143,36 +158,10 @@ function productError(x, y, product) {
     const hx = sx - (sx - x), hy = sy - (sy - y), lx = x - hx, ly = y - hy;
     return hx * hy - product + hx * ly + lx * hy + lx * ly;
 }
-/** Original dimensionless Lagrange-multiplier solve, restricted to surface/exterior. */
-function exteriorInverse(x, y, z, p, a, b, result) {
-    const ratio = b / a, ratio2 = ratio * ratio, horizontal = p / a, vertical = z / a;
-    const norm = Math.hypot(horizontal, z / b);
-    if (!Number.isFinite(norm) || norm < 1 - 8 * Number.EPSILON || !(ratio2 > 0))
-        return false;
-    let lower = 0, upper = Math.hypot(horizontal, vertical * ratio);
-    if (!Number.isFinite(upper))
-        return false;
-    let multiplier = norm <= 1 + 8 * Number.EPSILON ? 0 : upper / 2;
-    for (let i = 0; i < 96; i++) {
-        const equatorial = 1 + multiplier, polar = ratio2 + multiplier;
-        const u = horizontal / equatorial, v = (vertical * ratio) / polar;
-        const length = Math.hypot(u, v);
-        const residual = length - 1;
-        if (Math.abs(residual) <= 4 * Number.EPSILON) {
-            const latitude = Math.atan2(vertical / polar, horizontal / equatorial);
-            const sin = Math.sin(latitude), cos = Math.cos(latitude);
-            const height = p * cos + z * sin - a * Math.hypot(cos, ratio * sin);
-            return commit(result, Math.atan2(y, x), latitude, height);
-        }
-        if (residual > 0)
-            lower = multiplier;
-        else
-            upper = multiplier;
-        const derivative = -((u * u) / equatorial + (v * v) / polar) / length;
-        const next = multiplier - residual / derivative;
-        multiplier = Number.isFinite(next) && next > lower && next < upper ? next : (lower + upper) / 2;
-    }
-    return false;
+/** Shared signed distance along the solved surface normal. */
+function commitNormal(result, longitude, latitude, p, z, a, ratio) {
+    const sin = Math.sin(latitude), cos = Math.cos(latitude);
+    return commit(result, longitude, latitude, p * cos + z * sin - a * Math.hypot(cos, ratio * sin));
 }
 /** Values are captured before invoking application setters; recursive writes are safe. */
 function commit(result, x, y, z) {

package/dist/spheroid.js.map

@@ -1 +1 @@
-{"version":3,"file":"spheroid.js","sourceRoot":"","sources":["../src/spheroid.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,qQAAqQ;AAOrQ,6FAA6F;AAC7F,MAAM,UAAU,mBAAmB,CACjC,KAAoB,EACpB,QAA0B,EAC1B,SAAwB,KAAK;IAE7B,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,EACvB,QAAQ,GAAG,KAAK,CAAC,CAAC,EAClB,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC;IACnB,MAAM,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,mBAAmB,EAAE,EAAE,EAAC,GAAG,QAAQ,CAAC;IAC/E,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IAC3B,8EAA8E;IAC9E,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,MAAM,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;IAC/F,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC;IACxE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACjC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CACjC,KAAoB,EACpB,QAA0B,EAC1B,SAAwB,KAAK;IAE7B,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EACf,CAAC,GAAG,KAAK,CAAC,CAAC,EACX,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACd,MAAM,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,mBAAmB,EAAE,EAAE,EAAC,GAAG,QAAQ,CAAC;IAC/E,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACxB,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,KAAK,CAAC;IACpD,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IAC9F,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrF,MAAM,UAAU,GAAG,CAAC,GAAG,CAAC,EACtB,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,IAAI,UAAU,GAAG,UAAU,GAAG,QAAQ,GAAG,QAAQ,GAAG,CAAC,GAAG,EAAE,GAAG,MAAM,CAAC,OAAO;QACzE,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACnD,IAAI,EAAE,IAAI,GAAG;QAAE,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAChE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,EACf,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,IAAI,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACpD,IAAI,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAC1B,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QACrC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAChD,MAAM,OAAO,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAChC,OAAO,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC,IAAI,KAAK,EAAE,CAAC;YACrD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;gBAAE,MAAM;YAC9F,OAAO,MAAM,CAAC,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QAC1F,CAAC;QACD,GAAG,GAAG,OAAO,CAAC;QACd,GAAG,GAAG,OAAO,CAAC;IAChB,CAAC;IACD,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;AACnD,CAAC;AAED;;GAEG;AACH,SAAS,eAAe,CACtB,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAqB;IAErB,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,EACjB,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;IACzB,MAAM,UAAU,GAAG,CAAC,GAAG,CAAC,EACtB,QAAQ,GAAG,CAAC,GAAG,CAAC,EAChB,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC;IACrB,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,KAAK,GAAG,CAAC;QAAE,OAAO,KAAK,CAAC;IAC7C,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,IAAI,UAAU,GAAG,KAAK;YAAE,OAAO,KAAK,CAAC;QACrC,OAAO,MAAM,CAAC,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,CAAC;IACD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,UAAU,GACd,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAC3B,CAAC,CAAC,CAAC,EAAE;YACD,MAAM,GAAG,EAAE;YACX,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,MAAM,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,YAAY,CAAC,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,EAAE,CAAC,CAAC;YACxC,EAAE;QACJ,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,EAAE,GAAG,CAAC,GAAG,CAAC,EACV,EAAE,GAAG,CAAC,GAAG,CAAC,EACV,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;IAChB,MAAM,WAAW,GACf,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAC3B,CAAC,CAAC,CAAC,GAAG;YACF,EAAE;YACF,QAAQ,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC;YACrB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;YACzB,CAAC,CAAC,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,EACtB,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;IACvB,MAAM,UAAU,GACd,WAAW,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC;IAC/F,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,KAAK,CAAC;IAC5C,IAAI,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IAClC,IAAI,KAAK,GAAG,CAAC,EACX,KAAK,GAAG,MAAM,CAAC;IACjB,iFAAiF;IACjF,2EAA2E;IAC3E,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,UAAU,GAAG,KAAK,CAAC,CAAC;IAC/C,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC;QAAE,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;IAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,UAAU,GAAG,KAAK,GAAG,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,UAAU,GAAG,UAAU,EAC/B,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC;QACnB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,UAAU,GAAG,IAAI,GAAG,EAAE,CAAC;QACzB,MAAM,MAAM,GACV,CAAC,UAAU,GAAG,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,UAAU,CAAC,GAAG,UAAU,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;YACtF,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9B,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9B,IACE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACnE,KAAK,GAAG,KAAK,IAAI,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,CAAC,EACvC,CAAC;YACD,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,CAAC,EAAE,UAAU,GAAG,UAAU,CAAC,CAAC;YACnE,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YAC3B,OAAO,MAAM,CACX,MAAM,EACN,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAChB,QAAQ,EACR,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CACrD,CAAC;QACJ,CAAC;QACD,IAAI,QAAQ,GAAG,CAAC;YAAE,KAAK,GAAG,CAAC,CAAC;;YACvB,KAAK,GAAG,CAAC,CAAC;QACf,MAAM,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7D,MAAM,IAAI,GAAG,CAAC,GAAG,QAAQ,GAAG,UAAU,CAAC;QACvC,MAAM,SAAS,GACb,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,KAAK,IAAI,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QAC7F,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;YACpB,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,CAAC,EAAE,UAAU,GAAG,UAAU,CAAC,CAAC;YACnE,OAAO,MAAM,CACX,MAAM,EACN,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAChB,QAAQ,EACR,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;gBACpB,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;gBACtB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CACjE,CAAC;QACJ,CAAC;QACD,CAAC,GAAG,SAAS,CAAC;IAChB,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,yFAAyF;AACzF,SAAS,QAAQ,CAAC,CAAS,EAAE,CAAS,EAAE,GAAW;IACjD,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;IACrB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;AACvC,CAAC;AACD,SAAS,YAAY,CAAC,CAAS,EAAE,CAAS,EAAE,OAAe;IACzD,MAAM,EAAE,GAAG,SAAS,GAAG,CAAC,EACtB,EAAE,GAAG,SAAS,GAAG,CAAC,CAAC;IACrB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,CAAC,CAAC;IAC3D,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EACtB,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EAClB,EAAE,GAAG,CAAC,GAAG,EAAE,EACX,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,OAAO,EAAE,GAAG,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AACzD,CAAC;AAED,wFAAwF;AACxF,SAAS,eAAe,CACtB,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAqB;IAErB,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,EACjB,MAAM,GAAG,KAAK,GAAG,KAAK,EACtB,UAAU,GAAG,CAAC,GAAG,CAAC,EAClB,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,OAAO,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IAC3F,IAAI,KAAK,GAAG,CAAC,EACX,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,EAAE,QAAQ,GAAG,KAAK,CAAC,CAAC;IACnD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC;IAC1C,IAAI,UAAU,GAAG,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAChE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,UAAU,GAAG,CAAC,GAAG,UAAU,EAC/B,KAAK,GAAG,MAAM,GAAG,UAAU,CAAC;QAC9B,MAAM,CAAC,GAAG,UAAU,GAAG,UAAU,EAC/B,CAAC,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,QAAQ,GAAG,MAAM,GAAG,CAAC,CAAC;QAC5B,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC;YAC7C,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,GAAG,KAAK,EAAE,UAAU,GAAG,UAAU,CAAC,CAAC;YACvE,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YAC3B,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;YACpE,OAAO,MAAM,CAAC,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;QAC5D,CAAC;QACD,IAAI,QAAQ,GAAG,CAAC;YAAE,KAAK,GAAG,UAAU,CAAC;;YAChC,KAAK,GAAG,UAAU,CAAC;QACxB,MAAM,UAAU,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC;QACtE,MAAM,IAAI,GAAG,UAAU,GAAG,QAAQ,GAAG,UAAU,CAAC;QAChD,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,KAAK,IAAI,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;IAClG,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,0FAA0F;AAC1F,SAAS,MAAM,CAAC,MAAqB,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;IACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IACpF,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,OAAO,IAAI,CAAC;AACd,CAAC"}
+{"version":3,"file":"spheroid.js","sourceRoot":"","sources":["../src/spheroid.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,uYAAuY;AAOvY,MAAM,OAAO,GAAG,MAAM,CAAC,OAAO,CAAC;AAE/B,6FAA6F;AAC7F,MAAM,UAAU,mBAAmB,CACjC,KAAoB,EACpB,QAA0B,EAC1B,SAAwB,KAAK;IAE7B,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,EACvB,QAAQ,GAAG,KAAK,CAAC,CAAC,EAClB,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC;IACnB,MAAM,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,mBAAmB,EAAE,EAAE,EAAC,GAAG,QAAQ,CAAC;IAC/E,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IACpE,8EAA8E;IAC9E,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,MAAM,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;IAC/F,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC;IACxE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACjC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,mBAAmB,CACjC,KAAoB,EACpB,QAA0B,EAC1B,SAAwB,KAAK;IAE7B,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EACf,CAAC,GAAG,KAAK,CAAC,CAAC,EACX,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACd,MAAM,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,mBAAmB,EAAE,EAAE,EAAC,GAAG,QAAQ,CAAC;IAC/E,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACnC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACxB,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3B,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,KAAK,CAAC;IACpD,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IACvF,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrF,MAAM,UAAU,GAAG,CAAC,GAAG,CAAC,EACtB,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,IAAI,UAAU,GAAG,UAAU,GAAG,QAAQ,GAAG,QAAQ,GAAG,CAAC,GAAG,EAAE,GAAG,OAAO;QAClE,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IACnD,IAAI,EAAE,IAAI,GAAG;QAAE,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAChE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,EACf,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,IAAI,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;IACpD,IAAI,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAC1B,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QACrC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAChD,MAAM,OAAO,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,EAChC,OAAO,GAAG,EAAE,GAAG,EAAE,CAAC;QACpB,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC,IAAI,KAAK,EAAE,CAAC;YACrD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;gBAAE,MAAM;YACpC,OAAO,MAAM,CAAC,MAAM,EAAE,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QACnF,CAAC;QACD,GAAG,GAAG,OAAO,CAAC;QACd,GAAG,GAAG,OAAO,CAAC;IAChB,CAAC;IACD,OAAO,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;AACnD,CAAC;AAED;;GAEG;AACH,SAAS,eAAe,CACtB,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAqB;IAErB,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,EACjB,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;IACzB,MAAM,UAAU,GAAG,CAAC,GAAG,CAAC,EACtB,QAAQ,GAAG,CAAC,GAAG,CAAC,EAChB,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC;IACrB,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,KAAK,GAAG,CAAC;QAAE,OAAO,KAAK,CAAC;IAC7C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,UAAU,GACd,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAC3B,CAAC,CAAC,CAAC,EAAE;YACD,MAAM,GAAG,EAAE;YACX,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,MAAM,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,YAAY,CAAC,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,EAAE,CAAC,CAAC;YACxC,EAAE;QACJ,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,EAAE,GAAG,CAAC,GAAG,CAAC,EACV,EAAE,GAAG,CAAC,GAAG,CAAC,EACV,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;IAChB,MAAM,WAAW,GACf,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC;QAC3B,CAAC,CAAC,CAAC,GAAG;YACF,EAAE;YACF,QAAQ,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC;YACrB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YACtB,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;YACzB,CAAC,CAAC,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,CAAC,CAAC;IACR,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,EACtB,QAAQ,GAAG,CAAC,GAAG,CAAC,EAChB,IAAI,GAAG,QAAQ,GAAG,KAAK,CAAC;IAC1B,MAAM,UAAU,GACd,WAAW;QACX,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC;QACzB,QAAQ,CAAC,QAAQ,EAAE,KAAK,EAAE,IAAI,CAAC;QAC/B,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,CAAC;QAC9B,CAAC,GAAG,UAAU,CAAC;IACjB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,wEAAwE;QACxE,2EAA2E;QAC3E,IAAI,IAAI,GAAG,UAAU,GAAG,CAAC;YAAE,OAAO,KAAK,CAAC;QACxC,OAAO,MAAM,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;IAChD,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,KAAK,CAAC;IAC5C,IAAI,KAAK,GAAG,CAAC,EACX,KAAK,GAAG,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,EAAE,QAAQ,GAAG,KAAK,CAAC,CAAC;IAC5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC;IAC1C,iFAAiF;IACjF,2EAA2E;IAC3E,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,UAAU,GAAG,KAAK,CAAC,CAAC;IAC/C,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC;QAAE,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;IAC7C,oEAAoE;IACpE,yEAAyE;IACzE,4EAA4E;IAC5E,MAAM,GAAG,GAAG,CAAC,CAAC,IAAI,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;IACrC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC;IAC/D,IAAI,QAAQ,GACV,GAAG,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK,IAAI,OAAO,GAAG,GAAG;QAC9D,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,KAAK,EAAE,UAAU,GAAG,KAAK,CAAC;QACnE,CAAC,CAAC,SAAS,CAAC;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC,EAAE,EAAE,CAAC;QACxD,MAAM,UAAU,GAAG,KAAK,GAAG,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,UAAU,GAAG,UAAU;QAC/B,gEAAgE;QAChE,oEAAoE;QACpE,UAAU,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EACtB,CAAC,GAAG,UAAU,GAAG,KAAK,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EACd,UAAU,GAAG,IAAI,GAAG,EAAE,CAAC;QACzB,MAAM,MAAM,GACV,CAAC,UAAU,GAAG,CAAC,QAAQ,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,UAAU,CAAC,GAAG,UAAU,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;YACtF,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9B,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YAAE,OAAO,KAAK,CAAC;QAC7C,IAAI,QAAQ,GAAG,CAAC;YAAE,KAAK,GAAG,CAAC,CAAC;;YACvB,KAAK,GAAG,CAAC,CAAC;QACf,MAAM,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7D,MAAM,IAAI,GAAG,CAAC,GAAG,QAAQ,GAAG,UAAU,CAAC;QACvC,MAAM,SAAS,GAAG,IAAI,GAAG,KAAK,IAAI,IAAI,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QACpF,IACE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC5D,KAAK,GAAG,KAAK,IAAI,CAAC,GAAG,OAAO,GAAG,CAAC,EAChC,CAAC;YACD,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;YACrC,MAAM;QACR,CAAC;QACD,CAAC,GAAG,SAAS,CAAC;IAChB,CAAC;IACD,OAAO,QAAQ,KAAK,SAAS,IAAI,YAAY,CAAC,MAAM,EAAE,SAAS,EAAE,QAAQ,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;AAC7F,CAAC;AAED,yFAAyF;AACzF,SAAS,QAAQ,CAAC,CAAS,EAAE,CAAS,EAAE,GAAW;IACjD,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;IACrB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;AACvC,CAAC;AACD,SAAS,YAAY,CAAC,CAAS,EAAE,CAAS,EAAE,OAAe;IACzD,MAAM,EAAE,GAAG,SAAS,GAAG,CAAC,EACtB,EAAE,GAAG,SAAS,GAAG,CAAC,CAAC;IACrB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QAAE,OAAO,CAAC,CAAC;IAC3D,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EACtB,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EAClB,EAAE,GAAG,CAAC,GAAG,EAAE,EACX,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,OAAO,EAAE,GAAG,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AACzD,CAAC;AAED,8DAA8D;AAC9D,SAAS,YAAY,CACnB,MAAqB,EACrB,SAAiB,EACjB,QAAgB,EAChB,CAAS,EACT,CAAS,EACT,CAAS,EACT,KAAa;IAEb,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAC5B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IAC3B,OAAO,MAAM,CAAC,MAAM,EAAE,SAAS,EAAE,QAAQ,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;AACnG,CAAC;AAED,0FAA0F;AAC1F,SAAS,MAAM,CAAC,MAAqB,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;IACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,OAAO,KAAK,CAAC;IACpF,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,OAAO,IAAI,CAAC;AACd,CAAC"}

package/dist/types/bigint.d.ts

@@ -14,7 +14,7 @@ export declare const BigInt64ArrayAvailable: boolean;
 export declare const BigUint64ArrayAvailable: boolean;
 declare const BigIntCtor: typeof BigIntUnavailable | BigIntConstructor;
 declare const BigInt64ArrayCtor: BigInt64ArrayConstructor | typeof BigInt64ArrayUnavailable;
-declare const BigUint64ArrayCtor: BigUint64ArrayConstructor[] | (typeof BigInt64ArrayUnavailable)[];
+declare const BigUint64ArrayCtor: BigUint64ArrayConstructor | typeof BigInt64ArrayUnavailable;
 export { BigIntCtor as BigInt };
 export { BigInt64ArrayCtor as BigInt64Array };
 export { BigUint64ArrayCtor as BigUint64Array };

package/dist/types/bigint.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bigint.d.ts","sourceRoot":"","sources":["../../src/types/bigint.ts"],"names":[],"mappings":"AAYA,iBAAS,iBAAiB,SAEzB;kBAFQ,iBAAiB;;;;AAU1B,cAAM,wBAAwB;IAC5B,MAAM,KAAK,iBAAiB,WAE3B;IACD,MAAM,CAAC,EAAE;IAGT,MAAM,CAAC,IAAI;;CAMZ;AAED,eAAO,MAAM,eAAe,SAAgC,CAAC;AAC7D,eAAO,MAAM,sBAAsB,SAAuC,CAAC;AAC3E,eAAO,MAAM,uBAAuB,SAAwC,CAAC;AAE7E,QAAA,MAAM,UAAU,8CAEZ,CAAC;AAEL,QAAA,MAAM,iBAAiB,4DAEnB,CAAC;AAEL,QAAA,MAAM,kBAAkB,mEAEpB,CAAC;AAEL,OAAO,EAAC,UAAU,IAAI,MAAM,EAAC,CAAC;AAC9B,OAAO,EAAC,iBAAiB,IAAI,aAAa,EAAC,CAAC;AAC5C,OAAO,EAAC,kBAAkB,IAAI,cAAc,EAAC,CAAC"}
+{"version":3,"file":"bigint.d.ts","sourceRoot":"","sources":["../../src/types/bigint.ts"],"names":[],"mappings":"AA8BA,iBAAS,iBAAiB,SAEzB;kBAFQ,iBAAiB;;;;AAU1B,cAAM,wBAAwB;IAC5B,MAAM,KAAK,iBAAiB,WAE3B;IACD,MAAM,CAAC,EAAE;IAGT,MAAM,CAAC,IAAI;;CAMZ;AAED,eAAO,MAAM,eAAe,SAAgC,CAAC;AAC7D,eAAO,MAAM,sBAAsB,SAAuC,CAAC;AAC3E,eAAO,MAAM,uBAAuB,SAAwC,CAAC;AAE7E,QAAA,MAAM,UAAU,8CAEZ,CAAC;AAEL,QAAA,MAAM,iBAAiB,4DAEnB,CAAC;AAEL,QAAA,MAAM,kBAAkB,6DAEpB,CAAC;AAEL,OAAO,EAAC,UAAU,IAAI,MAAM,EAAC,CAAC;AAC9B,OAAO,EAAC,iBAAiB,IAAI,aAAa,EAAC,CAAC;AAC5C,OAAO,EAAC,kBAAkB,IAAI,cAAc,EAAC,CAAC"}

package/dist/types/bigint.js

@@ -1,9 +1,26 @@
 // math.gl
-// SPDX-License-Identifier: MIT
+// SPDX-License-Identifier: MIT AND Apache-2.0
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileCopyrightText: Copyright 2016-2019 The Apache Software Foundation
+// SPDX-FileComment: BigInt availability/fallback structure adapted from Apache Arrow JS util/compat.ts; compared with apache-arrow-4.0.0. See ../../ARROW-NOTICE.txt and ../../APACHE-2.0-LICENSE.txt. Modified to separate availability flags, share errors and narrow exports.
+// Licensed to the Apache Software Foundation (ASF) under one
+// or more contributor license agreements.  See the NOTICE file
+// distributed with this work for additional information
+// regarding copyright ownership.  The ASF licenses this file
+// to you 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.
 // BigInt compatibility layer
-// Inspired by ArrowJS (under Apache2 license)
-// https://github.com/apache/arrow/blob/master/js/src/util/compat.ts
+// Source comparison: https://github.com/apache/arrow/blob/apache-arrow-4.0.0/js/src/util/compat.ts
 // Requires tsconfig.json: target=esnext or (lib: esnext.bigint)
 // Requires eslint: env: {es2020: true}
 const ERR_BIGINT_UNAVAILABLE = 'BigInt is not available in this environment';
@@ -40,7 +57,7 @@ const BigInt64ArrayCtor = (() => {
     return BigInt64ArrayAvailable ? BigInt64Array : BigInt64ArrayUnavailable;
 })();
 const BigUint64ArrayCtor = (() => {
-    return BigUint64ArrayAvailable ? [BigUint64Array] : [BigInt64ArrayUnavailable];
+    return BigUint64ArrayAvailable ? BigUint64Array : BigInt64ArrayUnavailable;
 })();
 export { BigIntCtor as BigInt };
 export { BigInt64ArrayCtor as BigInt64Array };

package/dist/types/bigint.js.map

@@ -1 +1 @@
-{"version":3,"file":"bigint.js","sourceRoot":"","sources":["../../src/types/bigint.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,6BAA6B;AAC7B,8CAA8C;AAC9C,oEAAoE;AACpE,gEAAgE;AAChE,uCAAuC;AAEvC,MAAM,sBAAsB,GAAG,6CAA6C,CAAC;AAE7E,SAAS,iBAAiB;IACxB,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC;AACD,iBAAiB,CAAC,MAAM,GAAG,GAAG,EAAE;IAC9B,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC,CAAC;AACF,iBAAiB,CAAC,OAAO,GAAG,GAAG,EAAE;IAC/B,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC,CAAC;AAEF,MAAM,wBAAwB;IAC5B,MAAM,KAAK,iBAAiB;QAC1B,OAAO,CAAC,CAAC;IACX,CAAC;IACD,MAAM,CAAC,EAAE;QACP,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;IACD,MAAM,CAAC,IAAI;QACT,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;IACD;QACE,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;CACF;AAED,MAAM,CAAC,MAAM,eAAe,GAAG,OAAO,MAAM,KAAK,WAAW,CAAC;AAC7D,MAAM,CAAC,MAAM,sBAAsB,GAAG,OAAO,aAAa,KAAK,WAAW,CAAC;AAC3E,MAAM,CAAC,MAAM,uBAAuB,GAAG,OAAO,cAAc,KAAK,WAAW,CAAC;AAE7E,MAAM,UAAU,GAAG,CAAC,GAAG,EAAE;IACvB,OAAO,eAAe,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,iBAAiB,CAAC;AACtD,CAAC,CAAC,EAAE,CAAC;AAEL,MAAM,iBAAiB,GAAG,CAAC,GAAG,EAAE;IAC9B,OAAO,sBAAsB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,wBAAwB,CAAC;AAC3E,CAAC,CAAC,EAAE,CAAC;AAEL,MAAM,kBAAkB,GAAG,CAAC,GAAG,EAAE;IAC/B,OAAO,uBAAuB,CAAC,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC,wBAAwB,CAAC,CAAC;AACjF,CAAC,CAAC,EAAE,CAAC;AAEL,OAAO,EAAC,UAAU,IAAI,MAAM,EAAC,CAAC;AAC9B,OAAO,EAAC,iBAAiB,IAAI,aAAa,EAAC,CAAC;AAC5C,OAAO,EAAC,kBAAkB,IAAI,cAAc,EAAC,CAAC"}
+{"version":3,"file":"bigint.js","sourceRoot":"","sources":["../../src/types/bigint.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,4DAA4D;AAC5D,6EAA6E;AAC7E,iRAAiR;AAEjR,6DAA6D;AAC7D,+DAA+D;AAC/D,wDAAwD;AACxD,6DAA6D;AAC7D,oDAAoD;AACpD,6DAA6D;AAC7D,6DAA6D;AAC7D,EAAE;AACF,+CAA+C;AAC/C,EAAE;AACF,6DAA6D;AAC7D,8DAA8D;AAC9D,yDAAyD;AACzD,4DAA4D;AAC5D,0DAA0D;AAC1D,qBAAqB;AAErB,6BAA6B;AAC7B,mGAAmG;AACnG,gEAAgE;AAChE,uCAAuC;AAEvC,MAAM,sBAAsB,GAAG,6CAA6C,CAAC;AAE7E,SAAS,iBAAiB;IACxB,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC;AACD,iBAAiB,CAAC,MAAM,GAAG,GAAG,EAAE;IAC9B,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC,CAAC;AACF,iBAAiB,CAAC,OAAO,GAAG,GAAG,EAAE;IAC/B,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;AAC1C,CAAC,CAAC;AAEF,MAAM,wBAAwB;IAC5B,MAAM,KAAK,iBAAiB;QAC1B,OAAO,CAAC,CAAC;IACX,CAAC;IACD,MAAM,CAAC,EAAE;QACP,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;IACD,MAAM,CAAC,IAAI;QACT,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;IACD;QACE,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAC1C,CAAC;CACF;AAED,MAAM,CAAC,MAAM,eAAe,GAAG,OAAO,MAAM,KAAK,WAAW,CAAC;AAC7D,MAAM,CAAC,MAAM,sBAAsB,GAAG,OAAO,aAAa,KAAK,WAAW,CAAC;AAC3E,MAAM,CAAC,MAAM,uBAAuB,GAAG,OAAO,cAAc,KAAK,WAAW,CAAC;AAE7E,MAAM,UAAU,GAAG,CAAC,GAAG,EAAE;IACvB,OAAO,eAAe,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,iBAAiB,CAAC;AACtD,CAAC,CAAC,EAAE,CAAC;AAEL,MAAM,iBAAiB,GAAG,CAAC,GAAG,EAAE;IAC9B,OAAO,sBAAsB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,wBAAwB,CAAC;AAC3E,CAAC,CAAC,EAAE,CAAC;AAEL,MAAM,kBAAkB,GAAG,CAAC,GAAG,EAAE;IAC/B,OAAO,uBAAuB,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,wBAAwB,CAAC;AAC7E,CAAC,CAAC,EAAE,CAAC;AAEL,OAAO,EAAC,UAAU,IAAI,MAAM,EAAC,CAAC;AAC9B,OAAO,EAAC,iBAAiB,IAAI,aAAa,EAAC,CAAC;AAC5C,OAAO,EAAC,kBAAkB,IAAI,cAAc,EAAC,CAAC"}

package/dist/types/bounds-types.js

@@ -3,7 +3,7 @@
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
 /** Checks if a `Bounds` is a `Bounds2D` */
 export function isBounds2D(bounds) {
-    return bounds.length === 2;
+    return bounds.length === 2 && bounds[0].length === 2 && bounds[1].length === 2;
 }
 /** Accepts 2D and 3D bounds and returns a 2D bound, truncating 3D dimension if necessary */
 export function getBounds2D(bounds) {

package/dist/types/bounds-types.js.map

@@ -1 +1 @@
-{"version":3,"file":"bounds-types.js","sourceRoot":"","sources":["../../src/types/bounds-types.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAa5D,2CAA2C;AAC3C,MAAM,UAAU,UAAU,CAAC,MAAc;IACvC,OAAO,MAAM,CAAC,MAAM,KAAK,CAAC,CAAC;AAC7B,CAAC;AAED,4FAA4F;AAC5F,MAAM,UAAU,WAAW,CAAC,MAAc;IACxC,OAAO,UAAU,CAAC,MAAM,CAAC;QACvB,CAAC,CAAC,MAAM;QACR,CAAC,CAAC;YACE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5B,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC7B,CAAC;AACR,CAAC"}
+{"version":3,"file":"bounds-types.js","sourceRoot":"","sources":["../../src/types/bounds-types.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAa5D,2CAA2C;AAC3C,MAAM,UAAU,UAAU,CAAC,MAAc;IACvC,OAAO,MAAM,CAAC,MAAM,KAAK,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC;AACjF,CAAC;AAED,4FAA4F;AAC5F,MAAM,UAAU,WAAW,CAAC,MAAc;IACxC,OAAO,UAAU,CAAC,MAAM,CAAC;QACvB,CAAC,CAAC,MAAM;QACR,CAAC,CAAC;YACE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5B,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC7B,CAAC;AACR,CAAC"}

package/package.json

@@ -6,7 +6,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "5.0.0-alpha.19",
+  "version": "5.0.0-alpha.20",
   "keywords": [
     "webgl",
     "javascript",
@@ -90,10 +90,14 @@
   "files": [
     "dist",
     "src",
-    "PROJ4-LICENSE.md"
+    "PROJ4-LICENSE.md",
+    "LICENSE",
+    "THIRD-PARTY-NOTICES.md",
+    "ARROW-NOTICE.txt",
+    "APACHE-2.0-LICENSE.txt"
   ],
   "sideEffects": [
     "./src/lib/common.js"
   ],
-  "gitHead": "f5700027934faf25c37230905c7ffa14d7f677a7"
+  "gitHead": "ad8c205a98e5e057d9055300b1a71f422b275c73"
 }

package/src/classes/euler.ts

@@ -2,6 +2,8 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
 // SPDX-FileCopyrightText: Copyright (c) 2017 Uber Technologies, Inc.
+// SPDX-FileCopyrightText: Copyright 2010-2017 three.js authors
+// SPDX-FileComment: Euler conversion methods adapted from three.js; see ../../LICENSE.
 
 import {MathArray} from './base/math-array';
 import {NumericArray} from '../types';

package/src/expressions/function-registry.ts

@@ -5,7 +5,15 @@
 import type {ExpressionFunction, ExpressionFunctionLibrary} from './function-libraries';
 
 const FUNCTION_NAME_PATTERN = /^[A-Za-z_$][A-Za-z0-9_$]*$/;
-const DISALLOWED_FUNCTION_NAMES = new Set(['__proto__', 'prototype', 'constructor']);
+const DISALLOWED_FUNCTION_NAMES = new Set([
+  '__proto__',
+  'prototype',
+  'constructor',
+  'true',
+  'false',
+  'null',
+  'this'
+]);
 
 /**
  * Options for registering functions with an {@link ExpressionFunctionRegistry}.

package/src/expressions/get.ts

@@ -24,7 +24,7 @@ function isObject(value: unknown): value is Record<string, unknown> {
   return value !== null && typeof value === 'object';
 }
 
-const keyMap: Record<string, string[]> = {};
+const keyMap: Record<string, string[]> = Object.create(null);
 
 /** Returns a cached list of property names for a dot-separated path. */
 function getKeys(compositeKey: string): string[] {