xref: /aosp_15_r20/external/eigen/bench/sparse_randomsetter.cpp (revision bf2c37156dfe67e5dfebd6d394bad8b2ab5804d4)
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