numo-libsvm 1.1.2 → 2.1.0

data/ext/numo/libsvm/converter.c

@@ -1,204 +0,0 @@
-
-#include "converter.h"
-
-VALUE int_vec_to_nary(int* const arr, int const size)
-{
-  int i;
-  size_t shape[1] = { size };
-  VALUE v = rb_narray_new(numo_cInt32, 1, shape);
-  int32_t* vp = (int32_t*)na_get_pointer_for_write(v);
-  for (i = 0; i < size; i++) { vp[i] = (int32_t)arr[i]; }
-  return v;
-}
-
-int* nary_to_int_vec(VALUE vec_val)
-{
-  int i;
-  int n_elements;
-  narray_t* vec_nary;
-  int32_t* vec_pt;
-  int* vec;
-
-  if (vec_val == Qnil) return NULL;
-
-  GetNArray(vec_val, vec_nary);
-  n_elements = (int)NA_SHAPE(vec_nary)[0];
-
-  vec = ALLOC_N(int, n_elements);
-  vec_pt = (int32_t*)na_get_pointer_for_read(vec_val);
-  for (i = 0; i < n_elements; i++) { vec[i] = (int)vec_pt[i]; }
-
-  RB_GC_GUARD(vec_val);
-
-  return vec;
-}
-
-VALUE dbl_vec_to_nary(double* const arr, int const size)
-{
-  int i;
-  size_t shape[1] = { size };
-  VALUE v = rb_narray_new(numo_cDFloat, 1, shape);
-  double* vp = (double*)na_get_pointer_for_write(v);
-  for (i = 0; i < size; i++) { vp[i] = arr[i]; }
-  return v;
-}
-
-double* nary_to_dbl_vec(VALUE vec_val)
-{
-  int n_elements;
-  narray_t* vec_nary;
-  double* vec_pt;
-  double* vec;
-
-  if (vec_val == Qnil) return NULL;
-
-  GetNArray(vec_val, vec_nary);
-  n_elements = (int)NA_SHAPE(vec_nary)[0];
-
-  vec = ALLOC_N(double, n_elements);
-  vec_pt = (double*)na_get_pointer_for_read(vec_val);
-  memcpy(vec, vec_pt, n_elements * sizeof(double));
-
-  RB_GC_GUARD(vec_val);
-
-  return vec;
-}
-
-VALUE dbl_mat_to_nary(double** const mat, int const n_rows, int const n_cols)
-{
-  int i, j;
-  size_t shape[2] = { n_rows, n_cols };
-  VALUE v = rb_narray_new(numo_cDFloat, 2, shape);
-  double* vp = (double*)na_get_pointer_for_write(v);
-
-  for (i = 0; i < n_rows; i++) {
-    for (j = 0; j < n_cols; j++) {
-      vp[i * n_cols + j] = mat[i][j];
-    }
-  }
-
-  return v;
-}
-
-double** nary_to_dbl_mat(VALUE mat_val)
-{
-  int i, j;
-  int n_rows, n_cols;
-  narray_t* mat_nary;
-  double* mat_pt;
-  double** mat;
-
-  if (mat_val == Qnil) return NULL;
-
-  GetNArray(mat_val, mat_nary);
-  n_rows = (int)NA_SHAPE(mat_nary)[0];
-  n_cols = (int)NA_SHAPE(mat_nary)[1];
-
-  mat_pt = (double*)na_get_pointer_for_read(mat_val);
-  mat = ALLOC_N(double*, n_rows);
-  for (i = 0; i < n_rows; i++) {
-    mat[i] = ALLOC_N(double, n_cols);
-    for (j = 0; j < n_cols; j++) {
-      mat[i][j] = mat_pt[i * n_cols + j];
-    }
-  }
-
-  RB_GC_GUARD(mat_val);
-
-  return mat;
-}
-
-VALUE svm_nodes_to_nary(struct svm_node** const support_vecs, const int n_support_vecs)
-{
-  int i, j;
-  int n_dimensions = 0;
-  size_t shape[2] = { n_support_vecs, 1 };
-  VALUE v;
-  double* vp;
-
-  for (i = 0; i < n_support_vecs; i++) {
-    for (j = 0; support_vecs[i][j].index != -1; j++) {
-      if (n_dimensions < support_vecs[i][j].index) {
-        n_dimensions = support_vecs[i][j].index;
-      }
-    }
-  }
-
-  shape[1] = n_dimensions;
-  v = rb_narray_new(numo_cDFloat, 2, shape);
-  vp = (double*)na_get_pointer_for_write(v);
-  memset(vp, 0, n_support_vecs * n_dimensions * sizeof(double));
-
-  for (i = 0; i < n_support_vecs; i++) {
-    for (j = 0; support_vecs[i][j].index != -1; j++) {
-      vp[i * n_dimensions + support_vecs[i][j].index - 1] = support_vecs[i][j].value;
-    }
-  }
-
-  return v;
-}
-
-struct svm_node** nary_to_svm_nodes(VALUE nary_val)
-{
-  int i, j, k;
-  int n_rows, n_cols, n_nonzero_cols;
-  narray_t* nary;
-  double* nary_pt;
-  struct svm_node** support_vecs;
-
-  if (nary_val == Qnil) return NULL;
-
-  GetNArray(nary_val, nary);
-  n_rows = (int)NA_SHAPE(nary)[0];
-  n_cols = (int)NA_SHAPE(nary)[1];
-
-  nary_pt = (double*)na_get_pointer_for_read(nary_val);
-  support_vecs = ALLOC_N(struct svm_node*, n_rows);
-  for (i = 0; i < n_rows; i++) {
-    n_nonzero_cols = 0;
-    for (j = 0; j < n_cols; j++) {
-      if (nary_pt[i * n_cols + j] != 0) {
-        n_nonzero_cols++;
-      }
-    }
-    support_vecs[i] = ALLOC_N(struct svm_node, n_nonzero_cols + 1);
-    for (j = 0, k = 0; j < n_cols; j++) {
-      if (nary_pt[i * n_cols + j] != 0) {
-        support_vecs[i][k].index = j + 1;
-        support_vecs[i][k].value = nary_pt[i * n_cols + j];
-        k++;
-      }
-    }
-    support_vecs[i][n_nonzero_cols].index = -1;
-    support_vecs[i][n_nonzero_cols].value = 0.0;
-  }
-
-  RB_GC_GUARD(nary_val);
-
-  return support_vecs;
-}
-
-struct svm_node* dbl_vec_to_svm_node(double* const arr, int const size)
-{
-  int i, j;
-  int n_nonzero_elements;
-  struct svm_node* node;
-
-  n_nonzero_elements = 0;
-  for (i = 0; i < size; i++) {
-    if (arr[i] != 0.0) n_nonzero_elements++;
-  }
-
-  node = ALLOC_N(struct svm_node, n_nonzero_elements + 1);
-  for (i = 0, j = 0; i < size; i++) {
-    if (arr[i] != 0.0) {
-      node[j].index = i + 1;
-      node[j].value = arr[i];
-      j++;
-    }
-  }
-  node[n_nonzero_elements].index = -1;
-  node[n_nonzero_elements].value = 0.0;
-
-  return node;
-}

