Displaying 1 result from an estimated 1 matches for "gsl_mb_mvnorm_logpdf".
2007 Jan 04
1
Parameter changes and segfault when calling C code through .Call
...gsl_vector_view_array(xAr,k);
gsl_vector_view muView = gsl_vector_view_array(muAr,k);
gsl_matrix_view sigmaView = gsl_matrix_view_array(sigmaAr,k,k);
gsl_vector * x = &xView.vector;
gsl_vector * mu = &muView.vector;
gsl_matrix * sigma = &sigmaView.matrix;
1: double logans = gsl_MB_mvnorm_logpdf(x, mu, sigma, k); // <-call logpdf here
PROTECT(res=allocVector(REALSXP,1));
REAL(res)[0] = logans;
UNPROTECT(1);
return(res);
}
The logpdf function is here
double gsl_MB_mvnorm_logpdf(gsl_vector * beta, gsl_vector * betaMean, gsl_matrix * sigma, int k) {
// computes density of mu...