rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/multiset.c

@@ -1,210 +0,0 @@
-#include "include/rb_gsl.h"
-
-VALUE cMultiset;
-
-VALUE rb_gsl_multiset_alloc(VALUE klass, VALUE nn, VALUE kk)
-{
-
-  gsl_multiset *m;
-  m = gsl_multiset_alloc(FIX2INT(nn), FIX2INT(kk));
-  return Data_Wrap_Struct(klass, 0, gsl_multiset_free, m);
-}
-
-VALUE rb_gsl_multiset_calloc(VALUE klass, VALUE nn, VALUE kk)
-{
-  gsl_multiset *m;
-  m = gsl_multiset_alloc(FIX2INT(nn), FIX2INT(kk));
-  return Data_Wrap_Struct(klass, 0, gsl_multiset_free, m);
-}
-
-VALUE rb_gsl_multiset_init_first(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  gsl_multiset_init_first(m);
-  return Qnil;
-}
-
-VALUE rb_gsl_multiset_init_last(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  gsl_multiset_init_last(m);
-  return Qnil;
-}
-
-VALUE rb_gsl_multiset_memcpy(VALUE klass, VALUE m1, VALUE m2)
-{
-  gsl_multiset *mm1, *mm2;
-  if (!rb_obj_is_kind_of(m1, klass)) {
-    rb_raise(rb_eTypeError, "Wrong type %s (GSL::Multiset expected)", rb_class2name(CLASS_OF(m1)));
-  }
-  if (!rb_obj_is_kind_of(m2, klass)) {
-    rb_raise(rb_eTypeError, "Wrong type %s (GSL::Multiset expected)", rb_class2name(CLASS_OF(m2)));
-  }
-  Data_Get_Struct(m1, gsl_multiset, mm1);
-  Data_Get_Struct(m2, gsl_multiset, mm2);
-  return FIX2INT(gsl_multiset_memcpy(mm1, mm2));
-}
-
-VALUE rb_gsl_multiset_get(VALUE mm, VALUE i)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_get(m, FIX2INT(i)));
-}
-
-VALUE rb_gsl_multiset_n(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_n(m));
-}
-
-VALUE rb_gsl_multiset_k(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_k(m));
-}
-
-VALUE rb_gsl_multiset_valid(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_valid(m));
-}
-
-VALUE rb_gsl_multiset_valid2(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  if (gsl_multiset_valid(m)) return Qtrue;
-  else return Qfalse;
-}
-
-VALUE rb_gsl_multiset_next(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_next(m));
-}
-
-VALUE rb_gsl_multiset_prev(VALUE mm)
-{
-  gsl_multiset *m;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  return INT2FIX(gsl_multiset_prev(m));
-}
-
-VALUE rb_gsl_multiset_fwrite(VALUE mm, VALUE name)
-{
-  gsl_multiset *m;
-  FILE *fp = NULL;
-  int ret;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  fp = fopen(STR2CSTR(name), "wb");
-  if (fp == NULL) {
-    rb_raise(rb_eIOError, "Cannot open file %s", STR2CSTR(name));
-  }
-  ret = gsl_multiset_fwrite(fp, m);
-  fclose(fp);
-  return INT2FIX(ret);
-}
-
-VALUE rb_gsl_multiset_fread(VALUE mm, VALUE name)
-{
-  gsl_multiset *m;
-  FILE *fp = NULL;
-  int ret;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  fp = fopen(STR2CSTR(name), "wb");
-  if (fp == NULL) {
-    rb_raise(rb_eIOError, "Cannot open file %s", STR2CSTR(name));
-  }
-  ret = gsl_multiset_fread(fp, m);
-  fclose(fp);
-  return INT2FIX(ret);
-}
-
-VALUE rb_gsl_multiset_fprintf(VALUE mm, VALUE name, VALUE format)
-{
-  gsl_multiset *m;
-  FILE *fp = NULL;
-  int ret;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  fp = fopen(STR2CSTR(name), "w");
-  if (fp == NULL) {
-    rb_raise(rb_eIOError, "Cannot open file %s", STR2CSTR(name));
-  }
-  ret = gsl_multiset_fprintf(fp, m, STR2CSTR(format));
-  fclose(fp);
-  return INT2FIX(ret);
-}
-
-VALUE rb_gsl_multiset_fscanf(VALUE mm, VALUE name)
-{
-  gsl_multiset *m;
-  FILE *fp = NULL;
-  int ret;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  fp = fopen(STR2CSTR(name), "r");
-  if (fp == NULL) {
-    rb_raise(rb_eIOError, "Cannot open file %s", STR2CSTR(name));
-  }
-  ret = gsl_multiset_fscanf(fp, m);
-  fclose(fp);
-  return INT2FIX(ret);
-}
-
-VALUE rb_gsl_multiset_data(VALUE mm)
-{
-  gsl_multiset *m;
-  size_t *p;
-  gsl_vector_int *v;
-  size_t i;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  p = gsl_multiset_data(m);
-  v = gsl_vector_int_alloc(m->k);
-  for (i = 0; i < v->size; i++) gsl_vector_int_set(v, i, p[i]);
-  return Data_Wrap_Struct(cgsl_vector_int, 0, gsl_vector_int_free, v);
-}
-
-VALUE rb_gsl_multiset_data2(VALUE mm, VALUE i)
-{
-  gsl_multiset *m;
-  size_t *p;
-  Data_Get_Struct(mm, gsl_multiset, m);
-  p = gsl_multiset_data(m);
-  return INT2FIX(p[i]);
-}
-
-void Init_multiset(VALUE module)
-{
-  cMultiset = rb_define_class_under(module, "Multiset", cGSL_Object);
-  rb_define_singleton_method(cMultiset, "alloc", rb_gsl_multiset_alloc, 2);
-  rb_define_singleton_method(cMultiset, "calloc", rb_gsl_multiset_calloc, 2);
-  rb_define_singleton_method(cMultiset, "memcpy", rb_gsl_multiset_memcpy, 2);
-
-  rb_define_method(cMultiset, "init_first", rb_gsl_multiset_init_first, 0);
-  rb_define_method(cMultiset, "init_last", rb_gsl_multiset_init_last, 0);
-
-  rb_define_method(cMultiset, "get", rb_gsl_multiset_get, 1);
-  rb_define_alias(cMultiset, "[]", "get");
-
-  rb_define_method(cMultiset, "n", rb_gsl_multiset_n, 0);
-  rb_define_method(cMultiset, "k", rb_gsl_multiset_k, 0);
-  rb_define_method(cMultiset, "data", rb_gsl_multiset_data, 0);
-  rb_define_method(cMultiset, "data[]", rb_gsl_multiset_data2, 1);
-  rb_define_method(cMultiset, "valid", rb_gsl_multiset_valid, 0);
-  rb_define_method(cMultiset, "valid?", rb_gsl_multiset_valid2, 0);
-
-  rb_define_method(cMultiset, "next", rb_gsl_multiset_next, 0);
-  rb_define_method(cMultiset, "prev", rb_gsl_multiset_prev, 0);
-
-  rb_define_method(cMultiset, "fwrite", rb_gsl_multiset_fwrite, 1);
-  rb_define_method(cMultiset, "fread", rb_gsl_multiset_fread, 1);
-  rb_define_method(cMultiset, "fprintf", rb_gsl_multiset_fprintf, 2);
-  rb_define_method(cMultiset, "fscanf", rb_gsl_multiset_fscanf, 1);
-}
-

