passagemath-modules 10.6.31__cp314-cp314-musllinux_1_2_x86_64.whl

sage/libs/gsl/odeiv.pxd

@@ -0,0 +1,70 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+cdef extern from "gsl/gsl_odeiv.h":
+
+  ctypedef struct gsl_odeiv_system:
+    int (* function) (double t,  double y[], double dydt[], void * params)
+    int (* jacobian) (double t,  double y[], double * dfdy, double dfdt[], void * params)
+    size_t dimension
+    void * params
+
+  # no:
+  #define GSL_ODEIV_FN_EVAL(S,t,y,f)  (*((S)->function))(t,y,f,(S)->params)
+  #define GSL_ODEIV_JA_EVAL(S,t,y,dfdy,dfdt)  (*((S)->jacobian))(t,y,dfdy,dfdt,(S)->params)
+
+
+  ctypedef struct gsl_odeiv_step_type
+
+  ctypedef struct gsl_odeiv_step
+
+  gsl_odeiv_step_type *gsl_odeiv_step_rk2
+  gsl_odeiv_step_type *gsl_odeiv_step_rk4
+  gsl_odeiv_step_type *gsl_odeiv_step_rkf45
+  gsl_odeiv_step_type *gsl_odeiv_step_rkck
+  gsl_odeiv_step_type *gsl_odeiv_step_rk8pd
+  gsl_odeiv_step_type *gsl_odeiv_step_rk2imp
+  gsl_odeiv_step_type *gsl_odeiv_step_rk4imp
+  gsl_odeiv_step_type *gsl_odeiv_step_bsimp
+  gsl_odeiv_step_type *gsl_odeiv_step_gear1
+  gsl_odeiv_step_type *gsl_odeiv_step_gear2
+
+
+  gsl_odeiv_step * gsl_odeiv_step_alloc( gsl_odeiv_step_type * T, size_t dim)
+  int  gsl_odeiv_step_reset(gsl_odeiv_step * s)
+  void gsl_odeiv_step_free(gsl_odeiv_step * s)
+
+  char * gsl_odeiv_step_name( gsl_odeiv_step *)
+  unsigned int gsl_odeiv_step_order( gsl_odeiv_step * s)
+
+  int  gsl_odeiv_step_apply(gsl_odeiv_step *, double t, double h, double y[], double yerr[],  double dydt_in[], double dydt_out[],  gsl_odeiv_system * dydt)
+
+
+  ctypedef struct gsl_odeiv_control_type
+
+  ctypedef struct gsl_odeiv_control
+  cdef enum:
+    GSL_ODEIV_HADJ_DEC = -1
+    GSL_ODEIV_HADJ_NIL = 0
+    GSL_ODEIV_HADJ_INC = 1
+
+  gsl_odeiv_control * gsl_odeiv_control_alloc( gsl_odeiv_control_type * T)
+  int gsl_odeiv_control_init(gsl_odeiv_control * c, double eps_abs, double eps_rel, double a_y, double a_dydt)
+  void gsl_odeiv_control_free(gsl_odeiv_control * c)
+  int gsl_odeiv_control_hadjust (gsl_odeiv_control * c, gsl_odeiv_step * s,  double y0[],  double yerr[],  double dydt[], double * h)
+  char * gsl_odeiv_control_name( gsl_odeiv_control * c)
+
+
+  gsl_odeiv_control * gsl_odeiv_control_standard_new(double eps_abs, double eps_rel, double a_y, double a_dydt)
+  gsl_odeiv_control * gsl_odeiv_control_y_new(double eps_abs, double eps_rel)
+  gsl_odeiv_control * gsl_odeiv_control_yp_new(double eps_abs, double eps_rel)
+
+  gsl_odeiv_control * gsl_odeiv_control_scaled_new(double eps_abs, double eps_rel, double a_y, double a_dydt,  double scale_abs[], size_t dim)
+
+  ctypedef struct gsl_odeiv_evolve
+
+  gsl_odeiv_evolve * gsl_odeiv_evolve_alloc(size_t dim)
+  int gsl_odeiv_evolve_apply(gsl_odeiv_evolve *, gsl_odeiv_control * con, gsl_odeiv_step * step,  gsl_odeiv_system * dydt, double * t, double t1, double * h, double y[])
+  int gsl_odeiv_evolve_reset(gsl_odeiv_evolve *)
+  void gsl_odeiv_evolve_free(gsl_odeiv_evolve *)

