7 #include <RcppArmadillo.h> 10 int main(
int argc,
char *argv[]) {
14 std::string cmd =
"set.seed(42); matrix(rnorm(9),3,3)";
16 arma::mat m = Rcpp::as<arma::mat>(R.
parseEval(cmd));
17 arma::mat n = m.t() * m;
18 double nacc = arma::accu(n);
19 double nrnk = arma::rank(n);
21 m.print(
"Initial Matrix m");
22 n.print(
"Product n = m' * m");
23 std::cout <<
"accu(n) " << nacc <<
" " 24 <<
"rank(n) " << nrnk << std::endl;
int parseEval(const std::string &line, SEXP &ans)
int main(int argc, char *argv[])