1 #include <iostream>
2 #include <fstream>
3 #include <vector>
4 #include <string>
5 #include "eigen_src/Eigen/Core"
6 #include "../BenchTimer.h"
7 using namespace Eigen;
8
9 #ifndef SCALAR
10 #error SCALAR must be defined
11 #endif
12
13 typedef SCALAR Scalar;
14
15 typedef Matrix<Scalar,Dynamic,Dynamic> Mat;
16
17 template<typename Func>
18 EIGEN_DONT_INLINE
bench(long m,long n,long k,const Func & f)19 double bench(long m, long n, long k, const Func& f)
20 {
21 Mat A(m,k);
22 Mat B(k,n);
23 Mat C(m,n);
24 A.setRandom();
25 B.setRandom();
26 C.setZero();
27
28 BenchTimer t;
29
30 double up = 1e8*4/sizeof(Scalar);
31 double tm0 = 4, tm1 = 10;
32 if(NumTraits<Scalar>::IsComplex)
33 {
34 up /= 4;
35 tm0 = 2;
36 tm1 = 4;
37 }
38
39 double flops = 2. * m * n * k;
40 long rep = std::max(1., std::min(100., up/flops) );
41 long tries = std::max(tm0, std::min(tm1, up/flops) );
42
43 BENCH(t, tries, rep, f(A,B,C));
44
45 return 1e-9 * rep * flops / t.best();
46 }
47
48 template<typename Func>
main_gemm(int argc,char ** argv,const Func & f)49 int main_gemm(int argc, char **argv, const Func& f)
50 {
51 std::vector<double> results;
52
53 std::string filename = std::string("gemm_settings.txt");
54 if(argc>1)
55 filename = std::string(argv[1]);
56 std::ifstream settings(filename);
57 long m, n, k;
58 while(settings >> m >> n >> k)
59 {
60 //std::cerr << " Testing " << m << " " << n << " " << k << std::endl;
61 results.push_back( bench(m, n, k, f) );
62 }
63
64 std::cout << RowVectorXd::Map(results.data(), results.size());
65
66 return 0;
67 }
68