rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/histogram_find.c

@@ -1,99 +0,0 @@
-/* histogram/find.c
- *
- * Copyright (C) 1996, 1997, 1998, 1999, 2000 Brian Gough
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or (at
- * your option) any later version.
- *
- * This program 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.  See the GNU
- * General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
- */
-
-#include "include/rb_gsl_histogram.h"
-/* determines whether to optimize for linear ranges  */
-#define LINEAR_OPT 1
-
-int mygsl_find (const size_t n, const double range[], const double x, size_t * i)
-{
-  size_t i_linear, lower, upper, mid;
-  if (x < range[0]) return -1;
-  if (x >= range[n]) return +1;
-  /* optimize for linear case */
-
-#ifdef LINEAR_OPT
-  {
-    double u =  (x - range[0]) / (range[n] - range[0]);
-    i_linear = (size_t) (n*u);
-  }
-
-  if (x >= range[i_linear] && x < range[i_linear + 1])
-  {
-    *i = i_linear;
-    return 0;
-  }
-#endif
-
-  /* perform binary search */
-
-  upper = n;
-  lower = 0;
-
-  while (upper - lower > 1)
-  {
-    mid = (upper + lower) / 2;
-
-    if (x >= range[mid])
-    {
-      lower = mid;
-    }
-    else
-    {
-      upper = mid;
-    }
-  }
-
-  *i = lower;
-
-  /* sanity check the result */
-
-  if (x < range[lower] || x >= range[lower + 1])
-  {
-    GSL_ERROR ("x not found in range", GSL_ESANITY);
-  }
-  return 0;
-}
-
-int mygsl_find2d (const size_t nx, const double xrange[],
-                  const size_t ny, const double yrange[],
-                  const double x, const double y,
-                  size_t * i, size_t * j)
-{
-  int status = mygsl_find (nx, xrange, x, i);
-  if (status) return status;
-  status = mygsl_find (ny, yrange, y, j);
-  if (status) return status;
-  return 0;
-}
-
-int mygsl_find3d (const size_t nx, const double xrange[],
-                  const size_t ny, const double yrange[],
-                  const size_t nz, const double zrange[],
-                  const double x, const double y, const double z,
-                  size_t * i, size_t * j, size_t *k)
-{
-  int status = mygsl_find (nx, xrange, x, i);
-  if (status) return status;
-  status = mygsl_find (ny, yrange, y, j);
-  if (status) return status;
-  status = mygsl_find (nz, zrange, z, k);
-  if (status) return status;
-  return 0;
-}

data/ext/gsl_native/histogram_oper.c

@@ -1,150 +0,0 @@
-/* gsl_histogram_oper.c
- * Copyright (C) 2000  Simone Piccardi
- *
- * This library is free software; you can redistribute it and/or
- * modify it under the terms of the GNU General Public License as
- * published by the Free Software Foundation; either version 2 of the
- * License, or (at your option) any later version.
- *
- * This program 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.  See the GNU
- * General Public License for more details.
- *
- * You should have received a copy of the GNU General Public
- * License along with this library; if not, write to the
- * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
- * Boston, MA 02111-1307, USA.
- */
-/***************************************************************
- *
- * File gsl_histogram_oper.c:
- * Routine to make operation on histograms.
- * Need GSL library and header.
- * Contains the routines:
- * gsl_histogram_same_binning check if two histograms have the same binning
- * gsl_histogram_add          add two histograms
- * gsl_histogram_sub          subctract two histograms
- * gsl_histogram_mult         multiply two histograms
- * gsl_histogram_div          divide two histograms
- * gsl_histogram_scale        scale histogram contents
- *
- * Author: S. Piccardi
- * Jan. 2000
- *
- ***************************************************************/
-#include <stdlib.h>
-#include <gsl/gsl_errno.h>
-#include <gsl/gsl_math.h>
-#include <gsl/gsl_histogram.h>
-
-/*
- * gsl_histogram_same_binning:
- * control if two histograms have the
- * same binning
- */
-
-int
-mygsl_histogram_equal_bins_p (const gsl_histogram * h1, const gsl_histogram * h2)
-{
-  if (h1->n != h2->n)
-  {
-    return 0;
-  }
-  {
-    size_t i;
-    /* init ranges */
-
-    for (i = 0; i <= h1->n; i++)
-    {
-      if (gsl_fcmp(h1->range[i],h2->range[i], 1e-12))
-      {
-        return 0;
-      }
-    }
-  }
-
-  return 1;
-}
-
-/*
- * gsl_histogram_add:
- * add two histograms
- */
-int
-mygsl_histogram_add (gsl_histogram * h1, const gsl_histogram * h2)
-{
-  size_t i;
-
-  if (!mygsl_histogram_equal_bins_p (h1, h2))
-  {
-    GSL_ERROR ("histograms have different binning", GSL_EINVAL);
-  }
-  for (i = 0; i < h1->n; i++)
-  {
-    h1->bin[i] += h2->bin[i];
-  }
-  return GSL_SUCCESS;
-}
-
-/*
- * gsl_histogram_sub:
- * subtract two histograms
- */
-
-int
-mygsl_histogram_sub (gsl_histogram * h1, const gsl_histogram * h2)
-{
-  size_t i;
-
-  if (!mygsl_histogram_equal_bins_p (h1, h2))
-  {
-    GSL_ERROR ("histograms have different binning", GSL_EINVAL);
-  }
-  for (i = 0; i < h1->n; i++)
-  {
-    h1->bin[i] -= h2->bin[i];
-  }
-  return GSL_SUCCESS;
-
-}
-
-/*
- * gsl_histogram_mult:
- * multiply two histograms
- */
-
-int
-mygsl_histogram_mul (gsl_histogram * h1, const gsl_histogram * h2)
-{
-  size_t i;
-
-  if (!mygsl_histogram_equal_bins_p (h1, h2))
-  {
-    GSL_ERROR ("histograms have different binning", GSL_EINVAL);
-  }
-  for (i = 0; i < h1->n; i++)
-  {
-    h1->bin[i] *= h2->bin[i];
-  }
-  return GSL_SUCCESS;
-}
-/*
- * gsl_histogram_div:
- * divide two histograms
- */
-int
-mygsl_histogram_div (gsl_histogram * h1, const gsl_histogram * h2)
-{
-  size_t i;
-
-  if (!mygsl_histogram_equal_bins_p (h1, h2))
-  {
-    GSL_ERROR ("histograms have different binning", GSL_EINVAL);
-  }
-  for (i = 0; i < h1->n; i++)
-  {
-    h1->bin[i] /= h2->bin[i];
-  }
-  return GSL_SUCCESS;
-}