data/ext/numo/libsvm/converter.h

@@ -1,20 +0,0 @@
-#ifndef NUMO_LIBSVM_CONVERTER_H
-#define NUMO_LIBSVM_CONVERTER_H 1
-
-#include <string.h>
-#include <svm.h>
-#include <ruby.h>
-#include <numo/narray.h>
-#include <numo/template.h>
-
-VALUE int_vec_to_nary(int* const arr, int const size);
-int* nary_to_int_vec(VALUE vec_val);
-VALUE dbl_vec_to_nary(double* const arr, int const size);
-double* nary_to_dbl_vec(VALUE vec_val);
-VALUE dbl_mat_to_nary(double** const mat, int const n_rows, int const n_cols);
-double** nary_to_dbl_mat(VALUE mat_val);
-VALUE svm_nodes_to_nary(struct svm_node** const support_vecs, const int n_support_vecs);
-struct svm_node** nary_to_svm_nodes(VALUE nary_val);
-struct svm_node* dbl_vec_to_svm_node(double* const arr, int const size);
-
-#endif /* NUMO_LIBSVM_CONVERTER_H */

data/ext/numo/libsvm/kernel_type.c

@@ -1,22 +0,0 @@
-#include "kernel_type.h"
-
-RUBY_EXTERN VALUE mLibsvm;
-
-void rb_init_kernel_type_module()
-{
-  /**
-   * Document-module: Numo::Libsvm::KernelType
-   * The module consisting of constants for kernel type that used for parameter of LIBSVM.
-   */
-  VALUE mKernelType = rb_define_module_under(mLibsvm, "KernelType");
-  /* Linear kernel; u' * v */
-  rb_define_const(mKernelType, "LINEAR", INT2NUM(LINEAR));
-  /* Polynomial kernel; (gamma * u' * v + coef0)^degree */
-  rb_define_const(mKernelType, "POLY", INT2NUM(POLY));
-  /* RBF kernel; exp(-gamma * ||u - v||^2) */
-  rb_define_const(mKernelType, "RBF", INT2NUM(RBF));
-  /* Sigmoid kernel; tanh(gamma * u' * v + coef0) */
-  rb_define_const(mKernelType, "SIGMOID", INT2NUM(SIGMOID));
-  /* Precomputed kernel */
-  rb_define_const(mKernelType, "PRECOMPUTED", INT2NUM(PRECOMPUTED));
-}

data/ext/numo/libsvm/kernel_type.h

@@ -1,9 +0,0 @@
-#ifndef NUMO_LIBSVM_KERNEL_TYPE_H
-#define NUMO_LIBSVM_KERNEL_TYPE_H 1
-
-#include <svm.h>
-#include <ruby.h>
-
-void rb_init_kernel_type_module();
-
-#endif /* NUMO_LIBSVM_KERNEL_TYPE_H */