OSDN Git Service

Please enter the commit message for your changes. Lines starting
[eos/hostdependX86MAC64.git] / util / X86MAC64 / include / gsl / gsl_randist.h
index 6f4b0e3..d38ccb3 100644 (file)
@@ -20,6 +20,8 @@
 #ifndef __GSL_RANDIST_H__
 #define __GSL_RANDIST_H__
 #include <gsl/gsl_rng.h>
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_matrix.h>
 
 #undef __BEGIN_DECLS
 #undef __END_DECLS
@@ -93,6 +95,38 @@ double gsl_ran_ugaussian_tail_pdf (const double x, const double a);
 void gsl_ran_bivariate_gaussian (const gsl_rng * r, double sigma_x, double sigma_y, double rho, double *x, double *y);
 double gsl_ran_bivariate_gaussian_pdf (const double x, const double y, const double sigma_x, const double sigma_y, const double rho);
 
+int gsl_ran_multivariate_gaussian (const gsl_rng * r, const gsl_vector * mu, const gsl_matrix * L, gsl_vector * result);
+int gsl_ran_multivariate_gaussian_log_pdf (const gsl_vector * x,
+                                           const gsl_vector * mu,
+                                           const gsl_matrix * L,
+                                           double * result,
+                                           gsl_vector * work);
+int gsl_ran_multivariate_gaussian_pdf (const gsl_vector * x,
+                                       const gsl_vector * mu,
+                                       const gsl_matrix * L,
+                                       double * result,
+                                       gsl_vector * work);
+int gsl_ran_multivariate_gaussian_mean (const gsl_matrix * X, gsl_vector * mu_hat);
+int gsl_ran_multivariate_gaussian_vcov (const gsl_matrix * X, gsl_matrix * sigma_hat);
+
+int gsl_ran_wishart (const gsl_rng * r,
+                     const double df,
+                     const gsl_matrix * L,
+                     gsl_matrix * result,
+                     gsl_matrix * work);
+int gsl_ran_wishart_log_pdf (const gsl_matrix * X,
+                             const gsl_matrix * L_X,
+                             const double df,
+                             const gsl_matrix * L,
+                             double * result,
+                             gsl_matrix * work);
+int gsl_ran_wishart_pdf (const gsl_matrix * X,
+                         const gsl_matrix * L_X,
+                         const double df,
+                         const gsl_matrix * L,
+                         double * result,
+                         gsl_matrix * work);
+
 double gsl_ran_landau (const gsl_rng * r);
 double gsl_ran_landau_pdf (const double x);