5
4x4 크기의 정사각형 행렬을 선언하는 C 코드를 작성했습니다. 그리고 나서 R의 rgig
이라는 샘플링 함수로 샘플링을합니다. GeneralizedHyperbolic
은 gnu의 gsl 라이브러리를 사용하여 행렬을 반전하고 그 결과를 출력합니다.C에서 임베디드 R을 실행 중
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 SHLIB -lgsl -lgslcblas embedR_matinv.c
:이
#include <stdio.h>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <stddef.h>
// for gsl
#include <gsl/gsl_machine.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_cdf.h>
#include <gsl/gsl_cblas.h>
#include <gsl/gsl_sf_gamma.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_linalg.h>
// for R embedding in C
#include <Rinternals.h>
#include <Rdefines.h>
#include <Rembedded.h>
#include <R_ext/Parse.h>
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);
}
내가 사용하여 코드를 컴파일 C.에서 호출 R에있는 운동이다
R CMD embedR_matinv
그 오류를 제공합니다 :
/Library/Frameworks/R.framework/Resources/bin/Rcmd: line 62: exec: embedR_matinv: not found
내가 뭘 잘못하고 있니?
는 또한 test()
에 main()
을 변경하고 출력
R CMD SHLIB -lgsl -lgslcblas embedR_matinv.c -o embedR_matinv
로 공유 객체를 만든 : 해지하지 않고 내가 R Studio에서
dyn.load("embedR_matinv.so")
을 할 경우
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
, 코드가 실행 즉 멈춘다!
코드의 문제점에 대한 제안 사항이 있으십니까?
무엇 당신이'embedR_matinv.c'를 컴파일 할 때 생성됩니다
이 시도? 'R CMD embedR_matinv.o' 나'R CMD embedR_matinv.so'가 아닌'R CMD embedR_matinv'을 실행해야합니까? 어쩌면 생성 된 파일의 전체 경로를 지정해야합니까? – mangusta
코드가 중단 된 위치는 어디입니까? 내장 된 R 코드의 R 인터프리터에 붙어있는 것과 같은 냄새가납니다. –