carray-numo-narray 0.9.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 02abe0434b7e13bc5f1fe24ea8c4842dbcd14fca6f09df4d8b928fe580df94d0
+  data.tar.gz: 6a33460577eeb87854725041a5406c1eda8e72b38cf7d9d24a7a435f772fca5b
+SHA512:
+  metadata.gz: 4ceaa09d46044291c2bb6bc4b9163fff86726d6d0057158691f255b6c39cc1749454b2138d736f2209816cfa53fe5b39be92d4b5746dabe3790ee4e23e2114d5
+  data.tar.gz: 192ddac1ba18da317e19073dd1965556f27b230e77c284d3f1e94f0e7d23b702a9fedf784a3d5fefd99cc8cb2a926441a678a02b1968c681c68731c42d88df01

data/NOTE.ja.txt

@@ -0,0 +1,34 @@
+Numo::NArrayとのインターフェース
+
+Numo::NArrayは、数値の型でクラスが分かれている。
+
+CArrayを参照するNArrayオブジェクトを作成する方法
+
+Int32を例にとると、CArrayと同じ形のNArrayをnary_newを作成する。
+生成されたNarrayのクラスはNumo::Int32となっている。
+nary_newを使うと必要になるまでデータ領域はアロケートされない。
+そこで、CArrayのポインタをNarrayのポインタと共有してしまえばよいのだが、
+問題はオブジェクトをGCするときのメモリの解放である。
+
+CArrayとNArrayの両方が自分のデータ領域を解放しようとするためこれは問題。
+そのためNArrayのfree関数をデータ領域を解放しないものに変更してやる必要がある。
+
+そこで、かなりの荒技だが、NArrayの実体であるRTypedStructの中の
+typeメンバー（free関数を含む）を独自のものに置き換えてしまうというもの。
+
+typeメンバーは numo-narrayの実装ではstaticとして隠蔽されているため、
+コンパイル時には取得できないので、動的に取得することにした。
+
+それぞれの型(UINT8, UINT16, ...)に対して、
+
+static rb_data_type_t na_wrap_int32 = { 0, {0, 0, 0,}, 0, 0,  0, };
+
+というデータタイプを用意しておく。
+CArray#numoが呼ばれると、CArrayの型に応じて、たとえば、Numo::Int32のオブジェクト
+が生成される。na_wrap_int32が初期値のままの状態の場合には、
+生成したオブジェクトのtypeメンバをna_wrap_int32にコピーして、
+function.dfreeだけ独自のfree関数に変更する。
+na_wrap_int32は適切に設定されたはずなので、以降はそのまま使う。
+
+これで、CArrayと同じポインタを参照するNArrayができる。
+

data/Rakefile

@@ -0,0 +1,11 @@
+GEMSPEC = "carray-numo-narray.gemspec"
+
+task :install do
+  spec = eval File.read(GEMSPEC)
+  system %{
+    gem build #{GEMSPEC}; gem install #{spec.full_name}.gem
+  }
+end
+
+require 'rspec/core/rake_task'
+RSpec::Core::RakeTask.new

data/carray-numo-narray.gemspec

@@ -0,0 +1,25 @@
+Gem::Specification::new do |s|
+  version = "0.9.0"
+  files = Dir.glob("**/*") - [ 
+                               Dir.glob("carray-numo-narray-*.gem"), 
+                               Dir.glob("test/**/*"),
+                               Dir.glob("work/**/*"),
+                             ].flatten
+
+  s.platform    = Gem::Platform::RUBY
+  s.name        = "carray-numo-narray"
+  s.summary     = "CArray interface to Numo::NArray"
+  s.description = <<-HERE
+    CArray interface to Numo::NArray
+  HERE
+  s.version     = version
+  s.licenses    = ['MIT']
+  s.author      = "Hiroki Motoyoshi"
+  s.email       = ""
+  s.homepage    = 'https://github.com/himotoyoshi/carray-numo-narray'
+  s.files       = files
+  s.extensions  = [ "ext/extconf.rb" ]
+  s.required_ruby_version = ">= 1.8.1"
+  s.add_runtime_dependency 'carray', '~> 1.3'
+end
+

data/ext/ca_wrap_numo_narray.c

