gsl 1.12.108

data/include/rb_gsl_interp.h

@@ -0,0 +1,46 @@
+/*
+  rb_gsl_interp.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or INTERPOLATIONNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_INTERP_H___
+#define ___RB_GSL_INTERP_H___
+
+#include "ruby.h"
+#include "rb_gsl.h"
+#include <gsl/gsl_interp.h>
+#include <gsl/gsl_spline.h>
+#ifdef HAVE_NARRAY_H
+#include "narray.h"
+#endif
+
+
+typedef struct {
+  gsl_interp *p;
+  gsl_interp_accel *a;
+} rb_gsl_interp;
+
+typedef struct {
+  gsl_spline *s;
+  gsl_interp_accel *a;
+} rb_gsl_spline;
+
+enum {
+  GSL_INTERP_LINEAR,
+  GSL_INTERP_POLYNOMIAL,
+  GSL_INTERP_CSPLINE,
+  GSL_INTERP_CSPLINE_PERIODIC,
+  GSL_INTERP_AKIMA,
+  GSL_INTERP_AKIMA_PERIODIC,
+};
+
+const gsl_interp_type* get_interp_type(VALUE t);
+
+#endif

data/include/rb_gsl_linalg.h

@@ -0,0 +1,25 @@
+/*
+  rb_gsl_linalg.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY
+*/
+
+#ifndef ___RB_GSL_LINALG_H___
+#define ___RB_GSL_LINALG_H___
+
+#include "gsl/gsl_linalg.h"
+#include "gsl/gsl_math.h"
+
+#ifdef HAVE_NARRAY_H
+#include "rb_gsl_with_narray.h"
+#endif
+
+VALUE rb_gsl_linalg_complex_LU_decomp(int argc, VALUE *argv, VALUE obj);
+VALUE rb_gsl_linalg_complex_LU_decomp2(int argc, VALUE *argv, VALUE obj);
+
+#endif

data/include/rb_gsl_math.h

@@ -0,0 +1,26 @@
+/*
+  rb_gsl_math.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_MATH_H___
+#define ___RB_GSL_MATH_H___
+
+#include "ruby.h"
+#include <gsl/gsl_math.h>
+#include "rb_gsl.h"
+
+#ifndef GSL_1_3_LATER
+int gsl_fcmp (const double x1, const double x2, const double epsilon);
+#endif
+
+VALUE rb_gsl_math_complex_eval(gsl_complex (*func)(gsl_complex), VALUE obj);
+
+#endif

data/include/rb_gsl_odeiv.h

@@ -0,0 +1,21 @@
+/*
+  rb_gsl_odeiv.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_ODEIV_H___
+#define ___RB_GSL_ODEIV_H___
+
+#include "ruby.h"
+#include <gsl/gsl_odeiv.h>
+#include "rb_gsl.h"
+#include "rb_gsl_array.h"
+
+#endif

data/include/rb_gsl_poly.h

@@ -0,0 +1,71 @@
+/*
+  rb_gsl_poly.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or POLYNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_POLY_H___
+#define ___RB_GSL_POLY_H___
+
+#include "ruby.h"
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_poly.h>
+#include <gsl/gsl_complex.h>
+#include <gsl/gsl_complex_math.h>
+#include "rb_gsl_complex.h"
+#include "rb_gsl_array.h"
+
+EXTERN VALUE cgsl_poly;
+EXTERN VALUE cgsl_poly_int;
+EXTERN VALUE cgsl_poly_dd;
+EXTERN VALUE cgsl_poly_taylor;
+EXTERN VALUE cgsl_poly_workspace;
+EXTERN VALUE cgsl_rational;
+
+typedef gsl_vector gsl_poly;
+typedef gsl_vector_int gsl_poly_int;
+/*
+typedef struct ___gsl_rational 
+{
+  VALUE num, den;
+  gsl_poly *pnum;
+  gsl_poly *pden;
+} gsl_rational;
+*/
+
+int gsl_poly_conv(const double *a, size_t na, const double *b, size_t nb,
+		  double *c, size_t *nc);
+
+gsl_vector* gsl_poly_deconv_vector(const gsl_vector *c, const gsl_vector *a, gsl_vector **r);
+gsl_vector* gsl_poly_deriv(const gsl_vector *v);
+gsl_vector* gsl_poly_integ(const gsl_vector *v);
+gsl_vector* gsl_poly_reduce(const gsl_vector *v);
+VALUE rb_gsl_poly_complex_solve(int argc, VALUE *argv, VALUE obj);
+gsl_vector* gsl_poly_conv_vector(const gsl_vector *v1, const gsl_vector *v2);
+gsl_poly* gsl_poly_add(const gsl_poly *a, const gsl_poly *b);
+VALUE rb_gsl_poly_deconv(VALUE obj, VALUE bb);
+
+VALUE rb_gsl_poly_complex_solve2(int argc, VALUE *argv, VALUE obj);
+
+
+int gsl_poly_int_conv(const int *a, size_t na, const int *b, size_t nb,
+		  int *c, size_t *nc);
+
+gsl_vector_int* gsl_poly_int_deconv_vector(const gsl_vector_int *c, const gsl_vector_int *a, gsl_vector_int **r);
+gsl_vector_int* gsl_poly_int_deriv(const gsl_vector_int *v);
+gsl_vector_int* gsl_poly_int_integ(const gsl_vector_int *v);
+gsl_vector_int* gsl_poly_int_reduce(const gsl_vector_int *v);
+VALUE rb_gsl_poly_int_complex_solve(int argc, VALUE *argv, VALUE obj);
+gsl_vector_int* gsl_poly_int_conv_vector(const gsl_vector_int *v1, const gsl_vector_int *v2);
+gsl_poly_int* gsl_poly_int_add(const gsl_poly_int *a, const gsl_poly_int *b);
+VALUE rb_gsl_poly_int_deconv(VALUE obj, VALUE bb);
+
+gsl_poly* get_poly_get(VALUE obj, int *flag);
+gsl_poly_int* get_poly_int_get(VALUE obj, int *flag);
+#endif

