carray 1.3.6 → 1.5.3

data/spec/Features/feature_math_spec.rb

@@ -0,0 +1,97 @@
+
+require 'carray'
+require "rspec-power_assert"
+
+describe "TestMath " do
+
+  example "zerodiv" do
+    zero = CArray.int(3,3) { 0 }
+    one  = CArray.int(3,3) { 1 }
+    expect { 1/zero }.to raise_error(ZeroDivisionError)
+    expect { one/zero }.to raise_error(ZeroDivisionError) # div, div!
+    expect { one/0 }.to raise_error(ZeroDivisionError)
+    expect { one/zero }.to raise_error(ZeroDivisionError)
+    expect { one % 0 }.to raise_error(ZeroDivisionError)  # mod, mod!
+    expect { one % zero }.to raise_error(ZeroDivisionError)
+    expect { zero.rcp }.to raise_error(ZeroDivisionError) # rcp!
+    expect { zero.rcp_mul(1) }.to raise_error(ZeroDivisionError) # rcp_mul!
+    expect { zero.rcp_mul(one) }.to raise_error(ZeroDivisionError)
+    expect { 0 ** (-one) }.to raise_error(ZeroDivisionError) # pow, pow!
+    expect { zero ** (-one) }.to raise_error(ZeroDivisionError)
+  end
+
+  example "cmp" do
+    a = CA_INT([0,1,2,3,4])
+    b = CA_INT([4,3,2,1,0])
+
+    # ---
+    is_asserted_by { CA_INT8([-1, -1, 0, 1, 1]) == (a <=> 2) }
+    is_asserted_by { CA_INT8([1, 1, 0, -1, -1]) == (2 <=> a) }
+    is_asserted_by { CA_INT8([-1, -1, 0, 1, 1]) == (a <=> b) }
+  end
+
+  example "bit_op" do
+    a = CArray.int(10).seq!
+
+    # ---
+    is_asserted_by { a.convert {|x| ~x } == (~a) }
+
+    is_asserted_by { a.convert {|x| x & 1 } == (a & 1) }
+    is_asserted_by { (a & 1) == (1 & a) }
+
+    is_asserted_by { a.convert { |x| x | 1 } == (a | 1) }
+    is_asserted_by { (a | 1) == (1 | a) }
+
+    is_asserted_by { a.convert { |x| x ^ 1 } == (a ^ 1) }
+    is_asserted_by { (a ^ 1) == (1 ^ a) }
+
+    # ---
+    is_asserted_by { a.convert { |x| x << 1 } == (a << 1) }
+    is_asserted_by { a.convert { |x| 1 << x } == (1 << a) }
+
+    is_asserted_by { a.convert { |x| x >> 1 } == (a >> 1) }
+    is_asserted_by { a.convert { |x| 1 >> x } == (1 >> a) }
+  end
+
+  example "max_min" do
+    a = CArray.int(3,3).seq!
+    b = a.reverse
+
+    # ---
+    is_asserted_by { CA_INT32([[8,7,6],
+                               [5,4,5],
+                               [6,7,8]]) == a.pmax(b) }
+    is_asserted_by { CA_INT32([[0,1,2],
+                               [3,4,3],
+                               [2,1,0]]) == a.pmin(b) }
+    # ---
+    is_asserted_by { CA_INT32([[0,0,0],
+                               [0,0,0],
+                               [0,0,0]]) == a.pmin(0) }
+    is_asserted_by { CA_INT32([[8,8,8],
+                               [8,8,8],
+                               [8,8,8]]) == a.pmax(8) }
+    # ---
+    is_asserted_by { CA_INT32([[0,0,0],
+                               [0,0,0],
+                               [0,0,0]]) == CA_INT32(0).pmin(a) }
+    is_asserted_by { CA_INT32([[8,8,8],
+                               [8,8,8],
+                               [8,8,8]]) == CA_INT32(8).pmax(a) }
+    # ---
+    is_asserted_by { CA_INT32([[8,7,6],
+                               [5,4,5],
+                               [6,7,8]]) == CAMath.max(a, b) }
+    is_asserted_by { CA_INT32([[0,1,2],
+                               [3,4,3],
+                               [2,1,0]]) == CAMath.min(a, b) }
+    # ---
+    is_asserted_by { CA_INT32([[0,0,0],
+                               [0,0,0],
+                               [0,0,0]]) == CAMath.min(0,a,b,8) }
+    is_asserted_by { CA_INT32([[8,8,8],
+                               [8,8,8],
+                               [8,8,8]]) == CAMath.max(0,a,b,8) }
+  end
+
+end