@@ -0,0 +1,361 @@
+/* ---------------------------------------------------------------------------
+
+  carray/ca_wrap_narray.c
+
+  This file is part of Ruby/CArray extension library.
+  You can redistribute it and/or modify it under the terms of
+  the Ruby Licence.
+
+  Copyright (C) 2005 Hiroki Motoyoshi
+
+---------------------------------------------------------------------------- */
+
+#include "ruby.h"
+#include "carray.h"
+#include "numo/narray.h"
+
+#include <math.h>
+
+static char EMPTY_ARRAY_PTR;
+
+/* -------------------------------------------------------------------- */
+
+
+/* -------------------------------------------------------------------- */
+
+static int8_t
+na_class_to_ca_data_type (VALUE typecode)
+{
+  int8_t data_type;
+  
+  if ( typecode == numo_cUInt8 ) {
+    data_type = CA_UINT8;    
+  }
+  else if ( typecode == numo_cUInt16 ) {
+    data_type = CA_UINT16;
+  }
+  else if ( typecode == numo_cUInt32 ) {
+    data_type = CA_UINT32;
+  }
+  else if ( typecode == numo_cUInt64 ) {
+    data_type = CA_UINT64;
+  }
+  else if ( typecode == numo_cInt8 ) {
+    data_type = CA_INT8;    
+  }
+  else if ( typecode == numo_cInt16 ) {
+    data_type = CA_INT16;
+  }
+  else if ( typecode == numo_cInt32 ) {
+    data_type = CA_INT32;
+  }
+  else if ( typecode == numo_cInt64 ) {
+    data_type = CA_INT64;
+  }
+  else if ( typecode == numo_cSFloat ) {
+    data_type = CA_FLOAT32;
+  }
+  else if ( typecode == numo_cDFloat ) {
+    data_type = CA_FLOAT64;
+  }
+#ifdef HAVE_COMPLEX_H
+  else if ( typecode == numo_cSComplex ) {
+    data_type = CA_CMPLX64; 
+  }
+  else if ( typecode == numo_cDComplex ) {
+    data_type = CA_CMPLX128; 
+  }
+#endif
+  else if ( typecode == numo_cRObject ) {
+    data_type = CA_OBJECT; 
+  }
+  else {
+    rb_raise(rb_eRuntimeError, "invalid NArray class to convert to CArray");
+  }
+
+  return data_type;
+}
+
+static VALUE
+ca_data_type_to_na_class(int8_t data_type)
+{
+  volatile VALUE typecode;
+
+  switch ( data_type ) {
+  case CA_BOOLEAN:
+  case CA_UINT8:
+    typecode = numo_cUInt8; break;
+  case CA_UINT16:
+    typecode = numo_cUInt16; break;
+  case CA_UINT32:
+    typecode = numo_cUInt32; break;
+  case CA_UINT64:
+    typecode = numo_cUInt64; break;
+  case CA_INT8:
+    typecode = numo_cInt8; break;
+  case CA_INT16:
+    typecode = numo_cInt16; break;
+  case CA_INT32:
+    typecode = numo_cInt32; break;
+  case CA_INT64:
+    typecode = numo_cInt64; break;
+  case CA_FLOAT32:
+    typecode = numo_cSFloat; break;
+  case CA_FLOAT64:
+    typecode = numo_cDFloat; break;
+#ifdef HAVE_COMPLEX_H
+  case CA_CMPLX64:
+    typecode = numo_cSComplex; break;
+  case CA_CMPLX128:
+    typecode = numo_cDComplex; break;
+#endif
+  case CA_OBJECT:
+    typecode = numo_cRObject; break;
+  default:
+    rb_raise(rb_eRuntimeError, 
+             "no corresponding NArray typecode for CArray data type <%s>",
+             ca_type_name[data_type]);
+  }
+
+  return typecode;
+}
+
+/* -------------------------------------------------------------------- */
+
+static void
+na_wrap_free (void* ptr)
+{
+    narray_data_t *na = (narray_data_t*)ptr;
+
+    if (na->ptr != NULL) {
+        if (na->owned) {
+          /* xfree(na->ptr); */
+        }
+        na->ptr = NULL;
+    }
+    if (na->base.size > 0) {
+        if (na->base.shape != NULL && na->base.shape != &(na->base.size)) {
+            xfree(na->base.shape);
+            na->base.shape = NULL;
+        }
+    }
+    xfree(na);
+}
+
+static rb_data_type_t na_wrap_int8 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_int16 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_int32 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_int64 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_uint8 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_uint16 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_uint32 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_uint64 = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_sfloat = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_dfloat = { 0, {0, 0, 0,}, 0, 0,  0, };
+#ifdef HAVE_COMPLEX_H
+static rb_data_type_t na_wrap_scomplex = { 0, {0, 0, 0,}, 0, 0,  0, };
+static rb_data_type_t na_wrap_dcomplex = { 0, {0, 0, 0,}, 0, 0,  0, };
+#endif
+static rb_data_type_t na_wrap_robject = { 0, {0, 0, 0,}, 0, 0,  0, };
+
+static void
+na_set_wrap_type (VALUE obj)
+{
+  volatile VALUE typecode;
+  rb_data_type_t *data_type, *data_type_wrap;
+  
+  data_type = RTYPEDDATA_TYPE(obj);
+  typecode  = rb_obj_class(obj);
+  
+  if ( typecode == numo_cUInt8 ) {
+    data_type_wrap = &na_wrap_uint8;    
+  }
+  else if ( typecode == numo_cUInt16 ) {
+    data_type_wrap = &na_wrap_uint16;    
+  }
+  else if ( typecode == numo_cUInt32 ) {
+    data_type_wrap = &na_wrap_uint32;    
+  }
+  else if ( typecode == numo_cUInt64 ) {
+    data_type_wrap = &na_wrap_uint64;    
+  }
+  else if ( typecode == numo_cInt8 ) {
+    data_type_wrap = &na_wrap_int8;    
+  }
+  else if ( typecode == numo_cInt16 ) {
+    data_type_wrap = &na_wrap_int16;    
+  }
+  else if ( typecode == numo_cInt32 ) {
+    data_type_wrap = &na_wrap_int32;    
+  }
+  else if ( typecode == numo_cInt64 ) {
+    data_type_wrap = &na_wrap_int64;    
+  }
+  else if ( typecode == numo_cSFloat ) {
+    data_type_wrap = &na_wrap_sfloat;    
+  }
+  else if ( typecode == numo_cDFloat ) {
+    data_type_wrap = &na_wrap_dfloat;    
+  }
+#ifdef HAVE_COMPLEX_H
+  else if ( typecode == numo_cSComplex ) {
+    data_type_wrap = &na_wrap_scomplex;    
+  }
+  else if ( typecode == numo_cDComplex ) {
+    data_type_wrap = &na_wrap_dcomplex;    
+  }
+#endif
+  else if ( typecode == numo_cRObject ) {
+    data_type_wrap = &na_wrap_robject;    
+  }
+  else {
+    rb_raise(rb_eRuntimeError, "invalid NArray class to convert to CArray");
+  }
+
+  if ( ! ( data_type_wrap->wrap_struct_name ) ) {
+    *data_type_wrap = *data_type;
+    data_type_wrap->function.dfree = na_wrap_free;
+  }
+  
+  RTYPEDDATA_TYPE(obj) = data_type_wrap;
+}
+
+static VALUE
+rb_cary_na_ref_new (VALUE self)
+{
+  volatile VALUE obj, klass;
+  CArray *ca;
+  narray_t *na;
+
+  Data_Get_Struct(self, CArray, ca);
+
+  if ( ! ca_is_attached(ca) ) {
+    rb_raise(rb_eRuntimeError,
+             "cannot create NArray reference for not attached CArray");
+  }
+
+  klass = ca_data_type_to_na_class(ca->data_type);
+  obj = nary_new(klass, ca->rank, (size_t *)ca->dim);
+  na_set_wrap_type(obj);
+
+  GetNArray(obj, na);
+  NA_DATA_PTR(na)   = ca->ptr;
+  NA_DATA_OWNED(na) = TRUE;
+
+  rb_ivar_set(obj, rb_intern("@ref"), self);
+
+  return obj;
+}
+
+/* -------------------------------------------------------------------- */
+
+static VALUE
+rb_cary_to_na (VALUE self)
+{
+  volatile VALUE obj, type;
+  CArray *ca;
+  char *ptr;
+
+  Data_Get_Struct(self, CArray, ca);
+
+  type = ca_data_type_to_na_class(ca->data_type);
+  obj = nary_view_new(type, ca->rank, (size_t *)ca->dim);
+
+  ptr = na_get_pointer_for_write(obj);
+  ca_copy_data(ca, ptr);
+
+  return obj;
+}
+
+/* -------------------------------------------------------------------- */
+
+/* rdoc:
+  class NArray
+    def ca
+    end
+  end
+*/
+
+static VALUE
+rb_na_ca_ref_new (VALUE self)
+{
+  volatile VALUE obj = 0;
+  narray_t *na;
+  int8_t data_type;
+  char *ptr;
+
+  data_type = na_class_to_ca_data_type(rb_obj_class(self));
+
+  GetNArray(self, na);
+
+  if ( na->size == 0 ) {
+    ca_size_t zero = 0;
+    obj = rb_carray_wrap_ptr(data_type, 1, &zero, 0, NULL,
+                             &EMPTY_ARRAY_PTR, self); /* avoid ca->ptr == NULL */
+  }
+  else {
+    CA_CHECK_RANK(na->ndim);
+    ptr = na_get_pointer_for_read_write(self);
+    obj = rb_carray_wrap_ptr(data_type, na->ndim, (ca_size_t *)na->shape, 0, NULL, ptr, self);
+  }
+
+  return obj;
+}
+
+/* -------------------------------------------------------------------- */
+
+/* rdoc:
+  class NArray
+    def to_ca
+    end
+  end
+*/
+
+
+static VALUE
+rb_na_to_ca (VALUE self)
+{
+  volatile VALUE obj;
+  CArray *ca;
+  narray_t *na;
+  int8_t data_type;
+  char *ptr;
+
+  GetNArray(self, na);
+
+  data_type = na_class_to_ca_data_type(rb_obj_class(self));
+
+  if ( na->size == 0 ) {
+    ca_size_t zero = 0;
+    obj = rb_carray_new(data_type, 1, &zero, 0, NULL);
+  }
+  else {
+    CA_CHECK_RANK(na->ndim);
+
+    obj = rb_carray_new(data_type, na->ndim, (ca_size_t *) na->shape, 0, NULL);
+    Data_Get_Struct(obj, CArray, ca);
+    ptr = na_get_pointer_for_read(self);
+    ca_sync_data(ca, ptr);
+  }
+
+  return obj;
+}
+
+/* -------------------------------------------------------------------- */
+
+void
+Init_ca_wrap_numo_narray ()
+{
+  /* rb_require("narray"); */ /* "narray" should be loaded in config.rb */
+
+  rb_define_const(rb_cCArray, "HAVE_NUMO_NARRAY", Qtrue);
+
+  /* CArray -> NArray */
+  rb_define_method(rb_cCArray, "numo",     rb_cary_na_ref_new, 0);
+  rb_define_method(rb_cCArray, "to_numo",  rb_cary_to_na, 0);
+
+  /* NArray -> CArray */
+  rb_define_method(cNArray,    "ca",     rb_na_ca_ref_new, 0);
+  rb_define_method(cNArray,    "to_ca",  rb_na_to_ca, 0);
+}
+