data/include/rb_gsl_rational.h

@@ -0,0 +1,37 @@
+/*
+  rb_gsl_rational.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or POLYNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_RATIONAL_H___
+#define ___RB_GSL_RATIONAL_H___
+
+#include "ruby.h"
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_poly.h>
+#include <gsl/gsl_complex.h>
+#include <gsl/gsl_complex_math.h>
+#include "rb_gsl_complex.h"
+#include "rb_gsl_array.h"
+#include "rb_gsl_poly.h"
+
+typedef struct ___gsl_rational 
+{
+  VALUE num, den;
+  gsl_poly *pnum;
+  gsl_poly *pden;
+} gsl_rational;
+
+gsl_rational* gsl_rational_alloc();
+gsl_rational* gsl_rational_new(const gsl_poly *num, const gsl_poly *den);
+gsl_rational* gsl_rational_new2(const gsl_poly *num, const gsl_poly *den);
+void gsl_rational_free(gsl_rational *r);
+
+#endif

data/include/rb_gsl_rng.h

@@ -0,0 +1,21 @@
+/*
+  rb_gsl_rng.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY
+*/
+
+#ifndef ___RB_GSL_RNG_H___
+#define ___RB_GSL_RNG_H___
+
+#include "ruby.h"
+#include <gsl/gsl_rng.h>
+#include "rb_gsl.h"
+
+EXTERN VALUE cgsl_rng;
+
+#endif

data/include/rb_gsl_root.h

@@ -0,0 +1,22 @@
+/*
+  rb_gsl_roots.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY
+*/
+
+#ifndef ___RB_GSL_ROOT_H___
+#define ___RB_GSL_ROOT_H___
+
+#include <gsl/gsl_errno.h>
+#include <gsl/gsl_math.h>
+#include <gsl/gsl_roots.h>
+#include "rb_gsl.h"
+EXTERN VALUE cgsl_fsolver;
+EXTERN VALUE cgsl_fdfsolver;
+
+#endif

data/include/rb_gsl_sf.h