sage/libs/gsl/permutation.pxd

@@ -0,0 +1,78 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_permutation.h":
+  # Allocation
+  gsl_permutation *  gsl_permutation_alloc(size_t n)
+
+  gsl_permutation *  gsl_permutation_calloc(size_t n)
+
+  void  gsl_permutation_init(gsl_permutation * p)
+
+  void  gsl_permutation_free(gsl_permutation * p)
+
+  # Reading and writing permutations
+  int  gsl_permutation_fread(FILE * stream, gsl_permutation * p)
+
+  int  gsl_permutation_fwrite(FILE * stream, gsl_permutation * p)
+
+  int  gsl_permutation_fscanf(FILE * stream, gsl_permutation * p)
+
+  int  gsl_permutation_fprintf(FILE * stream, gsl_permutation * p, char *format)
+
+  # Permutation properties
+  size_t  gsl_permutation_size(gsl_permutation * p)
+
+  size_t *  gsl_permutation_data(gsl_permutation * p)
+
+  int  gsl_permutation_memcpy(gsl_permutation * dest, gsl_permutation * src)
+
+  size_t  gsl_permutation_get(gsl_permutation * p, size_t i)
+
+  int  gsl_permutation_swap(gsl_permutation * p, size_t i, size_t j)
+
+  int  gsl_permutation_valid(gsl_permutation * p)
+
+  # Permutation functions
+  void  gsl_permutation_reverse(gsl_permutation * p)
+
+  int  gsl_permutation_inverse(gsl_permutation * inv, gsl_permutation * p)
+
+  int  gsl_permutation_next(gsl_permutation * p)
+
+  int  gsl_permutation_prev(gsl_permutation * p)
+
+  int  gsl_permutation_mul(gsl_permutation * p, gsl_permutation * pa, gsl_permutation * pb)
+
+  # Permutations in Cyclic Form
+  int  gsl_permutation_linear_to_canonical(gsl_permutation * q, gsl_permutation * p)
+
+  int  gsl_permutation_canonical_to_linear(gsl_permutation * p, gsl_permutation * q)
+
+  size_t  gsl_permutation_inversions(gsl_permutation * p)
+
+  size_t  gsl_permutation_linear_cycles(gsl_permutation * p)
+
+  size_t  gsl_permutation_canonical_cycles(gsl_permutation * q)
+
+
+# Applying Permutations
+cdef extern from "gsl/gsl_permute_double.h":
+  int  gsl_permute(size_t * p, double * data, size_t stride, size_t n)
+
+  int  gsl_permute_inverse(size_t * p, double * data, size_t stride, size_t n)
+
+
+cdef extern from "gsl/gsl_permute_vector_double.h":
+  int  gsl_permute_vector(gsl_permutation * p, gsl_vector * v)
+
+  int  gsl_permute_vector_inverse(gsl_permutation * p, gsl_vector * v)
+
+cdef extern from "gsl/gsl_permute_vector_complex_double.h":
+  int gsl_permute_vector_complex ( gsl_permutation * p, gsl_vector_complex * v)
+
+  int gsl_permute_vector_complex_inverse ( gsl_permutation * p, gsl_vector_complex * v)
+

sage/libs/gsl/poly.pxd

