GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-09-17 20:03:17 +08:00
Code Issues Packages Projects Releases Wiki Activity

update this test to build against current eigen.

Browse Source 
remove the 'normal' path as it was not compiling anymore and I couldn't see the point of it (?)
This commit is contained in:
Benoit Jacob 2010-08-23 23:21:25 -04:00
parent 6261f4629f
commit 6924d4eec5
1 changed files with 3 additions and 27 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
bench/benchBlasGemm.cpp
View File

@ -1,5 +1,4 @@
// g++ -O3 -DNDEBUG -I.. -L /usr/lib64/atlas/ benchBlasGemm.cpp -o benchBlasGemm -lrt -lcblas
// g++-4.2 -O3 -DNDEBUG -I.. benchBlasGemm.cpp /usr/lib/libcblas.so.3 - o benchBlasGemm
// possible options: // possible options:
// -DEIGEN_DONT_VECTORIZE // -DEIGEN_DONT_VECTORIZE
// -msse2 // -msse2
@ -13,7 +12,9 @@
#include "BenchTimer.h" #include "BenchTimer.h"
// include the BLAS headers // include the BLAS headers
extern "C" {
#include <cblas.h> #include <cblas.h>
}
#include <string> #include <string>
#ifdef _FLOAT #ifdef _FLOAT
@ -27,7 +28,6 @@ typedef double Scalar;
typedef Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic> MyMatrix; typedef Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic> MyMatrix;
void bench_eigengemm(MyMatrix& mc, const MyMatrix& ma, const MyMatrix& mb, int nbloops); void bench_eigengemm(MyMatrix& mc, const MyMatrix& ma, const MyMatrix& mb, int nbloops);
void bench_eigengemm_normal(MyMatrix& mc, const MyMatrix& ma, const MyMatrix& mb, int nbloops);
void check_product(int M, int N, int K); void check_product(int M, int N, int K);
void check_product(void); void check_product(void);
@ -153,24 +153,6 @@ int main(int argc, char *argv[])
std::cout << M << " : " << timer.value() << " ; " << 1e-3*floor(1e-6*nbmad/timer.value()) << "\n"; std::cout << M << " : " << timer.value() << " ; " << 1e-3*floor(1e-6*nbmad/timer.value()) << "\n";
} }
// clear
ma = MyMatrix::Random(M,K);
mb = MyMatrix::Random(K,N);
mc = MyMatrix::Random(M,N);
// eigen normal
if (!(std::string(argv[1])=="auto"))
{
timer.reset();
for (uint k=0 ; k<nbtries ; ++k)
{
timer.start();
bench_eigengemm_normal(mc, ma, mb, nbloops);
timer.stop();
}
std::cout << "eigen : " << timer.value() << " (" << 1e-3*floor(1e-6*nbmad/timer.value()) << " GFlops/s)\n";
}
return 0; return 0;
} }
@ -182,12 +164,6 @@ void bench_eigengemm(MyMatrix& mc, const MyMatrix& ma, const MyMatrix& mb, int n
mc.noalias() += ma * mb; mc.noalias() += ma * mb;
} }
void bench_eigengemm_normal(MyMatrix& mc, const MyMatrix& ma, const MyMatrix& mb, int nbloops)
{
for (uint j=0 ; j<nbloops ; ++j)
mc.noalias() += GeneralProduct<MyMatrix,MyMatrix,UnrolledProduct>(ma,mb);
}
#define MYVERIFY(A,M) if (!(A)) { \ #define MYVERIFY(A,M) if (!(A)) { \
std::cout << "FAIL: " << M << "\n"; \ std::cout << "FAIL: " << M << "\n"; \
} }