@@ -0,0 +1,119 @@
+/*
+  rb_gsl_sf.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY
+*/
+#ifndef ___RB_GSL_SF_H___
+#define ___RB_GSL_SF_H___
+
+#include "ruby.h"
+#include <gsl/gsl_sf.h>
+#include <gsl/gsl_errno.h>
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_math.h>
+#include <gsl/gsl_complex.h>
+#include "rb_gsl_config.h"
+#include "rb_gsl.h"
+
+EXTERN VALUE cgsl_sf_result, cgsl_sf_result_e10;
+
+VALUE rb_gsl_sf_result_new(VALUE klass);
+
+VALUE rb_gsl_sf_eval1(double (*func)(double), VALUE argv);
+VALUE rb_gsl_sf_eval_int_double(double (*func)(int, double), VALUE jj, VALUE argv);
+VALUE rb_gsl_sf_eval_double_double(double (*func)(double, double), VALUE ff, VALUE argv);
+VALUE rb_gsl_sf_eval1_uint(double (*func)(unsigned int), VALUE argv);
+VALUE rb_gsl_sf_eval_double_m(double (*func)(double, gsl_mode_t), VALUE argv, VALUE m);;
+VALUE rb_gsl_sf_eval_int_double_double(double (*func)(int, double, double), VALUE jj, 
+				       VALUE ff, VALUE argv);
+VALUE rb_gsl_sf_eval_int_int_double(double (*func)(int, int, double), VALUE jj,
+				    VALUE jj2, VALUE argv);
+VALUE rb_gsl_sf_eval_double3(double (*func)(double, double, double), 
+					   VALUE ff, VALUE ff2, VALUE argv);
+VALUE rb_gsl_sf_eval_double4(double (*func)(double, double, double, double), 
+			     VALUE ff, VALUE ff2, VALUE ff3, VALUE argv);
+VALUE rb_gsl_sf_eval_double_int(double (*func)(double, int), VALUE argv, VALUE jj);
+VALUE rb_gsl_sf_eval1_int(double (*func)(int), VALUE argv);
+VALUE rb_gsl_sf_eval_e(int (*func)(double, gsl_sf_result*), VALUE x);
+VALUE rb_gsl_sf_eval_e_int_double(int (*func)(int, double, gsl_sf_result*), 
+				  VALUE n, VALUE x);
+VALUE rb_gsl_sf_eval_e_int_int_double(int (*func)(int, int, double, gsl_sf_result*), 
+				      VALUE n1, VALUE n2, VALUE x);
+VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*), 
+			       VALUE x1, VALUE x2);
+VALUE rb_gsl_sf_eval_e_uint(int (*func)(unsigned int, gsl_sf_result*), VALUE x);
+VALUE rb_gsl_sf_eval_e_int_uint(int (*func)(int, unsigned int, gsl_sf_result*), 
+				VALUE n, VALUE x);
+VALUE rb_gsl_sf_eval_e_double_uint(int (*func)(double, unsigned int, gsl_sf_result*), 
+				   VALUE y, VALUE x);
+VALUE rb_gsl_sf_eval_e_m(int (*func)(double, gsl_mode_t, gsl_sf_result*), 
+			 VALUE x, VALUE m);
+VALUE rb_gsl_sf_eval_e_double2_m(int (*func)(double, double, gsl_mode_t, gsl_sf_result*), 
+				 VALUE x1, VALUE x2, VALUE m);
+VALUE rb_gsl_sf_eval_e_double3_m(int (*func)(double, double, double, gsl_mode_t, gsl_sf_result*), 
+				 VALUE x1, VALUE x2, VALUE x3, VALUE m);
+VALUE rb_gsl_sf_eval_double2_m(double (*func)(double, double, gsl_mode_t), 
+			       VALUE argv, VALUE x2, VALUE m);
+VALUE rb_gsl_sf_eval_double3_m(double (*func)(double, double, double, gsl_mode_t), 
+			       VALUE argv, VALUE x2, VALUE x3, VALUE m);
+
+VALUE rb_gsl_sf_eval_e_double4_m(int (*func)(double, double, double, double, gsl_mode_t, gsl_sf_result*), 
+				 VALUE x1, VALUE x2, VALUE x3, VALUE x4, VALUE m);
+
+VALUE rb_gsl_sf_eval_e_int(int (*func)(int, gsl_sf_result*), VALUE x);
+VALUE rb_gsl_sf_eval_e_double3(int (*func)(double, double, double, gsl_sf_result*), 
+			      VALUE x1, VALUE x2, VALUE x3);
+VALUE rb_gsl_sf_eval_e_int_double2(int (*func)(int, double, double, gsl_sf_result*), 
+				   VALUE n, VALUE x1, VALUE x2);
+VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*), 
+			       VALUE x1, VALUE x2);
+
+VALUE eval_sf(double (*func)(double, gsl_mode_t), VALUE argv);
+
+VALUE rb_gsl_sf_eval_double4_m(double (*func)(double, double, double, double,
+					      gsl_mode_t), 
+			       VALUE argv, VALUE x2, VALUE x3, VALUE x4, VALUE m);
+
+VALUE rb_gsl_sf_eval_complex(double (*f)(double), VALUE obj);
+
+void Init_gsl_sf_airy(VALUE module);
+void Init_gsl_sf_bessel(VALUE module);
+void Init_gsl_sf_clausen(VALUE module);
+void Init_gsl_sf_coulomb(VALUE module);
+void Init_gsl_sf_coupling(VALUE module);
+void Init_gsl_sf_dawson(VALUE module);
+void Init_gsl_sf_debye(VALUE module);
+void Init_gsl_sf_dilog(VALUE module);
+void Init_gsl_sf_elementary(VALUE module);
+void Init_gsl_sf_ellint(VALUE module);
+void Init_gsl_sf_elljac(VALUE module);
+void Init_gsl_sf_erfc(VALUE module);
+void Init_gsl_sf_exp(VALUE module);
+void Init_gsl_sf_expint(VALUE module);
+void Init_gsl_sf_fermi_dirac(VALUE module);
+void Init_gsl_sf_gamma(VALUE module);
+void Init_gsl_sf_gegenbauer(VALUE module);
+void Init_gsl_sf_hyperg(VALUE module);
+void Init_gsl_sf_laguerre(VALUE module);
+void Init_gsl_sf_lambert(VALUE module);
+void Init_gsl_sf_legendre(VALUE module);
+void Init_gsl_sf_log(VALUE module);
+void Init_gsl_sf_power(VALUE module);
+void Init_gsl_sf_psi(VALUE module);
+void Init_gsl_sf_synchrotron(VALUE module);
+void Init_gsl_sf_transport(VALUE module);
+void Init_gsl_sf_trigonometric(VALUE module);
+void Init_gsl_sf_zeta(VALUE module);
+
+
+#ifdef GSL_1_9_LATER
+#include <gsl/gsl_sf_mathieu.h>
+void Init_sf_mathieu(VALUE module);
+#endif
+
+#endif

