passagemath-modules 10.6.31__cp314-cp314-macosx_13_0_x86_64.whl

sage/libs/gsl/bessel.pxd

@@ -0,0 +1,208 @@
+# 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_bessel.h":
+
+  double  gsl_sf_bessel_J0(double x)
+
+  int  gsl_sf_bessel_J0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_J1(double x)
+
+  int  gsl_sf_bessel_J1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Jn(int n, double x)
+
+  int  gsl_sf_bessel_Jn_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_Jn_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_Y0(double x)
+
+  int  gsl_sf_bessel_Y0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Y1(double x)
+
+  int  gsl_sf_bessel_Y1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Yn(int n, double x)
+
+  int  gsl_sf_bessel_Yn_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_Yn_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_I0(double x)
+
+  int  gsl_sf_bessel_I0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_I1(double x)
+
+  int  gsl_sf_bessel_I1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_In(int n, double x)
+
+  int  gsl_sf_bessel_In_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_In_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_I0_scaled(double x)
+
+  int  gsl_sf_bessel_I0_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_I1_scaled(double x)
+
+  int  gsl_sf_bessel_I1_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_In_scaled(int n, double x)
+
+  int  gsl_sf_bessel_In_scaled_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_In_scaled_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_K0(double x)
+
+  int  gsl_sf_bessel_K0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_K1(double x)
+
+  int  gsl_sf_bessel_K1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Kn(int n, double x)
+
+  int  gsl_sf_bessel_Kn_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_Kn_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_K0_scaled(double x)
+
+  int  gsl_sf_bessel_K0_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_K1_scaled(double x)
+
+  int  gsl_sf_bessel_K1_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Kn_scaled(int n, double x)
+
+  int  gsl_sf_bessel_Kn_scaled_e(int n, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_Kn_scaled_array(int nmin, int nmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_j0(double x)
+
+  int  gsl_sf_bessel_j0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_j1(double x)
+
+  int  gsl_sf_bessel_j1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_j2(double x)
+
+  int  gsl_sf_bessel_j2_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_jl(int l, double x)
+
+  int  gsl_sf_bessel_jl_e(int l, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_jl_array(int lmax, double x, double result_array[])
+
+  int  gsl_sf_bessel_jl_steed_array(int lmax, double x, double * jl_x_array)
+
+  double  gsl_sf_bessel_y0(double x)
+
+  int  gsl_sf_bessel_y0_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_y1(double x)
+
+  int  gsl_sf_bessel_y1_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_y2(double x)
+
+  int  gsl_sf_bessel_y2_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_yl(int l, double x)
+
+  int  gsl_sf_bessel_yl_e(int l, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_yl_array(int lmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_i0_scaled(double x)
+
+  int  gsl_sf_bessel_i0_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_i1_scaled(double x)
+
+  int  gsl_sf_bessel_i1_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_i2_scaled(double x)
+
+  int  gsl_sf_bessel_i2_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_il_scaled(int l, double x)
+
+  int  gsl_sf_bessel_il_scaled_e(int l, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_il_scaled_array(int lmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_k0_scaled(double x)
+
+  int  gsl_sf_bessel_k0_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_k1_scaled(double x)
+
+  int  gsl_sf_bessel_k1_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_k2_scaled(double x)
+
+  int  gsl_sf_bessel_k2_scaled_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_kl_scaled(int l, double x)
+
+  int  gsl_sf_bessel_kl_scaled_e(int l, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_kl_scaled_array(int lmax, double x, double result_array[])
+
+  double  gsl_sf_bessel_Jnu(double nu, double x)
+
+  int  gsl_sf_bessel_Jnu_e(double nu, double x, gsl_sf_result * result)
+
+  int  gsl_sf_bessel_sequence_Jnu_e(double nu, gsl_mode_t mode, size_t size, double v[])
+
+  double  gsl_sf_bessel_Ynu(double nu, double x)
+
+  int  gsl_sf_bessel_Ynu_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Inu(double nu, double x)
+
+  int  gsl_sf_bessel_Inu_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Inu_scaled(double nu, double x)
+
+  int  gsl_sf_bessel_Inu_scaled_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Knu(double nu, double x)
+
+  int  gsl_sf_bessel_Knu_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_lnKnu(double nu, double x)
+
+  int  gsl_sf_bessel_lnKnu_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_Knu_scaled(double nu, double x)
+
+  int  gsl_sf_bessel_Knu_scaled_e(double nu, double x, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_zero_J0(unsigned int s)
+
+  int  gsl_sf_bessel_zero_J0_e(unsigned int s, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_zero_J1(unsigned int s)
+
+  int  gsl_sf_bessel_zero_J1_e(unsigned int s, gsl_sf_result * result)
+
+  double  gsl_sf_bessel_zero_Jnu(double nu, unsigned int s)
+
+  int  gsl_sf_bessel_zero_Jnu_e(double nu, unsigned int s, gsl_sf_result * result)
+

sage/libs/gsl/blas.pxd

@@ -0,0 +1,116 @@
+# 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_blas.h":
+
+  # Level 1
+  int  gsl_blas_ddot(gsl_vector * x, gsl_vector * y, double * result)
+
+  int  gsl_blas_zdotu(gsl_vector_complex * x, gsl_vector_complex * y, gsl_complex * dotu)
+
+  int  gsl_blas_zdotc(gsl_vector_complex * x, gsl_vector_complex * y, gsl_complex * dotc)
+
+  double  gsl_blas_dnrm2(gsl_vector * x)
+
+  double  gsl_blas_dasum(gsl_vector * x)
+
+  double  gsl_blas_dznrm2(gsl_vector_complex * x)
+
+  double  gsl_blas_dzasum(gsl_vector_complex * x)
+
+  CBLAS_INDEX_t  gsl_blas_idamax(gsl_vector * x)
+
+  CBLAS_INDEX_t  gsl_blas_izamax(gsl_vector_complex * x)
+
+  int  gsl_blas_dswap(gsl_vector * x, gsl_vector * y)
+
+  int  gsl_blas_dcopy(gsl_vector * x, gsl_vector * y)
+
+  int  gsl_blas_daxpy(double alpha, gsl_vector * x, gsl_vector * y)
+
+  int  gsl_blas_zswap(gsl_vector_complex * x, gsl_vector_complex * y)
+
+  int  gsl_blas_zcopy(gsl_vector_complex * x, gsl_vector_complex * y)
+
+  int  gsl_blas_zaxpy(gsl_complex alpha, gsl_vector_complex * x, gsl_vector_complex * y)
+
+  int  gsl_blas_drotg(double a[], double b[], double c[], double s[])
+
+  int  gsl_blas_drotmg(double d1[], double d2[], double b1[], double b2, double P[])
+
+  int  gsl_blas_drot(gsl_vector * x, gsl_vector * y, double c, double s)
+
+  int  gsl_blas_drotm(gsl_vector * x, gsl_vector * y, double P[])
+
+  void  gsl_blas_dscal(double alpha, gsl_vector * x)
+
+  void  gsl_blas_zscal(gsl_complex alpha, gsl_vector_complex * x)
+
+  void  gsl_blas_zdscal(double alpha, gsl_vector_complex * x)
+
+
+  # Level 2
+  int  gsl_blas_dgemv(CBLAS_TRANSPOSE_t TransA, double alpha, gsl_matrix * A, gsl_vector * x, double beta, gsl_vector * y)
+
+  int  gsl_blas_dtrmv(CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, gsl_matrix * A, gsl_vector * x)
+
+  int  gsl_blas_dtrsv(CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, gsl_matrix * A, gsl_vector * x)
+
+  int  gsl_blas_zgemv (CBLAS_TRANSPOSE_t TransA,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_vector_complex * X,  gsl_complex beta, gsl_vector_complex * Y)
+
+  int  gsl_blas_ztrmv (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag,  gsl_matrix_complex * A, gsl_vector_complex * X)
+
+  int  gsl_blas_ztrsv (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag,  gsl_matrix_complex * A, gsl_vector_complex *X)
+
+  int  gsl_blas_dsymv(CBLAS_UPLO_t Uplo, double alpha, gsl_matrix * A, gsl_vector * x, double beta, gsl_vector * y)
+
+  int  gsl_blas_dger(double alpha, gsl_vector * x, gsl_vector * y, gsl_matrix * A)
+
+  int  gsl_blas_dsyr(CBLAS_UPLO_t Uplo, double alpha, gsl_vector * x, gsl_matrix * A)
+
+  int  gsl_blas_dsyr2(CBLAS_UPLO_t Uplo, double alpha, gsl_vector * x, gsl_vector * y, gsl_matrix * A)
+
+  int  gsl_blas_zhemv (CBLAS_UPLO_t Uplo,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_vector_complex * X,  gsl_complex beta, gsl_vector_complex * Y)
+
+  int  gsl_blas_zgeru ( gsl_complex alpha,  gsl_vector_complex * X,  gsl_vector_complex * Y, gsl_matrix_complex * A)
+
+  int  gsl_blas_zgerc ( gsl_complex alpha,  gsl_vector_complex * X,  gsl_vector_complex * Y, gsl_matrix_complex * A)
+
+  int  gsl_blas_zher (CBLAS_UPLO_t Uplo, double alpha,  gsl_vector_complex * X, gsl_matrix_complex * A)
+
+  int  gsl_blas_zher2 (CBLAS_UPLO_t Uplo,  gsl_complex alpha,  gsl_vector_complex * X,  gsl_vector_complex * Y, gsl_matrix_complex * A)
+
+
+  # Level 3
+  int  gsl_blas_dgemm(CBLAS_TRANSPOSE_t TransA, CBLAS_TRANSPOSE_t TransB, double alpha, gsl_matrix * A, gsl_matrix * B, double beta, gsl_matrix * C)
+
+  int  gsl_blas_dsymm(CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, double alpha, gsl_matrix * A, gsl_matrix * B, double beta, gsl_matrix * C)
+
+  int  gsl_blas_dsyrk(CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, double alpha, gsl_matrix * A, double beta, gsl_matrix * C)
+
+  int  gsl_blas_dsyr2k(CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, double alpha, gsl_matrix * A, gsl_matrix * B, double beta, gsl_matrix * C)
+
+  int  gsl_blas_dtrmm(CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, double alpha, gsl_matrix * A, gsl_matrix * B)
+
+  int  gsl_blas_dtrsm(CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, double alpha, gsl_matrix * A, gsl_matrix * B)
+
+  int  gsl_blas_zgemm (CBLAS_TRANSPOSE_t TransA, CBLAS_TRANSPOSE_t TransB,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_matrix_complex * B,  gsl_complex beta, gsl_matrix_complex * C)
+
+  int  gsl_blas_zsymm (CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_matrix_complex * B,  gsl_complex beta, gsl_matrix_complex * C)
+
+  int  gsl_blas_zsyrk (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_complex beta, gsl_matrix_complex * C)
+
+  int  gsl_blas_zsyr2k (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_matrix_complex * B,  gsl_complex beta, gsl_matrix_complex *C)
+
+  int  gsl_blas_ztrmm (CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag,  gsl_complex alpha,  gsl_matrix_complex * A, gsl_matrix_complex * B)
+
+  int  gsl_blas_ztrsm (CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag,  gsl_complex alpha,  gsl_matrix_complex * A, gsl_matrix_complex * B)
+
+  int  gsl_blas_zhemm (CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_matrix_complex * B,  gsl_complex beta, gsl_matrix_complex * C)
+
+  int  gsl_blas_zherk (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, double alpha,  gsl_matrix_complex * A, double beta, gsl_matrix_complex * C)
+
+  int  gsl_blas_zher2k (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans,  gsl_complex alpha,  gsl_matrix_complex * A,  gsl_matrix_complex * B, double beta, gsl_matrix_complex * C)

sage/libs/gsl/blas_types.pxd

@@ -0,0 +1,34 @@
+# sage_setup: distribution = sagemath-modules
+# distutils: include_dirs = GSL_INCDIR
+cdef extern from "gsl/gsl_cblas.h":
+  cdef enum CBLAS_ORDER:
+    CblasRowMajor=101
+    CblasColMajor=102
+
+  cdef enum CBLAS_TRANSPOSE:
+    CblasNoTrans=111
+    CblasTrans=112
+    CblasConjTrans=113
+
+  cdef enum CBLAS_UPLO:
+    CblasUpper=121
+    CblasLower=122
+
+  cdef enum CBLAS_DIAG:
+    CblasNonUnit=131
+    CblasUnit=132
+
+  cdef enum CBLAS_SIDE:
+    CblasLeft=141
+    CblasRight=142
+
+ctypedef size_t CBLAS_INDEX
+
+cdef extern from "gsl/gsl_blas_types.h":
+  ctypedef  CBLAS_INDEX  CBLAS_INDEX_t
+  ctypedef  CBLAS_ORDER       CBLAS_ORDER_t
+  ctypedef  CBLAS_TRANSPOSE   CBLAS_TRANSPOSE_t
+  ctypedef  CBLAS_UPLO        CBLAS_UPLO_t
+  ctypedef  CBLAS_DIAG        CBLAS_DIAG_t
+  ctypedef  CBLAS_SIDE        CBLAS_SIDE_t
+

sage/libs/gsl/block.pxd

@@ -0,0 +1,52 @@
+# 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_block_double.h":
+
+  ctypedef struct gsl_block:
+    size_t size
+    double * data
+
+  gsl_block *  gsl_block_alloc(size_t n)
+
+  gsl_block *  gsl_block_calloc(size_t n)
+
+  void  gsl_block_free(gsl_block * b)
+
+  int  gsl_block_fread(FILE * stream, gsl_block * b)
+
+  int  gsl_block_fwrite(FILE * stream, gsl_block * b)
+
+  int  gsl_block_fscanf(FILE * stream, gsl_block * b)
+
+  int  gsl_block_fprintf(FILE * stream, gsl_block * b, char * format)
+
+  size_t gsl_block_size (gsl_block * b)
+  double * gsl_block_data (gsl_block * b)
+
+
+cdef extern from "gsl/gsl_block_complex_double.h":
+
+  ctypedef struct gsl_block_complex:
+    size_t size
+    double * data
+
+  gsl_block_complex *  gsl_block_complex_alloc(size_t n)
+
+  gsl_block_complex *  gsl_block_complex_calloc(size_t n)
+
+  void  gsl_block_complex_free(gsl_block_complex * b)
+
+  int  gsl_block_complex_fread(FILE * stream, gsl_block_complex * b)
+
+  int  gsl_block_complex_fwrite(FILE * stream, gsl_block_complex * b)
+
+  int  gsl_block_complex_fscanf(FILE * stream, gsl_block_complex * b)
+
+  int  gsl_block_complex_fprintf(FILE * stream, gsl_block_complex * b, char * format)
+
+  size_t gsl_block_complex_size (gsl_block_complex * b)
+  double * gsl_block_complex_data (gsl_block_complex * b)

sage/libs/gsl/chebyshev.pxd

@@ -0,0 +1,37 @@
+# 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_chebyshev.h":
+
+  ctypedef struct gsl_cheb_series
+
+
+  gsl_cheb_series * gsl_cheb_alloc( size_t order)
+
+  void gsl_cheb_free(gsl_cheb_series * cs)
+
+  int gsl_cheb_init(gsl_cheb_series * cs,  gsl_function * func,
+                     double a,  double b)
+
+
+  double gsl_cheb_eval( gsl_cheb_series * cs,  double x)
+  int gsl_cheb_eval_err( gsl_cheb_series * cs,  double x,
+                        double * result, double * abserr)
+
+
+  double gsl_cheb_eval_n( gsl_cheb_series * cs,  size_t order, double x)
+  int gsl_cheb_eval_n_err( gsl_cheb_series * cs,  size_t order,
+                           double x, double * result, double * abserr)
+
+
+  double gsl_cheb_eval_mode( gsl_cheb_series * cs, double x, gsl_mode_t mode)
+  int gsl_cheb_eval_mode_e( gsl_cheb_series * cs,  double x, gsl_mode_t mode, double * result, double * abserr)
+
+
+  int gsl_cheb_calc_deriv(gsl_cheb_series * deriv,  gsl_cheb_series * cs)
+
+  int gsl_cheb_calc_integ(gsl_cheb_series * integ,  gsl_cheb_series * cs)
+

sage/libs/gsl/clausen.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_clausen.h":
+
+  double  gsl_sf_clausen(double x)
+
+  int  gsl_sf_clausen_e(double x, gsl_sf_result * result)
+

sage/libs/gsl/combination.pxd

@@ -0,0 +1,47 @@
+# 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_combination.h":
+
+  ctypedef struct gsl_combination:
+    size_t n
+    size_t k
+    size_t *data
+
+  gsl_combination *  gsl_combination_alloc(size_t n, size_t k)
+
+  gsl_combination *  gsl_combination_calloc(size_t n, size_t k)
+
+  void  gsl_combination_init_first(gsl_combination * c)
+
+  void  gsl_combination_init_last(gsl_combination * c)
+
+  void  gsl_combination_free(gsl_combination * c)
+
+  int  gsl_combination_memcpy(gsl_combination * dest, gsl_combination * src)
+
+  int  gsl_combination_fread(FILE * stream, gsl_combination * c)
+
+  int  gsl_combination_fwrite(FILE * stream, gsl_combination * c)
+
+  int  gsl_combination_fscanf(FILE * stream, gsl_combination * c)
+
+  int  gsl_combination_fprintf(FILE * stream, gsl_combination * c, char *format)
+
+  size_t  gsl_combination_n(gsl_combination * c)
+
+  size_t  gsl_combination_k(gsl_combination * c)
+
+  size_t *  gsl_combination_data(gsl_combination * c)
+
+  size_t  gsl_combination_get(gsl_combination * c, size_t i)
+
+  int  gsl_combination_valid(gsl_combination * c)
+
+  int  gsl_combination_next(gsl_combination * c)
+
+  int  gsl_combination_prev(gsl_combination * c)
+

sage/libs/gsl/complex.pxd

@@ -0,0 +1,151 @@
+# 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_complex.h":
+  double GSL_REAL(gsl_complex z)
+
+  double GSL_IMAG(gsl_complex z)
+
+  int GSL_COMPLEX_EQ(gsl_complex z1,gsl_complex z2)
+
+  void GSL_SET_COMPLEX(gsl_complex * zp, double  x, double  y)
+
+  void GSL_SET_REAL(gsl_complex * zp, double x)
+
+  void GSL_SET_IMAG(gsl_complex * zp, double y)
+
+
+cdef extern from "gsl/gsl_complex_math.h":
+
+# Complex numbers
+  gsl_complex  gsl_complex_rect(double x, double y)
+
+  gsl_complex  gsl_complex_polar(double r, double theta)
+
+  gsl_complex GSL_COMPLEX_ONE
+  gsl_complex GSL_COMPLEX_ZERO
+  gsl_complex GSL_COMPLEX_NEGONE
+
+  # Properties of complex numbers
+  double  gsl_complex_arg(gsl_complex z)
+
+  double  gsl_complex_abs(gsl_complex z)
+
+  double  gsl_complex_abs2(gsl_complex z)
+
+  double  gsl_complex_logabs(gsl_complex z)
+
+  #Complex arithmetic operators
+  gsl_complex  gsl_complex_add(gsl_complex a, gsl_complex b)
+
+  gsl_complex  gsl_complex_sub(gsl_complex a, gsl_complex b)
+
+  gsl_complex  gsl_complex_mul(gsl_complex a, gsl_complex b)
+
+  gsl_complex  gsl_complex_div(gsl_complex a, gsl_complex b)
+
+  gsl_complex  gsl_complex_add_real(gsl_complex a, double x)
+
+  gsl_complex  gsl_complex_sub_real(gsl_complex a, double x)
+
+  gsl_complex  gsl_complex_mul_real(gsl_complex a, double x)
+
+  gsl_complex  gsl_complex_div_real(gsl_complex a, double x)
+
+  gsl_complex  gsl_complex_add_imag(gsl_complex a, double y)
+
+  gsl_complex  gsl_complex_sub_imag(gsl_complex a, double y)
+
+  gsl_complex  gsl_complex_mul_imag(gsl_complex a, double y)
+
+  gsl_complex  gsl_complex_div_imag(gsl_complex a, double y)
+
+  gsl_complex  gsl_complex_conjugate(gsl_complex z)
+
+  gsl_complex  gsl_complex_inverse(gsl_complex z)
+
+  gsl_complex  gsl_complex_negative(gsl_complex z)
+
+  #Elementary Complex Functions
+  gsl_complex  gsl_complex_sqrt(gsl_complex z)
+
+  gsl_complex  gsl_complex_sqrt_real(double x)
+
+  gsl_complex  gsl_complex_pow(gsl_complex z, gsl_complex a)
+
+  gsl_complex  gsl_complex_pow_real(gsl_complex z, double x)
+
+  gsl_complex  gsl_complex_exp(gsl_complex z)
+
+  gsl_complex  gsl_complex_log(gsl_complex z)
+
+  gsl_complex  gsl_complex_log10(gsl_complex z)
+
+  gsl_complex  gsl_complex_log_b(gsl_complex z, gsl_complex b)
+
+  #Complex Trigonometric Functions
+  gsl_complex  gsl_complex_sin(gsl_complex z)
+
+  gsl_complex  gsl_complex_cos(gsl_complex z)
+
+  gsl_complex  gsl_complex_tan(gsl_complex z)
+
+  gsl_complex  gsl_complex_sec(gsl_complex z)
+
+  gsl_complex  gsl_complex_csc(gsl_complex z)
+
+  gsl_complex  gsl_complex_cot(gsl_complex z)
+
+  #Inverse Complex Trigonometric Functions
+  gsl_complex  gsl_complex_arcsin(gsl_complex z)
+
+  gsl_complex  gsl_complex_arcsin_real(double z)
+
+  gsl_complex  gsl_complex_arccos(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccos_real(double z)
+
+  gsl_complex  gsl_complex_arctan(gsl_complex z)
+
+  gsl_complex  gsl_complex_arcsec(gsl_complex z)
+
+  gsl_complex  gsl_complex_arcsec_real(double z)
+
+  gsl_complex  gsl_complex_arccsc(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccsc_real(double z)
+
+  gsl_complex  gsl_complex_arccot(gsl_complex z)
+
+  #Complex Hyperbolic Functions
+  gsl_complex  gsl_complex_sinh(gsl_complex z)
+
+  gsl_complex  gsl_complex_cosh(gsl_complex z)
+
+  gsl_complex  gsl_complex_tanh(gsl_complex z)
+
+  gsl_complex  gsl_complex_sech(gsl_complex z)
+
+  gsl_complex  gsl_complex_csch(gsl_complex z)
+
+  gsl_complex  gsl_complex_coth(gsl_complex z)
+
+  #Inverse Complex Hyperbolic Functions
+  gsl_complex  gsl_complex_arcsinh(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccosh(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccosh_real(double z)
+
+  gsl_complex  gsl_complex_arctanh(gsl_complex z)
+
+  gsl_complex  gsl_complex_arctanh_real(double z)
+
+  gsl_complex  gsl_complex_arcsech(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccsch(gsl_complex z)
+
+  gsl_complex  gsl_complex_arccoth(gsl_complex z)

sage/libs/gsl/coulomb.pxd

@@ -0,0 +1,30 @@
+# 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_coulomb.h":
+
+  double  gsl_sf_hydrogenicR_1(double Z, double r)
+
+  int  gsl_sf_hydrogenicR_1_e(double Z, double r, gsl_sf_result * result)
+
+  double  gsl_sf_hydrogenicR(int n, int l, double Z, double r)
+
+  int  gsl_sf_hydrogenicR_e(int n, int l, double Z, double r, gsl_sf_result * result)
+
+  int  gsl_sf_coulomb_wave_FG_e(double eta, double x, double L_F, int k, gsl_sf_result * F, gsl_sf_result * Fp, gsl_sf_result * G, gsl_sf_result * Gp, double * exp_F, double * exp_G)
+
+  int  gsl_sf_coulomb_wave_F_array(double L_min, int kmax, double eta, double x, double fc_array[], double * F_exponent)
+
+  int  gsl_sf_coulomb_wave_FG_array(double L_min, int kmax, double eta, double x, double fc_array[], double gc_array[], double * F_exponent, double * G_exponent)
+
+  int  gsl_sf_coulomb_wave_FGp_array(double L_min, int kmax, double eta, double x, double fc_array[], double fcp_array[], double gc_array[], double gcp_array[], double * F_exponent, double * G_exponent)
+
+  int  gsl_sf_coulomb_wave_sphF_array(double L_min, int kmax, double eta, double x, double fc_array[], double F_exponent[])
+
+  int  gsl_sf_coulomb_CL_e(double L, double eta, gsl_sf_result * result)
+
+  int  gsl_sf_coulomb_CL_array(double Lmin, int kmax, double eta, double cl[])
+