mirror of
https://github.com/josdejong/mathjs.git
synced 2026-01-25 15:07:57 +00:00
122 lines
5.6 KiB
JavaScript
122 lines
5.6 KiB
JavaScript
/**
|
|
* Benchmark
|
|
*
|
|
* Compare performance of basic matrix operations of a number of math libraries.
|
|
*
|
|
* These are some rough benchmarks to get an idea of the performance of math.js
|
|
* compared to other JavaScript libraries and to Octave (C++). They only give an
|
|
* _indication_ of the order of magnitude difference meant to see were math.js
|
|
* has room for improvements, it's not a fully fletched benchmark suite.
|
|
*/
|
|
|
|
import Benchmark from 'benchmark'
|
|
import det from 'ndarray-determinant'
|
|
import gemm from 'ndarray-gemm'
|
|
import ops from 'ndarray-ops'
|
|
import pack from 'ndarray-pack'
|
|
import numeric from 'numericjs'
|
|
import padRight from 'pad-right'
|
|
import sylvester from 'sylvester'
|
|
import zeros from 'zeros'
|
|
import { all, create } from '../../lib/esm/index.js'
|
|
|
|
const suite = new Benchmark.Suite()
|
|
const math = create(all)
|
|
|
|
function pad (text) {
|
|
return padRight(text, 40, ' ')
|
|
}
|
|
|
|
// fiedler matrix 25 x 25
|
|
const fiedler = [
|
|
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24],
|
|
[1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23],
|
|
[2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22],
|
|
[3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21],
|
|
[4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20],
|
|
[5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19],
|
|
[6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18],
|
|
[7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17],
|
|
[8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16],
|
|
[9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15],
|
|
[10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14],
|
|
[11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13],
|
|
[12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12],
|
|
[13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11],
|
|
[14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10],
|
|
[15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9],
|
|
[16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7, 8],
|
|
[17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6, 7],
|
|
[18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5, 6],
|
|
[19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4, 5],
|
|
[20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3, 4],
|
|
[21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2, 3],
|
|
[22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1, 2],
|
|
[23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, 1],
|
|
[24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0]
|
|
];
|
|
|
|
// mathjs
|
|
(function () {
|
|
const A = math.matrix(fiedler, 'dense', 'number')
|
|
|
|
suite.add(pad('matrix operations mathjs (number) A+A'), function () { return math.add(A, A) })
|
|
suite.add(pad('matrix operations mathjs (number) A*A'), function () { return math.multiply(A, A) })
|
|
suite.add(pad('matrix operations mathjs (number) A\''), function () { return math.transpose(A) })
|
|
suite.add(pad('matrix operations mathjs (number) det(A)'), function () { return math.det(A) })
|
|
})();
|
|
|
|
// mathjs
|
|
(function () {
|
|
const A = math.matrix(fiedler)
|
|
|
|
suite.add(pad('matrix operations mathjs (generic) A+A'), function () { return math.add(A, A) })
|
|
suite.add(pad('matrix operations mathjs (generic) A*A'), function () { return math.multiply(A, A) })
|
|
suite.add(pad('matrix operations mathjs (generic) A\''), function () { return math.transpose(A) })
|
|
suite.add(pad('matrix operations mathjs (generic) det(A)'), function () { return math.det(A) })
|
|
})();
|
|
|
|
// sylvester
|
|
(function () {
|
|
const A = sylvester.Matrix.create(fiedler)
|
|
|
|
suite.add(pad('matrix operations sylvester A+A'), function () { return A.add(A) })
|
|
suite.add(pad('matrix operations sylvester A*A'), function () { return A.multiply(A) })
|
|
suite.add(pad('matrix operations sylvester A\''), function () { return A.transpose() })
|
|
suite.add(pad('matrix operations sylvester det(A)'), function () { return A.det() })
|
|
})();
|
|
|
|
// numericjs
|
|
(function () {
|
|
const A = fiedler
|
|
|
|
suite.add(pad('matrix operations numericjs A+A'), function () { return numeric.add(A, A) })
|
|
suite.add(pad('matrix operations numericjs A*A'), function () { return numeric.dot(A, A) })
|
|
suite.add(pad('matrix operations numericjs A\''), function () { return numeric.transpose(A) })
|
|
suite.add(pad('matrix operations numericjs det(A)'), function () { return numeric.det(A) })
|
|
})();
|
|
|
|
// ndarray
|
|
(function () {
|
|
const A = pack(fiedler)
|
|
const B = zeros([25, 25])
|
|
|
|
suite.add(pad('matrix operations ndarray A+A'), function () { return ops.add(B, A, A) })
|
|
suite.add(pad('matrix operations ndarray A*A'), function () { return gemm(B, A, A) })
|
|
suite.add(pad('matrix operations ndarray A\''), function () { ops.assign(B, A); return B.transpose(1, 0) })
|
|
suite.add(pad('matrix operations ndarray det(A)'), function () { return det(A) })
|
|
})()
|
|
|
|
const durations = []
|
|
|
|
suite
|
|
.on('cycle', function (event) {
|
|
const benchmark = event.target
|
|
console.log(String(event.target))
|
|
durations.push(benchmark.name + ' avg duration per operation: ' + Math.round(benchmark.stats.mean * 1e6) + ' microseconds')
|
|
})
|
|
.run()
|
|
|
|
console.log()
|
|
console.log(durations.join('\n'))
|