@@ -0,0 +1,40 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_poly.h":
+
+  #  Evaluate polynomial
+  double  gsl_poly_eval(double c[], int len, double x)
+
+  # divided-difference polynomials
+  int  gsl_poly_dd_init(double dd[], double xa[], double ya[], size_t size)
+
+  double  gsl_poly_dd_eval(double dd[], double xa[], size_t size, double x)
+
+  int  gsl_poly_dd_taylor(double c[], double xp, double dd[], double xa[], size_t size, double w[])
+
+  # quadratic equation
+  int  gsl_poly_solve_quadratic(double a, double b, double c, double *x0, double *x1)
+
+  int  gsl_poly_complex_solve_quadratic(double a, double b, double c, gsl_complex *z0, gsl_complex *z1)
+
+  # cubic equation
+  int  gsl_poly_solve_cubic(double a, double b, double c, double *x0, double *x1, double *x2)
+
+  int  gsl_poly_complex_solve_cubic(double a, double b, double c, gsl_complex *z0, gsl_complex *z1, gsl_complex *z2)
+
+  # General Polynomial Equations
+  ctypedef struct gsl_poly_complex_workspace:
+    size_t nc
+    double * matrix
+
+  gsl_poly_complex_workspace *  gsl_poly_complex_workspace_alloc(size_t n)
+
+  void  gsl_poly_complex_workspace_free(gsl_poly_complex_workspace * w)
+
+  int  gsl_poly_complex_solve(double * a, size_t n, gsl_poly_complex_workspace * w, gsl_complex_packed_ptr z)
+
+

sage/libs/gsl/pow_int.pxd

@@ -0,0 +1,12 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_sf_pow_int.h":
+
+  double  gsl_sf_pow_int(double x, int n)
+
+  int  gsl_sf_pow_int_e(double x, int n, gsl_sf_result * result)
+

sage/libs/gsl/psi.pxd

@@ -0,0 +1,28 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_sf_psi.h":
+
+  double  gsl_sf_psi_int(int n)
+
+  int  gsl_sf_psi_int_e(int n, gsl_sf_result * result)
+
+  double  gsl_sf_psi(double x)
+
+  int  gsl_sf_psi_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_psi_1piy(double y)
+
+  int  gsl_sf_psi_1piy_e(double y, gsl_sf_result * result)
+
+  double  gsl_sf_psi_1_int(int n)
+
+  int  gsl_sf_psi_1_int_e(int n, gsl_sf_result * result)
+
+  double  gsl_sf_psi_n(int m, double x)
+
+  int  gsl_sf_psi_n_e(int m, double x, gsl_sf_result * result)
+

sage/libs/gsl/qrng.pxd

@@ -0,0 +1,29 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+cdef extern from "gsl/gsl_qrng.h":
+
+  ctypedef struct gsl_qrng
+  ctypedef struct gsl_qrng_type
+  gsl_qrng_type * gsl_qrng_niederreiter_2
+  gsl_qrng_type * gsl_qrng_sobol
+
+  gsl_qrng *  gsl_qrng_alloc(gsl_qrng_type * T, unsigned int d)
+
+  void  gsl_qrng_free(gsl_qrng * q)
+
+  void  gsl_qrng_init(gsl_qrng * q)
+
+  int  gsl_qrng_get(gsl_qrng * q, double x[])
+
+  char *  gsl_qrng_name(gsl_qrng * q)
+
+  size_t  gsl_qrng_size(gsl_qrng * q)
+
+  void *  gsl_qrng_state(gsl_qrng * q)
+
+  int  gsl_qrng_memcpy(gsl_qrng * dest, gsl_qrng * src)
+
+  gsl_qrng *  gsl_qrng_clone(gsl_qrng * q)
+

sage/libs/gsl/random.pxd

