extern "C" { SEXP myfun(SEXP S,SEXP A,SEXP B,SEXP C) { //will call this function from R. SEXP rate,dir,List,List_names; //declare variables PROTECT( rate = allocMatrix(REALSXP,10,2) ); //allocate 10x2 matrix of double type? PROTECT( dir = allocVector(INTSXP,10) ); //allocated vector(10) of int type? double* p_rate = REAL(rate); //why do I need pointers? int* p_dir = INTEGER(dir);
//这里我称之为C函数,它计算向量< vector< double> > someVal和vector< int> someVal2.
现在我想将这些值传递给rate和dir.
for(int i =0; i < 10; i++){ rate[i][0] = someVal1[i][0]; rate[i][1] = someVal1[i][1]; dir[i] = someVal2[i]; }
但编译器不喜欢它.在这种情况下,我不确定如何使用p_rate和p_dir指针.
另外,我不太明白为什么下面的代码是这样的.它似乎将字符串,rate和dir传递给List_names,并将值传递给rate和dir.为什么我们不能简单地做cout并打印所需的值?
char *tag[2] = {"rate","dir"}; PROTECT(List_names = allocVector(STRSXP,2)); SET_STRING_ELT(List_names,mkChar(tag[0])); SET_STRING_ELT(List_names,1,mkChar(tag[1])); PROTECT(List = allocVector(VECSXP,2)); SET_VECTOR_ELT(List,rate); SET_VECTOR_ELT(List,dir); setAttrib(List,R_namesSymbol,List_names); UNPROTECT(8); return (List);}
我没有得到在R中返回变量的基本想法,这看起来很混乱.如果有人能指出适当的资源解释这些东西,我真的很感激.
感谢您的帮助!
解决方法 如果您已经熟悉C,那么您可能会发现我们用于无缝R和C整合的 Rcpp项目非常有用.它包含大量文档和示例,而且CRAN现在有超过80个使用它的包,它提供了另一个大型的示例集.你可以从对应于我们JSS paper的“Introduction”vignette开始,也可以在这里阅读[rcpp]标签.
否则,你回到R API,这是更简单的基于C的.
总结以上是内存溢出为你收集整理的C和R接口,获得输出全部内容,希望文章能够帮你解决C和R接口,获得输出所遇到的程序开发问题。
如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)