Add bench of LU

Author: termoshtt

ndarray-linalg/Cargo.toml

@@ -58,5 +58,9 @@ harness = false
 name = "eigh"
 harness = false
 
+[[bench]]
+name = "solve"
+harness = false
+
 [package.metadata.docs.rs]
 rustdoc-args = ["--html-in-header", "katex-header.html"]

ndarray-linalg/benches/eigh.rs

@@ -5,13 +5,13 @@ use ndarray_linalg::*;
 fn eigh_small(c: &mut Criterion) {
     let mut group = c.benchmark_group("eigh");
     for &n in &[4, 8, 16, 32, 64, 128] {
-        group.bench_with_input(BenchmarkId::new("eigh", n), &n, |b, n| {
+        group.bench_with_input(BenchmarkId::new("C", n), &n, |b, n| {
             let a: Array2<f64> = random((*n, *n));
             b.iter(|| {
                 let (_e, _vecs) = a.eigh(UPLO::Upper).unwrap();
             })
         });
-        group.bench_with_input(BenchmarkId::new("eigh_t", n), &n, |b, n| {
+        group.bench_with_input(BenchmarkId::new("F", n), &n, |b, n| {
             let a: Array2<f64> = random((*n, *n).f());
             b.iter(|| {
                 let (_e, _vecs) = a.eigh(UPLO::Upper).unwrap();

ndarray-linalg/benches/solve.rs

@@ -0,0 +1,24 @@
+use criterion::*;
+use ndarray::*;
+use ndarray_linalg::*;
+
+fn factorize_small(c: &mut Criterion) {
+    let mut group = c.benchmark_group("factorize");
+    for &n in &[4, 8, 16, 32, 64, 128] {
+        group.bench_with_input(BenchmarkId::new("C", n), &n, |b, n| {
+            let a: Array2<f64> = random((*n, *n));
+            b.iter(|| {
+                let _lu = a.factorize().unwrap();
+            })
+        });
+        group.bench_with_input(BenchmarkId::new("F", n), &n, |b, n| {
+            let a: Array2<f64> = random((*n, *n).f());
+            b.iter(|| {
+                let _lu = a.factorize().unwrap();
+            })
+        });
+    }
+}
+
+criterion_group!(factorize, factorize_small);
+criterion_main!(factorize);