data/include/rb_gsl_statistics.h

@@ -0,0 +1,17 @@
+/*
+  rb_gsl_statistics.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2004 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY
+*/
+
+#ifndef ___RB_GSL_STATISTICS_H___
+#define ___RB_GSL_STATISTICS_H___
+
+#include "gsl/gsl_statistics.h"
+
+#endif

data/include/rb_gsl_tensor.h

@@ -0,0 +1,45 @@
+#ifndef ___RB_TENSOR_H___
+#define ___RB_TENSOR_H___
+
+#ifdef HAVE_TENSOR_TENSOR_H
+#include "rb_gsl.h"
+#include "rb_gsl_common.h"
+//#include <tensor/tensor.h>
+#include <tensor/tensor.h>
+
+EXTERN VALUE cgsl_tensor, cgsl_tensor_int;
+
+enum {
+  TENSOR_ADD,
+  TENSOR_SUB,
+  TENSOR_MUL_ELEMENTS,
+  TENSOR_DIV_ELEMENTS,
+  TENSOR_SCALE,
+  TENSOR_ADD_CONSTANT,
+  TENSOR_ADD_DIAGONAL,
+  TENSOR_PRODUCT,
+  TENSOR_CONTRACT,
+};
+
+typedef gsl_permutation tensor_indices;
+
+typedef struct __rbgsl_tensor {
+  /*  tensor *tensor;
+      tensor_indices *indices;*/
+  tensor *tensor;
+  tensor_indices *indices;
+} rbgsl_tensor;
+
+typedef struct __rbgsl_tensor_int {
+  /*  tensor_int *tensor;
+      tensor_indices *indices;*/
+  tensor_int *tensor;
+  tensor_indices *indices;
+} rbgsl_tensor_int;
+
+rbgsl_tensor* rbgsl_tensor_alloc(const unsigned int rank, const size_t dimension);
+rbgsl_tensor_int* rbgsl_tensor_int_alloc(const unsigned int rank, const size_t dimension);
+void rbgsl_tensor_free(rbgsl_tensor*);
+void rbgsl_tensor_int_free(rbgsl_tensor_int*);
+#endif
+#endif