data/ext/gsl_native/ieee.c

@@ -1,88 +0,0 @@
-/*
-  ieee.c
-  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
-    (C) Copyright 2001-2006 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.
-*/
-
-#include "include/rb_gsl.h"
-
-static VALUE rb_gsl_ieee_env_setup(VALUE obj)
-{
-  gsl_ieee_env_setup();
-  return obj;
-}
-
-static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
-{
-  rb_io_t *fptr = NULL;
-  FILE *fp = NULL;
-  int flag = 0;
-  VALUE vtmp;
-  double ftmp;
-
-  switch (argc) {
-  case 2:
-    switch (TYPE(argv[0])) {
-    case T_STRING:
-      fp = fopen(RSTRING_PTR(argv[0]), "w");
-      flag = 1;
-      break;
-    case T_FILE:
-      GetOpenFile(argv[0], fptr);
-      rb_io_check_writable(fptr);
-      fp = rb_io_stdio_file(fptr);
-      break;
-    default:
-      rb_raise(rb_eTypeError, "wrong type argument %s (IO or String expected)",
-               rb_class2name(CLASS_OF(argv[0])));
-    }
-    vtmp = argv[1];
-    break;
-  case 1:
-    vtmp = argv[0];
-    fp = stdout;
-    break;
-  default:
-    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
-  }
-  if (TYPE(vtmp) != T_FLOAT)
-    rb_raise(rb_eTypeError, "wrong argument type %s (Float expected)",
-             rb_class2name(CLASS_OF(vtmp)));
-  ftmp = RFLOAT_VALUE(vtmp);
-  gsl_ieee_fprintf_double(fp, &ftmp);
-  if (fp == stdout) fprintf(stdout, "\n");
-  if (flag == 1) fclose(fp);
-  return obj;
-}
-
-static VALUE rb_gsl_ieee_printf_double(VALUE obj, VALUE xx)
-{
-  double x;
-  x = NUM2DBL(xx);
-  gsl_ieee_printf_double(&x);
-  return xx;
-}
-
-void Init_gsl_ieee(VALUE module)
-{
-  VALUE mgsl_ieee;
-  mgsl_ieee = rb_define_module_under(module, "IEEE");
-
-  rb_define_singleton_method(mgsl_ieee, "env_setup",
-                             rb_gsl_ieee_env_setup, 0);
-  rb_define_module_function(module, "ieee_env_setup", rb_gsl_ieee_env_setup, 0);
-  rb_define_singleton_method(mgsl_ieee, "fprintf_double",
-                             rb_gsl_ieee_fprintf_double, -1);
-  rb_define_singleton_method(mgsl_ieee, "fprintf",
-                             rb_gsl_ieee_fprintf_double, -1);
-  rb_define_singleton_method(mgsl_ieee, "printf",
-                             rb_gsl_ieee_printf_double, -1);
-  rb_define_singleton_method(mgsl_ieee, "printf_double",
-                             rb_gsl_ieee_printf_double, -1);
-
-}