@@ -0,0 +1,257 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_randist.h":
+  unsigned int gsl_ran_bernoulli ( gsl_rng * r, double p)
+  double gsl_ran_bernoulli_pdf ( unsigned int k, double p)
+
+  double gsl_ran_beta ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_beta_pdf ( double x,  double a,  double b)
+
+  unsigned int gsl_ran_binomial ( gsl_rng * r, double p, unsigned int n)
+  unsigned int gsl_ran_binomial_tpe ( gsl_rng * r, double pp, unsigned int n)
+  double gsl_ran_binomial_pdf ( unsigned int k,  double p,  unsigned int n)
+
+  double gsl_ran_exponential ( gsl_rng * r,  double mu)
+  double gsl_ran_exponential_pdf ( double x,  double mu)
+
+  double gsl_ran_exppow ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_exppow_pdf ( double x,  double a,  double b)
+
+  double gsl_ran_cauchy ( gsl_rng * r,  double a)
+  double gsl_ran_cauchy_pdf ( double x,  double a)
+
+  double gsl_ran_chisq ( gsl_rng * r,  double nu)
+  double gsl_ran_chisq_pdf ( double x,  double nu)
+
+  void gsl_ran_dirichlet ( gsl_rng * r,  size_t K,  double alpha[], double theta[])
+  double gsl_ran_dirichlet_pdf ( size_t K,  double alpha[],  double theta[])
+  double gsl_ran_dirichlet_lnpdf ( size_t K,  double alpha[],  double theta[])
+
+  double gsl_ran_erlang ( gsl_rng * r,  double a,  double n)
+  double gsl_ran_erlang_pdf ( double x,  double a,  double n)
+
+  double gsl_ran_fdist ( gsl_rng * r,  double nu1,  double nu2)
+  double gsl_ran_fdist_pdf ( double x,  double nu1,  double nu2)
+
+  double gsl_ran_flat ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_flat_pdf (double x,  double a,  double b)
+
+  double gsl_ran_gamma ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_gamma_int ( gsl_rng * r,  unsigned int a)
+  double gsl_ran_gamma_pdf ( double x,  double a,  double b)
+
+  double gsl_ran_gaussian ( gsl_rng * r,  double sigma)
+  double gsl_ran_gaussian_ratio_method ( gsl_rng * r,  double sigma)
+  double gsl_ran_gaussian_pdf ( double x,  double sigma)
+
+  double gsl_ran_ugaussian ( gsl_rng * r)
+  double gsl_ran_ugaussian_ratio_method ( gsl_rng * r)
+  double gsl_ran_ugaussian_pdf ( double x)
+
+  double gsl_ran_gaussian_tail ( gsl_rng * r,  double a,  double sigma)
+  double gsl_ran_gaussian_tail_pdf ( double x,  double a,  double sigma)
+
+  double gsl_ran_ugaussian_tail ( gsl_rng * r,  double a)
+  double gsl_ran_ugaussian_tail_pdf ( double x,  double a)
+
+  void gsl_ran_bivariate_gaussian ( gsl_rng * r, double sigma_x, double sigma_y, double rho, double *x, double *y)
+  double gsl_ran_bivariate_gaussian_pdf ( double x,  double y,  double sigma_x,  double sigma_y,  double rho)
+
+  double gsl_ran_landau ( gsl_rng * r)
+  double gsl_ran_landau_pdf ( double x)
+
+  unsigned int gsl_ran_geometric ( gsl_rng * r,  double p)
+  double gsl_ran_geometric_pdf ( unsigned int k,  double p)
+
+  unsigned int gsl_ran_hypergeometric ( gsl_rng * r, unsigned int n1, unsigned int n2, unsigned int t)
+  double gsl_ran_hypergeometric_pdf ( unsigned int k,  unsigned int n1,  unsigned int n2, unsigned int t)
+
+  double gsl_ran_gumbel1 ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_gumbel1_pdf ( double x,  double a,  double b)
+
+  double gsl_ran_gumbel2 ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_gumbel2_pdf ( double x,  double a,  double b)
+
+  double gsl_ran_logistic ( gsl_rng * r,  double a)
+  double gsl_ran_logistic_pdf ( double x,  double a)
+
+  double gsl_ran_lognormal ( gsl_rng * r,  double zeta,  double sigma)
+  double gsl_ran_lognormal_pdf ( double x,  double zeta,  double sigma)
+
+  unsigned int gsl_ran_logarithmic ( gsl_rng * r,  double p)
+  double gsl_ran_logarithmic_pdf ( unsigned int k,  double p)
+
+  void gsl_ran_multinomial ( gsl_rng * r,  size_t K,
+                             unsigned int N,  double p[],
+                            unsigned int n[] )
+  double gsl_ran_multinomial_pdf ( size_t K,
+                                   double p[],  unsigned int n[] )
+  double gsl_ran_multinomial_lnpdf ( size_t K,
+                              double p[],  unsigned int n[] )
+
+
+  unsigned int gsl_ran_negative_binomial ( gsl_rng * r, double p, double n)
+  double gsl_ran_negative_binomial_pdf ( unsigned int k,  double p, double n)
+
+  unsigned int gsl_ran_pascal ( gsl_rng * r, double p, unsigned int n)
+  double gsl_ran_pascal_pdf ( unsigned int k,  double p, unsigned int n)
+
+  double gsl_ran_pareto ( gsl_rng * r, double a,  double b)
+  double gsl_ran_pareto_pdf ( double x,  double a,  double b)
+
+  unsigned int gsl_ran_poisson ( gsl_rng * r, double mu)
+  void gsl_ran_poisson_array ( gsl_rng * r, size_t n, unsigned int array[],
+                              double mu)
+  double gsl_ran_poisson_pdf ( unsigned int k,  double mu)
+
+  double gsl_ran_rayleigh ( gsl_rng * r,  double sigma)
+  double gsl_ran_rayleigh_pdf ( double x,  double sigma)
+
+  double gsl_ran_rayleigh_tail ( gsl_rng * r,  double a,  double sigma)
+  double gsl_ran_rayleigh_tail_pdf ( double x,  double a,  double sigma)
+
+  double gsl_ran_tdist ( gsl_rng * r,  double nu)
+  double gsl_ran_tdist_pdf ( double x,  double nu)
+
+  double gsl_ran_fdist ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_fdist_pdf ( double x,  double a,  double b)
+
+  double gsl_ran_laplace ( gsl_rng * r,  double a)
+  double gsl_ran_laplace_pdf ( double x,  double a)
+
+  double gsl_ran_levy ( gsl_rng * r,  double c,  double alpha)
+  double gsl_ran_levy_skew ( gsl_rng * r,  double c,  double alpha,  double beta)
+
+  double gsl_ran_weibull ( gsl_rng * r,  double a,  double b)
+  double gsl_ran_weibull_pdf ( double x,  double a,  double b)
+
+  void gsl_ran_dir_2d ( gsl_rng * r, double * x, double * y)
+  void gsl_ran_dir_2d_trig_method ( gsl_rng * r, double * x, double * y)
+  void gsl_ran_dir_3d ( gsl_rng * r, double * x, double * y, double * z)
+  void gsl_ran_dir_nd ( gsl_rng * r, size_t n, double * x)
+
+  void gsl_ran_shuffle ( gsl_rng * r, void * base, size_t nmembm, size_t size)
+  int gsl_ran_choose ( gsl_rng * r, void * dest, size_t k, void * src, size_t n, size_t size)
+  void gsl_ran_sample ( gsl_rng * r, void * dest, size_t k, void * src, size_t n, size_t size)
+
+
+  ctypedef struct gsl_ran_discrete_t
+
+  gsl_ran_discrete_t * gsl_ran_discrete_preproc (size_t K,  double *P)
+  void gsl_ran_discrete_free(gsl_ran_discrete_t *g)
+  size_t gsl_ran_discrete ( gsl_rng *r,  gsl_ran_discrete_t *g)
+  double gsl_ran_discrete_pdf (size_t k,  gsl_ran_discrete_t *g)
+
+cdef extern from "gsl/gsl_cdf.h":
+  double gsl_cdf_ugaussian_P ( double x)
+  double gsl_cdf_ugaussian_Q ( double x)
+
+  double gsl_cdf_ugaussian_Pinv ( double P)
+  double gsl_cdf_ugaussian_Qinv ( double Q)
+
+  double gsl_cdf_gaussian_P ( double x,  double sigma)
+  double gsl_cdf_gaussian_Q ( double x,  double sigma)
+
+  double gsl_cdf_gaussian_Pinv ( double P,  double sigma)
+  double gsl_cdf_gaussian_Qinv ( double Q,  double sigma)
+
+  double gsl_cdf_gamma_P ( double x,  double a,  double b)
+  double gsl_cdf_gamma_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_gamma_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_gamma_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_cauchy_P ( double x,  double a)
+  double gsl_cdf_cauchy_Q ( double x,  double a)
+
+  double gsl_cdf_cauchy_Pinv ( double P,  double a)
+  double gsl_cdf_cauchy_Qinv ( double Q,  double a)
+
+  double gsl_cdf_laplace_P ( double x,  double a)
+  double gsl_cdf_laplace_Q ( double x,  double a)
+
+  double gsl_cdf_laplace_Pinv ( double P,  double a)
+  double gsl_cdf_laplace_Qinv ( double Q,  double a)
+
+  double gsl_cdf_rayleigh_P ( double x,  double sigma)
+  double gsl_cdf_rayleigh_Q ( double x,  double sigma)
+
+  double gsl_cdf_rayleigh_Pinv ( double P,  double sigma)
+  double gsl_cdf_rayleigh_Qinv ( double Q,  double sigma)
+
+  double gsl_cdf_chisq_P ( double x,  double nu)
+  double gsl_cdf_chisq_Q ( double x,  double nu)
+
+  double gsl_cdf_chisq_Pinv ( double P,  double nu)
+  double gsl_cdf_chisq_Qinv ( double Q,  double nu)
+
+  double gsl_cdf_exponential_P ( double x,  double mu)
+  double gsl_cdf_exponential_Q ( double x,  double mu)
+
+  double gsl_cdf_exppow_P(double x, double a,double b)
+
+  double gsl_cdf_exponential_Pinv ( double P,  double mu)
+  double gsl_cdf_exponential_Qinv ( double Q,  double mu)
+
+  double gsl_cdf_tdist_P ( double x,  double nu)
+  double gsl_cdf_tdist_Q ( double x,  double nu)
+
+  double gsl_cdf_tdist_Pinv ( double P,  double nu)
+  double gsl_cdf_tdist_Qinv ( double Q,  double nu)
+
+  double gsl_cdf_fdist_P ( double x,  double nu1,  double nu2)
+  double gsl_cdf_fdist_Q ( double x,  double nu1,  double nu2)
+
+  double gsl_cdf_fdist_Pinv(double x, double a, double b)
+  double gsl_cdf_fdist_Qinv(double x, double a, double b)
+
+  double gsl_cdf_beta_P ( double x,  double a,  double b)
+  double gsl_cdf_beta_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_beta_Pinv(double x, double a, double b)
+
+  double gsl_cdf_flat_P ( double x,  double a,  double b)
+  double gsl_cdf_flat_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_flat_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_flat_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_lognormal_P ( double x,  double zeta,  double sigma)
+  double gsl_cdf_lognormal_Q ( double x,  double zeta,  double sigma)
+
+  double gsl_cdf_lognormal_Pinv ( double P,  double zeta,  double sigma)
+  double gsl_cdf_lognormal_Qinv ( double Q,  double zeta,  double sigma)
+
+  double gsl_cdf_gumbel1_P ( double x,  double a,  double b)
+  double gsl_cdf_gumbel1_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_gumbel1_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_gumbel1_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_gumbel2_P ( double x,  double a,  double b)
+  double gsl_cdf_gumbel2_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_gumbel2_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_gumbel2_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_weibull_P ( double x,  double a,  double b)
+  double gsl_cdf_weibull_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_weibull_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_weibull_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_pareto_P ( double x,  double a,  double b)
+  double gsl_cdf_pareto_Q ( double x,  double a,  double b)
+
+  double gsl_cdf_pareto_Pinv ( double P,  double a,  double b)
+  double gsl_cdf_pareto_Qinv ( double Q,  double a,  double b)
+
+  double gsl_cdf_logistic_P ( double x,  double a)
+  double gsl_cdf_logistic_Q ( double x,  double a)
+
+  double gsl_cdf_logistic_Pinv ( double P,  double a)
+  double gsl_cdf_logistic_Qinv ( double Q,  double a)

