carray 1.4.0 → 1.5.5

data/spec/Features/feature_serialization_spec.rb

@@ -0,0 +1,125 @@
+require "carray"
+require "rspec-power_assert"
+
+describe CArray do
+
+  example "dump_io ()" do
+    begin
+      a = CArray.int32(3, 3).seq!
+      b = CArray.int32(3, 3)
+      open("bintest", "w") { |io| a.dump_binary(io) }
+      open("bintest")      { |io| b.load_binary(io) }
+      is_asserted_by { a == b }
+
+      a = CArray.int32(3, 3).seq!
+      b = CArray.int32(3, 3).seq!
+      open("bintest", "w") { |io| a[0..1, 0..1].dump_binary(io) }
+      open("bintest")      { |io| b[0..1, 0..1].load_binary(io) }
+      is_asserted_by { a == b }
+    ensure
+      File.unlink("bintest")
+    end
+  end
+
+  example "dump_str ()" do
+    a = CArray.int32(3, 3).seq!
+    b = CArray.int32(3, 3)
+    s = ""
+    a.dump_binary(s)
+    b.load_binary(s)
+    is_asserted_by { a == b }
+
+    a = CArray.int32(3, 3).seq!
+    b = CArray.int32(3, 3).seq!
+    s = ""
+    a[0..1, 0..1].dump_binary(s)
+    b[0..1, 0..1].load_binary(s)
+    is_asserted_by { a == b }
+  end
+
+
+  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
+
+    example "should equal to original" do
+      is_asserted_by { @it == @original }
+      is_asserted_by { @it.has_mask? == true }
+    end
+
+    example "should dump string of same contents with file from which it was loaded" do 
+      is_asserted_by { CArray.dump(@it) == 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
+
+    example "should equal to original" do
+      is_asserted_by { @it == @original }
+      is_asserted_by { @it.has_mask? == 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
+
+    example "should equal to original" do
+      is_asserted_by { @it == @original }
+      is_asserted_by { @it.has_mask? == 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
+
+    example "should equal to original" do
+      is_asserted_by { @it == @original }
+      is_asserted_by { @it.has_mask? == true }
+    end
+
+    after do
+      File.unlink("test.ca")
+    end
+
+  end
+
+end

data/spec/Features/feature_stat_spec.rb

@@ -0,0 +1,397 @@
+
+require 'carray'
+require "rspec-power_assert"
+
+describe "TestCArrayStat " do
+
+  example "min" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.min
+    #assert_instance_of(Integer, s)
+    is_asserted_by {  1 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.min
+    is_asserted_by { s.class == Float }
+    is_asserted_by { 1 == s }
+
+    # ---
+    if CArray::HAVE_COMPLEX
+      a = CArray.cmplx128(10)
+      a.real.seq!(1)     ### [1, 2, ..., 10]
+      a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      expect { a.min }.to raise_error(CArray::DataTypeError)
+    end
+  end
+
+  example "min_addr" do
+    # ---
+    a = CA_INT32([8,7,6,5,4,5,6,7,8])
+    is_asserted_by {  4 == a.min_addr }
+
+    # ---
+    a = CA_INT32([1,1,1,2,2,2,3,3,3])
+    is_asserted_by {  0 == a.min_addr }
+    is_asserted_by {  6 == a.reverse.min_addr }
+
+    # ---
+    _ = UNDEF
+    a = CA_INT32([_,2,3])
+    is_asserted_by {  1 == a.min_addr }
+    is_asserted_by {  1 == a.min_addr(mask_limit: 0) }
+    is_asserted_by {  UNDEF == a.min_addr(mask_limit: 1) }
+    is_asserted_by {  (-9999) == a.min_addr(mask_limit: 1, fill_value: -9999) }
+    is_asserted_by {  1 == a.min_addr(mask_limit: 2) }
+  end
+
+  example "max" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.max
+    #assert_instance_of(Integer, s)
+    is_asserted_by {  10 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.max
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  10 == s }
+
+    # ---
+    if CArray::HAVE_COMPLEX
+      a = CArray.cmplx128(10)
+      a.real.seq!(1)     ### [1, 2, ..., 10]
+      a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      expect { a.max }.to raise_error(CArray::DataTypeError)
+    end
+  end
+
+  example "max_addr" do
+    # ---
+    a = CA_INT32([0,1,2,3,4,3,2,1,0])
+    is_asserted_by {  4 == a.max_addr }
+
+    # ---
+    a = CA_INT32([1,1,1,2,2,2,3,3,3])
+    is_asserted_by {  6 == a.max_addr }
+    is_asserted_by {  0 == a.reverse.max_addr }
+
+    # ---
+    _ = UNDEF
+    a = CA_INT32([1,2,_])
+    is_asserted_by {  1 == a.max_addr }
+    is_asserted_by {  1 == a.max_addr(mask_limit: 0) }
+    is_asserted_by {  UNDEF == a.max_addr(mask_limit: 1) }
+    is_asserted_by {  (-9999) == a.max_addr(mask_limit: 1, fill_value: -9999) }
+    is_asserted_by {  1 == a.max_addr(mask_limit: 2) }
+  end
+
+#  example "max_and_min_addr" do
+#    a = CArray.float(100).span!(1..100).shuffle
+#    addr1 = a.max_addr
+#    addr2 = (-a).min_addr
+#    is_asserted_by { addr1, addr2)
+#    is_asserted_by { 100, a[addr1])
+#  end
+
+  example "sum" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.sum
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  55 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.sum
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  55 == s }
+
+  end
+
+  example "prod" do
+    # ---
+    a = CArray.int32(5).seq!(1) ### [1, 2, ..., 5]
+    s = a.prod
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  120 == s }
+
+    # ---
+    a = CArray.float64(5).seq!(1) ### [1, 2, ..., 5]
+    s = a.prod
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  120 == s }
+
+  end
+
+  example "mean" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.mean
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  5.5 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.mean
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  5.5 == s }
+
+  end
+
+  example "wmean" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    w = a.reverse
+    s = a.wmean(w)
+    is_asserted_by {  4 == s }
+
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    w = a.reverse
+    a[0] = UNDEF
+    w[9] = UNDEF
+    s = a.wmean(w)
+    is_asserted_by {  (200.0 / 44) == s }
+  end
+
+  example "cumsum" do
+    # ---
+    a = CArray.int32(10).seq!(1)
+    is_asserted_by {  CA_DOUBLE([1, 3, 6, 10, 15, 21, 28, 36, 45, 55]) == a.cumsum }
+  end
+
+  example "cumprod" do
+    # ---
+    a = CArray.int32(5).seq!(1)
+    is_asserted_by {  CA_DOUBLE([1, 2, 6, 24, 120]) == a.cumprod }
+  end
+
+  example "variancep" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.variancep
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  8.25 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.variancep
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  8.25 == s }
+
+    # ---
+    if CArray::HAVE_COMPLEX
+      a = CArray.cmplx128(10)
+      a.real.seq!(1)     ### [1, 2, ..., 10]
+      a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      expect { a.variancep }.to raise_error(CArray::DataTypeError)
+    end
+  end
+
+  example "variance" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.variance
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  true == ((9.1666667 - s).abs < 1.0e-05) }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.variance
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  true == ((9.1666667 - s).abs < 1.0e-05) }
+
+    # ---
+    if CArray::HAVE_COMPLEX
+      a = CArray.cmplx128(10)
+      a.real.seq!(1)     ### [1, 2, ..., 10]
+      a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      expect { s = a.variance }.to raise_error(CArray::DataTypeError)
+    end
+  end
+
+  example "wsum" do
+    # ---
+    a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    w = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.wsum(w)
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  385 == s }
+
+    # ---
+    a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    w = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
+    s = a.wsum(w)
+    is_asserted_by { s.class == Float }
+    is_asserted_by {  385 == s }
+
+    # ---
+    if CArray::HAVE_COMPLEX
+      a = CArray.cmplx128(10)
+      w = CArray.cmplx128(10)
+      a.real.seq!(1)     ### [1, 2, ..., 10]
+      a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      w.real.seq!(1)     ### [1, 2, ..., 10]
+      w.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
+      s = a.wsum(w.conj)
+      is_asserted_by { s.class == CComplex }
+      is_asserted_by {  770 == s.real }
+      is_asserted_by {  0 == s.imag }
+    end
+  end
+
+  example "accumulate" do
+    # ---
+    a = CArray.uint8(255) {1}
+    b = CArray.uint8(256) {1}
+    is_asserted_by {  255 == a.accumulate }
+    is_asserted_by {  0 == b.accumulate }
+
+    # ---
+    a = CArray.int8(127) {1}
+    b = CArray.int8(128) {1}
+    c = CArray.int8(255) {1}
+    d = CArray.int8(256) {1}
+    is_asserted_by {  127 == a.accumulate }
+    is_asserted_by {  (-128) == b.accumulate }
+    is_asserted_by {  (-1) == c.accumulate }
+    is_asserted_by {  0 == d.accumulate }
+  end
+
+  example "count_equal" do
+    # ---
+    a = CArray.int32(10)
+    a[0] = 0
+    a[1..3] = 1
+    a[4..7] = 2
+    a[8..9] = 3          ### [0,1,1,1,2,2,2,2,3,3]
+
+    is_asserted_by {  1 == a.count_equal(0) }
+    is_asserted_by {  3 == a.count_equal(1) }
+    is_asserted_by {  4 == a.count_equal(2) }
+    is_asserted_by {  2 == a.count_equal(3) }
+    is_asserted_by {  0 == a.count_equal(4) }
+  end
+
+  example "count_equiv" do
+    # ---
+    a = CArray.int32(10)
+    a[0] = 0
+    a[1..3] = 1
+    a[4..7] = 2
+    a[8..9] = 3          ### [0,1,1,1,2,2,2,2,3,3]
+
+    is_asserted_by {  1 == a.count_equiv(0, 0.001) }
+    is_asserted_by {  3 == a.count_equiv(1, 0.001) }
+    is_asserted_by {  4 == a.count_equiv(2, 0.001) }
+    is_asserted_by {  2 == a.count_equiv(3, 0.001) }
+    is_asserted_by {  0 == a.count_equiv(4, 0.001) }
+
+    # ---
+    a = CArray.float64(10)
+    a[0] = 0.00005
+    a[1..3] = 1.00005
+    a[4..7] = 2.0001
+    a[8..9] = 3.00015          ### [0,1,1,1,2,2,2,2,3,3]
+
+    is_asserted_by {  0 == a.count_equiv(0, 1.0e-05) }
+    is_asserted_by {  0 == a.count_equiv(1, 1.0e-05) }
+    is_asserted_by {  0 == a.count_equiv(2, 1.0e-05) }
+    is_asserted_by {  0 == a.count_equiv(3, 1.0e-05) }
+    is_asserted_by {  0 == a.count_equiv(4, 1.0e-05) }
+
+    is_asserted_by {  0 == a.count_equiv(0, 0.0001) }
+
+    is_asserted_by {  3 == a.count_equiv(1, 0.0001) }
+    is_asserted_by {  4 == a.count_equiv(2, 0.0001) }
+    is_asserted_by {  2 == a.count_equiv(3, 0.0001) }
+    is_asserted_by {  0 == a.count_equiv(4, 0.0001) }
+  end
+
+  example "count_close" do
+    # ---
+    a = CArray.int32(10)
+    a[0] = 0
+    a[1..3] = 1
+    a[4..7] = 2
+    a[8..9] = 3          ### [0,1,1,1,2,2,2,2,3,3]
+
+    is_asserted_by {  1 == a.count_close(0, 0.001) }
+    is_asserted_by {  3 == a.count_close(1, 0.001) }
+    is_asserted_by {  4 == a.count_close(2, 0.001) }
+    is_asserted_by {  2 == a.count_close(3, 0.001) }
+    is_asserted_by {  0 == a.count_close(4, 0.001) }
+
+    # ---
+    a = CArray.float64(10)
+    a[0] = 0.00005
+    a[1..3] = 1.00005
+    a[4..7] = 2.00005
+    a[8..9] = 3.00005          ### [0,1,1,1,2,2,2,2,3,3]
+
+    is_asserted_by {  0 == a.count_close(0, 1.0e-05) }
+    is_asserted_by {  0 == a.count_close(1, 1.0e-05) }
+    is_asserted_by {  0 == a.count_close(2, 1.0e-05) }
+    is_asserted_by {  0 == a.count_close(3, 1.0e-05) }
+    is_asserted_by {  0 == a.count_close(4, 1.0e-05) }
+
+    is_asserted_by {  1 == a.count_close(0, 0.0001) }
+    is_asserted_by {  3 == a.count_close(1, 0.0001) }
+    is_asserted_by {  4 == a.count_close(2, 0.0001) }
+    is_asserted_by {  2 == a.count_close(3, 0.0001) }
+    is_asserted_by {  0 == a.count_close(4, 0.0001) }
+  end
+
+  example "all_equal" do
+    #---
+    a = CArray.int32(3,3) {1}
+    is_asserted_by {  false == a.all_equal?(0) }
+    is_asserted_by {  true == a.all_equal?(1) }
+  end
+
+  example "all_equiv" do
+    #---
+    a = CArray.float64(3,3) {1.00005}
+    is_asserted_by {  false == a.all_equiv?(0, 0.0001) }
+    is_asserted_by {  true == a.all_equiv?(1, 0.0001) }
+    is_asserted_by {  false == a.all_equiv?(1, 1.0e-05) }
+  end
+
+  example "all_close" do
+    #---
+    a = CArray.float64(3,3) {1.0001}
+    is_asserted_by {  false == a.all_close?(0, 0.0001) }
+    is_asserted_by {  true == a.all_close?(1, 0.0001) }
+    is_asserted_by {  false == a.all_close?(1, 1.0e-05) }
+  end
+
+  example "any_equal" do
+    #---
+    a = CArray.int32(3,3)
+    a[1,1] = 1
+    is_asserted_by {  true == a.any_equal?(0) }
+    is_asserted_by {  true == a.any_equal?(1) }
+  end
+
+  example "any_equiv" do
+    #---
+    a = CArray.float64(3,3) { 0.00005 }
+    a[1,1] = 1.00005
+    is_asserted_by {  false == a.any_equiv?(0, 0.0001) }
+    is_asserted_by {  true == a.any_equiv?(1, 0.0001) }
+    is_asserted_by {  false == a.any_equiv?(1, 1.0e-05) }
+  end
+
+  example "any_close" do
+    #---
+    a = CArray.float64(3,3) { 0.0001 }
+    a[1,1] = 1.0001
+    is_asserted_by {  true == a.any_close?(0, 0.0001) }
+    is_asserted_by {  true == a.any_close?(1, 0.0001) }
+    is_asserted_by {  false == a.any_close?(1, 1.0e-05) }
+  end
+
+end