data/ext/carray_numo_narray.c

@@ -0,0 +1,21 @@
+/* ---------------------------------------------------------------------------
+
+  carray/carray_narray.c
+
+  This file is part of Ruby/CArray extension library.
+  You can redistribute it and/or modify it under the terms of
+  the Ruby Licence.
+
+  Copyright (C) 2005 Hiroki Motoyoshi
+
+---------------------------------------------------------------------------- */
+
+#include "ruby.h"
+
+void Init_ca_wrap_numo_narray();
+
+void
+Init_carray_numo_narray ()
+{
+  Init_ca_wrap_numo_narray();
+}

data/ext/extconf.rb

@@ -0,0 +1,15 @@
+require 'numo/narray'
+require 'mkmf'
+require 'carray/mkmf'
+
+$LOAD_PATH.each do |path|
+  if File.exist? File.join(path, 'numo/numo/narray.h')
+    dir_config("narray", File.join(path,"numo"), File.join(path,"numo"))      
+    break
+  end
+end
+
+if have_carray() && have_header("numo/narray.h")
+  create_makefile("carray_numo_narray")
+end
+

data/lib/carray-numo-narray.rb

@@ -0,0 +1,3 @@
+require "numo/narray"
+require "carray"
+require "carray_numo_narray"

metadata

@@ -0,0 +1,66 @@
+--- !ruby/object:Gem::Specification
+name: carray-numo-narray
+version: !ruby/object:Gem::Version
+  version: 0.9.0
+platform: ruby
+authors:
+- Hiroki Motoyoshi
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2020-06-16 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: carray
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+description: "    CArray interface to Numo::NArray\n"
+email: ''
+executables: []
+extensions:
+- ext/extconf.rb
+extra_rdoc_files: []
+files:
+- NOTE.ja.txt
+- README.md
+- Rakefile
+- carray-numo-narray.gemspec
+- ext/ca_wrap_numo_narray.c
+- ext/carray_numo_narray.c
+- ext/extconf.rb
+- lib/carray-numo-narray.rb
+homepage: https://github.com/himotoyoshi/carray-numo-narray
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 1.8.1
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.7.7
+signing_key: 
+specification_version: 4
+summary: CArray interface to Numo::NArray
+test_files: []