Add benchmarks for Cholesky::solve

AtheMathmo · Feb 17, 2017 · 25046d6 · 25046d6
1 parent 69479fd
commit 25046d6
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 0 deletions.
diff --git a/benches/lib.rs b/benches/lib.rs
@@ -1,5 +1,6 @@
 #![feature(test)]
 
+#[macro_use]
 extern crate rulinalg;
 extern crate num as libnum;
 extern crate test;

diff --git a/benches/linalg/cholesky.rs b/benches/linalg/cholesky.rs
@@ -47,3 +47,23 @@ fn cholesky_500x500(b: &mut Bencher) {
         x.cholesky().expect("Matrix is invertible")
     })
 }
+
+#[bench]
+fn cholesky_solve_1000x1000(b: &mut Bencher) {
+    let n = 1000;
+    let x = Matrix::identity(n);
+    let cholesky = Cholesky::decompose(x).unwrap();
+    b.iter(|| {
+        cholesky.solve(vector![0.0; n])
+    });
+}
+
+#[bench]
+fn cholesky_solve_100x100(b: &mut Bencher) {
+    let n = 100;
+    let x = Matrix::identity(n);
+    let cholesky = Cholesky::decompose(x).unwrap();
+    b.iter(|| {
+        cholesky.solve(vector![0.0; n])
+    });
+}