MatrixXdSpecializations.hpp

Go to the documentation of this file.

#ifndef SRC_TYPES_MATRIXXDSPECIALIZATIONS_HPP_
#define SRC_TYPES_MATRIXXDSPECIALIZATIONS_HPP_

#include <types/Matrix.hpp>

#ifdef OpenBLAS_FOUND
#include <cblas.h>
#endif

namespace mic {
namespace types {

template<>
inline Eigen::MatrixXd mic::types::Matrix<double>::operator *(const Eigen::MatrixXd & mat_) {

#ifdef OpenBLAS_FOUND
    // Using openBLAS boost!
    //printf("Using specialized (double) BLAS boost!\n");
    // Get dimensions.
    size_t M = this->rows();
    size_t K = this->cols();
    size_t N = mat_.cols();
    // Create temporary matrix.
    Eigen::MatrixXd c(M,N);

    cblas_dgemm( CblasColMajor, CblasNoTrans, CblasNoTrans, M, N, K, 1.0,
            this->data(), M,
            mat_.data(), K, 0.0, c.data(), M );
    return c;
#else
    // Calling base EIGEN operator *
    //printf("Calling specialized (double) base EIGEN operator *\n");
    return Eigen::MatrixXd::operator*(mat_);
#endif
}


typedef Eigen::VectorXd VectorXd;


typedef mic::types::Matrix<double> MatrixXd;


typedef std::shared_ptr< mic::types::MatrixXd > MatrixXdPtr;


typedef std::pair< MatrixXdPtr, std::shared_ptr<char> > MatrixXdCharPair;

typedef std::pair< MatrixXdPtr, std::shared_ptr<unsigned int> > MatrixXdUintPair;


}//: namespace types
}//: namespace mic



#endif /* SRC_TYPES_MATRIXXDSPECIALIZATIONS_HPP_ */