在C中运行嵌入式R.
我写了一段C代码，声明了一个大小为4×4的方阵。 然后它从R中的一个名为rgig的采样函数中进行采样。它使用来自gnu的gsl库反转矩阵并吐出结果。 这是从C调用R的练习。
#include #include #include #include #include // for gsl #include #include #include #include #include #include #include #include #include #include // for R embedding in C #include #include #include #include void gsl_square_matrix_inverse (gsl_matrix *, gsl_matrix *, int); SEXP get_rInvGauss(void); int main(void) { // Define the dimension n of the matrix // and the signum s (for LU decomposition) int s, i, j, n = 4; // Define all the used matrices gsl_matrix * m = gsl_matrix_alloc (n, n); gsl_matrix * inverse = gsl_matrix_alloc (n, n); // R embedding in C char *localArgs[] = {"R", "--no-save","--silent"}; SEXP rInvGauss; // init R embedding Rf_initEmbeddedR(3, localArgs); printf("\n Printing matrix m before set. size %d by %d... \n", n, n); for(i=0; i<n; i++){ printf("\n"); for(j=0; j<n; j++){ printf(" %f ", gsl_matrix_get(m, i, j)); } } // set diagonal elements of matrix m from Inverse Gaussian Random samples for(i=0; i<n; i++){ rInvGauss = get_rInvGauss(); gsl_matrix_set(m, i, i, *REAL(rInvGauss)); } Rf_endEmbeddedR(0); // end the R embedding in C printf("\n Printing matrix m ..... \n"); for(i=0; i<n; i++){ printf("\n"); for(j=0; j<n; j++){ printf(" %f ", gsl_matrix_get(m, i, j)); } } // inverse of matrix m gsl_square_matrix_inverse (m, inverse, n); printf("\n Printing inverse of matrix m ..... \n"); for(i=0; i<n; i++){ printf("\n"); for(j=0; j<n; j++){ printf(" %f", gsl_matrix_get(inverse, i, j)); } } return 0; } SEXP get_rInvGauss(void) { SEXP e, s, t, tmp, result; int errorOccurred, n=1; double chi=5, psi=4, lambda=0.5; // create and evaluate 'require(GeneralizedHyperbolic)' PROTECT(e = lang2(install("require"), mkString("GeneralizedHyperbolic"))); R_tryEval(e, R_GlobalEnv, &errorOccurred); if (errorOccurred) { // handle error printf("\n Error loading library GeneralizedHyperbolic:"); } UNPROTECT(1); // Create the R expressions using a paired list // rgig(n = 1, chi = 5, psi = 4, lambda = 0.5) with the R API. PROTECT(t = s = allocVector(LANGSXP, 5)); // could also be done by: PROTECT(t = s = allocList(5)); SET_TYPEOF(s, LANGSXP); tmp = findFun(install("rgig"), R_GlobalEnv); if(tmp == R_NilValue) { printf("No definition for function rgig.\n"); UNPROTECT(1); exit(1); } SETCAR(t, tmp); t = CDR(t); SETCAR(t, ScalarInteger(n)); SET_TAG(t, install("n")); t= CDR(t); SETCAR(t, ScalarReal(chi)); SET_TAG(t, install("chi")); t= CDR(t); SETCAR(t, ScalarReal(psi)); SET_TAG(t, install("psi")); t= CDR(t); SETCAR(t, ScalarReal(lambda)); SET_TAG(t, install("lambda")); t= CDR(t); PROTECT(result = R_tryEval(VECTOR_ELT(s, 0), R_GlobalEnv, NULL)); UNPROTECT(2); return(result); } void gsl_square_matrix_inverse (gsl_matrix *m, gsl_matrix *inverse, int n){ int s, i, j; gsl_permutation * perm = gsl_permutation_alloc (n); // Make LU decomposition of matrix m gsl_linalg_LU_decomp (m, perm, &s); // Invert the matrix m gsl_linalg_LU_invert (m, perm, inverse); }
我使用以下代码编译代码：
R CMD SHLIB -lgsl -lgslcblas embedR_matinv.c
输出：
gcc -arch x86_64 -std=gnu99 -I/Library/Frameworks/R.framework/Resources/include -I/Library/Frameworks/R.framework/Resources/include/x86_64 -DNDEBUG -I/usr/local/include -fPIC -g -O2 -c embedR_matinv.c -o embedR_matinv.o gcc -arch x86_64 -std=gnu99 -dynamiclib -Wl,-headerpad_max_install_names -undefined dynamic_lookup -single_module -multiply_defined suppress -L/usr/local/lib -o embedR_matinv.so embedR_matinv.o -lgsl -lgslcblas -F/Library/Frameworks/R.framework/.. -framework R -Wl,-framework -Wl,CoreFoundation
当我提交时：
R CMD embedR_matinv
它给出了错误：
/Library/Frameworks/R.framework/Resources/bin/Rcmd: line 62: exec: embedR_matinv: not found
我究竟做错了什么？
我还将main()更改为test()并将共享对象更改为
R CMD SHLIB -lgsl -lgslcblas embedR_matinv.c -o embedR_matinv
与输出：
gcc -arch x86_64 -std=gnu99 -I/Library/Frameworks/R.framework/Resources/include -I/Library/Frameworks/R.framework/Resources/include/x86_64 -DNDEBUG -I/usr/local/include -fPIC -g -O2 -c embedR_matinv.c -o embedR_matinv.o
如果我在R Studio中执行dyn.load("embedR_matinv.so") ，代码运行时没有任何终止，即它挂起！
关于代码中有什么问题的任何建议？
您使用的是Mac OSX系统吗？
试试这个：
sudo cp /Library/Frameworks/R.framework/Resources/library/Rserve/libs/x86_64/Rserve-bin.so\ /Library/Frameworks/R.framework/Resources/bin/Rserve