data/ext/gsl_native/include/rb_gsl.h

@@ -1,136 +0,0 @@
-/*
-  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
-*/
-
-#ifndef ___RB_GSL_H___
-#define ___RB_GSL_H___
-
-#include <string.h>
-#include <stdlib.h>
-#include <gsl/gsl_mode.h>
-#include "rb_gsl_poly.h"
-#include "rb_gsl_linalg.h"
-#include "rb_gsl_math.h"
-#include "rb_gsl_root.h"
-#include "rb_gsl_histogram.h"
-#include "rb_gsl_sf.h"
-#include "rb_gsl_eigen.h"
-#include "rb_gsl_fft.h"
-#include "rb_gsl_integration.h"
-#include "rb_gsl_rng.h"
-#include "rb_gsl_statistics.h"
-#include "rb_gsl_odeiv.h"
-#include "rb_gsl_interp.h"
-#include "rb_gsl_cheb.h"
-#include "rb_gsl_fit.h"
-#include "rb_gsl_const.h"
-
-void Init_gsl_error(VALUE module);
-void Init_gsl_math(VALUE module);
-void Init_gsl_complex(VALUE module);
-void Init_gsl_array(VALUE module);
-void Init_gsl_blas(VALUE module);
-void Init_gsl_sort(VALUE module);
-void Init_gsl_poly(VALUE module);
-void Init_gsl_poly_int(VALUE module);
-void Init_gsl_poly2(VALUE module);
-void Init_gsl_rational(VALUE module);
-void Init_gsl_sf(VALUE module);
-void Init_gsl_linalg(VALUE module);
-void Init_gsl_eigen(VALUE module);
-void Init_gsl_fft(VALUE module);
-void Init_gsl_signal(VALUE module);
-void Init_gsl_function(VALUE module);
-void Init_gsl_integration(VALUE module);
-
-void Init_gsl_rng(VALUE module);
-void Init_gsl_qrng(VALUE module);
-void Init_gsl_ran(VALUE module);
-void Init_gsl_cdf(VALUE module);
-void Init_gsl_stats(VALUE module);
-
-void Init_gsl_histogram(VALUE module);
-void Init_gsl_histogram2d(VALUE module);
-void Init_gsl_histogram3d(VALUE module);
-void Init_gsl_ntuple(VALUE module);
-void Init_gsl_monte(VALUE module);
-void Init_gsl_siman(VALUE module);
-
-void Init_gsl_odeiv(VALUE module);
-void Init_gsl_interp(VALUE module);
-void Init_gsl_spline(VALUE module);
-void Init_gsl_diff(VALUE module);
-void Init_gsl_deriv(VALUE module);
-
-void Init_gsl_cheb(VALUE module);
-void Init_gsl_sum(VALUE module);
-void Init_gsl_dht(VALUE module);
-
-void Init_gsl_root(VALUE module);
-void Init_gsl_multiroot(VALUE module);
-void Init_gsl_min(VALUE module);
-void Init_gsl_multimin(VALUE module);
-void Init_gsl_fit(VALUE module);
-void Init_gsl_multifit(VALUE module);
-
-void Init_gsl_const(VALUE module);
-
-void Init_gsl_ieee(VALUE module);
-
-#ifdef HAVE_NARRAY_H
-void Init_gsl_narray(VALUE module);
-#endif
-
-void Init_wavelet(VALUE module);
-
-void Init_gsl_graph(VALUE module);
-
-#ifdef HAVE_TENSOR_TENSOR_H
-void Init_tensor_init(VALUE module);
-void Init_tensor_int_init(VALUE module);
-#endif
-
-void Init_gsl_dirac(VALUE module);
-
-EXTERN VALUE cGSL_Object;
-
-void Init_tamu_anova(VALUE module);
-
-#ifdef HAVE_TAMU_ANOVA_TAMU_ANOVA_H
-#include "tamu_anova/tamu_anova.h"
-#endif
-
-#ifdef HAVE_OOL_OOL_VERSION_H
-void Init_ool(VALUE module);
-#endif
-
-#ifdef HAVE_JACOBI_H
-void Init_jacobi(VALUE module);
-#endif
-
-#ifdef HAVE_GSL_GSL_CQP_H
-void Init_cqp(VALUE module);
-#endif
-
-void Init_fresnel(VALUE module);
-
-void Init_bspline(VALUE module);
-
-#ifdef HAVE_ALF_ALF_H
-#include "alf/alf.h"
-#endif
-void Init_alf(VALUE module);
-void Init_geometry(VALUE module);
-
-#include <gsl/gsl_multiset.h>
-extern VALUE cMultiset;
-void Init_multiset(VALUE module);
-
-#endif