carray 1.1.4 → 1.1.5

data/ruby_float_func.c

@@ -0,0 +1,83 @@
+/* ---------------------------------------------------------------------------
+
+  ruby_numeric.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 <math.h>
+
+#define rb_num_func(name) \
+static VALUE rb_num_ ## name (VALUE self) \
+{ return rb_float_new(name(NUM2DBL(self))); }
+
+rb_num_func(sqrt);
+rb_num_func(exp);
+rb_num_func(log);
+rb_num_func(log10);
+
+rb_num_func(sin);
+rb_num_func(cos);
+rb_num_func(tan);
+rb_num_func(sinh);
+rb_num_func(cosh);
+rb_num_func(tanh);
+
+rb_num_func(asin);
+rb_num_func(acos);
+rb_num_func(atan);
+rb_num_func(asinh);
+rb_num_func(acosh);
+rb_num_func(atanh);
+
+static VALUE 
+rb_num_rad (VALUE self)
+{ return rb_float_new(0.0174532925199433*NUM2DBL(self)); }
+
+static VALUE 
+rb_num_deg (VALUE self)
+{ return rb_float_new(57.2957795130823*NUM2DBL(self)); }
+
+static VALUE 
+rb_num_distance (VALUE self, VALUE other)
+{ 
+  double fs = NUM2DBL(self), fo = NUM2DBL(other);
+  return rb_float_new(fabs(fs - fo)); 
+}
+
+void
+Init_numeric_float_function ()
+{
+  VALUE mod = rb_define_module_under(rb_cNumeric, "FloatFunction");
+
+  rb_define_method(mod, "sqrt", rb_num_sqrt, 0);
+  rb_define_method(mod, "exp", rb_num_exp, 0);
+  rb_define_method(mod, "log", rb_num_log, 0);
+  rb_define_method(mod, "log10", rb_num_log10, 0);
+  rb_define_method(mod, "sin", rb_num_sin, 0);
+  rb_define_method(mod, "cos", rb_num_cos, 0);
+  rb_define_method(mod, "tan", rb_num_tan, 0);  
+  rb_define_method(mod, "sinh", rb_num_sinh, 0);
+  rb_define_method(mod, "cosh", rb_num_cosh, 0);
+  rb_define_method(mod, "tanh", rb_num_tanh, 0);  
+  rb_define_method(mod, "asin", rb_num_asin, 0);
+  rb_define_method(mod, "acos", rb_num_acos, 0);
+  rb_define_method(mod, "atan", rb_num_atan, 0);  
+  rb_define_method(mod, "asinh", rb_num_asinh, 0);
+  rb_define_method(mod, "acosh", rb_num_acosh, 0);
+  rb_define_method(mod, "atanh", rb_num_atanh, 0);  
+  rb_define_method(mod, "rad", rb_num_rad, 0);  
+  rb_define_method(mod, "deg", rb_num_deg, 0);    
+  rb_define_method(mod, "distance", rb_num_distance, 1);    
+
+  rb_include_module(rb_cFloat, mod);
+  rb_include_module(rb_cFixnum, mod);
+  rb_include_module(rb_cBignum, mod);  
+}
+

data/spec/CABlockIterator/CABlockIterator_spec.rb

@@ -0,0 +1,113 @@
+require "carray"
+
+describe CABlockIterator do
+
+  describe "each method" do
+
+    before do
+      @original = CArray.object(4,4).seq!
+      @it = @original.blocks(0..1, 0..1)
+    end
+
+    it "should call 4-times its block" do
+      count = 0
+      @it.each do |blk|
+        count += 1
+      end
+      ( count ).should == 4
+    end
+
+    it "should have the block parameter of 2x2 CABlock object" do
+      @it.each do |blk|
+        ( blk ).should be_an_instance_of(CABlock)
+        ( blk.data_type ).should == CA_OBJECT
+        ( blk.dim ).should == [2,2]
+      end
+    end
+
+  end
+
+  describe "reference method []" do
+
+    before do
+      @original = CArray.object(4,4).seq!
+      @it = @original.blocks(0..1, 0..1)
+    end
+
+    it "should return a 2x2 CABlock object" do
+      blk = @it[0]
+      ( blk ).should be_an_instance_of(CABlock)
+      ( blk.data_type ).should == CA_OBJECT
+      ( blk ).should == CA_OBJECT([[0,1],[4,5]])
+      blk = @it[1,0]
+      ( blk ).should be_an_instance_of(CABlock)
+      ( blk.data_type ).should == CA_OBJECT
+      ( blk ).should == CA_OBJECT([[8,9],[12,13]])
+    end
+
+    it "should return a 2x2 CABlock object" do
+      ( @it[0] ).should == @it.kernel_at_addr(0, @it.reference)
+      ( @it[1,0] ).should == @it.kernel_at_index([1,0], @it.reference)
+    end
+
+  end
+
+  describe "store method []=" do
+
+    before do
+      @original = CArray.object(4,4).seq!
+      @it = @original.blocks(0..1, 0..1)
+    end
+
+    it "should return a 2x2 CABlock object" do
+      @it[0] = 1
+      @it[1,0] = 2
+      ( @it.reference ).should == CA_OBJECT( [ [ 1, 1, 2, 3 ],
+                                               [ 1, 1, 6, 7 ],
+                                               [ 2, 2, 10, 11 ],
+                                               [ 2, 2, 14, 15 ] ] )
+    end
+
+  end
+
+  describe "pick method" do
+
+    before do
+      @original = CArray.object(4,4).seq!
+      @it = @original.blocks(0..1, 0..1)
+    end
+
+    it "should return a 2x2 CABlock object" do
+      blk = @it.pick(0)
+      ( blk ).should be_an_instance_of(CArray)
+      ( blk.data_type ).should == CA_OBJECT
+      ( blk ).should == CA_OBJECT([[0,2],[8,10]])
+      blk = @it.pick(1,0)
+      ( blk ).should be_an_instance_of(CArray)
+      ( blk.data_type ).should == CA_OBJECT
+      ( blk ).should == CA_OBJECT([[4,6],[12,14]])
+    end
+
+  end
+
+  describe "put method" do
+
+    before do
+      @original = CArray.object(4,4).seq!
+      @it = @original.blocks(0..1, 0..1)
+    end
+
+    it "should store value via 2x2 CBlock object" do
+      @it.put(0,1)
+      @it.put(0,1,2)
+      ( @it.pick(0) ).should == CA_OBJECT([[1,1],[1,1]])
+      ( @it.pick(0,1) ).should == CA_OBJECT([[2,2],[2,2]])
+      ( @it.reference ).should == CA_OBJECT([[ 1, 2, 1, 2 ],
+                                             [ 4, 5, 6, 7 ],
+                                             [ 1, 2, 1, 2 ],
+                                             [ 12, 13, 14, 15 ]])
+    end
+
+  end
+
+end

data/spec/CArray/bug/store_spec.rb

@@ -0,0 +1,27 @@
+
+require "carray"
+
+describe "CArray#[]" do
+
+  it "should return one element array with [nil,-1]" do
+    ca = CArray.int(4,4).seq!
+    ( ca[[nil,-1]].elements ).should == 1
+    ( ca[[nil,-1]] ).should == ca[[0]]
+  end
+
+  it "should return self[[0],nil] with [nil,-1], nil" do
+    ca = CArray.int(4,4).seq!
+    ( ca[[nil,-1],nil] ).should == ca[[0],nil]
+  end
+
+  it "should return self[[0],nil] with nil, [nil,-1]" do
+    ca = CArray.int(4,4).seq!
+    ( ca[nil,[nil,-1]] ).should == ca[nil,[0]]
+  end
+
+  it "should return self[[0],nil] with [nil,-1], [nil,-1]" do
+    ca = CArray.int(4,4).seq!
+    ( ca[[nil,-1],[nil,-1]] ).should == ca[[0],[0]]
+  end
+
+end

data/spec/CArray/index/repeat_spec.rb

@@ -0,0 +1,10 @@
+
+require "carray"
+
+describe "CArray#[]" do
+  it "should return CARepeat object" do
+    ca = CArray.int(3).seq!
+    ca[:%,3] # =>
+    ca[3,:%] # => 
+  end
+end

data/spec/CArray/method/eq_spec.rb

@@ -0,0 +1,80 @@
+
+require "carray"
+
+describe "CArray#eq" do
+
+  it "returns 1 if self and other equal (int)" do
+    ca123 = CA_INT([123])
+
+    ( (ca123).eq(123  )[0] ).should == 1
+    ( (ca123).eq(ca123)[0] ).should == 1
+  end
+
+  it "returns 1 if self and other equal (object)" do
+    ca123 = CA_OBJECT(["123"])
+
+    ( (ca123).eq("123")[0] ).should == 1
+    ( (ca123).eq(ca123)[0] ).should == 1
+  end
+
+  it "returns 1 if self and other equal (fixlen)" do
+    ca123 = CA_FIXLEN(["123"], :bytes=>3)
+
+    ( (ca123).eq("123")[0] ).should == 1
+    ( (ca123).eq(ca123)[0] ).should == 1
+  end
+
+  it "returns 0 if self and other equal (int)" do
+    ca123 = CA_INT([123])
+    ca124 = CA_INT([124])
+
+    ( (ca123).eq(124  )[0] ).should == 0
+    ( (ca123).eq(ca124)[0] ).should == 0
+  end
+
+  it "returns 0 if self and other equal (object)" do
+    ca123 = CA_OBJECT(["123"])
+    ca124 = CA_OBJECT(["124"])
+
+    ( (ca123).eq("124")[0] ).should == 0
+    ( (ca123).eq(ca124)[0] ).should == 0
+  end
+
+  it "returns 0 if self and other equal (fixlen)" do
+    ca123 = CA_FIXLEN(["123"], :bytes=>3)
+    ca124 = CA_FIXLEN(["124"], :bytes=>3)
+    ca12  = CA_FIXLEN(["12"], :bytes=>3)
+
+    ( (ca123).eq("124")[0] ).should == 0
+    ( (ca123).eq(ca124)[0] ).should == 0
+
+    ( (ca123).eq("12")[0] ).should == 0
+    ( (ca123).eq(ca12)[0] ).should == 0
+  end
+
+  it "returns 0 if self is a NaN" do
+    nan = 0.0/0.0
+
+    ( (nan             ).eq(CA_DOUBLE([0.0]))[0] ).should == 0
+    ( (CA_DOUBLE([nan])).eq(0.0             )[0] ).should == 0
+    ( (CA_DOUBLE([nan])).eq(CA_DOUBLE([0.0]))[0] ).should == 0
+  end
+
+  it "returns 0 if self and other are NaN" do
+    nan = 0.0/0.0
+
+    ( (nan             ).eq(CA_DOUBLE([nan]))[0] ).should == 0
+    ( (CA_DOUBLE([nan])).eq(nan             )[0] ).should == 0
+    ( (CA_DOUBLE([nan])).eq(CA_DOUBLE([nan]))[0] ).should == 0
+  end
+
+  it "returns 0 if other is a NaN" do
+    nan = 0.0/0.0
+
+    ( (0.0             ).eq(CA_DOUBLE([nan]))[0] ).should == 0
+    ( (CA_DOUBLE([0.0])).eq(nan             )[0] ).should == 0
+    ( (CA_DOUBLE([0.0])).eq(CA_DOUBLE([nan]))[0] ).should == 0
+  end
+
+end
+

data/spec/CArray/method/is_nan_spec.rb

@@ -0,0 +1,12 @@
+
+require "carray"
+
+describe "CArray#is_nan" do
+  it "returns 1 if self is not a valid IEEE floating-point number" do
+    CA_DOUBLE(0.0).is_nan[0].should == 0
+    CA_DOUBLE(-1.5).is_nan[0].should == 0
+    (CA_DOUBLE(0.0)/CA_DOUBLE(0.0)).is_nan[0].should == 1
+    CA_DOUBLE(0.0/0.0).is_nan[0].should == 1
+  end
+end
+

data/spec/CArray/method/ne_spec.rb

@@ -0,0 +1,18 @@
+
+require "carray"
+
+describe "CArray#ne" do
+  it "returns 0 if self is not a valid IEEE floating-point number" do
+    nan = 0.0/0.0
+    0.0.ne(CA_DOUBLE(nan))[0].should == 1
+    nan.ne(CA_DOUBLE(nan))[0].should == 1
+    nan.ne(CA_DOUBLE(0.0))[0].should == 1
+    CA_DOUBLE(nan).ne(0.0)[0].should == 1
+    CA_DOUBLE(nan).ne(nan)[0].should == 1
+    CA_DOUBLE(0.0).ne(nan)[0].should == 1
+    CA_DOUBLE(nan).ne(CA_DOUBLE(0.0))[0].should == 1
+    CA_DOUBLE(nan).ne(CA_DOUBLE(nan))[0].should == 1
+    CA_DOUBLE(0.0).ne(CA_DOUBLE(nan))[0].should == 1
+  end
+end
+    

data/spec/CArray/method/round_spec.rb

@@ -0,0 +1,11 @@
+
+require "carray"
+
+describe "CArray#round" do
+  it "returns the nearest Integer" do
+    CA_DOUBLE(5.5).round[0].should == 6
+    CA_DOUBLE(0.4).round[0].should == 0
+    CA_DOUBLE(-2.8).round[0].should == -3
+    CA_DOUBLE(0.0).round[0].should == 0
+  end
+end

data/spec/CArray/object/_attribute_spec.rb

@@ -0,0 +1,32 @@
+
+require "carray"
+
+describe "CArray#elements" do
+
+  it "should return element number of array" do
+    ca = CArray.int32(2,3,4)
+    ( ca.elements ).should == 2*3*4
+    ( ca.length ).should == 2*3*4
+    ( ca.size ).should == 2*3*4
+  end
+
+end
+
+describe "CArray#data_type" do
+  
+  it "should return data type of array" do
+    ca = CArray.int32(2,3,4)
+    ( ca.data_type ).should == CA_INT32
+  end
+
+end
+
+describe "CArray#dim" do
+
+  it "should return a array contains the shape of array" do
+    ca = CArray.int32(2,3,4)
+    ( ca.dim ).should == [2,3,4]
+  end
+
+end
+

data/spec/CArray/object/s_new_spec.rb

@@ -0,0 +1,31 @@
+
+require "carray"
+
+describe "CArray.new" do
+
+  it "should return a CA_INT type carray if CA_INT is specified as data_type" do
+    ca = CArray.new(CA_INT, [2,3,4])
+    ( ca.class ).should == CArray
+    ( ca.data_type ).should == CA_INT
+  end
+
+  it "should return a CA_FIXLEN type carray if CA_FIXLEN is specified as data_type" do
+    ca = CArray.new(CA_FIXLEN, [2,3,4], :bytes=>3)
+    ( ca.class ).should == CArray
+    ( ca.data_type ).should == CA_FIXLEN
+    ( ca.bytes ).should == 3
+  end
+
+  it "should return a carray filled by 0 when data_type == CA_INT" do
+    ca = CArray.new(CA_INT, [2,3,4])
+    ( ca[0] ).should == 0
+    ( ca.all_equal?(0) ).should == true
+  end
+
+  it "should return a carray filled by string filled by 0 when data_type == CA_FIXLEN" do
+    ca = CArray.new(CA_FIXLEN, [2,3,4], :bytes=>3)
+    ( ca[0] ).should == "\0\0\0"
+    ( ca.all_equal?("\0\0\0") ).should == true
+  end
+
+end

data/spec/CArray/serialize/Serialization_spec.rb

@@ -0,0 +1,89 @@
+require "carray"
+
+describe CArray do
+
+  describe "loaded from binary format (int32 type)" do
+
+    before do
+      @original = CArray.int32(10,10) { 10 }
+      @original[5,5] = UNDEF
+      CArray.save(@original, "test.ca")
+      @it = CArray.load("test.ca")
+    end
+
+    it "should equal to original" do
+      ( @it ).should == @original
+      ( @it.has_mask? ).should == true
+    end
+
+    it "should dump string of same contents with file from which it was loaded" do 
+      ( CArray.dump(@it) ).should == File.read("test.ca")
+    end
+
+    after do
+      File.unlink("test.ca")
+    end
+
+  end
+
+  describe "loaded from binary format (object type)" do
+
+    before do
+      @original = CArray.object(10,10) { 3.times { Time.now } }
+      @original[5,5] = UNDEF
+      CArray.save(@original, "test.ca")
+      @it = CArray.load("test.ca")
+    end
+
+    it "should equal to original" do
+      ( @it ).should == @original
+      ( @it.has_mask? ).should == true
+    end
+
+    after do
+      File.unlink("test.ca")
+    end
+
+  end
+
+  describe "loaded by marshal (int32 type)" do
+    
+    before do
+      @original = CArray.int32(10,10) { 10 }
+      @original[5,5] = UNDEF
+      open("test.ca", "w") { |io| Marshal.dump(@original, io) }
+      @it = open("test.ca") { |io| Marshal.load(io) }
+    end
+
+    it "should equal to original" do
+      ( @it ).should == @original
+      ( @it.has_mask? ).should == true
+    end
+
+    after do
+      File.unlink("test.ca")
+    end
+
+  end
+
+  describe "loaded by marshal (object type)" do
+    
+    before do
+      @original = CArray.object(10,10) { 3.times { Time.now } }
+      @original[5,5] = UNDEF
+      open("test.ca", "w") { |io| Marshal.dump(@original, io) }
+      @it = open("test.ca") { |io| Marshal.load(io) }
+    end
+
+    it "should equal to original" do
+      ( @it ).should == @original
+      ( @it.has_mask? ).should == true
+    end
+
+    after do
+      File.unlink("test.ca")
+    end
+
+  end
+
+end