sage/libs/gsl/rng.pxd

@@ -0,0 +1,100 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_rng.h":
+  cdef gsl_rng_type *gsl_rng_borosh13
+  cdef gsl_rng_type *gsl_rng_coveyou
+  cdef gsl_rng_type *gsl_rng_cmrg
+  cdef gsl_rng_type *gsl_rng_fishman18
+  cdef gsl_rng_type *gsl_rng_fishman20
+  cdef gsl_rng_type *gsl_rng_fishman2x
+  cdef gsl_rng_type *gsl_rng_gfsr4
+  cdef gsl_rng_type *gsl_rng_knuthran
+  cdef gsl_rng_type *gsl_rng_knuthran2
+  cdef gsl_rng_type *gsl_rng_lecuyer21
+  cdef gsl_rng_type *gsl_rng_minstd
+  cdef gsl_rng_type *gsl_rng_mrg
+  cdef gsl_rng_type *gsl_rng_mt19937
+  cdef gsl_rng_type *gsl_rng_mt19937_1999
+  cdef gsl_rng_type *gsl_rng_mt19937_1998
+  cdef gsl_rng_type *gsl_rng_r250
+  cdef gsl_rng_type *gsl_rng_ran0
+  cdef gsl_rng_type *gsl_rng_ran1
+  cdef gsl_rng_type *gsl_rng_ran2
+  cdef gsl_rng_type *gsl_rng_ran3
+  cdef gsl_rng_type *gsl_rng_rand
+  cdef gsl_rng_type *gsl_rng_rand48
+  cdef gsl_rng_type *gsl_rng_random128_bsd
+  cdef gsl_rng_type *gsl_rng_random128_glibc2
+  cdef gsl_rng_type *gsl_rng_random128_libc5
+  cdef gsl_rng_type *gsl_rng_random256_bsd
+  cdef gsl_rng_type *gsl_rng_random256_glibc2
+  cdef gsl_rng_type *gsl_rng_random256_libc5
+  cdef gsl_rng_type *gsl_rng_random32_bsd
+  cdef gsl_rng_type *gsl_rng_random32_glibc2
+  cdef gsl_rng_type *gsl_rng_random32_libc5
+  cdef gsl_rng_type *gsl_rng_random64_bsd
+  cdef gsl_rng_type *gsl_rng_random64_glibc2
+  cdef gsl_rng_type *gsl_rng_random64_libc5
+  cdef gsl_rng_type *gsl_rng_random8_bsd
+  cdef gsl_rng_type *gsl_rng_random8_glibc2
+  cdef gsl_rng_type *gsl_rng_random8_libc5
+  cdef gsl_rng_type *gsl_rng_random_bsd
+  cdef gsl_rng_type *gsl_rng_random_glibc2
+  cdef gsl_rng_type *gsl_rng_random_libc5
+  cdef gsl_rng_type *gsl_rng_randu
+  cdef gsl_rng_type *gsl_rng_ranf
+  cdef gsl_rng_type *gsl_rng_ranlux
+  cdef gsl_rng_type *gsl_rng_ranlux389
+  cdef gsl_rng_type *gsl_rng_ranlxd1
+  cdef gsl_rng_type *gsl_rng_ranlxd2
+  cdef gsl_rng_type *gsl_rng_ranlxs0
+  cdef gsl_rng_type *gsl_rng_ranlxs1
+  cdef gsl_rng_type *gsl_rng_ranlxs2
+  cdef gsl_rng_type *gsl_rng_ranmar
+  cdef gsl_rng_type *gsl_rng_slatec
+  cdef gsl_rng_type *gsl_rng_taus
+  cdef gsl_rng_type *gsl_rng_taus2
+  cdef gsl_rng_type *gsl_rng_taus113
+  cdef gsl_rng_type *gsl_rng_transputer
+  cdef gsl_rng_type *gsl_rng_tt800
+  cdef gsl_rng_type *gsl_rng_uni
+  cdef gsl_rng_type *gsl_rng_uni32
+  cdef gsl_rng_type *gsl_rng_vax
+  cdef gsl_rng_type *gsl_rng_waterman14
+  cdef gsl_rng_type *gsl_rng_zuf
+
+
+  cdef gsl_rng_type *gsl_rng_default
+  unsigned long int gsl_rng_default_seed
+
+  gsl_rng *gsl_rng_alloc ( gsl_rng_type * T)
+  int gsl_rng_memcpy (gsl_rng * dest, gsl_rng * src)
+  gsl_rng *gsl_rng_clone ( gsl_rng * r)
+
+  void gsl_rng_free (gsl_rng * r)
+
+  void gsl_rng_set ( gsl_rng * r, unsigned long int seed)
+  unsigned long int gsl_rng_max ( gsl_rng * r)
+  unsigned long int gsl_rng_min ( gsl_rng * r)
+  char *gsl_rng_name ( gsl_rng * r)
+
+  int gsl_rng_fread (FILE * stream, gsl_rng * r)
+  int gsl_rng_fwrite (FILE * stream, gsl_rng * r)
+
+  size_t gsl_rng_size ( gsl_rng * r)
+  void * gsl_rng_state ( gsl_rng * r)
+
+  void gsl_rng_print_state ( gsl_rng * r)
+
+  gsl_rng_type * gsl_rng_env_setup ()
+
+  unsigned long int gsl_rng_get ( gsl_rng * r)
+  double gsl_rng_uniform ( gsl_rng * r)
+  double gsl_rng_uniform_pos ( gsl_rng * r)
+  unsigned long int gsl_rng_uniform_int ( gsl_rng * r, unsigned long int n)
+
+

