austince · austince · Dec 14, 2019 · Dec 14, 2019
diff --git a/README.md b/README.md
@@ -1,8 +1,6 @@
 # maptastic
 
-=========
-
-Javascript/CSS projection mapping utility.  Put your internets on things!
+Javascript/CSS projection mapping utility with zero dependencies.  Put your internets on things!
 
 ![maptastic animation](https://glowbox.github.io/maptasticjs/images/maptastic.gif "Maptastic JS")
 
@@ -31,8 +29,6 @@ When you include `maptastic.js` in your page, a new class `maptastic.Maptastic`
 
         <div id="so-simple">This is pretty simple.</div>
 
-        <!-- The only dependency is on `numeric`, which can be loaded from cdnjs -->
-        <script src="https://cdnjs.cloudflare.com/ajax/libs/numeric/1.2.6/numeric.js"></script>
         <script src="./dist/maptastic.js"></script>
         <script>
             new maptastic.Maptastic("so-simple");

diff --git a/example/index.html b/example/index.html
@@ -34,7 +34,6 @@
 	</pre>
 	<iframe src="//player.vimeo.com/video/108531008" width="500" height="281" frameborder="0" webkitallowfullscreen mozallowfullscreen allowfullscreen id="video-test"></iframe>
 
-	<script src="../node_modules/numeric/numeric-1.2.6.min.js"></script>
 	<script src="../dist/maptastic.js"></script>
 	<script>
 		new maptastic.Maptastic("video-test", "controls");

diff --git a/example/simple.html b/example/simple.html
@@ -12,7 +12,6 @@
 
 		<div id="so-simple">This is pretty simple.</div>
 
-		<script src="../node_modules/numeric/numeric-1.2.6.min.js"></script>
 		<script src="../dist/maptastic.js"></script>
 		<script>
 			new maptastic.Maptastic("so-simple");

diff --git a/example/three.html b/example/three.html
@@ -8,7 +8,6 @@
 	</head>
 	<body>
 
-		<script src="../node_modules/numeric/numeric-1.2.6.min.js"></script>
 		<script src="../dist/maptastic.js"></script>
 		<script src="https://cdnjs.cloudflare.com/ajax/libs/three.js/r69/three.js"></script>
 		<script>

diff --git a/package-lock.json b/package-lock.json
diff --git a/package.json b/package.json
@@ -25,9 +25,7 @@
     "rollup": "^1.27.9",
     "rollup-plugin-babel": "^4.3.3",
     "rollup-plugin-commonjs": "^10.1.0",
-    "rollup-plugin-uglify": "^6.0.3"
+    "rollup-plugin-uglify": "^6.0.4"
   },
-  "dependencies": {
-    "numeric": "^1.2.6"
-  }
+  "dependencies": {}
 }
diff --git a/src/maptastic.js b/src/maptastic.js
@@ -1,4 +1,4 @@
-import { solve } from 'numeric';
+import solve from './solve';
 
 /**
  * @typedef {Object} Config

diff --git a/src/solve.js b/src/solve.js
@@ -0,0 +1,164 @@
+/**
+ * @fileoverview A snippet from the {@link https://github.com/sloisel/numeric | numeric} project, as it is unmaintained and doesn't export correctly.
+ * @see {@link https://github.com/sloisel/numeric/blob/656fa1254be540f428710738ca9c1539625777f1/src/numeric.js#L3001-L3089 | Source}
+ */
+
+/**
+ * @see {@link https://github.com/sloisel/numeric/blob/master/src/numeric.js#L315-L329}
+ */
+function _foreach2(x, s, k, f) {
+  if (k === s.length - 1) {
+    return f(x);
+  }
+  var i, n = s[k], ret = Array(n);
+  for (i = n - 1; i >= 0; i--) {
+    ret[i] = _foreach2(x[i], s, k + 1, f);
+  }
+  return ret;
+}
+
+
+function _dim(x) {
+  var ret = [];
+  while (typeof x === "object") {
+    ret.push(x.length);
+    x = x[0];
+  }
+  return ret;
+}
+
+/**
+ * @see {@link https://github.com/sloisel/numeric/blob/656fa1254be540f428710738ca9c1539625777f1/src/numeric.js#L309-L329}
+ * @param {*} x 
+ */
+function dim(x) {
+  var y, z;
+  if (typeof x === "object") {
+    y = x[0];
+    if (typeof y === "object") {
+      z = y[0];
+      if (typeof z === "object") {
+        return _dim(x);
+      }
+      return [x.length, y.length];
+    }
+    return [x.length];
+  }
+  return [];
+}
+
+
+function cloneV(x) {
+  var _n = x.length;
+  var i, ret = Array(_n);
+
+  for (i = _n - 1; i !== -1; --i) {
+    ret[i] = (x[i]);
+  }
+  return ret;
+}
+
+function clone(x) {
+  if (typeof x !== "object") {
+    return (x);
+  }
+  var s = dim(x);
+  return _foreach2(x, s, 0, cloneV);
+}
+
+// 11. Ax = b
+/**
+ * @param {*} A 
+ * @param {boolean} [fast=false] 
+ * @returns {{ LU: *, P: * }}
+ */
+function LU(A, fast = false) {
+  var abs = Math.abs;
+  var i, j, k, absAjk, Akk, Ak, Pk, Ai;
+  var max;
+  var n = A.length, n1 = n - 1;
+  var P = new Array(n);
+  if (!fast) {
+    A = clone(A);
+  }
+
+  for (k = 0; k < n; ++k) {
+    Pk = k;
+    Ak = A[k];
+    max = abs(Ak[k]);
+    for (j = k + 1; j < n; ++j) {
+      absAjk = abs(A[j][k]);
+      if (max < absAjk) {
+        max = absAjk;
+        Pk = j;
+      }
+    }
+    P[k] = Pk;
+
+    if (Pk != k) {
+      A[k] = A[Pk];
+      A[Pk] = Ak;
+      Ak = A[k];
+    }
+
+    Akk = Ak[k];
+
+    for (i = k + 1; i < n; ++i) {
+      A[i][k] /= Akk;
+    }
+
+    for (i = k + 1; i < n; ++i) {
+      Ai = A[i];
+      for (j = k + 1; j < n1; ++j) {
+        Ai[j] -= Ai[k] * Ak[j];
+        ++j;
+        Ai[j] -= Ai[k] * Ak[j];
+      }
+      if (j === n1) Ai[j] -= Ai[k] * Ak[j];
+    }
+  }
+
+  return {
+    LU: A,
+    P: P
+  };
+}
+
+function LUsolve(LUP, b) {
+  var i, j;
+  var LU = LUP.LU;
+  var n = LU.length;
+  var x = clone(b);
+  var P = LUP.P;
+  var Pi, LUi, tmp;
+
+  for (i = n - 1; i !== -1; --i) x[i] = b[i];
+  for (i = 0; i < n; ++i) {
+    Pi = P[i];
+    if (P[i] !== i) {
+      tmp = x[i];
+      x[i] = x[Pi];
+      x[Pi] = tmp;
+    }
+
+    LUi = LU[i];
+    for (j = 0; j < i; ++j) {
+      x[i] -= x[j] * LUi[j];
+    }
+  }
+
+  for (i = n - 1; i >= 0; --i) {
+    LUi = LU[i];
+    for (j = i + 1; j < n; ++j) {
+      x[i] -= x[j] * LUi[j];
+    }
+
+    x[i] /= LUi[i];
+  }
+
+  return x;
+}
+
+export default function solve(A, b, fast) {
+  return LUsolve(LU(A, fast), b);
+}