特征-返回类型。CWISE产品?

我在RcppEigen中写了一个加权协方差的函数。在其中一个步骤中,我想取矩阵X的第I列和第j列,并计算cwiseProduct,它应该返回某种向量。cwiseProduct的输出将进入一个中间变量,该变量可以多次重用。从文件来看 cwiseProduct 返回a CwiseBinaryOp 它本身有两种类型。我的cwiseProduct在两个列向量上运行,因此我认为正确的返回类型应该是 Eigen::CwiseBinaryOp<Eigen::ColXpr, Eigen::ColXpr> 命名空间Eigen中没有名为ColXpr的成员

#include <RcppEigen.h>
// [[Rcpp::depends(RcppEigen)]]

Rcpp::List Crossprod_sparse(Eigen::MappedSparseMatrix<double> X, Eigen::Map<Eigen::MatrixXd> W) {
  int K = W.cols();
  int p = X.cols();

  Rcpp::List crossprods(W.cols());

  for (int i = 0; i < p; i++) {
    for (int j = i; j < p; j++) {
      Eigen::CwiseBinaryOp<Eigen::ColXpr, Eigen::ColXpr> prod = X.col(i).cwiseProduct(X.col(j));
      for (int k = 0; k < K; k++) {
        //double out = prod.dot(W.col(k));
      }
    }
  }
  return crossprods;
}

我也试着把钱存到SparseVector里

Eigen::SparseVector<double> prod = X.col(i).cwiseProduct(X.col(j));

X.col(i).cwiseProduct(X.col(j));