//testRcppMatrix.cpp #include RcppExport SEXP plusonediag_RcppMatrix(SEXP matrix) { Rcpp::NumericMatrix orig(matrix); Rcpp::NumericMatrix plus(matrix); int n = orig.rows(); int k = orig.cols(); for (int i=0; i