SparseMatrix kronecker产品崩溃

时间:2018-08-21 10:53:31

标签: eigen

我正在针对本征3.3.5编译以下内容:

#include <Eigen/Dense>
#include <Eigen/Sparse>
#include <unsupported/Eigen/KroneckerProduct>

template <typename eT, int SIDE>
inline
auto
xLeft(const Eigen::Matrix<eT, SIDE, SIDE> & M)
{
  Eigen::SparseMatrix<double> Id(SIDE, SIDE);
  Id.setIdentity();
  return kroneckerProduct(M.transpose(), Id);
}

void fail()
{
  Eigen::Matrix<double,2,2> M = Eigen::Matrix<double,2,2>::Random();
  std::cout << "die:" << xLeft(M) << "\n";
}

int main() { fail(); return 0; }

崩溃了,其他人都看到了,这是本征错误还是我做错了我看不见的东西?

还:

g++ (GCC) 7.3.1 20180712 (Red Hat 7.3.1-6)

0 个答案:

没有答案