关于R和C的混合编程

因为R是解释类型语言，所以考虑到效率问题我们多采用主要部分用R、迭代部分用C。但是我不知道如何将R和C进行混合编程。
譬如我们求1+2^3+3^3+4^3+....+n^3 (当然这个用R的向量计算也相当方便)，我们可以很容易地写出用R循环来计算该问题的代码。
sum=function(n){
for(i in 1:n)
循环部分
}

下面我们写出用C来计算循环
#include<stdio.h>
main()
{int sum=0,n,i;
scanf("%d",&n);
for(i=1;i<=n;i++)
sum+=i*i*i;
printf("\n%d",sum);
getch();
}
我们将上面的c代码命名为sum.c
我电脑上只安装了win-TC，上述代码运行通过（我输入3，返回36），并且产生了三个文件，分别是：
sum.c和sum.exe及sum.obj
好了，现在我的问题就是如何将c嵌入到R中去代替R的循环计算。
希望高手解答，以便以后碰到更复杂的循环可以节省大量时间，R有时候还是稍显慢。

1.你可选择自行安装MinGW (GCC) Compiler
   http://www.mingw.org/wiki/InstallationHOWTOforMinGW
  也可选择安装 Rtools (建议)
   http://www.murdoch-sutherland.com/Rtools/index.html


2. C:/Users/lenovo/Documents

建立文件夹convolve,放入文件convolve.c

然后在R运行:
############
setwd('C:/Users/lenovo/Documents/convolve')
system("R CMD SHLIB convolve.c")  #产生两个文件convolve.dll,convolve.o
                                  #windows 用到convolve.dll

dyn.load("convolve.dll")
conv <- function(a, b){
.C("convolve",
as.double(a), as.integer(length(a)),
as.double(b), as.integer(length(b)),
ab=double(length(a)+length(b)-1))$ab
}
a=3
b=4
conv(a,b)

##########
convolve.c, convolve.dll
convolve.rar (3.2 KB)

2011-6-30 09:40:40 上传

R中自带的两个函数.c()  .call()可以用来处理LZ想运行的C函数，两者的使用略有差别。比如lz可以另加一个wrapper函数，其中的变量类型用C和R接口的变量类型（具体参见CRAN manual），用这个函数来调用你写好的C函数。在R中运行的时候，通过.call()在R中调用wrapper函数。wrapper的作用主要就是让R和C函数中的变量可以相互转换使用，相当于一个接口。头文件可能需要稍做改动。

另外，更简单的方式是用package Rcpp 。 这个包绝对是一个非常出色的产品，LZ可以下载来看看，可以用于C或C++。

关于你所说的三个文件，在用上述方式调用的时候，不会产生exe等文件，你要做的就是有源文件即可。编译的工作由R调用你系统下的C编译器运行，结果直接返回R。我不用windows，不清楚在windows下面跑对编译器有没有什么特殊要求，但是常用的编译器应该可以正常工作。

只能用GCC。

版主正解，Rcpp有PDF帮助文件，看完就会了。

首先感谢版主ltx5151 。但是我给出具体例子的初衷还是想请教高人结合具体例子给出详细操作，我也有看过不少资料，其中李东风书上给出了一个计算卷积的例子，算是写的比较详细的。但是我不知道“C程序编译成了一个DLL库testdll.dll放在当前工作目录中”这句话，因为后面好像要用到testdll.dll。我用Win-TC编译C源程序得到了sum.c、sum.exe和sum.obj三个文件（就上面例子而言），并没有什么.dll文件啊。所以，我还是等等看。谨再谢！

我把它们贴出来吧
Ｃ计算卷积的源程序为：
void convolve(double *a, int *na,
double *b, int *nb, double *ab){
int i, j, nab = *na + *nb - 1;
for(i = 0; i < nab; i++)
ab[i] = 0.0;
for(i = 0; i < *na; i++)
for(j = 0; j < *nb; j++)
ab[i + j] += a[i] * b[j];
}
为了在Ｒ中调用这个Ｃ函数需要如下步骤：
conv <- function(a, b){
dll.filename <-
file.path(paste("testdll",
.Platform$dynlib.ext, sep=""))
if(!is.loaded(symbol.C("convolve"),
PACKAGE="testdll")){
dyn.load(dll.filename)
}
.C("convolve",
as.double(a), as.integer(length(a)),
as.double(b), as.integer(length(b)),
ab=double(length(a)+length(b)-1))$ab
}

是不是要用什么特别的Ｃ编辑器才能产生ｄｌｌ文件？

把R里面运行的结果也贴出来吧：
> conv <- function(a, b){
+ dll.filename <-
+ file.path(paste("testdll",
+ .Platform$dynlib.ext, sep=""))
+ if(!is.loaded(symbol.C("convolve"),
+ PACKAGE="testdll")){
+ dyn.load(dll.filename)
+ }
+ .C("convolve",
+ as.double(a), as.integer(length(a)),
+ as.double(b), as.integer(length(b)),
+ ab=double(length(a)+length(b)-1))$ab
+ }
> a=3
> b=4
> conv(a,b)
symbol.C("convolve")里有警告:
  'symbol.C' is not needed: please remove it
错误于inDL(x, as.logical(local), as.logical(now), ...) :
  无法载入共享目标对象‘C:/Users/lenovo/Documents/testdll.dll’：:
  LoadLibrary failure:  找不到指定的模块。
>

1.你可选择自行安装MinGW (GCC) Compiler
   http://www.mingw.org/wiki/InstallationHOWTOforMinGW
  也可选择安装 Rtools (建议)
   http://www.murdoch-sutherland.com/Rtools/index.html


2. C:/Users/lenovo/Documents

建立文件夹convolve,放入文件convolve.c

然后在R运行:
###########
setwd('C:/Users/lenovo/Documents/convolve')
system("R CMD SHLIB convolve.c")  #产生两个文件convolve.dll,convolve.o
                                  #windows 用到convolve.dll
dyn.load("convolve.dll")
conv <- function(a, b){
.C("convolve",
as.double(a), as.integer(length(a)),
as.double(b), as.integer(length(b)),
ab=double(length(a)+length(b)-1))$ab
}
a=3
b=4
conv(a,b)
###########