diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2010-07-15 20:45:45 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2010-07-15 20:45:45 +0200 |
| commit | 6a370f50c72a320a4a5a2740e881b61601d2a917 (patch) | |
| tree | 79499ae4a74752dc97e121e41d5997a34bc83aa9 /bench/bench_gemm.cpp | |
| parent | b08c26aefaa30d7e3a8d218909c93a36df25ce6e (diff) | |
MPRealSupport was missing
Diffstat (limited to 'bench/bench_gemm.cpp')
| -rw-r--r-- | bench/bench_gemm.cpp | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/bench/bench_gemm.cpp b/bench/bench_gemm.cpp index 06a124f8f..d04c24f44 100644 --- a/bench/bench_gemm.cpp +++ b/bench/bench_gemm.cpp @@ -10,8 +10,8 @@ using namespace std; using namespace Eigen; #ifndef SCALAR -#define SCALAR std::complex<double> -// #define SCALAR double +// #define SCALAR std::complex<double> +#define SCALAR float #endif typedef SCALAR Scalar; @@ -87,7 +87,9 @@ void blas_gemm(const MatrixXd& a, const MatrixXd& b, MatrixXd& c) template<typename M> EIGEN_DONT_INLINE void gemm(const M& a, const M& b, M& c) { - c.noalias() += a * b; + //c.noalias() += a * b; + for(int i=0;i<b.cols();++i) + c.col(i).noalias() += a * b.col(i); } int main(int argc, char ** argv) |