Félicitations !

Optimiser le code R avec Rcpp

Romain François

Consulting Datactive, ThinkR

evalCpp et cppFunction

Évaluer de simples instructions C++

evalCpp( "40+2" )

Créer une fonction C++ depuis la console R

cppFunction( "double add( double x, double y){
    return x + y ;
}")
add( 40, 2 )

Boucles for

for( init ; condition ; increment ){
    body
}

init : ce qui se passe au début
condition : faut-il poursuivre la boucle
increment : après chaque itération
body : ce que fait la boucle

Boucles for

for( int i=0; i<n; i++){
    // do something with i
}

Indexation de vecteurs

NumericVector x = ... ; 
int n = x.size() ;

// first value
x[0]

// second value
x[1]

// last value
x[n-1]

Fichiers C++ avec Rcpp

#include <Rcpp.h>

using namespace Rcpp ;

// [[Rcpp::export]]
double add( double x, double y){
    return x + y ;
}

// [[Rcpp::export]]
double twice( double x){
    return 2.0 * x;
}

Fonction Rcpp type

#include <Rcpp.h>
using namespace Rcpp ;
// [[Rcpp::export]]
double fun( NumericVector x ){

    // extract data from input and prepare outputs
    int n = x.size() ;
    double res = 0.0 ;

    // loop around input and/or output
    for(int i=0; i<n; i++){
        // do something with x[i]
    }

    // return output
    return res ;
}

Félicitations !

Optimiser le code R avec Rcpp