1*bf2c3715SXin Li
2*bf2c3715SXin Li #define NOGMM
3*bf2c3715SXin Li #define NOMTL
4*bf2c3715SXin Li #define EIGEN_GOOGLEHASH_SUPPORT 1
5*bf2c3715SXin Li
6*bf2c3715SXin Li #include <map>
7*bf2c3715SXin Li #include <ext/hash_map>
8*bf2c3715SXin Li #include <google/dense_hash_map>
9*bf2c3715SXin Li #include <google/sparse_hash_map>
10*bf2c3715SXin Li
11*bf2c3715SXin Li #ifndef SIZE
12*bf2c3715SXin Li #define SIZE 10000
13*bf2c3715SXin Li #endif
14*bf2c3715SXin Li
15*bf2c3715SXin Li #ifndef DENSITY
16*bf2c3715SXin Li #define DENSITY 0.01
17*bf2c3715SXin Li #endif
18*bf2c3715SXin Li
19*bf2c3715SXin Li #ifndef REPEAT
20*bf2c3715SXin Li #define REPEAT 1
21*bf2c3715SXin Li #endif
22*bf2c3715SXin Li
23*bf2c3715SXin Li #include "BenchSparseUtil.h"
24*bf2c3715SXin Li
25*bf2c3715SXin Li #ifndef MINDENSITY
26*bf2c3715SXin Li #define MINDENSITY 0.0004
27*bf2c3715SXin Li #endif
28*bf2c3715SXin Li
29*bf2c3715SXin Li #ifndef NBTRIES
30*bf2c3715SXin Li #define NBTRIES 10
31*bf2c3715SXin Li #endif
32*bf2c3715SXin Li
33*bf2c3715SXin Li #define BENCH(X) \
34*bf2c3715SXin Li timer.reset(); \
35*bf2c3715SXin Li for (int _j=0; _j<NBTRIES; ++_j) { \
36*bf2c3715SXin Li timer.start(); \
37*bf2c3715SXin Li for (int _k=0; _k<REPEAT; ++_k) { \
38*bf2c3715SXin Li X \
39*bf2c3715SXin Li } timer.stop(); }
40*bf2c3715SXin Li
41*bf2c3715SXin Li
42*bf2c3715SXin Li static double rtime;
43*bf2c3715SXin Li static double nentries;
44*bf2c3715SXin Li
45*bf2c3715SXin Li template<typename SetterType>
dostuff(const char * name,EigenSparseMatrix & sm1)46*bf2c3715SXin Li void dostuff(const char* name, EigenSparseMatrix& sm1)
47*bf2c3715SXin Li {
48*bf2c3715SXin Li int rows = sm1.rows();
49*bf2c3715SXin Li int cols = sm1.cols();
50*bf2c3715SXin Li sm1.setZero();
51*bf2c3715SXin Li BenchTimer t;
52*bf2c3715SXin Li SetterType* set1 = new SetterType(sm1);
53*bf2c3715SXin Li t.reset(); t.start();
54*bf2c3715SXin Li for (int k=0; k<nentries; ++k)
55*bf2c3715SXin Li (*set1)(internal::random<int>(0,rows-1),internal::random<int>(0,cols-1)) += 1;
56*bf2c3715SXin Li t.stop();
57*bf2c3715SXin Li std::cout << "std::map => \t" << t.value()-rtime
58*bf2c3715SXin Li << " nnz=" << set1->nonZeros() << std::flush;
59*bf2c3715SXin Li
60*bf2c3715SXin Li // getchar();
61*bf2c3715SXin Li
62*bf2c3715SXin Li t.reset(); t.start(); delete set1; t.stop();
63*bf2c3715SXin Li std::cout << " back: \t" << t.value() << "\n";
64*bf2c3715SXin Li }
65*bf2c3715SXin Li
main(int argc,char * argv[])66*bf2c3715SXin Li int main(int argc, char *argv[])
67*bf2c3715SXin Li {
68*bf2c3715SXin Li int rows = SIZE;
69*bf2c3715SXin Li int cols = SIZE;
70*bf2c3715SXin Li float density = DENSITY;
71*bf2c3715SXin Li
72*bf2c3715SXin Li EigenSparseMatrix sm1(rows,cols), sm2(rows,cols);
73*bf2c3715SXin Li
74*bf2c3715SXin Li
75*bf2c3715SXin Li nentries = rows*cols*density;
76*bf2c3715SXin Li std::cout << "n = " << nentries << "\n";
77*bf2c3715SXin Li int dummy;
78*bf2c3715SXin Li BenchTimer t;
79*bf2c3715SXin Li
80*bf2c3715SXin Li t.reset(); t.start();
81*bf2c3715SXin Li for (int k=0; k<nentries; ++k)
82*bf2c3715SXin Li dummy = internal::random<int>(0,rows-1) + internal::random<int>(0,cols-1);
83*bf2c3715SXin Li t.stop();
84*bf2c3715SXin Li rtime = t.value();
85*bf2c3715SXin Li std::cout << "rtime = " << rtime << " (" << dummy << ")\n\n";
86*bf2c3715SXin Li const int Bits = 6;
87*bf2c3715SXin Li for (;;)
88*bf2c3715SXin Li {
89*bf2c3715SXin Li dostuff<RandomSetter<EigenSparseMatrix,StdMapTraits,Bits> >("std::map ", sm1);
90*bf2c3715SXin Li dostuff<RandomSetter<EigenSparseMatrix,GnuHashMapTraits,Bits> >("gnu::hash_map", sm1);
91*bf2c3715SXin Li dostuff<RandomSetter<EigenSparseMatrix,GoogleDenseHashMapTraits,Bits> >("google::dense", sm1);
92*bf2c3715SXin Li dostuff<RandomSetter<EigenSparseMatrix,GoogleSparseHashMapTraits,Bits> >("google::sparse", sm1);
93*bf2c3715SXin Li
94*bf2c3715SXin Li // {
95*bf2c3715SXin Li // RandomSetter<EigenSparseMatrix,GnuHashMapTraits,Bits> set1(sm1);
96*bf2c3715SXin Li // t.reset(); t.start();
97*bf2c3715SXin Li // for (int k=0; k<n; ++k)
98*bf2c3715SXin Li // set1(internal::random<int>(0,rows-1),internal::random<int>(0,cols-1)) += 1;
99*bf2c3715SXin Li // t.stop();
100*bf2c3715SXin Li // std::cout << "gnu::hash_map => \t" << t.value()-rtime
101*bf2c3715SXin Li // << " nnz=" << set1.nonZeros() << "\n";getchar();
102*bf2c3715SXin Li // }
103*bf2c3715SXin Li // {
104*bf2c3715SXin Li // RandomSetter<EigenSparseMatrix,GoogleDenseHashMapTraits,Bits> set1(sm1);
105*bf2c3715SXin Li // t.reset(); t.start();
106*bf2c3715SXin Li // for (int k=0; k<n; ++k)
107*bf2c3715SXin Li // set1(internal::random<int>(0,rows-1),internal::random<int>(0,cols-1)) += 1;
108*bf2c3715SXin Li // t.stop();
109*bf2c3715SXin Li // std::cout << "google::dense => \t" << t.value()-rtime
110*bf2c3715SXin Li // << " nnz=" << set1.nonZeros() << "\n";getchar();
111*bf2c3715SXin Li // }
112*bf2c3715SXin Li // {
113*bf2c3715SXin Li // RandomSetter<EigenSparseMatrix,GoogleSparseHashMapTraits,Bits> set1(sm1);
114*bf2c3715SXin Li // t.reset(); t.start();
115*bf2c3715SXin Li // for (int k=0; k<n; ++k)
116*bf2c3715SXin Li // set1(internal::random<int>(0,rows-1),internal::random<int>(0,cols-1)) += 1;
117*bf2c3715SXin Li // t.stop();
118*bf2c3715SXin Li // std::cout << "google::sparse => \t" << t.value()-rtime
119*bf2c3715SXin Li // << " nnz=" << set1.nonZeros() << "\n";getchar();
120*bf2c3715SXin Li // }
121*bf2c3715SXin Li std::cout << "\n\n";
122*bf2c3715SXin Li }
123*bf2c3715SXin Li
124*bf2c3715SXin Li return 0;
125*bf2c3715SXin Li }
126*bf2c3715SXin Li
127