diff options
-rw-r--r-- | bench/bench_gemm.cpp | 8 | ||||
-rw-r--r-- | unsupported/Eigen/CMakeLists.txt | 2 |
2 files changed, 6 insertions, 4 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) 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 |