Run results for: Solving System of Linear Equations with Matrix

Run details:

User agent: Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/120.0.0.0 Safari/537.36

Browser: Chrome 120

Operating system: Linux

Device Platform: Desktop

Date tested: 2 years ago

Benchmark engine: Benchmark.js

jStat

1.4M/s

Numeric JS (LU Decomposed first)

879K/s

Sparse - math.js (LU Decomposed first)

623K/s

math.js (LU Decomposed first)

617K/s

Numeric JS (uses LUsolve under the hood)

523K/s

math.js

373K/s

Sparse - Numeric JS (LU Decomposed first)

180K/s

Sparse - jStat

85K/s

Sparse - Numeric JS (uses LUsolve under the hood)

62K/s

Sparse - math.js

10K/s

View exact numbers

Test name	Executions per second
✓ jStat	1,443,914 Ops/sec
Numeric JS (LU Decomposed first)	879,450 Ops/sec
Sparse - math.js (LU Decomposed first)	622,829 Ops/sec
math.js (LU Decomposed first)	616,685 Ops/sec
Numeric JS (uses LUsolve under the hood)	522,868 Ops/sec
math.js	373,476 Ops/sec
Sparse - Numeric JS (LU Decomposed first)	179,594 Ops/sec
Sparse - jStat	85,040 Ops/sec
Sparse - Numeric JS (uses LUsolve under the hood)	62,372 Ops/sec
Sparse - math.js	9,731 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 A = [[11, -1, 0], [-1, 12, -1], [0, -1, 11]],
    B = [[87], [12] , [0]];

var Numeric_LUa = numeric.LU(A),
    math_LUa = math.lup(A);

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(A);

Tests:

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