d3js-plugins-rails 0.0.2 → 0.0.3

data/lib/d3js-plugins-rails/version.rb

@@ -1,7 +1,8 @@
 module D3js
   module Plugins
     module Rails
-      VERSION = "0.0.2"
+      VERSION = "0.0.3"
+      D3_PLUGINS_VERSION = "c690548cb4173f5edf312133f364e949482c3821"
     end
   end
 end

data/vendor/assets/javascripts/d3/plugins/Makefile

@@ -12,6 +12,7 @@ d3.geo.projection.js: \
 	geo/projection/parallel1.js \
 	geo/projection/parallel2.js \
 	geo/projection/interrupt.js \
+	geo/projection/elliptic.js \
 	geo/projection/aitoff.js \
 	geo/projection/guyou.js \
 	geo/projection/mollweide.js \
@@ -37,6 +38,7 @@ d3.geo.projection.js: \
 	geo/projection/eckert5.js \
 	geo/projection/eckert6.js \
 	geo/projection/eisenlohr.js \
+	geo/projection/fahey.js \
 	geo/projection/gringorten.js \
 	geo/projection/hammer-retroazimuthal.js \
 	geo/projection/hammer.js \
@@ -51,6 +53,7 @@ d3.geo.projection.js: \
 	geo/projection/littrow.js \
 	geo/projection/loximuthal.js \
 	geo/projection/miller.js \
+	geo/projection/modified-stereographic.js \
 	geo/projection/mt-flat-polar-parabolic.js \
 	geo/projection/mt-flat-polar-quartic.js \
 	geo/projection/mt-flat-polar-sinusoidal.js \
@@ -58,10 +61,15 @@ d3.geo.projection.js: \
 	geo/projection/nell-hammer.js \
 	geo/projection/peirce-quincuncial.js \
 	geo/projection/polyconic.js \
+	geo/projection/rectangular-polyconic.js \
 	geo/projection/robinson.js \
 	geo/projection/satellite.js \
+	geo/projection/times.js \
+	geo/projection/two-point-azimuthal.js \
 	geo/projection/two-point-equidistant.js \
 	geo/projection/van-der-grinten.js \
+	geo/projection/van-der-grinten2.js \
+	geo/projection/van-der-grinten3.js \
 	geo/projection/van-der-grinten4.js \
 	geo/projection/wagner4.js \
 	geo/projection/wagner6.js \

data/vendor/assets/javascripts/d3/plugins/geo/projection/README.md