data/ext/gsl_native/ndlinear.c

@@ -1,320 +0,0 @@
-#include "include/rb_gsl.h"
-
-#ifdef HAVE_NDLINEAR_GSL_MULTIFIT_NDLINEAR_H
-#include <gsl/gsl_math.h>
-#include <gsl/gsl_matrix.h>
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_multifit.h>
-#include <ndlinear/gsl_multifit_ndlinear.h>
-
-static VALUE cWorkspace;
-
-enum Index_Ndlinear {
-  INDEX_NDIM = 0,
-  INDEX_N = 1,
-  INDEX_PROCS = 2,
-  INDEX_PARAMS = 3,
-  INDEX_FUNCS = 4,
-  INDEX_NDIM_I = 5,
-
-  NDLINEAR_ARY_SIZE = 6,
-};
-
-static void multifit_ndlinear_mark(gsl_multifit_ndlinear_workspace *w)
-{
-  rb_gc_mark((VALUE) w->params);
-}
-
-typedef int (*UFUNC)(double, double[], void*);
-typedef struct ufunc_struct
-{
-  UFUNC *fptr;
-} ufunc_struct;
-
-static VALUE cUFunc;
-static ufunc_struct* ufunc_struct_alloc(size_t n_dim) {
-  ufunc_struct *p;
-  p = (ufunc_struct*) malloc(sizeof(ufunc_struct));
-  p->fptr = malloc(sizeof(UFUNC)*n_dim);
-  return p;
-}
-static void ufunc_struct_free(ufunc_struct *p)
-{
-  free(p->fptr);
-  free(p);
-}
-
-static int func_u(double x, double y[], void *data);
-static VALUE rb_gsl_multifit_ndlinear_alloc(int argc, VALUE *argv, VALUE klass)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  int istart = 0;
-  size_t n_dim = 0, *N, i;
-  struct ufunc_struct *p;
-  VALUE params, wspace, pp;
-  switch (argc) {
-  case 4:
-    istart = 1;
-    CHECK_FIXNUM(argv[0]);
-    n_dim = FIX2INT(argv[0]);
-  /* no break */
-  case 3:
-    if (TYPE(argv[istart]) != T_ARRAY) {
-      rb_raise(rb_eTypeError, "Wrong argument type %s (Array expected)",
-               rb_class2name(CLASS_OF(argv[istart])));
-    }
-    if (TYPE(argv[istart+1]) != T_ARRAY) {
-      rb_raise(rb_eTypeError, "Wrong argument type %s (Array expected)",
-               rb_class2name(CLASS_OF(argv[istart+1])));
-    }
-    //    n_dim = RARRAY(argv[istart])->len;
-    n_dim = RARRAY_LEN(argv[istart]);
-    N = (size_t*) malloc(sizeof(size_t)*n_dim);
-    break;
-  default:
-    rb_raise(rb_eArgError, "Wrong number of arguments (%d for 3 or 4)", argc);
-  }
-  for (i = 0; i < n_dim; i++) {
-    N[i] = FIX2INT(rb_ary_entry(argv[istart], i));
-  }
-  params = rb_ary_new2(NDLINEAR_ARY_SIZE);
-  rb_ary_store(params, INDEX_NDIM, INT2FIX((int) n_dim));
-  rb_ary_store(params, INDEX_N, argv[istart]);   /* N */
-  rb_ary_store(params, INDEX_PROCS, argv[istart+1]); /* procs */
-  rb_ary_store(params, INDEX_PARAMS, argv[istart+2]); /* params */
-  rb_ary_store(params, INDEX_NDIM_I, INT2FIX(0)); /* for the first parameter */
-
-  p = ufunc_struct_alloc(n_dim);
-  for (i = 0; i < n_dim; i++) p->fptr[i] = func_u;
-  pp = Data_Wrap_Struct(cUFunc, 0, ufunc_struct_free, p);
-  rb_ary_store(params, INDEX_FUNCS, pp);
-
-  w = gsl_multifit_ndlinear_alloc(n_dim, N, p->fptr, (void*) params);
-
-  free((size_t*) N);
-
-  wspace = Data_Wrap_Struct(cWorkspace, multifit_ndlinear_mark, gsl_multifit_ndlinear_free, w);
-
-  return wspace;
-}
-
-static int func_u(double x, double y[], void *data)
-{
-  VALUE ary, vN, procs, proc, vy, params;
-  gsl_vector_view ytmp;
-  size_t i, n_dim;
-  int rslt;
-  ary = (VALUE) data;
-  n_dim = FIX2INT(rb_ary_entry(ary, INDEX_NDIM));
-  vN = rb_ary_entry(ary, INDEX_N);
-  procs = rb_ary_entry(ary, INDEX_PROCS);
-  params = rb_ary_entry(ary, INDEX_PARAMS);
-  i = FIX2INT(rb_ary_entry(ary, INDEX_NDIM_I));
-  proc = rb_ary_entry(procs, i);
-
-  ytmp.vector.data = (double*) y;
-  ytmp.vector.stride = 1;
-  ytmp.vector.size = FIX2INT(rb_ary_entry(vN, i));
-  vy = Data_Wrap_Struct(cgsl_vector_view, 0, NULL, &ytmp);
-
-  rslt = rb_funcall((VALUE) proc, RBGSL_ID_call, 3, rb_float_new(x), vy, params);
-
-  /* for the next parameter */
-  i += 1;
-  if (i == n_dim) i = 0;
-  rb_ary_store(ary, INDEX_NDIM_I, INT2FIX(i));
-
-  return GSL_SUCCESS;
-}
-
-static VALUE rb_gsl_multifit_ndlinear_design(int argc, VALUE *argv, VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  gsl_matrix *vars = NULL, *X = NULL;
-  int argc2, flag = 0, ret;
-  switch (TYPE(obj)) {
-  case T_MODULE:
-  case T_CLASS:
-  case T_OBJECT:
-    if (!rb_obj_is_kind_of(argv[argc-1], cWorkspace)) {
-      rb_raise(rb_eTypeError, "Wrong argument type %s (GSL::MultiFit::Ndlinear::Workspace expected)",
-               rb_class2name(CLASS_OF(argv[argc-1])));
-    }
-    Data_Get_Struct(argv[argc-1], gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc-1;
-    break;
-  default:
-    Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc;
-  }
-  switch (argc2) {
-  case 1:
-    CHECK_MATRIX(argv[0]);
-    Data_Get_Struct(argv[0], gsl_matrix, vars);
-    X = gsl_matrix_alloc(vars->size1, w->n_coeffs);
-    flag = 1;
-    break;
-  case 2:
-    CHECK_MATRIX(argv[0]);
-    CHECK_MATRIX(argv[1]);
-    Data_Get_Struct(argv[0], gsl_matrix, vars);
-    Data_Get_Struct(argv[1], gsl_matrix, X);
-    break;
-  default:
-    rb_raise(rb_eArgError, "Wrong number of arguments.");
-  }
-  ret = gsl_multifit_ndlinear_design(vars, X, w);
-
-  if (flag == 1) {
-    return Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, X);
-  } else {
-    return INT2FIX(ret);
-  }
-}
-
-static VALUE rb_gsl_multifit_ndlinear_est(int argc, VALUE *argv, VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  gsl_vector *x = NULL, *c = NULL;
-  gsl_matrix *cov = NULL;
-  double y, yerr;
-  int argc2;
-  switch (TYPE(obj)) {
-  case T_MODULE:
-  case T_CLASS:
-  case T_OBJECT:
-    if (!rb_obj_is_kind_of(argv[argc-1], cWorkspace)) {
-      rb_raise(rb_eTypeError, "Wrong argument type %s (GSL::MultiFit::Ndlinear::Workspace expected)",
-               rb_class2name(CLASS_OF(argv[argc-1])));
-    }
-    Data_Get_Struct(argv[argc-1], gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc-1;
-    break;
-  default:
-    Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc;
-  }
-  switch (argc2) {
-  case 3:
-    CHECK_VECTOR(argv[0]);
-    CHECK_VECTOR(argv[1]);
-    CHECK_MATRIX(argv[2]);
-    Data_Get_Struct(argv[0], gsl_vector, x);
-    Data_Get_Struct(argv[1], gsl_vector, c);
-    Data_Get_Struct(argv[2], gsl_matrix, cov);
-    break;
-  default:
-    rb_raise(rb_eArgError, "Wrong number of arguments.");
-  }
-  gsl_multifit_ndlinear_est(x, c, cov, &y, &yerr, w);
-  return rb_ary_new3(2, rb_float_new(y), rb_float_new(yerr));
-}
-
-static VALUE rb_gsl_multifit_ndlinear_calc(int argc, VALUE *argv, VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  gsl_vector *x = NULL, *c = NULL;
-  double val;
-  int argc2;
-  switch (TYPE(obj)) {
-  case T_MODULE:
-  case T_CLASS:
-  case T_OBJECT:
-    if (!rb_obj_is_kind_of(argv[argc-1], cWorkspace)) {
-      rb_raise(rb_eTypeError,
-               "Wrong argument type %s (GSL::MultiFit::Ndlinear::Workspace expected)",
-               rb_class2name(CLASS_OF(argv[argc-1])));
-    }
-    Data_Get_Struct(argv[argc-1], gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc-1;
-    break;
-  default:
-    Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-    argc2 = argc;
-  }
-  switch (argc2) {
-  case 2:
-    CHECK_VECTOR(argv[0]);
-    CHECK_VECTOR(argv[1]);
-    Data_Get_Struct(argv[0], gsl_vector, x);
-    Data_Get_Struct(argv[1], gsl_vector, c);
-    break;
-  default:
-    rb_raise(rb_eArgError, "Wrong number of arguments.");
-  }
-  val = gsl_multifit_ndlinear_calc(x, c, w);
-  return rb_float_new(val);
-}
-
-static VALUE rb_gsl_multifit_ndlinear_n_coeffs(VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-  return INT2FIX(w->n_coeffs);
-}
-
-static VALUE rb_gsl_multifit_ndlinear_n_dim(VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-  return INT2FIX(w->n_dim);
-}
-
-static VALUE rb_gsl_multifit_ndlinear_N(VALUE obj)
-{
-  gsl_multifit_ndlinear_workspace *w;
-  VALUE ary;
-  Data_Get_Struct(obj, gsl_multifit_ndlinear_workspace, w);
-  ary = (VALUE) w->params;
-  return rb_ary_entry(ary, INDEX_N);
-}
-/*
-static VALUE rb_gsl_multifit_linear_Rsq(VALUE module, VALUE vy, VALUE vchisq)
-{
-  gsl_vector *y;
-  double chisq, Rsq;
-  CHECK_VECTOR(vy);
-  Data_Get_Struct(vy, gsl_vector, y);
-  chisq = NUM2DBL(vchisq);
-  gsl_multifit_linear_Rsq(y, chisq, &Rsq);
-  return rb_float_new(Rsq);
-}
-*/
-void Init_ndlinear(VALUE module)
-{
-  VALUE mNdlinear;
-  mNdlinear = rb_define_module_under(module, "Ndlinear");
-  cUFunc = rb_define_class_under(mNdlinear, "UFunc", rb_cObject);
-  cWorkspace = rb_define_class_under(mNdlinear, "Workspace", cGSL_Object);
-
-  rb_define_singleton_method(mNdlinear, "alloc",
-                             rb_gsl_multifit_ndlinear_alloc, -1);
-  rb_define_singleton_method(cWorkspace, "alloc",
-                             rb_gsl_multifit_ndlinear_alloc, -1);
-
-  rb_define_singleton_method(mNdlinear, "design",
-                             rb_gsl_multifit_ndlinear_design, -1);
-  rb_define_singleton_method(cWorkspace, "design",
-                             rb_gsl_multifit_ndlinear_design, -1);
-  rb_define_method(cWorkspace, "design",rb_gsl_multifit_ndlinear_est, -1);
-  rb_define_singleton_method(mNdlinear, "est",
-                             rb_gsl_multifit_ndlinear_est, -1);
-  rb_define_singleton_method(cWorkspace, "est",
-                             rb_gsl_multifit_ndlinear_est, -1);
-  rb_define_method(cWorkspace, "est",rb_gsl_multifit_ndlinear_est, -1);
-
-  rb_define_singleton_method(mNdlinear, "calc",
-                             rb_gsl_multifit_ndlinear_calc, -1);
-  rb_define_singleton_method(cWorkspace, "calc",
-                             rb_gsl_multifit_ndlinear_calc, -1);
-  rb_define_method(cWorkspace, "calc",rb_gsl_multifit_ndlinear_calc, -1);
-
-  rb_define_method(cWorkspace, "n_coeffs",rb_gsl_multifit_ndlinear_n_coeffs, 0);
-  rb_define_method(cWorkspace, "n_dim",rb_gsl_multifit_ndlinear_n_dim, 0);
-  rb_define_method(cWorkspace, "N",rb_gsl_multifit_ndlinear_N, 0);
-
-  //  rb_define_module_function(module, "linear_Rsq", rb_gsl_multifit_linear_Rsq, 2);
-}
-
-#endif
-