rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/include/rb_gsl_array.h

@@ -1,214 +0,0 @@
-/*
-  rb_gsl_array.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_ARRAY_H___
-#define ___RB_GSL_ARRAY_H___
-
-#include <math.h>
-#include <gsl/gsl_block.h>
-#include <gsl/gsl_vector_complex.h>
-#include <gsl/gsl_vector_complex_double.h>
-#include <gsl/gsl_matrix_complex_double.h>
-#include <gsl/gsl_combination.h>
-#include <gsl/gsl_sort.h>
-#include <gsl/gsl_sort_vector.h>
-#include <gsl/gsl_linalg.h>
-#include <gsl/gsl_eigen.h>
-#include <gsl/gsl_odeiv.h>
-#include <gsl/gsl_blas.h>
-#include "rb_gsl_common.h"
-
-typedef gsl_permutation gsl_index;
-
-#ifdef HAVE_NARRAY_H
-EXTERN VALUE cNArray;
-#endif
-
-EXTERN VALUE cgsl_block, cgsl_block_int;
-EXTERN VALUE cgsl_block_uchar;
-EXTERN VALUE cgsl_block_complex;
-EXTERN VALUE cgsl_vector, cgsl_vector_complex;
-EXTERN VALUE cgsl_vector_col;
-EXTERN VALUE cgsl_vector_col_view;
-EXTERN VALUE cgsl_vector_complex_col;
-EXTERN VALUE cgsl_vector_complex_col_view;
-EXTERN VALUE cgsl_vector_view, cgsl_vector_complex_view;
-EXTERN VALUE cgsl_vector_view_ro, cgsl_vector_col_view_ro;
-EXTERN VALUE cgsl_vector_complex_view_ro;
-
-EXTERN VALUE cgsl_vector_int, cgsl_vector_int_col;
-EXTERN VALUE cgsl_vector_int_view, cgsl_vector_int_col_view;
-EXTERN VALUE cgsl_vector_int_view_ro, cgsl_vector_int_col_view_ro;
-
-EXTERN VALUE cgsl_matrix, cgsl_matrix_complex;
-EXTERN VALUE cgsl_matrix_view_ro;
-EXTERN VALUE cgsl_matrix_complex_view_ro;
-EXTERN VALUE cgsl_matrix_view, cgsl_matrix_complex_view;
-EXTERN VALUE cgsl_matrix_int, cgsl_matrix_int_view;
-EXTERN VALUE cgsl_matrix_int_view_ro;
-EXTERN VALUE cgsl_permutation;
-EXTERN VALUE cgsl_index;
-EXTERN VALUE cgsl_function;
-EXTERN VALUE mgsl_narray;
-
-EXTERN VALUE mDirac;
-
-gsl_matrix_view* gsl_matrix_view_alloc();
-void gsl_matrix_view_free(gsl_matrix_view * mv);
-gsl_vector_view* gsl_vector_view_alloc();
-void gsl_vector_view_free(gsl_vector_view * v);
-gsl_vector_complex_view* gsl_vector_complex_view_alloc();
-void gsl_vector_complex_view_free(gsl_vector_view * vv);
-gsl_matrix_complex_view* gsl_matrix_complex_view_alloc();
-void gsl_matrix_complex_view_free(gsl_matrix_view * vv);
-
-VALUE rb_gsl_vector_new(int argc, VALUE *argv, VALUE klass);
-
-gsl_vector* get_cvector(VALUE v);
-VALUE make_rarray_from_cvector(const gsl_vector *v);
-VALUE make_rarray_from_cpermutation(const gsl_permutation *v);
-gsl_vector* make_cvector_from_rarrays(VALUE a);
-gsl_vector* make_cvector_from_rarray(VALUE a);
-void cvector_set_from_carray(gsl_vector *v, const double *a);
-void cvector_set_from_rarrays(gsl_vector *v, VALUE ary);
-void cvector_set_from_rarray(gsl_vector *v, VALUE ary);
-void carray_set_from_cvector(double *a, const gsl_vector *v);
-void carray_set_from_rarrays(double *a, VALUE ary);
-void carray_set_from_rarray(double *a, VALUE ary);
-int is_vector_p(VALUE obj);
-void check_vector(VALUE obj);
-int is_vector_complex_p(VALUE obj);
-void check_vector_complex(VALUE obj);
-int is_matrix_p(VALUE obj);
-void check_matrix(VALUE obj);
-int is_matrix_complex_p(VALUE obj);
-void check_matrix_complex(VALUE obj);
-gsl_complex* make_complex(double re, double im);
-int is_permutation_p(VALUE obj);
-void check_permutation(VALUE obj);
-int is_combination_p(VALUE obj);
-void check_combination(VALUE obj);
-gsl_vector* get_vector(VALUE ary);
-gsl_matrix* make_matrix_clone(const gsl_matrix *m);
-gsl_matrix_int* make_matrix_int_clone(const gsl_matrix_int *m);
-VALUE make_matrix_clone2(VALUE vm);
-gsl_matrix_complex* make_matrix_complex_clone(const gsl_matrix_complex *m);
-int is_matrix_complex_p(VALUE obj);
-void check_matrix_complex(VALUE obj);
-
-double* get_vector_ptr(VALUE ary, size_t *stride, size_t *n);
-
-
-gsl_matrix_complex* matrix_to_complex(const gsl_matrix *m);
-
-void gsl_matrix_complex_mul(gsl_matrix_complex *mnew, const gsl_matrix_complex *m,
-                            const gsl_matrix_complex *mb);
-void gsl_matrix_mul(gsl_matrix *mnew, gsl_matrix *m, gsl_matrix *b);
-void gsl_matrix_complex_mul_vector(gsl_vector_complex *vnew,
-                                   const gsl_matrix_complex *m,
-                                   const gsl_vector_complex *v);
-void gsl_matrix_mul_vector(gsl_vector *vnew,
-                           const gsl_matrix *m, const gsl_vector *v);
-gsl_vector_complex* vector_to_complex(const gsl_vector *v);
-
-gsl_vector* make_vector_clone(const gsl_vector *v);
-gsl_vector_complex* make_vector_complex_clone(const gsl_vector_complex *v);
-int gsl_vector_complex_add(gsl_vector_complex *cv, const gsl_vector_complex *cv2);
-int gsl_vector_complex_sub(gsl_vector_complex *cv, const gsl_vector_complex *cv2);
-int gsl_vector_complex_mul(gsl_vector_complex *cv, const gsl_vector_complex *cv2);
-int gsl_vector_complex_div(gsl_vector_complex *cv, const gsl_vector_complex *cv2);
-int gsl_vector_complex_add_constant(gsl_vector_complex *cv, gsl_complex b);
-int gsl_vector_complex_scale(gsl_vector_complex *cv, gsl_complex b);
-gsl_vector_view* rb_gsl_make_vector_view(double *data, size_t size, size_t stride);
-gsl_vector_int_view* rb_gsl_make_vector_int_view(int *data, size_t size, size_t stride);
-
-void Init_gsl_array_complex(VALUE module);
-void Init_gsl_vector(VALUE module);
-void Init_gsl_vector_complex(VALUE module);
-void Init_gsl_matrix(VALUE module);
-void Init_gsl_matrix_complex(VALUE module);
-void Init_gsl_matrix(VALUE module);
-void Init_gsl_permutation(VALUE module);
-void Init_gsl_combination(VALUE module);
-void Init_gsl_matrix_int(VALUE module);
-
-VALUE rb_gsl_range2ary(VALUE obj);
-
-void Init_gsl_vector_int();
-gsl_vector_int_view* rb_gsl_vector_int_view_alloc(size_t n);
-void rb_gsl_vector_int_view_free(gsl_vector_int_view *v);
-gsl_vector_int* make_vector_int_clone(const gsl_vector_int *v);
-gsl_matrix_int_view* rb_gsl_matrix_int_view_alloc();
-void rb_gsl_matrix_int_view_free(gsl_matrix_int_view *v);
-VALUE rb_gsl_matrix_to_i(VALUE obj);
-VALUE rb_gsl_matrix_int_to_f(VALUE obj);
-void gsl_matrix_int_mul_vector(gsl_vector_int *vnew,
-                               const gsl_matrix_int *m, const gsl_vector_int *v);
-VALUE rb_gsl_vector_to_i(VALUE obj);
-VALUE make_rarray_from_cvector_int(const gsl_vector_int *v);
-VALUE rb_gsl_vector_int_to_f(VALUE obj);
-VALUE rb_gsl_vector_uminus(VALUE obj);
-VALUE rb_gsl_vector_print(VALUE obj);
-void gsl_vector_print(const gsl_vector *v, VALUE klass);
-int rbgsl_vector_equal(const gsl_vector *v1, const gsl_vector *v2, double eps);
-
-void Init_gsl_vector_init(VALUE module);
-void Init_gsl_vector_int_init(VALUE module);
-void Init_gsl_matrix_init(VALUE module);
-void Init_gsl_matrix_int_init(VALUE module);
-
-gsl_matrix* gsl_matrix_alloc_from_array_sizes(VALUE ary,
-                                              VALUE nn1, VALUE nn2);
-gsl_matrix* gsl_matrix_alloc_from_arrays(int argc, VALUE *argv);
-
-gsl_matrix* gsl_matrix_alloc_from_vector_sizes(VALUE ary,
-                                               VALUE nn1, VALUE nn2);
-gsl_matrix* gsl_matrix_alloc_from_vectors(int argc, VALUE *argv);
-gsl_matrix_int* gsl_matrix_int_alloc_from_array_sizes(VALUE ary,
-                                                      VALUE nn1, VALUE nn2);
-gsl_matrix_int* gsl_matrix_int_alloc_from_arrays(int argc, VALUE *argv);
-
-gsl_matrix_int* gsl_matrix_int_alloc_from_vector_sizes(VALUE ary,
-                                                       VALUE nn1, VALUE nn2);
-gsl_matrix_int* gsl_matrix_int_alloc_from_vectors(int argc, VALUE *argv);
-
-VALUE rb_gsl_matrix_do_something(VALUE obj, void (*f)(gsl_matrix *));
-VALUE rb_gsl_matrix_int_do_something(VALUE obj, void (*f)(gsl_matrix_int *));
-
-VALUE rb_gsl_matrix_power(VALUE a, VALUE b);
-VALUE rb_gsl_matrix_int_power(VALUE a, VALUE b);
-
-void mygsl_vector_shift(gsl_vector *p, size_t n);
-void mygsl_vector_int_shift(gsl_vector_int *p, size_t n);
-void mygsl_vector_shift_scale2(gsl_vector *p, size_t n);
-void mygsl_vector_int_shift_scale2(gsl_vector_int *p, size_t n);
-
-VALUE rb_gsl_vector_to_s(VALUE obj);
-VALUE rb_gsl_vector_int_to_s(VALUE obj);
-
-VALUE rb_gsl_vector_add_constant(VALUE obj, VALUE x);
-VALUE rb_gsl_vector_int_add_constant(VALUE obj, VALUE x);
-VALUE rb_gsl_vector_scale(VALUE obj, VALUE x);
-VALUE rb_gsl_vector_scale_bang(VALUE obj, VALUE x);
-VALUE rb_gsl_vector_int_scale(VALUE obj, VALUE x);
-VALUE rb_gsl_vector_int_scale_bang(VALUE obj, VALUE x);
-gsl_vector_int* make_cvector_int_from_rarray(VALUE ary);
-void cvector_int_set_from_rarray(gsl_vector_int *v, VALUE ary);
-VALUE rb_gsl_range2vector(VALUE obj);
-VALUE rb_gsl_range2vector_int(VALUE obj);
-
-void Init_gsl_block_init(VALUE module);
-void Init_gsl_block_int_init(VALUE module);
-void Init_gsl_block_uchar_init(VALUE module);
-
-void Init_gsl_matrix_nmf(void);
-
-#endif

data/ext/gsl_native/include/rb_gsl_cheb.h

@@ -1,19 +0,0 @@
-/*
-  rb_gsl_cheb.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_CHEB_H___
-#define ___RB_GSL_CHEB_H___
-
-#include "rb_gsl_function.h"
-#include <gsl/gsl_chebyshev.h>
-
-#endif

data/ext/gsl_native/include/rb_gsl_common.h

@@ -1,348 +0,0 @@
-/*
-  rb_gsl_common.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_COMMON_H___
-#define ___RB_GSL_COMMON_H___
-
-#include <ruby.h>
-#include <ruby/io.h>
-
-#include <ctype.h>
-#include <gsl/gsl_errno.h>
-#include <gsl/gsl_version.h>
-#include <gsl/gsl_math.h>
-#include <gsl/gsl_complex.h>
-#include <gsl/gsl_matrix.h>
-#include <gsl/gsl_permutation.h>
-#include <gsl/gsl_ieee_utils.h>
-#include "rb_gsl_with_narray.h"
-
-EXTERN ID rb_gsl_id_beg, rb_gsl_id_end, rb_gsl_id_excl, rb_gsl_id_to_a;
-
-#ifndef CHECK_FIXNUM
-#define CHECK_FIXNUM(x) if(!FIXNUM_P(x)) rb_raise(rb_eTypeError,"Fixnum expected");
-#endif
-
-#ifndef Need_Float
-#define Need_Float(x) (x) = rb_Float(x)
-#endif
-
-#ifndef Need_Float2
-#define Need_Float2(x,y) do { \
-    Need_Float(x); \
-    Need_Float(y); } while (0)
-#endif
-
-#ifndef COMPLEX_P
-#define COMPLEX_P(x) (rb_obj_is_kind_of(x,cgsl_complex))
-#endif
-
-#ifndef CHECK_RNG
-#define CHECK_RNG(x) if(!rb_obj_is_kind_of(x,cgsl_rng)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Rng expected)");
-#endif
-
-#ifndef CHECK_COMPLEX
-#define CHECK_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_complex)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Complex expected)");
-#endif
-
-#ifndef POLY_P
-#define POLY_P(x) (rb_obj_is_kind_of(x,cgsl_poly))
-#endif
-
-#ifndef CHECK_POLY
-#define CHECK_POLY(x) if(!rb_obj_is_kind_of(x,cgsl_poly)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Poly expected)");
-#endif
-
-#ifndef CHECK_POLY_INT
-#define CHECK_POLY_INT(x) if(!rb_obj_is_kind_of(x,cgsl_poly_int)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Poly::Int expected)");
-#endif
-
-#ifndef RATIONAL_P
-#define RATIONAL_P(x) (rb_obj_is_kind_of(x,cgsl_rational))
-#endif
-
-#ifndef CHECK_RATIONAL
-#define CHECK_RATIONAL(x) if(!rb_obj_is_kind_of(x,cgsl_rational)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Rational expected)");
-#endif
-
-/*****/
-#ifndef BLOCK_P
-#define BLOCK_P(x) (rb_obj_is_kind_of(x,cgsl_block))
-#endif
-
-#ifndef CHECK_BLOCK
-#define CHECK_BLOCK(x) if(!rb_obj_is_kind_of(x,cgsl_block)) \
-    rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block expected)", rb_class2name(CLASS_OF(x)));
-#endif
-#ifndef BLOCK_INT_P
-#define BLOCK_INT_P(x) (rb_obj_is_kind_of(x,cgsl_block_int))
-#endif
-
-#ifndef CHECK_BLOCK_INT
-#define CHECK_BLOCK_INT(x) if(!rb_obj_is_kind_of(x,cgsl_block_int)) \
-    rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block::Int expected)", rb_class2name(CLASS_OF(x)));
-#endif
-#ifndef BLOCK_UCHAR_P
-#define BLOCK_UCHAR_P(x) (rb_obj_is_kind_of(x,cgsl_block_uchar))
-#endif
-
-#ifndef CHECK_BLOCK_UCHAR
-#define CHECK_BLOCK_UCHAR(x) if(!rb_obj_is_kind_of(x,cgsl_block_uchar)) \
-    rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Block::Byte expected)", rb_class2name(CLASS_OF(x)));
-#endif
-/*****/
-
-#ifndef VECTOR_P
-#define VECTOR_P(x) (rb_obj_is_kind_of(x,cgsl_vector))
-#endif
-
-#ifndef VECTOR_VIEW_P
-#define VECTOR_VIEW_P(x) ((CLASS_OF(x)==cgsl_vector_view||CLASS_OF(x)==cgsl_vector_col_view||CLASS_OF(x)==cgsl_vector_view_ro||CLASS_OF(x)==cgsl_vector_col_view_ro))
-#endif
-
-#ifndef VECTOR_ROW_P
-#define VECTOR_ROW_P(x) ((CLASS_OF(x)==cgsl_vector||CLASS_OF(x)==cgsl_vector_view||CLASS_OF(x)==cgsl_vector_view_ro))
-#endif
-
-#ifndef VECTOR_COL_P
-#define VECTOR_COL_P(x) ((CLASS_OF(x)==cgsl_vector_col||CLASS_OF(x)==cgsl_vector_col_view||CLASS_OF(x)==cgsl_vector_col_view_ro))
-#endif
-
-#ifndef VECTOR_ROW_COL
-#define VECTOR_ROW_COL(x) ((rb_obj_is_kind_of(x,cgsl_vector_col)||rb_obj_is_kind_of(x,cgsl_vector_int_col)) ? cgsl_vector_col : cgsl_vector)
-#endif
-
-#ifndef CHECK_VECTOR
-#define CHECK_VECTOR(x) if(!rb_obj_is_kind_of(x,cgsl_vector)) \
-    rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Vector expected)", rb_class2name(CLASS_OF(x)));
-#endif
-
-#ifdef HAVE_NARRAY_H
-#define Data_Get_Vector(obj,sval) do { \
-    if (NA_IsNArray(obj)) { \
-      /* Convert obj to GSL::Vector::View */ \
-      obj = rb_gsl_na_to_gsl_vector_view_method(obj); \
-    } \
-    CHECK_VECTOR(obj); \
-    Data_Get_Struct(obj,gsl_vector,sval); \
-} while (0)
-#else
-#define Data_Get_Vector(obj,sval) do { \
-    CHECK_VECTOR(obj); \
-    Data_Get_Struct(obj,gsl_vector,sval); \
-} while (0)
-#endif
-
-/******/
-#ifndef VECTOR_INT_P
-#define VECTOR_INT_P(x) (rb_obj_is_kind_of(x,cgsl_vector_int))
-#endif
-
-#ifndef VECTOR_INT_VIEW_P
-#define VECTOR_INT_VIEW_P(x) ((CLASS_OF(x)==cgsl_vector_int_view||CLASS_OF(x)==cgsl_vector_int_col_view||CLASS_OF(x)==cgsl_vector_int_view_ro||CLASS_OF(x)==cgsl_vector_int_col_view_ro))
-#endif
-
-#ifndef VECTOR_INT_ROW_P
-#define VECTOR_INT_ROW_P(x) ((CLASS_OF(x)==cgsl_vector_int||CLASS_OF(x)==cgsl_vector_int_view||CLASS_OF(x)==cgsl_vector_int_view_ro))
-#endif
-
-#ifndef VECTOR_INT_COL_P
-#define VECTOR_INT_COL_P(x) ((CLASS_OF(x)==cgsl_vector_int_col||CLASS_OF(x)==cgsl_vector_int_col_view||CLASS_OF(x)==cgsl_vector_int_col_view_ro))
-#endif
-
-#ifndef VECTOR_INT_ROW_COL
-#define VECTOR_INT_ROW_COL(x) (VECTOR_INT_ROW_P(x) ? cgsl_vector_int : cgsl_vector_int_col)
-#endif
-
-#ifndef CHECK_VECTOR_INT
-#define CHECK_VECTOR_INT(x) if(!rb_obj_is_kind_of(x,cgsl_vector_int)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Vector::Int expected)");
-#endif
-
-/******/
-#ifndef VECTOR_COMPLEX_P
-#define VECTOR_COMPLEX_P(x) (rb_obj_is_kind_of(x,cgsl_vector_complex))
-#endif
-
-#ifndef VECTOR_COMPLEX_ROW_P
-#define VECTOR_COMPLEX_ROW_P(x) ((CLASS_OF(x)==cgsl_vector_complex||CLASS_OF(x)==cgsl_vector_complex_view))
-#endif
-
-#ifndef VECTOR_COMPLEX_COL_P
-#define VECTOR_COMPLEX_COL_P(x) ((CLASS_OF(x)==cgsl_vector_complex_col||CLASS_OF(x)==cgsl_vector_complex_col_view))
-#endif
-
-#ifndef VECTOR_COMPLEX_ROW_COL
-#define VECTOR_COMPLEX_ROW_COL(x) (VECTOR_COMPLEX_ROW_P(x) ? cgsl_vector_complex : cgsl_vector_complex_col)
-#endif
-
-#ifndef CHECK_VECTOR_COMPLEX
-#define CHECK_VECTOR_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_vector_complex)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Vector::Complex expected)");
-#endif
-
-#ifndef MATRIX_P
-#define MATRIX_P(x) (rb_obj_is_kind_of(x,cgsl_matrix))
-#endif
-
-#ifndef CHECK_MATRIX
-#define CHECK_MATRIX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix expected)");
-#endif
-
-#define Data_Get_Matrix(obj,sval) do { \
-    CHECK_MATRIX(obj); \
-    Data_Get_Struct(obj,gsl_matrix,sval); \
-} while (0)
-
-
-#ifndef MATRIX_INT_P
-#define MATRIX_INT_P(x) (rb_obj_is_kind_of(x,cgsl_matrix_int))
-#endif
-
-#ifndef CHECK_MATRIX_INT
-#define CHECK_MATRIX_INT(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_int)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix::Int expected)");
-#endif
-
-#ifndef MATRIX_COMPLEX_P
-#define MATRIX_COMPLEX_P(x) (rb_obj_is_kind_of(x,cgsl_matrix_complex))
-#endif
-
-#ifndef CHECK_MATRIX_COMPLEX
-#define CHECK_MATRIX_COMPLEX(x) if(!rb_obj_is_kind_of(x,cgsl_matrix_complex)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Matrix::Complex expected)");
-#endif
-
-#ifndef TENSOR_P
-#define TENSOR_P(x) ((CLASS_OF(x)==cgsl_tensor))
-#endif
-
-#ifndef CHECK_TENSOR
-#define CHECK_TENSOR(x) if(CLASS_OF(x)!=cgsl_tensor) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Tensor expected)");
-#endif
-
-#ifndef TENSOR_INT_P
-#define TENSOR_INT_P(x) ((CLASS_OF(x)==cgsl_tensor_int))
-#endif
-
-#ifndef CHECK_TENSOR_INT
-#define CHECK_TENSOR_INT(x) if(CLASS_OF(x)!=cgsl_tensor_int) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Tensor::Int expected)");
-#endif
-
-#ifndef PERMUTATION_P
-#define PERMUTATION_P(x) (rb_obj_is_kind_of(x,cgsl_permutation))
-#endif
-
-#ifndef CHECK_PERMUTATION
-#define CHECK_PERMUTATION(x) if(!rb_obj_is_kind_of(x,cgsl_permutation)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Permutation expected)");
-#endif
-
-#ifndef PROC_P
-#define PROC_P(x) (rb_obj_is_kind_of(x,rb_cProc))
-#endif
-
-#ifndef CHECK_PROC
-#define CHECK_PROC(x) if(!rb_obj_is_kind_of(x,rb_cProc)) \
-    rb_raise(rb_eTypeError, "wrong argument type (Proc expected)");
-#endif
-
-#ifndef FUNCTION_P
-#define FUNCTION_P(x) (rb_obj_is_kind_of(x,cgsl_function))
-#endif
-
-#ifndef CHECK_FUNCTION
-#define CHECK_FUNCTION(x) if(!rb_obj_is_kind_of(x,cgsl_function)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Function expected)");
-#endif
-
-#ifndef FUNCTION_FDF_P
-#define FUNCTION_FDF_P(x) (rb_obj_is_kind_of(x,cgsl_function_fdf))
-#endif
-
-#ifndef CHECK_FUNCTION_FDF
-#define CHECK_FUNCTION_FDF(x) if(!rb_obj_is_kind_of(x,cgsl_function_fdf)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Function_fdf expected)");
-#endif
-
-#ifndef HISTOGRAM_P
-#define HISTOGRAM_P(x) (rb_obj_is_kind_of(x,cgsl_histogram))
-#endif
-
-#ifndef CHECK_HISTOGRAM
-#define CHECK_HISTOGRAM(x) if(!rb_obj_is_kind_of(x,cgsl_histogram)) \
-    rb_raise(rb_eTypeError, "wrong argument type (GSL::Histogram expected)");
-#endif
-
-#ifndef RBGSL_SET_CLASS
-#ifdef RB_OBJ_WRITE
-#define RBGSL_SET_CLASS0(obj0, cls) RB_OBJ_WRITE(obj0, &(RBASIC_CLASS(obj0)), cls)
-#else
-#define RBGSL_SET_CLASS0(obj0, cls) RBASIC(obj0)->klass = cls
-#endif
-#define RBGSL_SET_CLASS(obj, cls) do { \
-    VALUE _obj_ = (obj); \
-    RBGSL_SET_CLASS0(_obj_, cls); \
-} while (0)
-#endif
-
-void rb_gsl_error_handler(const char *reason, const char *file,
-                          int line, int gsl_errno);
-
-FILE* rb_gsl_open_writefile(VALUE io, int *flag);
-FILE* rb_gsl_open_readfile(VALUE io, int *flag);
-
-VALUE rb_gsl_obj_read_only(int argc, VALUE *argv, VALUE obj);
-
-int str_tail_grep(const char *s0, const char *s1);
-int str_head_grep(const char *s0, const char *s1);
-
-#ifndef STR2CSTR
-#define STR2CSTR StringValuePtr
-#endif
-
-#ifndef RCLASS_SUPER
-#define RCLASS_SUPER(cls) RCLASS(cls)->super
-#endif
-
-void make_graphcommand(char *command, VALUE hash);
-int rbgsl_complex_equal(const gsl_complex *z1, const gsl_complex *z2, double eps);
-
-gsl_vector* mygsl_vector_down(gsl_vector *p);
-void mygsl_vector_up2(gsl_vector *pnew, gsl_vector *p);
-gsl_vector* mygsl_vector_up(gsl_vector *p);
-
-gsl_vector_int* mygsl_vector_int_down(gsl_vector_int *p);
-void mygsl_vector_int_up2(gsl_vector_int *pnew, gsl_vector_int *p);
-gsl_vector_int* mygsl_vector_int_up(gsl_vector_int *p);
-
-size_t count_columns(const char *str);
-char* str_scan_double(const char *str, double *val);
-char* str_scan_int(const char *str, int *val);
-double* get_ptr_double3(VALUE obj, size_t *size, size_t *stride, int *flag);
-gsl_complex ary2complex(VALUE obj);
-VALUE vector_eval_create(VALUE obj, double (*func)(double));
-VALUE matrix_eval_create(VALUE obj, double (*func)(double));
-VALUE rb_gsl_ary_eval1(VALUE ary, double (*f)(double));
-#ifdef HAVE_NARRAY_H
-VALUE rb_gsl_nary_eval1(VALUE ary, double (*f)(double));
-#endif
-
-EXTERN VALUE cGSL_Object;
-#endif