From 6a370f50c72a320a4a5a2740e881b61601d2a917 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 15 Jul 2010 20:45:45 +0200
Subject: MPRealSupport was missing

---
 bench/bench_gemm.cpp | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

(limited to 'bench/bench_gemm.cpp')

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)
-- 
cgit v1.2.3