test: add some basic code example for benchmark

Author: lpatiny

benchmark/denseMatrix.mjs

@@ -0,0 +1,19 @@
+import { Matrix } from 'ml-matrix';
+
+let size = 200;
+
+console.time('dense');
+const matrix = new Matrix(size, size).fill(1);
+const matrix2 = new Matrix(size, size).fill(1);
+
+const product = matrix.mmul(matrix2);
+// sum of all the elements
+let sum = 0;
+for (let i = 0; i < size; i++) {
+  for (let j = 0; j < size; j++) {
+    sum += product.get(i, j);
+  }
+}
+console.log(sum)
+console.timeEnd('dense');
+

benchmark/sparseMatrix.byhand.mjs

@@ -0,0 +1,41 @@
+
+
+import { SparseMatrix } from '../lib/index.js';
+
+let size = 200;
+
+
+
+console.time('dense');
+
+const matrix1 = new SparseMatrix(size, size);
+fillSparseMatrix(matrix1)
+
+const matrix2 = new SparseMatrix(size, size);
+fillSparseMatrix(matrix2)
+
+const product = new SparseMatrix(size, size);
+for (let i = 0; i < size; i++) {
+  for (let j = 0; j < size; j++) {
+    for (let k = 0; k < size; k++) {
+      product.set(i, j, product.get(i, j) + matrix1.get(i, k) * matrix2.get(k, j));
+    }
+  }
+}
+
+let sum = 0;
+for (let i = 0; i < size; i++) {
+  for (let j = 0; j < size; j++) {
+    sum += product.get(i, j);
+  }
+}
+console.timeEnd('dense');
+console.log(sum)
+
+function fillSparseMatrix(matrix) {
+  for (let row = 0; row < matrix.rows; row++) {
+    for (let column = 0; column < matrix.columns; column++) {
+      matrix.set(row, column, 1);
+    }
+  }
+}

benchmark/sparseMatrix.mjs

@@ -0,0 +1,30 @@
+import { SparseMatrix } from '../lib/index.js';
+
+let size = 200;
+
+console.time('dense');
+
+const matrix = new SparseMatrix(size, size);
+fillSparseMatrix(matrix)
+
+const matrix2 = new SparseMatrix(size, size);
+fillSparseMatrix(matrix2)
+
+const product = matrix.mmul(matrix2);
+// sum of all the elements
+let sum = 0;
+for (let i = 0; i < size; i++) {
+  for (let j = 0; j < size; j++) {
+    sum += product.get(i, j);
+  }
+}
+console.timeEnd('dense');
+console.log(sum)
+
+function fillSparseMatrix(matrix) {
+  for (let row = 0; row < matrix.rows; row++) {
+    for (let column = 0; column < matrix.columns; column++) {
+      matrix.set(row, column, 1);
+    }
+  }
+}

package.json

@@ -37,6 +37,7 @@
     "@vitest/coverage-v8": "^2.0.5",
     "eslint": "^8.10.0",
     "eslint-config-cheminfo": "^11.1.1",
+    "ml-matrix": "^6.11.1",
     "prettier": "^3.3.3",
     "rollup": "^4.21.0",
     "vitest": "^2.0.5"