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)
diff --git a/unsupported/Eigen/CMakeLists.txt b/unsupported/Eigen/CMakeLists.txt
index d01f95d71..a858e5113 100644
--- a/unsupported/Eigen/CMakeLists.txt
+++ b/unsupported/Eigen/CMakeLists.txt
@@ -1,6 +1,6 @@
 set(Eigen_HEADERS AdolcForward BVH IterativeSolvers MatrixFunctions MoreVectorization AutoDiff AlignedVector3 Polynomials
                   CholmodSupport FFT NonLinearOptimization SparseExtra SuperLUSupport UmfPackSupport IterativeSolvers
-                  NumericalDiff Skyline TaucsSupport
+                  NumericalDiff Skyline TaucsSupport MPRealSupport
    )
 
 install(FILES