From 0cfb8be29eb6229f1893a0737e76fc1711483d2c Mon Sep 17 00:00:00 2001
From: Olivier Parcollet <olivier.parcollet@cea.fr>
Date: Fri, 10 Jan 2014 15:29:18 +0100
Subject: [PATCH] add one speed test ...

machine dependant, to get an idea...
---
 test/speed/mat_vec_mat.cpp | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100644 test/speed/mat_vec_mat.cpp

diff --git a/test/speed/mat_vec_mat.cpp b/test/speed/mat_vec_mat.cpp
new file mode 100644
index 00000000..4a1d25ed
--- /dev/null
+++ b/test/speed/mat_vec_mat.cpp
@@ -0,0 +1,34 @@
+#include <triqs/arrays.hpp>
+
+using namespace triqs::arrays;
+using namespace triqs;
+const int N1= 25; //, N2 = 300;
+
+struct mat_vec {
+ void operator()() {
+  matrix<double> a (N1,N1), V(N1, N1), f(N1, N1);
+  for (int u =0; u<5000; ++u)
+   for (int j=0; j<N1; ++j)
+  triqs::arrays::blas::gemv(1.0, a, V(j,range()), 0.0,f(j,range()));
+ }
+};
+
+struct mat_mat {
+ void operator()() {
+  matrix<double> a (N1,N1), V(N1, N1), f(N1, N1);
+  for (int u =0; u<5000; ++u)
+    triqs::arrays::blas::gemm(1.0, a, V, 0.0,f);
+ }
+};
+
+
+#include "./speed_tester.hpp"
+int main() {
+ const int l = 200;
+
+ speed_tester<mat_mat> (l);
+ speed_tester<mat_vec> (l);
+
+ return 0;
+}
+