MatrixXfSpecializations.hpp

Go to the documentation of this file.

#ifndef SRC_TYPES_MATRIXXFSPECIALIZATIONS_HPP_
#define SRC_TYPES_MATRIXXFSPECIALIZATIONS_HPP_

#include <types/Matrix.hpp>

//#include <types/Pair.hpp>


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

namespace mic {
namespace types {

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

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

    cblas_sgemm( 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 (float) base EIGEN operator *\n");
    return Eigen::MatrixXf::operator*(mat_);
#endif
}


typedef Eigen::VectorXf VectorXf;


typedef std::shared_ptr< Eigen::VectorXf >VectorXfPtr;


typedef mic::types::Matrix<float> MatrixXf;


typedef std::shared_ptr< mic::types::MatrixXf > MatrixXfPtr;


typedef std::pair< MatrixXfPtr, std::shared_ptr<char> > MatrixXfCharPair;

typedef std::pair< MatrixXfPtr, std::shared_ptr<unsigned int> > MatrixXfUintPair;

typedef std::pair< MatrixXfPtr, MatrixXfPtr > MatrixXfMatrixXfPair;


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



#endif /* SRC_TYPES_MATRIXXFSPECIALIZATIONS_HPP_ */