#include #include #include // Choose either one: //typedef Eigen::MatrixXd Mat; typedef Eigen::MatrixXf Mat; EIGEN_DONT_INLINE void trivial_product(Mat& res, const Mat& A, const Mat& B) { typedef Mat::Scalar Scalar; typedef Mat::Index Index; Index rows = A.rows(), cols = B.cols(), inner = A.cols(); res.resize(rows, cols); for(Index i=0; i