passagemath-modules 10.6.31rc3__cp314-cp314-musllinux_1_2_x86_64.whl

sage/libs/gsl/linalg.pxd

@@ -0,0 +1,185 @@
+# 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_linalg.h":
+
+  cdef enum gsl_linalg_matrix_mod_t:
+    GSL_LINALG_MOD_NONE = 0
+    GSL_LINALG_MOD_TRANSPOSE = 1
+    GSL_LINALG_MOD_CONJUGATE = 2
+
+  int gsl_linalg_matmult ( gsl_matrix * A,  gsl_matrix * B, gsl_matrix * C)
+
+  int gsl_linalg_matmult_mod ( gsl_matrix * A, gsl_linalg_matrix_mod_t modA,  gsl_matrix * B, gsl_linalg_matrix_mod_t modB, gsl_matrix * C)
+
+  int gsl_linalg_exponential_ss(gsl_matrix * A, gsl_matrix * eA, gsl_mode_t mode)
+
+  # Householder Transformations
+  double gsl_linalg_householder_transform (gsl_vector * v)
+
+  gsl_complex gsl_linalg_complex_householder_transform (gsl_vector_complex * v)
+
+  int gsl_linalg_householder_hm (double tau,  gsl_vector * v, gsl_matrix * A)
+
+  int gsl_linalg_householder_mh (double tau,  gsl_vector * v, gsl_matrix * A)
+
+  int gsl_linalg_householder_hv (double tau,  gsl_vector * v, gsl_vector * w)
+
+  int gsl_linalg_householder_hm1 (double tau, gsl_matrix * A)
+
+  int gsl_linalg_complex_householder_hm (gsl_complex tau,  gsl_vector_complex * v, gsl_matrix_complex * A)
+
+  int gsl_linalg_complex_householder_hv (gsl_complex tau,  gsl_vector_complex * v, gsl_vector_complex * w)
+
+  # Singular Value Decomposition
+  int  gsl_linalg_SV_decomp(gsl_matrix * A, gsl_matrix * V, gsl_vector * S, gsl_vector * work)
+
+  int  gsl_linalg_SV_decomp_mod(gsl_matrix * A, gsl_matrix * X, gsl_matrix * V, gsl_vector * S, gsl_vector * work)
+
+  int  gsl_linalg_SV_decomp_jacobi(gsl_matrix * A, gsl_matrix * V, gsl_vector * S)
+
+  int  gsl_linalg_SV_solve(gsl_matrix * U, gsl_matrix * V, gsl_vector * S, gsl_vector * b, gsl_vector * x)
+
+
+  # LU Decomposition, Gaussian elimination with partial pivoting
+  int  gsl_linalg_LU_decomp(gsl_matrix * A, gsl_permutation * p, int *signum)
+
+  int  gsl_linalg_LU_solve(gsl_matrix * LU, gsl_permutation * p, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_LU_svx(gsl_matrix * LU, gsl_permutation * p, gsl_vector * x)
+
+  int  gsl_linalg_LU_refine(gsl_matrix * A, gsl_matrix * LU, gsl_permutation * p, gsl_vector * b, gsl_vector * x, gsl_vector * residual)
+
+  int  gsl_linalg_LU_invert(gsl_matrix * LU, gsl_permutation * p, gsl_matrix * inverse)
+
+  double  gsl_linalg_LU_det(gsl_matrix * LU, int signum)
+
+  double  gsl_linalg_LU_lndet(gsl_matrix * LU)
+
+  int  gsl_linalg_LU_sgndet(gsl_matrix * LU, int signum)
+
+
+  # Complex LU Decomposition
+
+  int  gsl_linalg_complex_LU_decomp(gsl_matrix_complex * A, gsl_permutation * p, int *signum)
+
+  int  gsl_linalg_complex_LU_solve(gsl_matrix_complex * LU, gsl_permutation * p, gsl_vector_complex * b, gsl_vector_complex * x)
+
+  int  gsl_linalg_complex_LU_svx(gsl_matrix_complex * LU, gsl_permutation * p, gsl_vector_complex * x)
+
+  int  gsl_linalg_complex_LU_refine(gsl_matrix_complex * A, gsl_matrix_complex * LU, gsl_permutation * p, gsl_vector_complex * b, gsl_vector_complex * x, gsl_vector_complex * residual)
+
+#  int  gsl_complex_linalg_LU_invert(gsl_matrix_complex * LU, gsl_permutation * p, gsl_matrix_complex * inverse)
+  int  gsl_linalg_complex_LU_invert(gsl_matrix_complex * LU, gsl_permutation * p, gsl_matrix_complex * inverse)
+
+  gsl_complex  gsl_linalg_complex_LU_det(gsl_matrix_complex * LU, int signum)
+
+  double  gsl_linalg_complex_LU_lndet(gsl_matrix_complex * LU)
+
+  gsl_complex  gsl_linalg_complex_LU_sgndet(gsl_matrix_complex * LU, int signum)
+
+
+  # QR decomposition
+
+  int  gsl_linalg_QR_decomp(gsl_matrix * A, gsl_vector * tau)
+
+  int  gsl_linalg_QR_solve(gsl_matrix * QR, gsl_vector * tau, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QR_svx(gsl_matrix * QR, gsl_vector * tau, gsl_vector * x)
+
+  int  gsl_linalg_QR_lssolve(gsl_matrix * QR, gsl_vector * tau, gsl_vector * b, gsl_vector * x, gsl_vector * residual)
+
+  int  gsl_linalg_QR_QRsolve(gsl_matrix * Q, gsl_matrix * R, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QR_Rsolve(gsl_matrix * QR, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QR_Rsvx(gsl_matrix * QR, gsl_vector * x)
+
+  int  gsl_linalg_QR_update(gsl_matrix * Q, gsl_matrix * R, gsl_vector * w, gsl_vector * v)
+
+  int  gsl_linalg_QR_QTvec(gsl_matrix * QR, gsl_vector * tau, gsl_vector * v)
+
+  int  gsl_linalg_QR_Qvec(gsl_matrix * QR, gsl_vector * tau, gsl_vector * v)
+
+  int  gsl_linalg_QR_unpack(gsl_matrix * QR, gsl_vector * tau, gsl_matrix * Q, gsl_matrix * R)
+
+  int  gsl_linalg_R_solve(gsl_matrix * R, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_R_svx(gsl_matrix * R, gsl_vector * x)
+
+
+  # Q R P^T decomposition
+
+  int  gsl_linalg_QRPT_decomp(gsl_matrix * A, gsl_vector * tau, gsl_permutation * p, int *signum, gsl_vector * norm)
+
+  int  gsl_linalg_QRPT_decomp2(gsl_matrix * A, gsl_matrix * q, gsl_matrix * r, gsl_vector * tau, gsl_permutation * p, int *signum, gsl_vector * norm)
+
+  int  gsl_linalg_QRPT_solve(gsl_matrix * QR, gsl_vector * tau, gsl_permutation * p, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QRPT_svx(gsl_matrix * QR, gsl_vector * tau, gsl_permutation * p, gsl_vector * x)
+
+  int  gsl_linalg_QRPT_QRsolve(gsl_matrix * Q, gsl_matrix * R, gsl_permutation * p, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QRPT_Rsolve(gsl_matrix * QR, gsl_permutation * p, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_QRPT_Rsvx(gsl_matrix * QR, gsl_permutation * p, gsl_vector * x)
+
+  int  gsl_linalg_QRPT_update(gsl_matrix * Q, gsl_matrix * R, gsl_permutation * p, gsl_vector * u, gsl_vector * v)
+
+
+  # Cholesky Decomposition
+  int  gsl_linalg_cholesky_decomp(gsl_matrix * A)
+
+  int  gsl_linalg_cholesky_solve(gsl_matrix * cholesky, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_cholesky_svx(gsl_matrix * cholesky, gsl_vector * x)
+
+  # Symmetric to symmetric tridiagonal decomposition
+  int  gsl_linalg_symmtd_decomp(gsl_matrix * A, gsl_vector * tau)
+
+  int  gsl_linalg_symmtd_unpack(gsl_matrix * A, gsl_vector * tau, gsl_matrix * Q, gsl_vector * diag, gsl_vector * subdiag)
+
+  int  gsl_linalg_symmtd_unpack_T(gsl_matrix * A, gsl_vector * diag, gsl_vector * subdiag)
+
+
+  # Hermitian to symmetric tridiagonal decomposition
+  int  gsl_linalg_hermtd_decomp(gsl_matrix_complex * A, gsl_vector_complex * tau)
+
+  int  gsl_linalg_hermtd_unpack(gsl_matrix_complex * A, gsl_vector_complex * tau, gsl_matrix_complex * Q, gsl_vector * diag, gsl_vector * subdiag)
+
+  int  gsl_linalg_hermtd_unpack_T(gsl_matrix_complex * A, gsl_vector * diag, gsl_vector * subdiag)
+
+  # Linear Solve Using Householder Transformations
+  int  gsl_linalg_HH_solve(gsl_matrix * A, gsl_vector * b, gsl_vector * x)
+
+  int  gsl_linalg_HH_svx(gsl_matrix * A, gsl_vector * x)
+
+  # Linear solve for a symmetric tridiagonal system
+  int  gsl_linalg_solve_symm_tridiag(gsl_vector * diag, gsl_vector * e, gsl_vector * b, gsl_vector * x)
+
+  # Linear solve for a nonsymmetric tridiagonal system
+  int gsl_linalg_solve_tridiag ( gsl_vector * diag,  gsl_vector * abovediag,  gsl_vector * belowdiag,  gsl_vector * b, gsl_vector * x)
+
+  # Linear solve for a symmetric cyclic tridiagonal system
+  int  gsl_linalg_solve_symm_cyc_tridiag(gsl_vector * diag, gsl_vector * e, gsl_vector * b, gsl_vector * x)
+
+
+  # Linear solve for a nonsymmetric cyclic tridiagonal system
+  int gsl_linalg_solve_cyc_tridiag ( gsl_vector * diag,  gsl_vector * abovediag,  gsl_vector * belowdiag,  gsl_vector * b, gsl_vector * x)
+
+  # Bidiagonal decomposition
+  int  gsl_linalg_bidiag_decomp(gsl_matrix * A, gsl_vector * tau_U, gsl_vector * tau_V)
+
+  int  gsl_linalg_bidiag_unpack(gsl_matrix * A, gsl_vector * tau_U, gsl_matrix * U, gsl_vector * tau_V, gsl_matrix * V, gsl_vector * diag, gsl_vector * superdiag)
+
+  int  gsl_linalg_bidiag_unpack2(gsl_matrix * A, gsl_vector * tau_U, gsl_vector * tau_V, gsl_matrix * V)
+
+  int  gsl_linalg_bidiag_unpack_B(gsl_matrix * A, gsl_vector * diag, gsl_vector * superdiag)
+
+
+  # Balancing
+  int gsl_linalg_balance_columns (gsl_matrix * A, gsl_vector * D)
+

sage/libs/gsl/log.pxd

@@ -0,0 +1,26 @@
+# 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_log.h":
+
+  double  gsl_sf_log(double x)
+
+  int  gsl_sf_log_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_log_abs(double x)
+
+  int  gsl_sf_log_abs_e(double x, gsl_sf_result * result)
+
+  int  gsl_sf_complex_log_e(double zr, double zi, gsl_sf_result * lnr, gsl_sf_result * theta)
+
+  double  gsl_sf_log_1plusx(double x)
+
+  int  gsl_sf_log_1plusx_e(double x, gsl_sf_result * result)
+
+  double  gsl_sf_log_1plusx_mx(double x)
+
+  int  gsl_sf_log_1plusx_mx_e(double x, gsl_sf_result * result)
+

sage/libs/gsl/math.pxd

@@ -0,0 +1,43 @@
+# 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_math.h":
+    bint gsl_isnan(double x)
+    bint gsl_isinf(double x)
+    bint gsl_finite(double x)
+    double gsl_log1p(double x)
+    double gsl_expm1(double x)
+    double gsl_hypot(double x, double y)
+    double gsl_acosh(double x)
+    double gsl_asinh(double x)
+    double gsl_atanh(double x)
+    double gsl_ldexp(double x, int e)
+    double gsl_frexp(double x, int * e)
+    double gsl_pow_int(double x, int n)
+    double gsl_pow_2(double x)
+    double gsl_pow_3(double x)
+    double gsl_pow_4(double x)
+    double gsl_pow_5(double x)
+    double gsl_pow_6(double x)
+    double gsl_pow_7(double x)
+    double gsl_pow_8(double x)
+    double gsl_pow_9(double x)
+    int GSL_SIGN(double x)
+    bint GSL_IS_ODD(int n)
+    bint GSL_IS_EVEN(int n)
+    double GSL_MAX(double a, double b)
+    double GSL_MIN(double a, double b)
+    double GSL_MAX_DBL(double a, double b)
+    double GSL_MIN_DBL(double a, double b)
+    int GSL_MAX_INT(int a, int b)
+    int GSL_MIN_INT(int a, int b)
+    long double GSL_MAX_LDBL(long double a, long double b)
+    long double GSL_MIN_LDBL(long double a, long double b)
+    int gsl_fcmp(double x, double y, double epsilon)
+    double GSL_FN_EVAL(gsl_function * F, double x)
+    double GSL_FN_FDF_EVAL_F(gsl_function_fdf * FDF, double x)
+    GSL_FN_FDF_EVAL_DF(gsl_function_fdf * FDF, double x)
+    GSL_FN_FDF_EVAL_F_DF(gsl_function_fdf * FDF, double x, double y, double dy)

sage/libs/gsl/matrix.pxd

@@ -0,0 +1,143 @@
+# 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_matrix_double.h":
+  # Allocation
+  gsl_matrix *  gsl_matrix_alloc(size_t n1, size_t n2)
+
+  gsl_matrix *  gsl_matrix_calloc(size_t n1, size_t n2)
+
+  gsl_matrix *  gsl_matrix_alloc_from_block(gsl_block * b,
+    size_t offset, size_t n1, size_t n2, size_t d2)
+
+  gsl_matrix * gsl_matrix_alloc_from_matrix (gsl_matrix * m,  size_t k1,  size_t k2,  size_t n1,  size_t n2)
+
+  gsl_vector * gsl_vector_alloc_row_from_matrix (gsl_matrix * m,  size_t i)
+
+  gsl_vector * gsl_vector_alloc_col_from_matrix (gsl_matrix * m,  size_t j)
+
+  void  gsl_matrix_free(gsl_matrix * m)
+
+  # Views
+  gsl_matrix_view  gsl_matrix_submatrix(gsl_matrix * m, size_t k1, size_t k2, size_t n1, size_t n2)
+
+  gsl_vector_view  gsl_matrix_row(gsl_matrix * m, size_t i)
+
+  gsl_vector_view  gsl_matrix_column(gsl_matrix * m, size_t j)
+
+  gsl_vector_view  gsl_matrix_diagonal(gsl_matrix * m)
+
+  gsl_vector_view  gsl_matrix_subdiagonal(gsl_matrix * m, size_t k)
+
+  gsl_vector_view  gsl_matrix_superdiagonal(gsl_matrix * m, size_t k)
+
+  gsl_matrix_view  gsl_matrix_view_array(double * base, size_t n1, size_t n2)
+
+  gsl_matrix_view  gsl_matrix_view_array_with_tda(double * base, size_t n1, size_t n2, size_t tda)
+
+  gsl_matrix_view  gsl_matrix_view_vector(gsl_vector * v, size_t n1, size_t n2)
+
+  gsl_matrix_view  gsl_matrix_view_vector_with_tda(gsl_vector * v, size_t n1, size_t n2, size_t tda)
+
+  gsl_matrix_const_view  gsl_matrix_const_submatrix(gsl_matrix * m, size_t k1, size_t k2, size_t n1, size_t n2)
+
+  gsl_vector_const_view  gsl_matrix_const_row(gsl_matrix * m, size_t i)
+
+  gsl_vector_const_view  gsl_matrix_const_column(gsl_matrix * m, size_t j)
+
+  gsl_vector_const_view  gsl_matrix_const_diagonal(gsl_matrix * m)
+
+  gsl_vector_const_view  gsl_matrix_const_subdiagonal(gsl_matrix * m, size_t k)
+
+  gsl_vector_const_view  gsl_matrix_const_superdiagonal(gsl_matrix * m, size_t k)
+
+  gsl_matrix_const_view  gsl_matrix_const_view_array(double * base, size_t n1, size_t n2)
+
+  gsl_matrix_const_view  gsl_matrix_const_view_array_with_tda(double * base, size_t n1, size_t n2, size_t tda)
+
+  gsl_matrix_const_view  gsl_matrix_const_view_vector(gsl_vector * v, size_t n1, size_t n2)
+
+  gsl_matrix_const_view  gsl_matrix_const_view_vector_with_tda(gsl_vector * v, size_t n1, size_t n2, size_t tda)
+
+
+  # Operations
+  double  gsl_matrix_get(gsl_matrix * m, size_t i, size_t j)
+
+  void  gsl_matrix_set(gsl_matrix * m, size_t i, size_t j, double x)
+
+  double *  gsl_matrix_ptr(gsl_matrix * m, size_t i, size_t j)
+
+  double *  gsl_matrix_const_ptr(gsl_matrix * m, size_t i, size_t j)
+
+  void  gsl_matrix_set_zero(gsl_matrix * m)
+
+  void  gsl_matrix_set_identity(gsl_matrix * m)
+
+  void  gsl_matrix_set_all(gsl_matrix * m, double x)
+
+  # Reading and writing matrices
+  int  gsl_matrix_fread(FILE * stream, gsl_matrix * m)
+
+  int  gsl_matrix_fwrite(FILE * stream, gsl_matrix * m)
+
+  int  gsl_matrix_fscanf(FILE * stream, gsl_matrix * m)
+
+  int  gsl_matrix_fprintf(FILE * stream, gsl_matrix * m, char * format)
+
+  # Copying or exchanging elements
+  int  gsl_matrix_memcpy(gsl_matrix * dest, gsl_matrix * src)
+
+  int  gsl_matrix_swap(gsl_matrix * m1, gsl_matrix * m2)
+
+  int  gsl_matrix_swap_rows(gsl_matrix * m, size_t i, size_t j)
+
+  int  gsl_matrix_swap_columns(gsl_matrix * m, size_t i, size_t j)
+
+  int  gsl_matrix_swap_rowcol(gsl_matrix * m, size_t i, size_t j)
+
+  int  gsl_matrix_transpose(gsl_matrix * m)
+
+  int  gsl_matrix_transpose_memcpy(gsl_matrix * dest, gsl_matrix * src)
+
+  # Finding maximum and minimum elements of matrices
+  double  gsl_matrix_max(gsl_matrix * m)
+
+  double  gsl_matrix_min(gsl_matrix * m)
+
+  void  gsl_matrix_minmax(gsl_matrix * m, double * min_out, double * max_out)
+
+  void  gsl_matrix_max_index(gsl_matrix * m, size_t * imax, size_t * jmax)
+
+  void  gsl_matrix_min_index(gsl_matrix * m, size_t * imax, size_t * jmax)
+
+  void  gsl_matrix_minmax_index(gsl_matrix * m, size_t * imin, size_t * jmin, size_t * imax, size_t * jmax)
+
+  int  gsl_matrix_isnull(gsl_matrix * m)
+
+  # Matrix operations
+  int  gsl_matrix_add(gsl_matrix * a, gsl_matrix * b)
+
+  int  gsl_matrix_sub(gsl_matrix * a, gsl_matrix * b)
+
+  int  gsl_matrix_mul_elements(gsl_matrix * a, gsl_matrix * b)
+
+  int  gsl_matrix_div_elements(gsl_matrix * a, gsl_matrix * b)
+
+  int  gsl_matrix_scale(gsl_matrix * a, double x)
+
+  int  gsl_matrix_add_constant(gsl_matrix * a, double x)
+
+  int gsl_matrix_add_diagonal (gsl_matrix * a,  double x)
+
+  # The functions below are obsolete
+  int  gsl_matrix_get_row(gsl_vector * v, gsl_matrix * m, size_t i)
+
+  int  gsl_matrix_get_col(gsl_vector * v, gsl_matrix * m, size_t j)
+
+  int  gsl_matrix_set_row(gsl_matrix * m, size_t i, gsl_vector * v)
+
+  int  gsl_matrix_set_col(gsl_matrix * m, size_t j, gsl_vector * v)
+

sage/libs/gsl/matrix_complex.pxd

@@ -0,0 +1,130 @@
+# 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_matrix_complex_double.h":
+  # Allocation
+  gsl_matrix_complex *  gsl_matrix_complex_alloc(size_t n1, size_t n2)
+
+  gsl_matrix_complex *  gsl_matrix_complex_calloc(size_t n1, size_t n2)
+
+  gsl_matrix_complex * gsl_matrix_complex_alloc_from_block (gsl_block_complex * b,  size_t offset,  size_t n1,  size_t n2,  size_t d2)
+
+  gsl_matrix_complex * gsl_matrix_complex_alloc_from_matrix (gsl_matrix_complex * b,  size_t k1,  size_t k2,  size_t n1,  size_t n2)
+
+  gsl_vector_complex * gsl_vector_complex_alloc_row_from_matrix (gsl_matrix_complex * m,  size_t i)
+
+  gsl_vector_complex * gsl_vector_complex_alloc_col_from_matrix (gsl_matrix_complex * m,  size_t j)
+
+  void  gsl_matrix_complex_free(gsl_matrix_complex * m)
+
+  # Views
+  gsl_matrix_complex_view  gsl_matrix_complex_submatrix(gsl_matrix_complex * m, size_t k1, size_t k2, size_t n1, size_t n2)
+
+  gsl_vector_complex_view  gsl_matrix_complex_row(gsl_matrix_complex * m, size_t i)
+
+  gsl_vector_complex_view  gsl_matrix_complex_column(gsl_matrix_complex * m, size_t j)
+
+  gsl_vector_complex_view  gsl_matrix_complex_diagonal(gsl_matrix_complex * m)
+
+  gsl_vector_complex_view  gsl_matrix_complex_subdiagonal(gsl_matrix_complex * m, size_t k)
+
+  gsl_vector_complex_view  gsl_matrix_complex_superdiagonal(gsl_matrix_complex * m, size_t k)
+
+  gsl_matrix_complex_view gsl_matrix_complex_view_array (double * base,  size_t n1,  size_t n2)
+
+  gsl_matrix_complex_view gsl_matrix_complex_view_array_with_tda (double * base,  size_t n1,  size_t n2,  size_t tda)
+
+  gsl_matrix_complex_view  gsl_matrix_complex_view_vector(gsl_vector_complex * v, size_t n1, size_t n2)
+
+  gsl_matrix_complex_view  gsl_matrix_complex_view_vector_with_tda(gsl_vector_complex* v, size_t n1, size_t n2, size_t tda)
+
+  gsl_matrix_complex_const_view  gsl_matrix_complex_const_submatrix(gsl_matrix_complex * m, size_t k1, size_t k2, size_t n1, size_t n2)
+
+  gsl_vector_complex_const_view  gsl_matrix_complex_const_row(gsl_matrix_complex * m, size_t i)
+
+  gsl_vector_complex_const_view  gsl_matrix_complex_const_column(gsl_matrix_complex * m, size_t j)
+
+  gsl_vector_complex_const_view  gsl_matrix_complex_const_diagonal(gsl_matrix_complex * m)
+
+  gsl_vector_complex_const_view  gsl_matrix_complex_const_subdiagonal(gsl_matrix_complex * m, size_t k)
+
+  gsl_vector_complex_const_view  gsl_matrix_complex_const_superdiagonal(gsl_matrix_complex * m, size_t k)
+
+  gsl_matrix_complex_const_view gsl_matrix_complex_const_view_array ( double * base,  size_t n1,  size_t n2)
+
+  gsl_matrix_complex_const_view gsl_matrix_complex_const_view_array_with_tda ( double * base,  size_t n1,  size_t n2,  size_t tda)
+
+  gsl_matrix_complex_const_view  gsl_matrix_complex_const_view_vector(gsl_vector_complex * v, size_t n1, size_t n2)
+
+  gsl_matrix_complex_const_view  gsl_matrix_complex_const_view_vector_with_tda(gsl_vector_complex * v, size_t n1, size_t n2, size_t tda)
+
+  # Operations
+  gsl_complex  gsl_matrix_complex_get(gsl_matrix_complex * m, size_t i, size_t j)
+
+  void  gsl_matrix_complex_set(gsl_matrix_complex * m, size_t i, size_t j, gsl_complex x)
+
+  gsl_complex *  gsl_matrix_complex_ptr(gsl_matrix_complex * m, size_t i, size_t j)
+
+  gsl_complex *  gsl_matrix_complex_const_ptr(gsl_matrix_complex * m, size_t i, size_t j)
+
+  void  gsl_matrix_complex_set_zero(gsl_matrix_complex * m)
+
+  void  gsl_matrix_complex_set_identity(gsl_matrix_complex * m)
+
+  void  gsl_matrix_complex_set_all(gsl_matrix_complex * m, gsl_complex x)
+
+  # Reading and writing matrices
+  int  gsl_matrix_complex_fread(FILE * stream, gsl_matrix_complex * m)
+
+  int  gsl_matrix_complex_fwrite(FILE * stream, gsl_matrix_complex * m)
+
+  int  gsl_matrix_complex_fscanf(FILE * stream, gsl_matrix_complex * m)
+
+  int  gsl_matrix_complex_fprintf(FILE * stream, gsl_matrix_complex * m, char * format)
+
+
+  # Copying or exchanging elements
+  int  gsl_matrix_complex_memcpy(gsl_matrix_complex * dest, gsl_matrix_complex * src)
+
+  int  gsl_matrix_complex_swap(gsl_matrix_complex * m1, gsl_matrix_complex * m2)
+
+  int  gsl_matrix_complex_swap_rows(gsl_matrix_complex * m, size_t i, size_t j)
+
+  int  gsl_matrix_complex_swap_columns(gsl_matrix_complex * m, size_t i, size_t j)
+
+  int  gsl_matrix_complex_swap_rowcol(gsl_matrix_complex * m, size_t i, size_t j)
+
+  int  gsl_matrix_complex_transpose(gsl_matrix_complex * m)
+
+  int  gsl_matrix_complex_transpose_memcpy(gsl_matrix_complex * dest, gsl_matrix_complex * src)
+
+  int  gsl_matrix_complex_isnull(gsl_matrix_complex * m)
+
+  # Matrix operations
+  int  gsl_matrix_complex_add(gsl_matrix_complex * a, gsl_matrix_complex * b)
+
+  int  gsl_matrix_complex_sub(gsl_matrix_complex * a, gsl_matrix_complex * b)
+
+  int  gsl_matrix_complex_mul_elements(gsl_matrix_complex * a, gsl_matrix_complex * b)
+
+  int  gsl_matrix_complex_div_elements(gsl_matrix_complex * a, gsl_matrix_complex * b)
+
+  int  gsl_matrix_complex_scale(gsl_matrix_complex * a, gsl_complex x)
+
+  int  gsl_matrix_complex_add_constant(gsl_matrix_complex * a, gsl_complex x)
+
+  int gsl_matrix_complex_add_diagonal (gsl_matrix_complex * a,  gsl_complex x)
+
+  # The functions below are obsolete
+  int  gsl_matrix_complex_get_row(gsl_vector * v, gsl_matrix_complex * m, size_t i)
+
+  int  gsl_matrix_complex_get_col(gsl_vector * v, gsl_matrix_complex * m, size_t j)
+
+  int  gsl_matrix_complex_set_row(gsl_matrix_complex * m, size_t i, gsl_vector * v)
+
+  int  gsl_matrix_complex_set_col(gsl_matrix_complex * m, size_t j, gsl_vector * v)
+
+

sage/libs/gsl/min.pxd

@@ -0,0 +1,67 @@
+# 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_min.h":
+
+  ctypedef struct gsl_min_fminimizer_type:
+      char *name
+      size_t size
+      int (*set) (void *state, gsl_function * f, double x_minimum, double f_minimum, double x_lower, double f_lower, double x_upper, double f_upper)
+      int (*iterate) (void *state, gsl_function * f, double * x_minimum, double * f_minimum, double * x_lower, double * f_lower, double * x_upper, double * f_upper)
+
+  ctypedef struct gsl_min_fminimizer:
+      gsl_min_fminimizer_type * type
+      gsl_function * function
+      double x_minimum
+      double x_lower
+      double x_upper
+      double f_minimum, f_lower, f_upper
+      void *state
+
+  gsl_min_fminimizer * gsl_min_fminimizer_alloc ( gsl_min_fminimizer_type * T)
+
+  void gsl_min_fminimizer_free (gsl_min_fminimizer * s)
+
+  int gsl_min_fminimizer_set (gsl_min_fminimizer * s,
+                              gsl_function * f, double x_minimum,
+                              double x_lower, double x_upper)
+
+  int gsl_min_fminimizer_set_with_values (gsl_min_fminimizer * s,
+                                          gsl_function * f,
+                                          double x_minimum, double f_minimum,
+                                          double x_lower, double f_lower,
+                                          double x_upper, double f_upper)
+
+  int gsl_min_fminimizer_iterate (gsl_min_fminimizer * s)
+
+  char * gsl_min_fminimizer_name ( gsl_min_fminimizer * s)
+
+  double gsl_min_fminimizer_x_minimum ( gsl_min_fminimizer * s)
+  double gsl_min_fminimizer_x_lower ( gsl_min_fminimizer * s)
+  double gsl_min_fminimizer_x_upper ( gsl_min_fminimizer * s)
+  double gsl_min_fminimizer_f_minimum ( gsl_min_fminimizer * s)
+  double gsl_min_fminimizer_f_lower ( gsl_min_fminimizer * s)
+  double gsl_min_fminimizer_f_upper ( gsl_min_fminimizer * s)
+
+  #Deprecated, use x_minimum instead
+  #double gsl_min_fminimizer_minimum ( gsl_min_fminimizer * s)
+
+  int gsl_min_test_interval (double x_lower, double x_upper, double epsabs, double epsrel)
+
+  gsl_min_fminimizer_type  * gsl_min_fminimizer_goldensection
+  gsl_min_fminimizer_type  * gsl_min_fminimizer_brent
+
+  ctypedef int (*gsl_min_bracketing_function)(gsl_function *f,
+                                     double *x_minimum,double * f_minimum,
+                                     double *x_lower, double * f_lower,
+                                     double *x_upper, double * f_upper,
+                                     size_t eval_max)
+
+  int gsl_min_find_bracket(gsl_function *f,double *x_minimum,double * f_minimum,
+                       double *x_lower, double * f_lower,
+                       double *x_upper, double * f_upper,
+                       size_t eval_max)
+

sage/libs/gsl/monte.pxd

@@ -0,0 +1,56 @@
+# 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_monte.h":
+  ctypedef struct gsl_monte_function:
+    double (*f)(double * x_array, size_t dim, void * params)
+    size_t dim
+    void * params
+
+cdef extern from "gsl/gsl_monte_plain.h":
+  ctypedef struct gsl_monte_plain_state
+
+  int gsl_monte_plain_integrate ( gsl_monte_function * f, \
+                              double xl[],  double xu[], \
+                              size_t dim, size_t calls, gsl_rng * r,  \
+                             gsl_monte_plain_state * state, \
+                             double *result, double *abserr)
+
+  gsl_monte_plain_state* gsl_monte_plain_alloc(size_t dim)
+
+  int gsl_monte_plain_init(gsl_monte_plain_state* state)
+
+  void gsl_monte_plain_free (gsl_monte_plain_state* state)
+
+cdef extern from "gsl/gsl_monte_miser.h":
+
+  ctypedef struct gsl_monte_miser_state
+
+  int gsl_monte_miser_integrate(gsl_monte_function * f, double xl[],  double xh[],  size_t dim, size_t calls, gsl_rng *r,  gsl_monte_miser_state* state, double *result, double *abserr)
+
+  gsl_monte_miser_state* gsl_monte_miser_alloc(size_t dim)
+
+  int gsl_monte_miser_init(gsl_monte_miser_state* state)
+
+  void gsl_monte_miser_free(gsl_monte_miser_state* state)
+
+  # vegas
+cdef extern from "gsl/gsl_monte_vegas.h":
+  cdef enum:
+    GSL_VEGAS_MODE_IMPORTANCE = 1
+    GSL_VEGAS_MODE_IMPORTANCE_ONLY = 0
+    GSL_VEGAS_MODE_STRATIFIED = -1
+
+  ctypedef struct gsl_monte_vegas_state:
+    double chisq
+
+  int gsl_monte_vegas_integrate(gsl_monte_function * f, double xl[], double xu[], size_t dim, size_t calls, gsl_rng * r, gsl_monte_vegas_state *state, double* result, double* abserr)
+
+  gsl_monte_vegas_state* gsl_monte_vegas_alloc(size_t dim)
+
+  int gsl_monte_vegas_init(gsl_monte_vegas_state* state)
+
+  void gsl_monte_vegas_free (gsl_monte_vegas_state* state)

sage/libs/gsl/ntuple.pxd

@@ -0,0 +1,32 @@
+# 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_ntuple.h":
+  ctypedef struct gsl_ntuple
+
+  ctypedef struct gsl_ntuple_select_fn:
+    int (* function) (void * ntuple_data, void * params)
+    void * params
+
+  ctypedef struct gsl_ntuple_value_fn:
+    double (* function) (void * ntuple_data, void * params)
+    void * params
+
+  gsl_ntuple * gsl_ntuple_open (char * filename, void * ntuple_data, size_t size)
+
+  gsl_ntuple * gsl_ntuple_create (char * filename, void * ntuple_data, size_t size)
+
+  int gsl_ntuple_write (gsl_ntuple * ntuple)
+  int gsl_ntuple_read (gsl_ntuple * ntuple)
+
+  int gsl_ntuple_bookdata (gsl_ntuple * ntuple)
+
+  int gsl_ntuple_project (gsl_histogram * h, gsl_ntuple * ntuple,
+                          gsl_ntuple_value_fn *value_func,
+                          gsl_ntuple_select_fn *select_func)
+
+  int gsl_ntuple_close (gsl_ntuple * ntuple)
+