@@ -27,6 +27,7 @@ Projections:
 * d3.geo.eckert6 - [Eckert VI](http://bl.ocks.org/3734329)
 * d3.geo.eisenlohr - [Eisenlohr conformal](http://bl.ocks.org/3797585)
 * d3.geo.equirectangular - [Equirectangular (Plate Carrée)](http://bl.ocks.org/3757119)
+* d3.geo.fahey - [Fahey](http://bl.ocks.org/4731228)
 * d3.geo.gnomonic - [gnomonic](http://bl.ocks.org/3757349)
 * d3.geo.gringorten - [Gringorten](http://bl.ocks.org/3796882)
 * d3.geo.guyou - [Guyou hemisphere-in-a-square](http://bl.ocks.org/3763867)
@@ -43,6 +44,7 @@ Projections:
 * d3.geo.loximuthal - [loximuthal](http://bl.ocks.org/3867220)
 * d3.geo.mercator - [Mercator](http://bl.ocks.org/3757132)
 * d3.geo.miller - [Miller](http://bl.ocks.org/3734333)
+* d3.geo.modifiedStereographic - Modified Stereographic
 * d3.geo.mollweide - [Mollweide](http://bl.ocks.org/3734336), [Atlantis](http://bl.ocks.org/4519975)
 * d3.geo.mtFlatPolarParabolic - [McBryde–Thomas flat-polar parabolic](http://bl.ocks.org/4465130)
 * d3.geo.mtFlatPolarQuartic - [McBryde–Thomas flat-polar quartic](http://bl.ocks.org/4465137)
@@ -52,12 +54,18 @@ Projections:
 * d3.geo.orthographic - [orthographic](http://bl.ocks.org/3757125)
 * d3.geo.peirceQuincuncial - [Pierce quincuncial](http://bl.ocks.org/4310087)
 * d3.geo.polyconic - [Polyconic](http://bl.ocks.org/3734343)
+* d3.geo.rectangularPolyconic - Rectangular Polyconic
 * d3.geo.robinson - [Robinson](http://bl.ocks.org/3710566)
 * d3.geo.satellite - [satellite (tilted perpsective)](http://bl.ocks.org/3790444)
 * d3.geo.sinusoidal - [sinusoidal](http://bl.ocks.org/3712399)
 * d3.geo.sinuMollweide - [Sinu-Mollweide](http://bl.ocks.org/4319903)
 * d3.geo.stereographic - [stereographic](http://bl.ocks.org/3757137)
+* d3.geo.times - Times
+* d3.geo.twoPointAzimuthal - Two-Point Azimuthal
+* d3.geo.twoPointEquidistant - Two-Point Equidistant
 * d3.geo.vanDerGrinten - [Van der Grinten](http://bl.ocks.org/3796831)
+* d3.geo.vanDerGrinten2 - Van der Grinten II
+* d3.geo.vanDerGrinten3 - Van der Grinten III
 * d3.geo.vanDerGrinten4 - [Van der Grinten IV](http://bl.ocks.org/4489365)
 * d3.geo.wagner4 - [Wagner IV, Putniṇš P2´](http://bl.ocks.org/4487674)
 * d3.geo.wagner6 - [Wagner VI](http://bl.ocks.org/3710148)

data/vendor/assets/javascripts/d3/plugins/geo/projection/august.js

@@ -1,6 +1,6 @@
 function august(λ, φ) {
   var tanφ = Math.tan(φ / 2),
-      k = 1 - tanφ * tanφ,
+      k = asqrt(1 - tanφ * tanφ),
       c = 1 + k * Math.cos(λ /= 2),
       x = Math.sin(λ) * k / c,
       y = tanφ / c,
@@ -12,4 +12,22 @@ function august(λ, φ) {
   ];
 }
 
+august.invert = function(x, y) {
+  x *= 3 / 8, y *= 3 / 8;
+  if (!x && Math.abs(y) > 1) return null;
+  var x2 = x * x,
+      y2 = y * y,
+      s = 1 + x2 + y2,
+      sin3η = Math.sqrt(.5 * (s - Math.sqrt(s * s - 4 * y * y))),
+      η = asin(sin3η) / 3,
+      ξ = sin3η ? arcosh(Math.abs(y / sin3η)) / 3 : arsinh(Math.abs(x)) / 3,
+      cosη = Math.cos(η),
+      coshξ = cosh(ξ),
+      d = coshξ * coshξ - cosη * cosη;
+  return [
+    sgn(x) * 2 * Math.atan2(sinh(ξ) * cosη, .25 - d),
+    sgn(y) * 2 * Math.atan2(coshξ * Math.sin(η), .25 + d)
+  ];
+};
+
 (d3.geo.august = function() { return projection(august); }).raw = august;

data/vendor/assets/javascripts/d3/plugins/geo/projection/berghaus.js

@@ -17,6 +17,20 @@ function berghaus(n) {
     return p;
   }
 
+  forward.invert = function(x, y) {
+    var r = Math.sqrt(x * x + y * y);
+    if (r > π / 2) {
+      var θ = Math.atan2(y, x),
+          θ0 = k * Math.round((θ - π / 2) / k) + π / 2,
+          s = θ > θ0 ? -1 : 1,
+          A = r * Math.cos(θ0 - θ)
+          cotα = 1 / Math.tan(s * Math.acos((A - π) / Math.sqrt(π * (π - 2 * A) + r * r)));
+      θ = θ0 + 2 * Math.atan((cotα + s * Math.sqrt(cotα * cotα - 3)) / 3);
+      x = r * Math.cos(θ), y = r * Math.sin(θ);
+    }
+    return berghausAzimuthalEquidistant.invert(x, y);
+  };
+
   return forward;
 }
 

data/vendor/assets/javascripts/d3/plugins/geo/projection/bonne.js

@@ -7,7 +7,7 @@ function bonne(φ0) {
 
   function forward(λ, φ) {
     var ρ = cotφ0 + φ0 - φ,
-        E = λ * Math.cos(φ) / ρ;
+        E = ρ ? λ * Math.cos(φ) / ρ : ρ;
     return [
       ρ * Math.sin(E),
       cotφ0 - ρ * Math.cos(E)

data/vendor/assets/javascripts/d3/plugins/geo/projection/conic-conformal.js

@@ -3,9 +3,11 @@
 function conicConformal(φ0, φ1) {
   var cosφ0 = Math.cos(φ0),
       t = function(φ) { return Math.tan(π / 4 + φ / 2); },
-      n = Math.log(cosφ0 / Math.cos(φ1)) / Math.log(t(φ1) / t(φ0)),
+      n = φ0 === φ1 ? Math.sin(φ0) : Math.log(cosφ0 / Math.cos(φ1)) / Math.log(t(φ1) / t(φ0)),
       F = cosφ0 * Math.pow(t(φ0), n) / n;
 
+  if (!n) return conicConformalMercator;
+
   function forward(λ, φ) {
     var ρ = Math.abs(Math.abs(φ) - π / 2) < ε ? 0 : F / Math.pow(t(φ), n);
     return [
@@ -26,4 +28,12 @@ function conicConformal(φ0, φ1) {
   return forward;
 }
 
+function conicConformalMercator(λ, φ) {
+  return [λ, Math.log(Math.tan(π / 4 + φ / 2))];
+}
+
+conicConformalMercator.invert = function(x, y) {
+  return [x, 2 * Math.atan(Math.exp(y)) - π / 2];
+};
+
 (d3.geo.conicConformal = function() { return parallel2Projection(conicConformal); }).raw = conicConformal;

data/vendor/assets/javascripts/d3/plugins/geo/projection/conic-equidistant.js

@@ -2,9 +2,11 @@
 
 function conicEquidistant(φ0, φ1) {
   var cosφ0 = Math.cos(φ0),
-      n = (cosφ0 - Math.cos(φ1)) / (φ1 - φ0),
+      n = φ0 === φ1 ? Math.sin(φ0) : (cosφ0 - Math.cos(φ1)) / (φ1 - φ0),
       G = cosφ0 / n + φ0;
 
+  if (Math.abs(n) < ε) return d3.geo.equirectangular.raw;
+
   function forward(λ, φ) {
     var ρ = G - φ;
     return [

data/vendor/assets/javascripts/d3/plugins/geo/projection/craig.js

@@ -10,13 +10,19 @@ function craig(φ0) {
     ];
   }
 
-  // TODO
-  // forward.invert = function(x, y) {
-  //   return [
-  //     x,
-  //     asin(y * (x ? Math.tan(x) / x : 1))
-  //   ];
-  // };
+  forward.invert = tanφ0 ? function(x, y) {
+    if (x) y *= Math.sin(x) / x;
+    var cosλ = Math.cos(x);
+    return [
+      x,
+      2 * Math.atan2(Math.sqrt(cosλ * cosλ + tanφ0 * tanφ0 - y * y) - cosλ, tanφ0 - y)
+    ];
+  } : function(x, y) {
+    return [
+      x,
+      asin(x ? y * Math.tan(x) / x : y)
+    ];
+  };
 
   return forward;
 }

data/vendor/assets/javascripts/d3/plugins/geo/projection/eisenlohr.js

@@ -1,16 +1,77 @@
+// @import august
+
 function eisenlohr(λ, φ) {
-  var f = 3 + Math.sqrt(8),
-      s1 = Math.sin(λ /= 2),
-      c1 = Math.cos(λ),
-      k = Math.sqrt(Math.cos(φ) / 2),
-      cosφ2 = Math.cos(φ /= 2),
-      t = Math.sin(φ) / (cosφ2 + 2 * c1 * k),
+  var s0 = Math.sin(λ /= 2),
+      c0 = Math.cos(λ),
+      k = Math.sqrt(Math.cos(φ)),
+      c1 = Math.cos(φ /= 2),
+      t = Math.sin(φ) / (c1 + Math.SQRT2 * c0 * k),
       c = Math.sqrt(2 / (1 + t * t)),
-      v = Math.sqrt((cosφ2 + (c1 + s1) * k) / (cosφ2 + (c1 - s1) * k));
+      v = Math.sqrt((Math.SQRT2 * c1 + (c0 + s0) * k) / (Math.SQRT2 * c1 + (c0 - s0) * k));
   return [
-    f * (c * (v - 1 / v) - 2 * Math.log(v)),
-    f * (c * t * (v + 1 / v) - 2 * Math.atan(t))
+    eisenlohrK * (c * (v - 1 / v) - 2 * Math.log(v)),
+    eisenlohrK * (c * t * (v + 1 / v) - 2 * Math.atan(t))
   ];
 }
 
+eisenlohr.invert = function(x, y) {
+  var p = d3.geo.august.raw.invert(x / 1.2, y * 1.065);
+  if (!p) return null;
+  var λ = p[0],
+      φ = p[1],
+      i = 20;
+  x /= eisenlohrK, y /= eisenlohrK;
+  do {
+    var _0 = λ / 2,
+        _1 = φ / 2,
+        s0 = Math.sin(_0),
+        c0 = Math.cos(_0),
+        s1 = Math.sin(_1),
+        c1 = Math.cos(_1),
+        cos1 = Math.cos(φ),
+
+        k = Math.sqrt(cos1),
+        t = s1 / (c1 + Math.SQRT2 * c0 * k),
+        t2 = t * t,
+        c = Math.sqrt(2 / (1 + t2)),
+        v0 = (Math.SQRT2 * c1 + (c0 + s0) * k),
+        v1 = (Math.SQRT2 * c1 + (c0 - s0) * k),
+        v2 = v0 / v1,
+        v = Math.sqrt(v2),
+        vm1v = v - 1 / v,
+        vp1v = v + 1 / v,
+
+        fx = c * vm1v - 2 * Math.log(v) - x,
+        fy = c * t * vp1v - 2 * Math.atan(t) - y,
+
+        δtδλ = s1 && Math.SQRT1_2 * k * s0 * t2 / s1,
+        δtδφ = (Math.SQRT2 * c0 * c1 + k) / (2 * (c1 + Math.SQRT2 * c0 * k) * (c1 + Math.SQRT2 * c0 * k) * k),
+
+        δcδt = -.5 * t * c * c * c,
+        δcδλ = δcδt * δtδλ,
+        δcδφ = δcδt * δtδφ,
+
+        A = (A = 2 * c1 + Math.SQRT2 * k * (c0 - s0)) * A * v,
+        δvδλ = (Math.SQRT2 * c0 * c1 * k + cos1) / A,
+        δvδφ = -(Math.SQRT2 * s0 * s1) / (k * A),
+
+        δxδλ = vm1v * δcδλ - 2 * δvδλ / v + c * (δvδλ + δvδλ / v2),
+        δxδφ = vm1v * δcδφ - 2 * δvδφ / v + c * (δvδφ + δvδφ / v2),
+        δyδλ = t * vp1v * δcδλ - 2 * δtδλ / (1 + t2) + c * vp1v * δtδλ + c * t * (δvδλ - δvδλ / v2),
+        δyδφ = t * vp1v * δcδφ - 2 * δtδφ / (1 + t2) + c * vp1v * δtδφ + c * t * (δvδφ - δvδφ / v2),
+
+        denominator = δxδφ * δyδλ - δyδφ * δxδλ;
+    if (!denominator) break;
+    var δλ = (fy * δxδφ - fx * δyδφ) / denominator,
+        δφ = (fx * δyδλ - fy * δxδλ) / denominator;
+    λ -= δλ;
+    φ = Math.max(-π / 2, Math.min(π / 2, φ - δφ));
+  } while ((Math.abs(δλ) > ε || Math.abs(δφ) > ε) && --i > 0);
+  return Math.abs(Math.abs(φ) - π / 2) < ε
+      ? [0, φ]
+      : i && [λ, φ];
+};
+
+var eisenlohrK = 3 + 2 * Math.SQRT2;
+
 (d3.geo.eisenlohr = function() { return projection(eisenlohr); }).raw = eisenlohr;

data/vendor/assets/javascripts/d3/plugins/geo/projection/elliptic.js

@@ -0,0 +1,115 @@
+// Returns [sn, cn, dn](u + iv|m).
+function ellipticJi(u, v, m) {
+  if (!u) {
+    var b = ellipticJ(v, 1 - m);
+    return [
+      [0, b[0] / b[1]],
+      [1 / b[1], 0],
+      [b[2] / b[1], 0]
+    ];
+  }
+  var a = ellipticJ(u, m);
+  if (!v) return [[a[0], 0], [a[1], 0], [a[2], 0]];
+  var b = ellipticJ(v, 1 - m),
+      denominator = b[1] * b[1] + m * a[0] * a[0] * b[0] * b[0];
+  return [
+    [a[0] * b[2] / denominator, a[1] * a[2] * b[0] * b[1] / denominator],
+    [a[1] * b[1] / denominator, -a[0] * a[2] * b[0] * b[2] / denominator],
+    [a[2] * b[1] * b[2] / denominator, -m * a[0] * a[1] * b[0] / denominator]
+  ];
+}
+
+// Returns [sn, cn, dn, ph](u|m).
+function ellipticJ(u, m) {
+  var ai, b, φ, t, twon;
+  if (m < ε) {
+    t = Math.sin(u);
+    b = Math.cos(u);
+    ai = .25 * m * (u - t * b);
+    return [
+      t - ai * b,
+      b + ai * t,
+      1 - .5 * m * t * t,
+      u - ai
+    ];
+  }
+  if (m >= 1 - ε) {
+    ai = .25 * (1 - m);
+    b = cosh(u);
+    t = tanh(u);
+    φ = 1 / b;
+    twon = b * sinh(u);
+    return [
+      t + ai * (twon - u) / (b * b),
+      φ - ai * t * φ * (twon - u),
+      φ + ai * t * φ * (twon + u),
+      2 * Math.atan(Math.exp(u)) - π / 2 + ai * (twon - u) / b
+    ];
+  }
+
+  var a = [1, 0, 0, 0, 0, 0, 0, 0, 0],
+      c = [Math.sqrt(m), 0, 0, 0, 0, 0, 0, 0, 0],
+      i = 0;
+  b = Math.sqrt(1 - m);
+  twon = 1;
+
+  while (Math.abs(c[i] / a[i]) > ε && i < 8) {
+    ai = a[i++];
+    c[i] = .5 * (ai - b);
+    a[i] = .5 * (ai + b);
+    b = asqrt(ai * b);
+    twon *= 2;
+  }
+
+  φ = twon * a[i] * u;
+  do {
+    t = c[i] * Math.sin(b = φ) / a[i];
+    φ = .5 * (asin(t) + φ);
+  } while (--i);
+
+  return [Math.sin(φ), t = Math.cos(φ), t / Math.cos(φ - b), φ];
+}
+
+// Calculate F(φ+iψ|m).
+// See Abramowitz and Stegun, 17.4.11.
+function ellipticFi(φ, ψ, m) {
+  var r = Math.abs(φ),
+      i = Math.abs(ψ),
+      sinhψ = sinh(i);
+  if (r) {
+    var cscφ = 1 / Math.sin(r),
+        cotφ2 = 1 / (Math.tan(r) * Math.tan(r)),
+        b = -(cotφ2 + m * (sinhψ * sinhψ * cscφ * cscφ) - 1 + m),
+        c = (m - 1) * cotφ2,
+        cotλ2 = .5 * (-b + Math.sqrt(b * b - 4 * c));
+    return [
+      ellipticF(Math.atan(1 / Math.sqrt(cotλ2)), m) * sgn(φ),
+      ellipticF(Math.atan(asqrt((cotλ2 / cotφ2 - 1) / m)), 1 - m) * sgn(ψ)
+    ];
+  }
+  return [
+    0,
+    ellipticF(Math.atan(sinhψ), 1 - m) * sgn(ψ)
+  ];
+}
+
+// Calculate F(φ|m) where m = k² = sin²α.
+// See Abramowitz and Stegun, 17.6.7.
+function ellipticF(φ, m) {
+  if (!m) return φ;
+  if (m === 1) return Math.log(Math.tan(φ / 2 + π / 4));
+  var a = 1,
+      b = Math.sqrt(1 - m),
+      c = Math.sqrt(m);
+  for (var i = 0; Math.abs(c) > ε; i++) {
+    if (φ % π) {
+      var dφ = Math.atan(b * Math.tan(φ) / a);
+      if (dφ < 0) dφ += π;
+      φ += dφ + ~~(φ / π) * π;
+    } else φ += φ;
+    c = (a + b) / 2;
+    b = Math.sqrt(a * b);
+    c = ((a = c) - b) / 2;
+  }
+  return φ / (Math.pow(2, i) * a);
+}

data/vendor/assets/javascripts/d3/plugins/geo/projection/fahey.js

@@ -0,0 +1,19 @@
+function fahey(λ, φ) {
+  var t = Math.tan(φ / 2);
+  return [
+    λ * faheyK * asqrt(1 - t * t),
+    (1 + faheyK) * t
+  ];
+}
+
+fahey.invert = function(x, y) {
+  var t = y / (1 + faheyK);
+  return [
+    x ? x / (faheyK * asqrt(1 - t * t)) : 0,
+    2 * Math.atan(t)
+  ];
+};
+
+var faheyK = Math.cos(35 * radians);
+
+(d3.geo.fahey = function() { return projection(fahey); }).raw = fahey;

data/vendor/assets/javascripts/d3/plugins/geo/projection/guyou.js

@@ -1,46 +1,60 @@
+// @import elliptic
+
+// √k' tn(½K - w) = exp(-ζ).
 function guyou(λ, φ) {
-  return guyouEllipticFi(λ, sgn(φ) * Math.log(Math.tan(.5 * (Math.abs(φ) + π / 2))), .5);
+
+  var k_ = (Math.SQRT2 - 1) / (Math.SQRT2 + 1),
+      k = Math.sqrt(1 - k_ * k_),
+      K = ellipticF(π / 2, k * k),
+      f = -1;
+
+  var s = λ > 0 ? -1 : 1;
+  λ += s * π / 2;
+
+  var ψ = Math.log(Math.tan(π / 4 + Math.abs(φ) / 2)),
+      r = Math.exp(f * ψ) / Math.sqrt(k_),
+      at = guyouComplexAtan(r * Math.cos(f * λ), r * Math.sin(f * λ)),
+      t = ellipticFi(at[0], at[1], k * k);
+
+  return [-s * .5 * K - t[1], sgn(φ) * (.5 * K - t[0])];
 }
 
-// Calculate F(φ+iψ|m).
-// See Abramowitz and Stegun, 17.4.11.
-function guyouEllipticFi(φ, ψ, m) {
-  var r = Math.abs(φ),
-      i = Math.abs(ψ),
-      sinhψ = .5 * ((sinhψ = Math.exp(i)) - 1 / sinhψ);
-  if (r) {
-    var cscφ = 1 / Math.sin(r),
-        cotφ2 = (cotφ2 = Math.cos(r) * cscφ) * cotφ2,
-        b = -(cotφ2 + m * (sinhψ * sinhψ * cscφ * cscφ + 1) - 1),
-        cotλ2 = .5 * (-b + Math.sqrt(b * b - 4 * (m - 1) * cotφ2));
-    return [
-      guyouEllipticF(Math.atan(1 / Math.sqrt(cotλ2)), m) * sgn(φ),
-      guyouEllipticF(Math.atan(Math.sqrt(cotλ2 / cotφ2 - 1) / m), 1 - m) * sgn(ψ)
-    ];
-  }
+function guyouComplexAtan(x, y) {
+  var x2 = x * x,
+      y_1 = y + 1,
+      t = 1 - x2 - y * y;
   return [
-    0,
-    guyouEllipticF(Math.atan(sinhψ), 1 - m) * sgn(ψ)
+    sgn(x) * π / 4 - .5 * Math.atan2(t, 2 * x),
+    -.25 * Math.log(t * t + 4 * x2) + .5 * Math.log(y_1 * y_1 + x2)
   ];
 }
 
-// Calculate F(φ|m) where m = k² = sin²α.
-// See Abramowitz and Stegun, 17.6.7.
-function guyouEllipticF(φ, m) {
-  var a = 1,
-      b = Math.sqrt(1 - m),
-      c = Math.sqrt(m);
-  for (var i = 0; Math.abs(c) > ε; i++) {
-    if (φ % π) {
-      var dφ = Math.atan(b * Math.tan(φ) / a);
-      if (dφ < 0) dφ += π;
-      φ += dφ + ~~(φ / π) * π;
-    } else φ += φ;
-    c = (a + b) / 2;
-    b = Math.sqrt(a * b);
-    c = ((a = c) - b) / 2;
-  }
-  return φ / (Math.pow(2, i) * a);
+function guyouComplexDivide(a, b) {
+  var denominator = b[0] * b[0] + b[1] * b[1];
+  return [
+    (a[0] * b[0] + a[1] * b[1]) / denominator,
+    (a[1] * b[0] - a[0] * b[1]) / denominator
+  ];
 }
 
+guyou.invert = function(x, y) {
+  var k_ = (Math.SQRT2 - 1) / (Math.SQRT2 + 1),
+      k = Math.sqrt(1 - k_ * k_),
+      K = ellipticF(π / 2, k * k),
+      f = -1;
+
+  var s = x > 0 ? -1 : 1;
+
+  var j = ellipticJi(.5 * K - y, -s * .5 * K - x, k * k),
+      tn = guyouComplexDivide(j[0], j[1]),
+      λ = Math.atan2(tn[1], tn[0]) / f - s * π / 2;
+  if (λ < -π) λ += 2 * π;
+  else if (λ > π) λ -= 2 * π;
+
+  return [
+    λ,
+    2 * Math.atan(Math.exp(.5 / f * Math.log(k_ * tn[0] * tn[0] + k_ * tn[1] * tn[1]))) - π / 2
+  ];
+};
+
 (d3.geo.guyou = function() { return projection(guyou); }).raw = guyou;