data/spec/Features/feature_order_spec.rb

@@ -0,0 +1,146 @@
+
+require 'carray'
+require "rspec-power_assert"
+
+describe "TestCArrayOrder " do
+
+  example "where" do
+    # ---
+    a = CArray.int(3,3).seq!
+    c = (a % 2).eq(0)
+    is_asserted_by { a[c] == c.where.int }
+  end
+
+  example "index" do
+    # ---
+    a = CArray.float(3,3).seq!
+    is_asserted_by { CA_INT32([[0,1,2],
+                               [3,4,5],
+                               [6,7,8]]) == a.address }
+    is_asserted_by { CA_INT32([0, 1, 2]) == a.index(0) }
+    is_asserted_by { CA_INT32([0, 1, 2]) == a.index(1) }
+  end
+
+  example "indices" do
+    # ---
+    a = CArray.float(3,3).seq!
+    a.indices {|x, y|
+      is_asserted_by { CA_INT32([[0,0,0],
+                                 [1,1,1],
+                                 [2,2,2]]) == x }
+      is_asserted_by { CA_INT32([[0,1,2],
+                                 [0,1,2],
+                                 [0,1,2]]) == y }
+    }
+  end
+
+  example "reverse" do
+    a = CArray.int32(3,3).seq!
+    b = CArray.int32(3,3).seq!(8,-1)
+
+    # ---
+    is_asserted_by { b == a.reverse }
+    is_asserted_by { b == a.reversed }
+
+    c = a.to_ca.reverse!
+    is_asserted_by { b == c }
+  end
+
+  example "sort" do
+    a = CA_INT32([[5,4,3],
+                  [0,1,2],
+                  [8,7,6]])
+    s = CA_INT32([[0,1,2],
+                  [3,4,5],
+                  [6,7,8]])
+    # ---
+    is_asserted_by { s == a.sort }
+    b = a.to_ca
+    b.sort!
+    is_asserted_by { s == b }
+  end
+
+  example "sort_addr" do
+    a = CA_INT32([[5,4,3],
+                  [0,1,2],
+                  [8,7,6]])
+    # ---
+    is_asserted_by { CA_SIZE([[3,4,5],
+                              [2,1,0],
+                              [8,7,6]]) == a.sort_addr }
+    is_asserted_by { CA_INT32([[5,4,3],
+                               [0,1,2],
+                               [8,7,6]]) == a.order }
+    # ---
+    is_asserted_by { CA_SIZE([[6,7,8],
+                              [0,1,2],
+                              [5,4,3]]) == a.sort_addr.reverse }
+    is_asserted_by { CA_INT32([[3,4,5],
+                               [8,7,6],
+                               [0,1,2]]) == a.order(-1) }
+  end
+
+#  example "uniq" do
+#    a = CArray.int(100).random!(10)
+#
+#    # ---
+#    is_asserted_by { a.to_a.uniq, a.uniq.to_a)
+#
+#  end
+
+  example "search" do
+    # ---
+    a = CArray.int(9,9).seq!
+    is_asserted_by { 50 == a.bsearch(50) }
+    is_asserted_by { 50 == a.search(50) }
+    is_asserted_by { [5, 5] == a.bsearch_index(50) }
+    is_asserted_by { [5, 5] == a.search_index(50) }
+    # ---
+    a = CArray.int(9,9).seq!(0,2)
+    is_asserted_by { nil == a.bsearch(51) }
+    is_asserted_by { nil == a.search(51) }
+    is_asserted_by { nil == a.bsearch_index(51) }
+    is_asserted_by { nil == a.search_index(51) }
+    # ---
+    a = CArray.float(9,9).span!(0..80)
+    is_asserted_by { 50 == a.bsearch(50) }
+    is_asserted_by { 50 == a.search(50) }
+    is_asserted_by { [5, 5] == a.bsearch_index(50) }
+    is_asserted_by { [5, 5] == a.search_index(50) }
+
+    is_asserted_by { nil == a.bsearch(50.0001) }
+    is_asserted_by { nil == a.search(50.0001) }
+    is_asserted_by { nil == a.bsearch_index(50.0001) }
+    is_asserted_by { nil == a.search_index(50.0001) }
+    # ---
+    a = CArray.float(9,9).span!(0..80)    
+    is_asserted_by { 50 == a.search(50.0001, 0.0001) }
+    is_asserted_by { [5, 5] == a.search_index(50.0001, 0.0001) }
+
+    is_asserted_by { nil == a.search(50.0001, 1.0e-05) }
+    is_asserted_by { nil == a.search_index(50.0001, 1.0e-05) }
+  end
+  
+  example "search_nearest   " do
+    # ---
+    a = CArray.int(9,9).seq!
+    is_asserted_by { 50 == a.search_nearest(50) }
+    is_asserted_by { [5, 5] == a.search_nearest_index(50) }
+    # ---
+    a = CArray.int(9,9).seq!(0,2)
+    is_asserted_by { 25 == a.search_nearest(51) }
+    is_asserted_by { [2, 7] == a.search_nearest_index(51) }
+    # ---
+    a = CArray.float(9,9).span!(0..80)
+    is_asserted_by { 50 == a.search_nearest(50) }
+    is_asserted_by { [5, 5] == a.search_nearest_index(50) }
+
+    is_asserted_by { 50 == a.search_nearest(50.5) }
+    is_asserted_by { [5, 5] == a.search_nearest_index(50.5) }
+
+    # first detected
+    is_asserted_by { 49 == a.search_nearest(49.5) }
+    is_asserted_by { [5, 4] == a.search_nearest_index(49.5) }
+    
+  end
+end

