rb-gsl 1.15.3.1

data/include/rb_gsl_complex.h

@@ -0,0 +1,25 @@
+/*
+  rb_gsl_complex.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_COMPLEX_H___
+#define ___RB_GSL_COMPLEX_H___
+
+#include <math.h>
+#include "ruby.h"
+#include <gsl/gsl_complex.h>
+#include <gsl/gsl_complex_math.h>
+
+EXTERN VALUE cgsl_complex;
+VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj);
+VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj);
+
+#endif

data/include/rb_gsl_config.h

@@ -0,0 +1,62 @@
+#ifndef ___RB_GSL_CONFIG_H___
+#define ___RB_GSL_CONFIG_H___
+
+#ifndef GSL_VERSION
+#define GSL_VERSION "1.15"
+#endif
+#ifndef GSL_0_9_4_LATER
+#define GSL_0_9_4_LATER
+#endif
+#ifndef GSL_1_0_LATER
+#define GSL_1_0_LATER
+#endif
+#ifndef GSL_1_1_LATER
+#define GSL_1_1_LATER
+#endif
+#ifndef GSL_1_1_1_LATER
+#define GSL_1_1_1_LATER
+#endif
+#ifndef GSL_1_2_LATER
+#define GSL_1_2_LATER
+#endif
+#ifndef GSL_1_3_LATER
+#define GSL_1_3_LATER
+#endif
+#ifndef GSL_1_4_LATER
+#define GSL_1_4_LATER
+#endif
+#ifndef GSL_1_4_9_LATER
+#define GSL_1_4_9_LATER
+#endif
+#ifndef GSL_1_6_LATER
+#define GSL_1_6_LATER
+#endif
+#ifndef GSL_1_8_LATER
+#define GSL_1_8_LATER
+#endif
+#ifndef GSL_1_9_LATER
+#define GSL_1_9_LATER
+#endif
+#ifndef GSL_1_10_LATER
+#define GSL_1_10_LATER
+#endif
+#ifndef GSL_1_11_LATER
+#define GSL_1_11_LATER
+#endif
+#ifndef GSL_1_13_LATER
+#define GSL_1_13_LATER
+#endif
+#ifndef GSL_1_14_LATER
+#define GSL_1_14_LATER
+#endif
+#ifndef GSL_1_15_LATER
+#define GSL_1_15_LATER
+#endif
+#ifndef RUBY_GSL_VERSION
+#define RUBY_GSL_VERSION "1.15.0"
+#endif
+#ifndef RUBY_1_8_LATER
+#define RUBY_1_8_LATER
+#endif
+
+#endif

data/include/rb_gsl_const.h

@@ -0,0 +1,29 @@
+/*
+  rb_gsl.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_CONST_H___
+#define ___RB_GSL_CONST_H___
+
+#include "ruby.h"
+#include "rb_gsl.h"
+#ifdef GSL_CONST_OLD
+#include <gsl/gsl_const_mks.h>
+#include <gsl/gsl_const_cgs.h>
+#else
+#include <gsl/gsl_const_mksa.h>
+#include <gsl/gsl_const_cgsm.h>
+#endif
+#include <gsl/gsl_const_num.h>
+
+EXTERN VALUE mgsl_const_mks, mgsl_const_cgs;
+
+#endif

data/include/rb_gsl_dirac.h

@@ -0,0 +1,13 @@
+#ifndef ___RB_GSL_DIRAC_H___
+#define ___RB_GSL_DIRAC_H___
+
+#include "rb_gsl.h"
+#include "rb_gsl_common.h"
+#include "rb_gsl_complex.h"
+#include "rb_gsl_array.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <math.h>
+
+#endif

data/include/rb_gsl_eigen.h

@@ -0,0 +1,17 @@
+/*
+  rb_gsl_eigen.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_EIGEN_H___
+#define ___RB_GSL_EIGEN_H___
+
+#include "gsl/gsl_eigen.h"
+
+#endif

data/include/rb_gsl_fft.h

@@ -0,0 +1,62 @@
+/*
+  rb_gsl_fft.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 FFTNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_FFT_H___
+#define ___RB_GSL_FFT_H___
+
+
+#include <stdio.h>
+#include <math.h>
+#include "ruby.h"
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_errno.h>
+#include <gsl/gsl_fft.h>
+#include <gsl/gsl_fft_complex.h>
+#include <gsl/gsl_fft_real.h>
+#include <gsl/gsl_fft_halfcomplex.h>
+#include "rb_gsl_array.h"
+#include "rb_gsl.h"
+
+#include "rb_gsl_complex.h"
+#include "rb_gsl_array.h"
+#ifdef HAVE_NARRAY_H
+#include "narray.h"
+#endif
+
+typedef struct
+{
+  size_t n;
+  size_t nf;
+  size_t factor[64];
+  gsl_complex *twiddle[64];
+  gsl_complex *trig;
+} GSL_FFT_Wavetable;
+
+typedef struct
+{
+  size_t n;
+  double *scratch;
+} GSL_FFT_Workspace;
+
+enum {
+  RB_GSL_FFT_INPLACE,
+  RB_GSL_FFT_COPY,
+};
+
+EXTERN VALUE mgsl_fft;
+EXTERN VALUE cgsl_fft_wavetable;
+EXTERN VALUE cgsl_fft_wavetable_factor;
+EXTERN VALUE cgsl_fft_complex_wavetable, cgsl_fft_complex_workspace;
+EXTERN VALUE cgsl_fft_real_wavetable, cgsl_fft_halfcomplex_wavetable;
+EXTERN VALUE cgsl_fft_real_workspace;
+
+#endif

data/include/rb_gsl_fit.h

@@ -0,0 +1,25 @@
+/*
+  rb_gsl_fit.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_FIT_H___
+#define ___RB_GSL_FIT_H___
+
+#include "ruby.h"
+#include <gsl/gsl_blas.h>
+#include <gsl/gsl_fit.h>
+#include <gsl/gsl_multifit.h>
+#include <gsl/gsl_multifit_nlin.h>
+#include "rb_gsl_array.h"
+
+EXTERN VALUE mgsl_multifit;
+
+#endif

data/include/rb_gsl_function.h

@@ -0,0 +1,27 @@
+/*
+  rb_gsl_function.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_FUNCTION_H___
+#define ___RB_GSL_FUNCTION_H___
+
+#include "ruby.h"
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_integration.h>
+#include <gsl/gsl_math.h>
+#include "rb_gsl.h"
+
+EXTERN VALUE cgsl_function;
+EXTERN VALUE cgsl_function_fdf;
+extern ID RBGSL_ID_call, RBGSL_ID_arity;
+void gsl_function_mark(gsl_function *f);
+void gsl_function_free(gsl_function *f);
+#endif

data/include/rb_gsl_graph.h

@@ -0,0 +1,70 @@
+/*
+  rb_gsl_graph.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 CHEBNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#ifndef ___RB_GSL_GRAPH_H___
+#define ___RB_GSL_GRAPH_H___
+
+#include "rb_gsl.h"
+#include "rb_gsl_array.h"
+#include "rb_gsl_common.h"
+
+typedef struct __rb_gsl_plot {
+  VALUE xdata, ydata;
+  VALUE T;
+  VALUE E;
+  VALUE f;
+  VALUE F;
+  VALUE g;
+  VALUE h;
+  VALUE k;
+  VALUE K;
+  VALUE l;
+  VALUE L;
+  VALUE N;
+  VALUE r;
+  VALUE R;
+  VALUE s;
+  VALUE t;
+  VALUE u;
+  VALUE w;
+  VALUE x;
+  VALUE y;
+  VALUE bg;
+  VALUE bitmap_size;
+  VALUE frame;
+  VALUE frame_line_width;
+  VALUE max_line_length;
+  VALUE page_size;
+  VALUE pen_colors;
+  VALUE rotation;
+  VALUE title_font_name;
+  VALUE title_font_size;
+  VALUE rotate_y_label;
+  VALUE I;
+  VALUE B;
+  VALUE m;
+  VALUE S;
+  VALUE W;
+  VALUE q;
+  VALUE C;
+  VALUE symbol_font_name;
+  VALUE reposition;
+  VALUE blankout;
+  VALUE O;
+  VALUE X, Y; 
+  
+} gsl_graph;
+
+gsl_graph* gsl_graph_new();
+void gsl_graph_free(gsl_graph *g);
+
+#endif

data/include/rb_gsl_histogram.h

@@ -0,0 +1,63 @@
+/*
+  rb_gsl_histogram.h
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001 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_HISTOGRAM_H___
+#define ___RB_GSL_HISTOGRAM_H___
+
+#include <gsl/gsl_vector.h>
+#include <gsl/gsl_histogram.h>
+#include <gsl/gsl_histogram2d.h>
+#include "rb_gsl.h"
+
+EXTERN VALUE cgsl_histogram;
+EXTERN VALUE cgsl_histogram_range;
+EXTERN VALUE cgsl_histogram_bin;
+EXTERN VALUE cgsl_histogram2d;
+EXTERN VALUE cgsl_histogram2d_view;
+
+typedef struct {
+  gsl_histogram h;
+} mygsl_histogram2d_view;
+
+#ifndef HISTOGRAM2D_P
+#define HISTOGRAM2D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram2d)?1:0)
+#endif
+
+#ifndef CHECK_HISTOGRAM2D
+#define CHECK_HISTOGRAM2D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram2d))\
+    rb_raise(rb_eTypeError, "wrong type (Histogram2d expected)");
+#endif
+
+
+#ifndef HISTOGRAM3D_P
+#define HISTOGRAM3D_P(x) (rb_obj_is_kind_of(x,cgsl_histogram3d)?1:0)
+#endif
+
+#ifndef CHECK_HISTOGRAM3D
+#define CHECK_HISTOGRAM3D(x) if(!rb_obj_is_kind_of(x,cgsl_histogram3d))\
+    rb_raise(rb_eTypeError, "wrong type (Histogram3d expected)");
+#endif
+
+#include "rb_gsl_histogram3d.h"
+
+int
+mygsl_histogram_equal_bins_p (const gsl_histogram * h1, const gsl_histogram * h2);
+int 
+mygsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2);
+int 
+mygsl_histogram_sub (gsl_histogram * h1, const gsl_histogram * h2);
+int 
+mygsl_histogram_mul (gsl_histogram * h1, const gsl_histogram * h2);
+int 
+mygsl_histogram_div (gsl_histogram * h1, const gsl_histogram * h2);
+
+#endif

data/include/rb_gsl_histogram3d.h

@@ -0,0 +1,97 @@
+#ifndef ___RB_GSL_HISTOGRAM3D_H___
+#define ___RB_GSL_HISTOGRAM3D_H___
+
+typedef struct {
+  size_t nx, ny, nz;
+  double * xrange ;
+  double * yrange ;
+  double * zrange ;
+  double * bin ;
+} mygsl_histogram3d ;
+
+typedef struct {
+  gsl_histogram2d h;
+} mygsl_histogram3d_view;
+
+
+mygsl_histogram3d* mygsl_histogram3d_alloc(const size_t nx, const size_t ny,
+					   const size_t nz);
+void mygsl_histogram3d_free (mygsl_histogram3d * h);
+mygsl_histogram3d* mygsl_histogram3d_calloc_uniform(const size_t nx, 
+						    const size_t ny,
+						    const size_t nz,
+						    const double xmin,
+						    const double xmax,
+						    const double ymin,
+						    const double ymax,
+						    const double zmin,
+						    const double zmax);
+mygsl_histogram3d* mygsl_histogram3d_calloc(const size_t nx, 
+					    const size_t ny,
+					    const size_t nz);
+int mygsl_histogram3d_set_ranges_uniform (mygsl_histogram3d * h, 
+					  double xmin, double xmax,
+					  double ymin, double ymax,
+					  double zmin, double zmax);
+int mygsl_histogram3d_set_ranges (mygsl_histogram3d * h, 
+				  const double xrange[], size_t xsize,
+				  const double yrange[], size_t ysize,
+				  const double zrange[], size_t zsize);
+int mygsl_histogram3d_memcpy(mygsl_histogram3d * dest, const mygsl_histogram3d * src);
+mygsl_histogram3d* mygsl_histogram3d_clone(const mygsl_histogram3d * src);
+int mygsl_histogram3d_fread(FILE * stream, mygsl_histogram3d * h);
+int mygsl_histogram3d_fwrite(FILE * stream, const mygsl_histogram3d * h);
+int mygsl_histogram3d_increment(mygsl_histogram3d * h, double x, double y, double z);
+int mygsl_histogram3d_accumulate (mygsl_histogram3d * h, 
+				  double x, double y, double z, double weight);
+int mygsl_histogram3d_increment2(mygsl_histogram3d * h, 
+				 double x, double y, double z);
+int mygsl_histogram3d_accumulate2(mygsl_histogram3d * h, 
+				  double x, double y, double z, double weight);
+double mygsl_histogram3d_get (const mygsl_histogram3d * h, const size_t i, 
+			      const size_t j, const size_t k);
+int mygsl_histogram3d_get_xrange(const mygsl_histogram3d * h, const size_t i,  
+				 double *xlower, double *xupper);
+int mygsl_histogram3d_get_yrange(const mygsl_histogram3d * h, const size_t j,  
+				 double *ylower, double *yupper);
+int mygsl_histogram3d_get_zrange(const mygsl_histogram3d * h, const size_t k,  
+				 double *zlower, double *zupper);
+int mygsl_histogram3d_find (const mygsl_histogram3d * h,
+			    const double x, const double y, const double z,
+			    size_t * i, size_t * j, size_t *k);
+gsl_histogram2d* mygsl_histogram3d_xyproject(const mygsl_histogram3d * h3,
+					     size_t kstart, size_t kend);
+gsl_histogram2d* mygsl_histogram3d_xzproject(const mygsl_histogram3d * h3,
+					     size_t jstart, size_t jend);
+gsl_histogram2d* mygsl_histogram3d_yzproject(const mygsl_histogram3d * h3,
+					     size_t istart, size_t iend);
+int mygsl_histogram3d_scale(mygsl_histogram3d * h, double scale);
+int mygsl_histogram3d_shift(mygsl_histogram3d * h, double shift);
+double mygsl_histogram3d_xmax(const mygsl_histogram3d * h);
+double mygsl_histogram3d_xmin(const mygsl_histogram3d * h);
+double mygsl_histogram3d_ymax(const mygsl_histogram3d * h);
+double mygsl_histogram3d_ymin(const mygsl_histogram3d * h);
+double mygsl_histogram3d_zmax(const mygsl_histogram3d * h);
+double mygsl_histogram3d_zmin(const mygsl_histogram3d * h);
+double mygsl_histogram3d_max_val(const mygsl_histogram3d * h);
+void mygsl_histogram3d_max_bin(const mygsl_histogram3d * h, 
+			       size_t *imax_out, size_t *jmax_out, size_t *kmax_out);
+double mygsl_histogram3d_min_val(const mygsl_histogram3d * h);
+void mygsl_histogram3d_min_bin(const mygsl_histogram3d * h, 
+			       size_t *imin_out, size_t *jmin_out, size_t *kmin_out);
+double mygsl_histogram3d_sum (const mygsl_histogram3d * h);
+double mygsl_histogram3d_xmean (const mygsl_histogram3d * h);
+double mygsl_histogram3d_ymean (const mygsl_histogram3d * h);
+double mygsl_histogram3d_zmean (const mygsl_histogram3d * h);
+double mygsl_histogram3d_xsigma(const mygsl_histogram3d * h);
+double mygsl_histogram3d_ysigma(const mygsl_histogram3d * h);
+double mygsl_histogram3d_zsigma(const mygsl_histogram3d * h);
+void mygsl_histogram3d_reset(mygsl_histogram3d * h);
+int mygsl_histogram3d_equal_bins_p(const mygsl_histogram3d * h1,
+				   const mygsl_histogram3d * h2);
+int mygsl_histogram3d_add(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
+int mygsl_histogram3d_sub(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
+int mygsl_histogram3d_mul(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
+int mygsl_histogram3d_div(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
+
+#endif