如何构造一个从使用.Call（）调用的C函数返回的命名列表（SEXP）？

Question

我通过.Call("foo", <args>)调用C代码，其中foo调用其他C函数，计算结果并返回它。 结果是长度为3的列表，我想列出这个列表。 为此，foo做到了这一点：

/* Construct result list from variables containing the results */
SEXP res = PROTECT(allocVector(VECSXP, 3)); /* list of length 3 */
SET_VECTOR_ELT(res, 0, ScalarReal(a)); /* numeric(1) */ 
SET_VECTOR_ELT(res, 1, somenumericvector); /* numeric(<some length>) */
SET_VECTOR_ELT(res, 2, ScalarInteger(i)); /* integer(1) */

/* Name components and return */
SEXP nms = PROTECT(allocVector(STRSXP, 3)); /* names as SEXP */
char *nms_ = CHAR(STRING_ELT(nms, 0)); /* pointer to names */
char *names[3] = {"result_numeric", "result_numeric_vector", "result_integer"};
for(i = 0; i < 3; i++) nms_[i] = names[i]; 
setAttrib(res, R_NamesSymbol, nms);
UNPROTECT(1); 
return res;

这是构建长度为3的命名列表的正确方法吗？

C函数确实返回R但是一旦我将输出分配给R中的变量，我立即得到分段错误。 可能有什么问题？ 我可以把'debug statement''（简单的printf("...\\n")放在上面的'return res;'之前，它们被执行得很好。有没有方便的方法来调试从R调用的C代码？

Answer 1

到BrodieG的回答另一种方法是使用mkNamed从Rinlinedfuns.h（其中包含如何使用一个例子mkNamed ）。

/* Construct named result list from variables containing the results */
const char *names[] = {"result_numeric", "result_numeric_vector",
    "result_integer", ""};                   /* note the null string */
SEXP res = PROTECT(mkNamed(VECSXP, names));  /* list of length 3 */
SET_VECTOR_ELT(res, 0, ScalarReal(a));       /* numeric(1) */ 
SET_VECTOR_ELT(res, 1, somenumericvector);   /* numeric(<some length>) */
SET_VECTOR_ELT(res, 2, ScalarInteger(i));    /* integer(1) */
UNPROTECT(1);
return res;

Answer 2

由于你要求普通的vanilla方式，你需要使用mkChar和SET_STRING_ELT从C字符串创建一个R字符向量：

for(i = 0; i < 3; i++) SET_STRING_ELT(nms, i, mkChar(names[i]));  

现在，您正在尝试将原始C字符串用作R对象，这将无法正常工作。

Re：调试，可以在C代码中使用PrintValue打印出R对象。

所有这些都说，除非你有一个非常具体的理由想要普通香草，你应该考虑Rcpp 。

Answer 3

Per @ nrussell的好建议，单一陈述中的解决方案（为了易读性而分为四行）

R> cppFunction('List marius(double a, NumericVector b, int c) \
       { return List::create(Named("resnum")=a,\ 
                             Named("resvec")=b, \
                             Named("resint")=c); }')
R> marius(1.2, sqrt(3:5), 42L)
$resnum
[1] 1.2

$resvec
[1] 1.73205 2.00000 2.23607

$resint
[1] 42

R>