[Rcpp-devel] Seamless Rcpp gives errors with RcppArmadillo

Wed, Jun 26, 2013 3:01 AM

Hello,

I don't have a copy of the book at hand. Below are some insights on what 
you are doing wrong.

Hope this helps.

Romain

Le 26/06/13 11:11, John Swan a ?crit :

This returns a colvec. See armadillo documentation:
http://arma.sourceforge.net/docs.html#sum

You need the "RcppArmadillo" plugin as you are using armadillo

You are calling rowSumsRcppArmadillo without argument, when it in fact 
takes an argument. This can not work.

If I try to execute this in R console with the source command, I get:

 > source("/Users/me/Desktop/R/somefile.R")

file29967eb06d23.cpp:20: error: ?arma? has not been declared
file29967eb06d23.cpp:20: error: expected constructor, destructor, or
type conversion before ?rowSumsRcppArmadillo?
file29967eb06d23.cpp: In function ?SEXPREC* file29967eb06d23()?:
file29967eb06d23.cpp:35: error: ?rowSumsRcppArmadillo? was not declared
in this scope
make: *** [file29967eb06d23.o] Error 1

ERROR(s) during compilation: source code errors or compiler
configuration errors!

Program source:
   1:
   2: // includes from the plugin
   3:
   4: #include <Rcpp.h>
   5:
   6:
   7: #ifndef BEGIN_RCPP
   8: #define BEGIN_RCPP
   9: #endif
  10:
  11: #ifndef END_RCPP
  12: #define END_RCPP
  13: #endif
  14:
  15: using namespace Rcpp;
  16:
  17:
  18: // user includes
  19:
  20: arma::rowvec rowSumsRcppArmadillo(NumericMatrix x){
  21:     arma::mat X = arma::mat(x.begin(), x.nrow(), x.ncol(), false);
  22:     return arma::sum(X, 1);
  23: }
  24:
  25: // declarations
  26: extern "C" {
  27: SEXP file29967eb06d23( ) ;
  28: }
  29:
  30: // definition
  31:
  32: SEXP file29967eb06d23(  ){
  33: BEGIN_RCPP
  34:
  35: return Rcpp::wrap(rowSumsRcppArmadillo());
  36:
  37: END_RCPP
  38: }
  39:
  40:
Error in compileCode(f, code, language = language, verbose = verbose) :
   Compilation ERROR, function(s)/method(s) not created!
file29967eb06d23.cpp:20: error: ?arma? has not been declared
file29967eb06d23.cpp:20: error: expected constructor, destructor, or
type conversion before ?rowSumsRcppArmadillo?
file29967eb06d23.cpp: In function ?SEXPREC* file29967eb06d23()?:
file29967eb06d23.cpp:35: error: ?rowSumsRcppArmadillo? was not declared
in this scope
make: *** [file29967eb06d23.o] Error 1
In addition: Warning message:
running command '/Library/Frameworks/R.framework/Resources/bin/R CMD
SHLIB file29967eb06d23.cpp 2> file29967eb06d23.cpp.err.txt' had status 1

I added these:
20: #include <RcppArmadillo.h>
  21: // [[Rcpp::depends(RcppArmadillo)]]
  22:
  23: // [[Rcpp::export]]
  24:
  25: arma::rowvec rowSumsRcppArmadillo(NumericMatrix x){
  26:     arma::mat X = arma::mat(x.begin(), x.nrow(), x.ncol(), false);
  27:     return arma::sum(X, 1);
  28: }

This works for me, put this code in a cpp file:

#include <RcppArmadillo.h>
using namespace Rcpp ;
// [[Rcpp::depends(RcppArmadillo)]]

// [[Rcpp::export]]
arma::colvec rowSumsRcppArmadillo(NumericMatrix x){
     arma::mat X = arma::mat(x.begin(), x.nrow(), x.ncol(), false);
     return arma::sum(X, 1);
}


And then sourceCpp() it.

 > sourceCpp( "arma.cpp" )

 > x <- matrix(rep(1, 20), nrow = 4)

 > rowSumsRcppArmadillo(x)
      [,1]
[1,]    5
[2,]    5
[3,]    5
[4,]    5

Do you have RcppArmadillo installed ?

sourceCpp is for sourcing c++ files, not R files.

Romain Francois
Professional R Enthusiast
+33(0) 6 28 91 30 30

R Graph Gallery: http://gallery.r-enthusiasts.com

blog:            http://blog.r-enthusiasts.com
|- http://bit.ly/13SrjxO : highlight 0.4.2
`- http://bit.ly/10X94UM : Mobile version of the graph gallery

[Rcpp-devel] Seamless Rcpp gives errors with RcppArmadillo

Thread (3 messages)