Run results for: Solving Sparse System of Linear Equations with Matrix

Run details:

User agent: Mozilla/5.0 (iPad; CPU OS 13_4 like Mac OS X) AppleWebKit/605.1.15 (KHTML, like Gecko) CriOS/93.0.4577.78 Mobile/15E148 Safari/604.1

Browser: Chrome Mobile iOS 93

Operating system: iOS 13.4

Device Platform: Tablet

Date tested: 2 years ago

Test name	Executions per second
jStat	49118.0 Ops/sec
Numeric JS (uses LUsolve under the hood)	31665.6 Ops/sec
Numeric JS (LU Decomposed first)	42579.6 Ops/sec
math.js	10265.2 Ops/sec
math.js (LU Decomposed first)	32204.5 Ops/sec

HTML Preparation code:

<!-- JStat -->
<script src="https://cdn.jsdelivr.net/jstat/latest/jstat.min.js"></script>
<!-- NumericJS -->
<script src="https://cdn.rawgit.com/sloisel/numeric/master/src/numeric.js"></script>
<!-- Math.js -->
<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjs/3.9.0/math.min.js"></script>

Script Preparation code:

var sparseA = [
  [1,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
  [0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
  [0,0,4,0,0,0,0,6,0,0,0,0,6,0,0,0,0,0,0,0],
  [0,1,0,2,0,0,0,0,0,0,0,0,0,6,0,6,0,0,0,0],
  [0,0,4,0,6,0,0,0,0,0,4,0,0,0,0,0,0,5,0,0],
  [0,0,0,2,0,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0],
  [0,0,0,0,6,0,23,4,0,0,0,0,0,0,6,0,0,0,0,0],
  [0,3,0,0,0,5,0,6,5,6,0,0,0,0,0,0,0,0,0,0],
  [0,0,0,0,0,2,6,45,66,0,0,0,0,0,0,0,6,0,0,0],
  [0,0,0,0,0,0,0,9,0,12,6,0,0,0,0,0,0,0,6,0],
  [0,0,3,0,0,0,0,0,23,0,45,0,0,5,0,0,0,0,0,0],
  [0,0,0,0,0,0,0,0,0,13,0,6,0,0,0,0,0,0,0,0],
  [3,0,0,0,3,0,56,0,0,0,5,0,5,0,0,0,0,6,0,0],
  [0,0,0,0,0,0,0,0,0,0,6,6,0,23,0,0,0,0,0,0],
  [0,3,0,0,0,0,0,0,0,0,0,0,8,6,67,0,0,0,0,0],
  [0,0,3,0,0,0,0,0,0,0,0,0,0,1,0,2,0,0,0,0],
  [0,0,0,0,0,0,0,5,0,0,0,4,0,0,5,0,43,6,0,0],
  [0,3,0,0,0,0,0,0,0,0,0,6,0,0,0,45,6,6,6,0],
  [0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,4,0,7,0],
  [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
],
sparseB = [[0],[0],[4],[0],[10],[0],[5],[0],[0],[0],[-5],[6],[0],[-6],[0],[0],[0],[10],[20],[4]];

var Numeric_sparse_LUa = numeric.LU(sparseA),
    math_sparse_LUa = math.lup(sparseA);

Tests:

jStat
jStat.gauss_elimination(sparseA, sparseB);
Numeric JS (uses LUsolve under the hood)
numeric.solve(sparseA, sparseB);
Numeric JS (LU Decomposed first)
numeric.LUsolve(Numeric_sparse_LUa, sparseB);
math.js
math.lusolve(sparseA, sparseB);
math.js (LU Decomposed first)
math.lusolve(math_sparse_LUa, sparseB);

Run results for: Solving Sparse System of Linear Equations with Matrix

jStat

Numeric JS (uses LUsolve under the hood)

Numeric JS (LU Decomposed first)

math.js

math.js (LU Decomposed first)