1

我一直在研究我用 C 语言编写的一个程序，该程序使用 GSL 中的函数运行 MCMC。我已经阅读了很多 GNU 文档和编写 R 扩展，并且我已经阅读了很多关于使用 RcppGSL 的内容，但在我看来，用 C 编写然后将函数动态加载到 R 中更容易。我见过许多资源描述了如何在 Windows 上构建该功能，但我看到的关于如何在基于 Unix 的系统上使用它的一切都是它“相对简单”和“简单”，但我无法让它工作。

我渴望有人开始工作的 C 脚本是一个简单的脚本。它只是取一个具有给定行数和列数的数组，将其转换为 gsl_matrix，然后将其转换回数组（这基本上是我的程序所做的，除了整个改变条目是什么）。C脚本是

#include <R.h>
#include <stdio.h>
#include <gsl/gsl_linalg.h>

void simple( int *n, int *rows, int *cols)
{
  int r,c;
  int Cols = *cols; //This step, and the step below it are unnecessary except for readability
  int Rows = *rows;
  gsl_matrix * m = gsl_matrix_alloc (Rows, Cols); // Declares a gsl_matrix m of size Rows x Cols

  for( r = 0; r< Rows; r++)
    for (c=0; c< Cols; c++)
      gsl_matrix_set(m,r,c, *(n+sizeof(int)*(Cols*r+c))); // The array is organized by rows, sets matrix values

  for( r=0; r< Rows; r++)
    for (c=0; c < Cols; c++)
      *(n+sizeof(int)*(Cols*r+c))=gsl_matrix_get(m,r,c); // This return matrix values to the array (should be the same as before)
}

我所做的是确保我的终端（Mac OS X）和 R 共享保存 simple.c 的同一目录。我编译上面的打字

R CMD SHLIB simple.c

进入我的终端创建对应的simple.so文件，也保存到桌面。然后在 R i 可以做

dyn.load("simple.so")

这是我收到错误的地方

Symbol not found: _gsl_matrix_alloc 
Expected in: flat namespace

我对这个错误感到困惑，因为编译器似乎识别了 make 文件< gsl/gsl_linalg.h >，所以我认为它无法识别函数的原因是因为 R 没有连接到库，但我不知道如何解决这个问题。

此时，如果函数都可以识别，那么我可以在 R 中执行该函数

x=.C("simple", c(as.integer(c(1,4,7,2,5,8,3,6,9)),as.integer(3),as.integer(3)))

如果我的函数动态加载正确，它将起作用，并且我将返回 X 正是我作为 .C 的第二个输入输入的内容

任何见解都会非常有帮助，无论是否有人成功地将 GSL 库链接到 Mac OS X 上的 R 都可以给我发消息或评论将不胜感激。我在帮助论坛上唯一能找到的是它“直截了当”或“相对简单”，但我不知道该怎么做！请帮忙！

4

1 回答 1

3

您需要指定 gsl 函数的位置：

R CMD SHLIB simple.c -lgsl -lgslcblas
于 2012-07-19T07:07:35.397 回答