내가 ++ C이 함수를 작성한

Question

는 R 함수 (C++)에 배열 송신 :

extern "C" 
{ 

    void add(int* first, int* second, int *n , int* sum) 
    { 
     for (int i = 0; i< *n; i++) 
     { 
      sum[i] = first[i] + second[i]; 
     } 


    } 
} 

이 드라이버 :

add <- function(twoColumn) 
{ 
    if(!is.data.frame(twoColumn))stop("twoColumn should be data.frame") 

    first <- twoColumn[1] 
    second <- twoColumn[2] 
    n <- length(first) 
    .C("add",first = as.integer(unlist(first)),second = as.integer(unlist(second)), n = as.integer(n),sum = as.integer(rep(0,n)))$sum 


} 

그러나 R내는 것은 단지 수이고 데이터 프레임의 첫 번째 행의 합계.

Answer 1

문제점은 무엇인지 모르겠지만 (다소 야심적인) 코드의 작동 버전을 제공 할 수 있습니다. 이 기능을 패키지로 당신을 필요로하며, 다음과 같이 작동

src/add.cpp :

SEXP add(SEXP Rx, SEXP Ry){ 
    SEXP xy_sum; 
    int i; 
    PROTECT(xy_sum = allocVector(REALSXP, 1)); 
    for(i = 0; i < length(Rx); i++){ 
     REAL(xy_sum)[i] = REAL(Rx)[i] + REAL(Ry)[i]; 
    } 
    UNPROTECT(1); 
    return xy_sum; 
} 

scr/add.h :

#ifndef _add_ADD_H 
#define _add_ADD_H 

#include <R.h> 
#include <Rdefines.h> 

extern "C" SEXP add(SEXP Rx, SEXP Ry); 

#endif 

R/add.R :

add <- function(x, y){ 
    if(length(x) != length(y)) stop("Vectors must be of the same length.") 
    # coerce into numeric in case the user supplied something silly 
    .Call("add", PACKAGE="add", as.numeric(x), as.numeric(y)) 
} 

R/zzz.R :

.onLoad <- function(lib, pkg){ 
    library.dynam("add", pkg, lib) 
}