10 int main(
int argc,
char *argv[]) {
14 std::string cmd =
"set.seed(42); matrix(rnorm(9),3,3)";
16 const Eigen::Map<Eigen::MatrixXd> m =
17 Rcpp::as<Eigen::Map<Eigen::MatrixXd> >(R.
parseEval(cmd));
18 Eigen::MatrixXd n = m.transpose() * m;
19 Eigen::ColPivHouseholderQR<Eigen::MatrixXd> nqr(n);
21 std::cout <<
"Initial Matrix m\n" << m << std::endl;
22 std::cout <<
"Product n = m' * m\n" << n << std::endl;
23 std::cout <<
"n.sum() " << n.sum() << std::endl;
24 std::cout <<
"nrq.rank() " << nqr.rank() << std::endl;
int main(int argc, char *argv[])
int parseEval(const std::string &line, SEXP &ans)