sage/libs/gsl/roots.pxd

@@ -0,0 +1,72 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: libraries = GSL_LIBRARIES
+# distutils: library_dirs = GSL_LIBDIR
+# distutils: include_dirs = GSL_INCDIR
+from sage.libs.gsl.types cimport *
+
+cdef extern from "gsl/gsl_roots.h":
+
+  ctypedef struct gsl_root_fsolver_type:
+      char *name
+      size_t size
+      int (*set) (void *state, gsl_function * f, double * root, double x_lower, double x_upper)
+      int (*iterate) (void *state, gsl_function * f, double * root, double * x_lower, double * x_upper)
+
+  ctypedef struct gsl_root_fsolver:
+      gsl_root_fsolver_type * type
+      gsl_function * function
+      double root
+      double x_lower
+      double x_upper
+      void *state
+
+  ctypedef struct gsl_root_fdfsolver_type:
+      char *name
+      size_t size
+      int (*set) (void *state, gsl_function_fdf * f, double * root)
+      int (*iterate) (void *state, gsl_function_fdf * f, double * root)
+
+  ctypedef struct gsl_root_fdfsolver:
+      gsl_root_fdfsolver_type * type
+      gsl_function_fdf * fdf
+      double root
+      void *state
+
+  gsl_root_fsolver * gsl_root_fsolver_alloc ( gsl_root_fsolver_type * T)
+  void gsl_root_fsolver_free (gsl_root_fsolver * s)
+
+  int gsl_root_fsolver_set (gsl_root_fsolver * s, gsl_function * f,
+                            double x_lower, double x_upper)
+
+  int gsl_root_fsolver_iterate (gsl_root_fsolver * s)
+
+  char * gsl_root_fsolver_name ( gsl_root_fsolver * s)
+  double gsl_root_fsolver_root ( gsl_root_fsolver * s)
+  double gsl_root_fsolver_x_lower ( gsl_root_fsolver * s)
+  double gsl_root_fsolver_x_upper ( gsl_root_fsolver * s)
+
+
+  gsl_root_fdfsolver * gsl_root_fdfsolver_alloc ( gsl_root_fdfsolver_type * T)
+
+  int gsl_root_fdfsolver_set (gsl_root_fdfsolver * s,
+                           gsl_function_fdf * fdf, double root)
+
+  int gsl_root_fdfsolver_iterate (gsl_root_fdfsolver * s)
+
+  void gsl_root_fdfsolver_free (gsl_root_fdfsolver * s)
+
+  char * gsl_root_fdfsolver_name ( gsl_root_fdfsolver * s)
+  double gsl_root_fdfsolver_root ( gsl_root_fdfsolver * s)
+
+  int gsl_root_test_interval (double x_lower, double x_upper, double epsabs, double epsrel)
+
+  int gsl_root_test_residual (double f, double epsabs)
+
+  int gsl_root_test_delta (double x1, double x0, double epsabs, double epsrel)
+
+  gsl_root_fsolver_type  * gsl_root_fsolver_bisection
+  gsl_root_fsolver_type  * gsl_root_fsolver_brent
+  gsl_root_fsolver_type  * gsl_root_fsolver_falsepos
+  gsl_root_fdfsolver_type  * gsl_root_fdfsolver_newton
+  gsl_root_fdfsolver_type  * gsl_root_fdfsolver_secant
+  gsl_root_fdfsolver_type  * gsl_root_fdfsolver_steffenson