data/spec/Features/feature_ref_store_spec.rb

@@ -0,0 +1,209 @@
+
+require 'carray'
+require "rspec-power_assert"
+
+describe "TestCArrayRefStore " do
+
+  example "access_point" do
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  0 == a[0, 0] }
+    is_asserted_by {  4 == a[1, 1] }
+    is_asserted_by {  8 == a[-1, -1] }
+    expect { a[0,3] }.to raise_error(IndexError)
+    expect { a[-4,0] }.to raise_error(IndexError)
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[0,0] = -1
+    a[1,1] = -2
+    a[-1,-1] = -3
+    is_asserted_by {  (-1) == a[0, 0] }
+    is_asserted_by {  (-2) == a[1, 1] }
+    is_asserted_by {  (-3) == a[-1, -1] }
+    expect { a[0,3] = -4 }.to raise_error(IndexError)
+    expect { a[-4,0] = -5 }.to raise_error(IndexError)
+  end
+
+  example "access_block" do
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  CA_INT([1, 4, 7]) == a[0..-1, 1] }
+    is_asserted_by {  CA_INT([[1], [4], [7]]) == a[0..-1, [1]] }
+    is_asserted_by {  CABlock == a[0..2, 0..2].class }
+    is_asserted_by {  a == a[0..2, 0..2] }
+    is_asserted_by {  CABlock == a[nil, nil].class }
+    is_asserted_by {  a == a[nil, nil] }
+    is_asserted_by {  a == a[-3..-1, -3..-1] }
+    is_asserted_by {  CA_INT([[0, 1, 2], [3, 4, 5]]) == a[0..1, nil] }
+    is_asserted_by {  CA_INT([[0, 2], [3, 5]]) == a[0..1, [nil, 2]] }
+    expect { a[0,0..3]  }.to raise_error(IndexError)
+    expect { a[-4..-1,0] }.to raise_error(IndexError)
+    expect { a[0,[0,3,2]] }.to raise_error(IndexError)
+    expect { a[0,[0,4]] }.to raise_error(IndexError)
+
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  CA_INT([[0, 1, 2], [3, 4, 5], [6, 7, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[0..-1, 1] = 9
+    is_asserted_by {  CA_INT([[0, 9, 2], [3, 9, 5], [6, 9, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[0..-1, [1]] = 9
+    is_asserted_by {  CA_INT([[0, 9, 2], [3, 9, 5], [6, 9, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[nil, nil] = 9
+    is_asserted_by {  CA_INT([[9, 9, 9], [9, 9, 9], [9, 9, 9]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[0..1, [nil,2]] = 9
+    is_asserted_by {  CA_INT([[9, 1, 9], [9, 4, 9], [6, 7, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    expect { a[0,0..3] = -3  }.to raise_error(IndexError)
+    expect { a[-4..-1,0] = -4 }.to raise_error(IndexError)
+    expect { a[0,[0,3,2]] = -5  }.to raise_error(IndexError)
+    expect { a[0,[0,4]] = -6  }.to raise_error(IndexError)
+  end
+
+  example "access_selection" do
+    # ---
+    a = CArray.int(3,3).seq!;
+    is_asserted_by {  CA_INT([0, 1, 2, 3, 4]) == a[a < 5] }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[a < 5] = 1
+    is_asserted_by {  CA_INT([[1, 1, 1], [1, 1, 5], [6, 7, 8]]) == a }
+  end
+
+#  example "access_mapping" do
+#    # ---
+#    a = CArray.int(3,3).seq!;
+#    m = a.shuffle
+#    is_asserted_by { 10*m, (10*a)[m])
+
+    # ---
+#    a = CArray.int(3,3).seq!
+#    b = 10*a
+#    m = a.shuffle
+#    a[m] = (10*m)
+#    is_asserted_by { b, a)
+#  end
+
+  example "access_grid" do
+    # ---
+    a = CArray.int(3,3).seq!;
+    is_asserted_by {  CA_INT([[1, 2], [4, 5]]) == a[CA_INT([0, 1]), 1..2] }
+    is_asserted_by {  CA_INT([[0, 2], [6, 8]]) == a[+[0, 2], +[0, 2]] }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[+[0,2],+[0,2]] = 9
+    is_asserted_by {  CA_INT([[9, 1, 9], [3, 4, 5], [9, 7, 9]]) == a }
+  end
+
+  example "access_repeat" do
+    # ---
+    # See test_CARepeat
+  end
+
+  example "access_method" do
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  CA_INT([0, 1, 2, 3, 4]) == a[:lt, 5] }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[:lt,5] = 1
+    is_asserted_by {  CA_INT([[1, 1, 1], [1, 1, 5], [6, 7, 8]]) == a }
+  end
+
+  example "access_address" do
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  0 == a[0] }
+    is_asserted_by {  4 == a[4] }
+    is_asserted_by {  8 == a[-1] }
+    expect { a[9] }.to raise_error(IndexError)
+    expect { a[-10] }.to raise_error(IndexError)
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[0] = -1
+    a[4] = -2
+    a[-1] = -3
+    is_asserted_by {  (-1) == a[0, 0] }
+    is_asserted_by {  (-2) == a[1, 1] }
+    is_asserted_by {  (-3) == a[-1, -1] }
+    expect { a[9] = -4 }.to raise_error(IndexError)
+    expect { a[-10] = -5  }.to raise_error(IndexError)
+  end
+
+  example "access_address_block" do
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  CA_INT([4]) == a[[4]] }
+    is_asserted_by {  CA_INT([0, 1, 2, 3]) == a[[(0..3)]] }
+    is_asserted_by {  CARefer == a[nil].class }
+    is_asserted_by {  CA_INT([0, 1, 2, 3, 4, 5, 6, 7, 8]) == a[nil] }
+    expect { a[0..9]  }.to raise_error(IndexError)
+    expect {  a[-10..-1] }.to raise_error(IndexError)
+
+    # ---
+    a = CArray.int(3,3).seq!
+    is_asserted_by {  CA_INT([[0, 1, 2], [3, 4, 5], [6, 7, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[[1..-1,3]] = 9
+    is_asserted_by {  CA_INT([[0, 9, 2], [3, 9, 5], [6, 9, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[[1,3,3]] = 9
+    is_asserted_by {  CA_INT([[0, 9, 2], [3, 9, 5], [6, 9, 8]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[nil] = 9
+    is_asserted_by {  CA_INT([[9, 9, 9], [9, 9, 9], [9, 9, 9]]) == a }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    expect { a[0..9] = -3  }.to raise_error(IndexError)
+    expect { a[-10..-1] = -4  }.to raise_error(IndexError)
+  end
+
+  example "access_address_selection" do
+    # ---
+    a = CArray.int(3,3).seq!;
+    a1 = a.flatten
+    is_asserted_by {  CA_INT([0, 1, 2, 3, 4]) == a[a1 < 5] }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[a1 < 5] = 1
+    is_asserted_by {  CA_INT([[1, 1, 1], [1, 1, 5], [6, 7, 8]]) == a }
+  end
+
+  example "access_address_grid" do
+    # ---
+    a = CArray.int(3,3).seq!;
+    is_asserted_by {  CA_INT([0, 2, 4, 6, 8]) == a[+[0, 2, 4, 6, 8]] }
+
+    # ---
+    a = CArray.int(3,3).seq!
+    a[+[0,2,4,6,8]] = 9
+    is_asserted_by {  CA_INT([[9, 1, 9], [3, 9, 5], [9, 7, 9]]) == a }
+  end
+
+end