carray 1.4.0 → 1.5.5

Sign up to get free protection for your applications and to get access to all the features.
Files changed (239) hide show
  1. checksums.yaml +4 -4
  2. data/.yardopts +32 -0
  3. data/LICENSE +22 -0
  4. data/NEWS.md +83 -0
  5. data/README.md +38 -18
  6. data/Rakefile +2 -2
  7. data/TODO.md +17 -0
  8. data/carray.gemspec +13 -18
  9. data/{ca_iter_block.c → ext/ca_iter_block.c} +16 -18
  10. data/{ca_iter_dimension.c → ext/ca_iter_dimension.c} +20 -21
  11. data/{ca_iter_window.c → ext/ca_iter_window.c} +12 -14
  12. data/{ca_obj_array.c → ext/ca_obj_array.c} +451 -176
  13. data/{ca_obj_bitarray.c → ext/ca_obj_bitarray.c} +18 -23
  14. data/{ca_obj_bitfield.c → ext/ca_obj_bitfield.c} +12 -16
  15. data/{ca_obj_block.c → ext/ca_obj_block.c} +47 -54
  16. data/{ca_obj_fake.c → ext/ca_obj_fake.c} +10 -12
  17. data/{ca_obj_farray.c → ext/ca_obj_farray.c} +21 -23
  18. data/{ca_obj_field.c → ext/ca_obj_field.c} +30 -32
  19. data/{ca_obj_grid.c → ext/ca_obj_grid.c} +32 -33
  20. data/{ca_obj_mapping.c → ext/ca_obj_mapping.c} +11 -13
  21. data/{ca_obj_object.c → ext/ca_obj_object.c} +40 -42
  22. data/{ca_obj_reduce.c → ext/ca_obj_reduce.c} +3 -5
  23. data/{ca_obj_refer.c → ext/ca_obj_refer.c} +44 -48
  24. data/{ca_obj_repeat.c → ext/ca_obj_repeat.c} +45 -47
  25. data/{ca_obj_select.c → ext/ca_obj_select.c} +4 -6
  26. data/{ca_obj_shift.c → ext/ca_obj_shift.c} +26 -28
  27. data/{ca_obj_transpose.c → ext/ca_obj_transpose.c} +26 -28
  28. data/{ca_obj_unbound_repeat.c → ext/ca_obj_unbound_repeat.c} +106 -160
  29. data/{ca_obj_window.c → ext/ca_obj_window.c} +33 -35
  30. data/{carray.h → ext/carray.h} +86 -59
  31. data/{carray_access.c → ext/carray_access.c} +194 -101
  32. data/{carray_attribute.c → ext/carray_attribute.c} +161 -207
  33. data/{carray_call_cfunc.c → ext/carray_call_cfunc.c} +1 -3
  34. data/{carray_cast.c → ext/carray_cast.c} +351 -150
  35. data/{carray_cast_func.rb → ext/carray_cast_func.rb} +1 -2
  36. data/{carray_class.c → ext/carray_class.c} +28 -36
  37. data/{carray_conversion.c → ext/carray_conversion.c} +63 -68
  38. data/{carray_copy.c → ext/carray_copy.c} +34 -50
  39. data/{carray_core.c → ext/carray_core.c} +75 -62
  40. data/ext/carray_data_type.c +66 -0
  41. data/{carray_element.c → ext/carray_element.c} +34 -53
  42. data/{carray_generate.c → ext/carray_generate.c} +71 -50
  43. data/{carray_iterator.c → ext/carray_iterator.c} +53 -53
  44. data/{carray_loop.c → ext/carray_loop.c} +77 -106
  45. data/{carray_mask.c → ext/carray_mask.c} +105 -114
  46. data/{carray_math.rb → ext/carray_math.rb} +29 -13
  47. data/ext/{mathfunc/carray_mathfunc.c → carray_mathfunc.c} +1 -3
  48. data/{carray_numeric.c → ext/carray_numeric.c} +43 -46
  49. data/{carray_operator.c → ext/carray_operator.c} +49 -36
  50. data/{carray_order.c → ext/carray_order.c} +232 -217
  51. data/{carray_sort_addr.c → ext/carray_sort_addr.c} +14 -21
  52. data/{carray_stat.c → ext/carray_stat.c} +6 -8
  53. data/{carray_stat_proc.rb → ext/carray_stat_proc.rb} +25 -27
  54. data/{carray_test.c → ext/carray_test.c} +63 -51
  55. data/{carray_undef.c → ext/carray_undef.c} +1 -11
  56. data/{carray_utils.c → ext/carray_utils.c} +12 -4
  57. data/{extconf.rb → ext/extconf.rb} +10 -7
  58. data/{mkmath.rb → ext/mkmath.rb} +2 -2
  59. data/{ruby_carray.c → ext/ruby_carray.c} +20 -8
  60. data/{ruby_ccomplex.c → ext/ruby_ccomplex.c} +2 -4
  61. data/{ruby_float_func.c → ext/ruby_float_func.c} +1 -3
  62. data/ext/version.h +16 -0
  63. data/{version.rb → ext/version.rb} +0 -0
  64. data/lib/carray.rb +51 -40
  65. data/lib/carray/{base/autoload.rb → autoload.rb} +1 -3
  66. data/lib/carray/autoload/autoload_base.rb +1 -1
  67. data/lib/carray/autoload/autoload_gem_numo_narray.rb +7 -6
  68. data/lib/carray/autoload/autoload_gem_random.rb +8 -0
  69. data/lib/carray/basic.rb +191 -0
  70. data/lib/carray/broadcast.rb +101 -0
  71. data/lib/carray/compose.rb +315 -0
  72. data/lib/carray/construct.rb +484 -0
  73. data/lib/carray/convert.rb +115 -0
  74. data/lib/carray/info.rb +1 -3
  75. data/lib/carray/{base/inspect.rb → inspect.rb} +9 -11
  76. data/lib/carray/io/imagemagick.rb +2 -4
  77. data/lib/carray/{base/iterator.rb → iterator.rb} +6 -6
  78. data/lib/carray/mask.rb +102 -0
  79. data/lib/carray/{base/math.rb → math.rb} +20 -52
  80. data/lib/carray/math/histogram.rb +8 -10
  81. data/lib/carray/math/recurrence.rb +1 -3
  82. data/lib/carray/mkmf.rb +9 -3
  83. data/lib/carray/object/ca_obj_iterator.rb +1 -3
  84. data/lib/carray/object/ca_obj_link.rb +1 -3
  85. data/lib/carray/object/ca_obj_pack.rb +9 -11
  86. data/lib/carray/obsolete.rb +256 -0
  87. data/lib/carray/ordering.rb +181 -0
  88. data/lib/carray/{base/serialize.rb → serialize.rb} +60 -76
  89. data/lib/carray/{base/string.rb → string.rb} +10 -64
  90. data/lib/carray/{base/struct.rb → struct.rb} +19 -21
  91. data/lib/carray/{io/table.rb → table.rb} +1 -10
  92. data/lib/carray/testing.rb +51 -0
  93. data/lib/carray/time.rb +76 -0
  94. data/lib/carray/transform.rb +109 -0
  95. data/misc/Methods.ja.md +182 -0
  96. data/{NOTE → misc/NOTE} +16 -38
  97. data/spec/Classes/CABitfield_spec.rb +58 -0
  98. data/spec/Classes/CABlockIterator_spec.rb +114 -0
  99. data/spec/Classes/CABlock_spec.rb +205 -0
  100. data/spec/Classes/CAField_spec.rb +39 -0
  101. data/spec/Classes/CAGrid_spec.rb +75 -0
  102. data/spec/Classes/CAMap_spec.rb +0 -0
  103. data/{test/test_CAMapping.rb → spec/Classes/CAMapping_spec.rb} +35 -36
  104. data/spec/Classes/CAObject_attribute_spec.rb +33 -0
  105. data/spec/Classes/CAObject_spec.rb +33 -0
  106. data/spec/Classes/CARefer_spec.rb +93 -0
  107. data/spec/Classes/CARepeat_spec.rb +65 -0
  108. data/spec/Classes/CASelect_spec.rb +22 -0
  109. data/spec/Classes/CAShift_spec.rb +16 -0
  110. data/spec/Classes/CAStruct_spec.rb +71 -0
  111. data/{test/test_CATranspose.rb → spec/Classes/CATranspose_spec.rb} +20 -21
  112. data/spec/Classes/CAUnboudRepeat_spec.rb +102 -0
  113. data/spec/Classes/CAWindow_spec.rb +54 -0
  114. data/spec/Classes/CAWrap_spec.rb +8 -0
  115. data/{test/test_CArray.rb → spec/Classes/CArray_spec.rb} +48 -92
  116. data/spec/Classes/CScalar_spec.rb +55 -0
  117. data/spec/Features/feature_130_spec.rb +19 -0
  118. data/spec/Features/feature_attributes_spec.rb +280 -0
  119. data/spec/Features/feature_boolean_spec.rb +98 -0
  120. data/spec/Features/feature_broadcast.rb +116 -0
  121. data/spec/Features/feature_cast_function.rb +19 -0
  122. data/spec/Features/feature_cast_spec.rb +33 -0
  123. data/spec/Features/feature_class_spec.rb +84 -0
  124. data/spec/Features/feature_complex_spec.rb +42 -0
  125. data/{test/test_composite.rb → spec/Features/feature_composite_spec.rb} +17 -18
  126. data/spec/Features/feature_convert_spec.rb +46 -0
  127. data/spec/Features/feature_copy_spec.rb +123 -0
  128. data/spec/Features/feature_creation_spec.rb +84 -0
  129. data/spec/Features/feature_element_spec.rb +144 -0
  130. data/spec/Features/feature_extream_spec.rb +54 -0
  131. data/spec/Features/feature_generate_spec.rb +74 -0
  132. data/spec/Features/feature_index_spec.rb +69 -0
  133. data/spec/Features/feature_mask_spec.rb +580 -0
  134. data/spec/Features/feature_math_spec.rb +97 -0
  135. data/spec/Features/feature_order_spec.rb +146 -0
  136. data/spec/Features/feature_ref_store_spec.rb +209 -0
  137. data/spec/Features/feature_serialization_spec.rb +125 -0
  138. data/spec/Features/feature_stat_spec.rb +397 -0
  139. data/spec/Features/feature_virtual_spec.rb +48 -0
  140. data/spec/Features/method_eq_spec.rb +81 -0
  141. data/spec/Features/method_is_nan_spec.rb +12 -0
  142. data/spec/Features/method_map_spec.rb +54 -0
  143. data/spec/Features/method_max_with.rb +20 -0
  144. data/spec/Features/method_min_with.rb +19 -0
  145. data/spec/Features/method_ne_spec.rb +18 -0
  146. data/spec/Features/method_project_spec.rb +188 -0
  147. data/spec/Features/method_ref_spec.rb +27 -0
  148. data/spec/Features/method_round_spec.rb +11 -0
  149. data/spec/Features/method_s_linspace_spec.rb +48 -0
  150. data/spec/Features/method_s_span_spec.rb +14 -0
  151. data/spec/Features/method_seq_spec.rb +47 -0
  152. data/spec/Features/method_sort_with.rb +43 -0
  153. data/spec/Features/method_sorted_with.rb +29 -0
  154. data/spec/Features/method_span_spec.rb +42 -0
  155. data/spec/Features/method_wrap_readonly_spec.rb +43 -0
  156. data/{test → spec/UnitTest}/test_CAVirtual.rb +0 -0
  157. data/spec/spec_all.rb +0 -1
  158. data/utils/convert_test.rb +73 -0
  159. data/utils/{extract_rdoc.rb → extract_yard.rb} +7 -12
  160. data/{devel → utils}/guess_shape.rb +0 -0
  161. data/utils/{diff_method.rb → monkey_patch_methods.rb} +17 -7
  162. metadata +159 -206
  163. data/COPYING +0 -56
  164. data/GPL +0 -340
  165. data/LEGAL +0 -50
  166. data/TODO +0 -5
  167. data/carray_random.c +0 -531
  168. data/devel/im2col.rb +0 -17
  169. data/ext/calculus/carray_calculus.c +0 -931
  170. data/ext/calculus/carray_interp.c +0 -358
  171. data/ext/calculus/extconf.rb +0 -12
  172. data/ext/calculus/lib/math/calculus.rb +0 -119
  173. data/ext/calculus/lib/math/interp/adapter_interp1d.rb +0 -31
  174. data/ext/mathfunc/extconf.rb +0 -18
  175. data/ext/mathfunc/test/test_hypot.rb +0 -5
  176. data/ext/mathfunc/test/test_j0.rb +0 -22
  177. data/ext/mathfunc/test/test_jn.rb +0 -8
  178. data/ext/mathfunc/test/test_sph.rb +0 -9
  179. data/lib/carray/autoload/autoload_io_table.rb +0 -1
  180. data/lib/carray/autoload/autoload_math_interp.rb +0 -4
  181. data/lib/carray/base/basic.rb +0 -1146
  182. data/lib/carray/base/obsolete.rb +0 -131
  183. data/lib/carray/math/interp.rb +0 -57
  184. data/lib/carray/math/interp/adapter_gsl_spline.rb +0 -47
  185. data/mt19937ar.c +0 -182
  186. data/mt19937ar.h +0 -86
  187. data/rdoc_main.rb +0 -27
  188. data/rdoc_math.rb +0 -5
  189. data/rdoc_stat.rb +0 -31
  190. data/spec/CABlockIterator/CABlockIterator_spec.rb +0 -113
  191. data/spec/CArray/bug/store_spec.rb +0 -27
  192. data/spec/CArray/index/repeat_spec.rb +0 -10
  193. data/spec/CArray/method/eq_spec.rb +0 -80
  194. data/spec/CArray/method/is_nan_spec.rb +0 -12
  195. data/spec/CArray/method/ne_spec.rb +0 -18
  196. data/spec/CArray/method/round_spec.rb +0 -11
  197. data/spec/CArray/object/_attribute_spec.rb +0 -32
  198. data/spec/CArray/object/s_new_spec.rb +0 -31
  199. data/spec/CArray/serialize/Serialization_spec.rb +0 -89
  200. data/test/test_130.rb +0 -23
  201. data/test/test_ALL.rb +0 -51
  202. data/test/test_CABitfield.rb +0 -59
  203. data/test/test_CABlock.rb +0 -208
  204. data/test/test_CAField.rb +0 -40
  205. data/test/test_CAGrid.rb +0 -76
  206. data/test/test_CAMmap.rb +0 -11
  207. data/test/test_CARefer.rb +0 -94
  208. data/test/test_CARepeat.rb +0 -66
  209. data/test/test_CASelect.rb +0 -23
  210. data/test/test_CAShift.rb +0 -17
  211. data/test/test_CAWindow.rb +0 -55
  212. data/test/test_CAWrap.rb +0 -9
  213. data/test/test_CComplex.rb +0 -83
  214. data/test/test_CScalar.rb +0 -91
  215. data/test/test_attribute.rb +0 -281
  216. data/test/test_block_iterator.rb +0 -17
  217. data/test/test_boolean.rb +0 -99
  218. data/test/test_cast.rb +0 -33
  219. data/test/test_class.rb +0 -85
  220. data/test/test_complex.rb +0 -43
  221. data/test/test_convert.rb +0 -79
  222. data/test/test_copy.rb +0 -141
  223. data/test/test_creation.rb +0 -85
  224. data/test/test_element.rb +0 -146
  225. data/test/test_extream.rb +0 -55
  226. data/test/test_generate.rb +0 -75
  227. data/test/test_index.rb +0 -71
  228. data/test/test_mask.rb +0 -578
  229. data/test/test_math.rb +0 -98
  230. data/test/test_narray.rb +0 -64
  231. data/test/test_order.rb +0 -147
  232. data/test/test_random.rb +0 -15
  233. data/test/test_ref_store.rb +0 -211
  234. data/test/test_stat.rb +0 -414
  235. data/test/test_struct.rb +0 -72
  236. data/test/test_virtual.rb +0 -49
  237. data/utils/create_rdoc.sh +0 -9
  238. data/utils/make_tgz.sh +0 -3
  239. data/version.h +0 -18
@@ -0,0 +1,125 @@
1
+ require "carray"
2
+ require "rspec-power_assert"
3
+
4
+ describe CArray do
5
+
6
+ example "dump_io ()" do
7
+ begin
8
+ a = CArray.int32(3, 3).seq!
9
+ b = CArray.int32(3, 3)
10
+ open("bintest", "w") { |io| a.dump_binary(io) }
11
+ open("bintest") { |io| b.load_binary(io) }
12
+ is_asserted_by { a == b }
13
+
14
+ a = CArray.int32(3, 3).seq!
15
+ b = CArray.int32(3, 3).seq!
16
+ open("bintest", "w") { |io| a[0..1, 0..1].dump_binary(io) }
17
+ open("bintest") { |io| b[0..1, 0..1].load_binary(io) }
18
+ is_asserted_by { a == b }
19
+ ensure
20
+ File.unlink("bintest")
21
+ end
22
+ end
23
+
24
+ example "dump_str ()" do
25
+ a = CArray.int32(3, 3).seq!
26
+ b = CArray.int32(3, 3)
27
+ s = ""
28
+ a.dump_binary(s)
29
+ b.load_binary(s)
30
+ is_asserted_by { a == b }
31
+
32
+ a = CArray.int32(3, 3).seq!
33
+ b = CArray.int32(3, 3).seq!
34
+ s = ""
35
+ a[0..1, 0..1].dump_binary(s)
36
+ b[0..1, 0..1].load_binary(s)
37
+ is_asserted_by { a == b }
38
+ end
39
+
40
+
41
+ describe "loaded from binary format (int32 type)" do
42
+
43
+ before do
44
+ @original = CArray.int32(10,10) { 10 }
45
+ @original[5,5] = UNDEF
46
+ CArray.save(@original, "test.ca")
47
+ @it = CArray.load("test.ca")
48
+ end
49
+
50
+ example "should equal to original" do
51
+ is_asserted_by { @it == @original }
52
+ is_asserted_by { @it.has_mask? == true }
53
+ end
54
+
55
+ example "should dump string of same contents with file from which it was loaded" do
56
+ is_asserted_by { CArray.dump(@it) == File.read("test.ca") }
57
+ end
58
+
59
+ after do
60
+ File.unlink("test.ca")
61
+ end
62
+
63
+ end
64
+
65
+ describe "loaded from binary format (object type)" do
66
+
67
+ before do
68
+ @original = CArray.object(10,10) { 3.times { Time.now } }
69
+ @original[5,5] = UNDEF
70
+ CArray.save(@original, "test.ca")
71
+ @it = CArray.load("test.ca")
72
+ end
73
+
74
+ example "should equal to original" do
75
+ is_asserted_by { @it == @original }
76
+ is_asserted_by { @it.has_mask? == true }
77
+ end
78
+
79
+ after do
80
+ File.unlink("test.ca")
81
+ end
82
+
83
+ end
84
+
85
+ describe "loaded by marshal (int32 type)" do
86
+
87
+ before do
88
+ @original = CArray.int32(10,10) { 10 }
89
+ @original[5,5] = UNDEF
90
+ open("test.ca", "w") { |io| Marshal.dump(@original, io) }
91
+ @it = open("test.ca") { |io| Marshal.load(io) }
92
+ end
93
+
94
+ example "should equal to original" do
95
+ is_asserted_by { @it == @original }
96
+ is_asserted_by { @it.has_mask? == true }
97
+ end
98
+
99
+ after do
100
+ File.unlink("test.ca")
101
+ end
102
+
103
+ end
104
+
105
+ describe "loaded by marshal (object type)" do
106
+
107
+ before do
108
+ @original = CArray.object(10,10) { 3.times { Time.now } }
109
+ @original[5,5] = UNDEF
110
+ open("test.ca", "w") { |io| Marshal.dump(@original, io) }
111
+ @it = open("test.ca") { |io| Marshal.load(io) }
112
+ end
113
+
114
+ example "should equal to original" do
115
+ is_asserted_by { @it == @original }
116
+ is_asserted_by { @it.has_mask? == true }
117
+ end
118
+
119
+ after do
120
+ File.unlink("test.ca")
121
+ end
122
+
123
+ end
124
+
125
+ end
@@ -0,0 +1,397 @@
1
+
2
+ require 'carray'
3
+ require "rspec-power_assert"
4
+
5
+ describe "TestCArrayStat " do
6
+
7
+ example "min" do
8
+ # ---
9
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
10
+ s = a.min
11
+ #assert_instance_of(Integer, s)
12
+ is_asserted_by { 1 == s }
13
+
14
+ # ---
15
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
16
+ s = a.min
17
+ is_asserted_by { s.class == Float }
18
+ is_asserted_by { 1 == s }
19
+
20
+ # ---
21
+ if CArray::HAVE_COMPLEX
22
+ a = CArray.cmplx128(10)
23
+ a.real.seq!(1) ### [1, 2, ..., 10]
24
+ a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
25
+ expect { a.min }.to raise_error(CArray::DataTypeError)
26
+ end
27
+ end
28
+
29
+ example "min_addr" do
30
+ # ---
31
+ a = CA_INT32([8,7,6,5,4,5,6,7,8])
32
+ is_asserted_by { 4 == a.min_addr }
33
+
34
+ # ---
35
+ a = CA_INT32([1,1,1,2,2,2,3,3,3])
36
+ is_asserted_by { 0 == a.min_addr }
37
+ is_asserted_by { 6 == a.reverse.min_addr }
38
+
39
+ # ---
40
+ _ = UNDEF
41
+ a = CA_INT32([_,2,3])
42
+ is_asserted_by { 1 == a.min_addr }
43
+ is_asserted_by { 1 == a.min_addr(mask_limit: 0) }
44
+ is_asserted_by { UNDEF == a.min_addr(mask_limit: 1) }
45
+ is_asserted_by { (-9999) == a.min_addr(mask_limit: 1, fill_value: -9999) }
46
+ is_asserted_by { 1 == a.min_addr(mask_limit: 2) }
47
+ end
48
+
49
+ example "max" do
50
+ # ---
51
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
52
+ s = a.max
53
+ #assert_instance_of(Integer, s)
54
+ is_asserted_by { 10 == s }
55
+
56
+ # ---
57
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
58
+ s = a.max
59
+ is_asserted_by { s.class == Float }
60
+ is_asserted_by { 10 == s }
61
+
62
+ # ---
63
+ if CArray::HAVE_COMPLEX
64
+ a = CArray.cmplx128(10)
65
+ a.real.seq!(1) ### [1, 2, ..., 10]
66
+ a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
67
+ expect { a.max }.to raise_error(CArray::DataTypeError)
68
+ end
69
+ end
70
+
71
+ example "max_addr" do
72
+ # ---
73
+ a = CA_INT32([0,1,2,3,4,3,2,1,0])
74
+ is_asserted_by { 4 == a.max_addr }
75
+
76
+ # ---
77
+ a = CA_INT32([1,1,1,2,2,2,3,3,3])
78
+ is_asserted_by { 6 == a.max_addr }
79
+ is_asserted_by { 0 == a.reverse.max_addr }
80
+
81
+ # ---
82
+ _ = UNDEF
83
+ a = CA_INT32([1,2,_])
84
+ is_asserted_by { 1 == a.max_addr }
85
+ is_asserted_by { 1 == a.max_addr(mask_limit: 0) }
86
+ is_asserted_by { UNDEF == a.max_addr(mask_limit: 1) }
87
+ is_asserted_by { (-9999) == a.max_addr(mask_limit: 1, fill_value: -9999) }
88
+ is_asserted_by { 1 == a.max_addr(mask_limit: 2) }
89
+ end
90
+
91
+ # example "max_and_min_addr" do
92
+ # a = CArray.float(100).span!(1..100).shuffle
93
+ # addr1 = a.max_addr
94
+ # addr2 = (-a).min_addr
95
+ # is_asserted_by { addr1, addr2)
96
+ # is_asserted_by { 100, a[addr1])
97
+ # end
98
+
99
+ example "sum" do
100
+ # ---
101
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
102
+ s = a.sum
103
+ is_asserted_by { s.class == Float }
104
+ is_asserted_by { 55 == s }
105
+
106
+ # ---
107
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
108
+ s = a.sum
109
+ is_asserted_by { s.class == Float }
110
+ is_asserted_by { 55 == s }
111
+
112
+ end
113
+
114
+ example "prod" do
115
+ # ---
116
+ a = CArray.int32(5).seq!(1) ### [1, 2, ..., 5]
117
+ s = a.prod
118
+ is_asserted_by { s.class == Float }
119
+ is_asserted_by { 120 == s }
120
+
121
+ # ---
122
+ a = CArray.float64(5).seq!(1) ### [1, 2, ..., 5]
123
+ s = a.prod
124
+ is_asserted_by { s.class == Float }
125
+ is_asserted_by { 120 == s }
126
+
127
+ end
128
+
129
+ example "mean" do
130
+ # ---
131
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
132
+ s = a.mean
133
+ is_asserted_by { s.class == Float }
134
+ is_asserted_by { 5.5 == s }
135
+
136
+ # ---
137
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
138
+ s = a.mean
139
+ is_asserted_by { s.class == Float }
140
+ is_asserted_by { 5.5 == s }
141
+
142
+ end
143
+
144
+ example "wmean" do
145
+ # ---
146
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
147
+ w = a.reverse
148
+ s = a.wmean(w)
149
+ is_asserted_by { 4 == s }
150
+
151
+ # ---
152
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
153
+ w = a.reverse
154
+ a[0] = UNDEF
155
+ w[9] = UNDEF
156
+ s = a.wmean(w)
157
+ is_asserted_by { (200.0 / 44) == s }
158
+ end
159
+
160
+ example "cumsum" do
161
+ # ---
162
+ a = CArray.int32(10).seq!(1)
163
+ is_asserted_by { CA_DOUBLE([1, 3, 6, 10, 15, 21, 28, 36, 45, 55]) == a.cumsum }
164
+ end
165
+
166
+ example "cumprod" do
167
+ # ---
168
+ a = CArray.int32(5).seq!(1)
169
+ is_asserted_by { CA_DOUBLE([1, 2, 6, 24, 120]) == a.cumprod }
170
+ end
171
+
172
+ example "variancep" do
173
+ # ---
174
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
175
+ s = a.variancep
176
+ is_asserted_by { s.class == Float }
177
+ is_asserted_by { 8.25 == s }
178
+
179
+ # ---
180
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
181
+ s = a.variancep
182
+ is_asserted_by { s.class == Float }
183
+ is_asserted_by { 8.25 == s }
184
+
185
+ # ---
186
+ if CArray::HAVE_COMPLEX
187
+ a = CArray.cmplx128(10)
188
+ a.real.seq!(1) ### [1, 2, ..., 10]
189
+ a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
190
+ expect { a.variancep }.to raise_error(CArray::DataTypeError)
191
+ end
192
+ end
193
+
194
+ example "variance" do
195
+ # ---
196
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
197
+ s = a.variance
198
+ is_asserted_by { s.class == Float }
199
+ is_asserted_by { true == ((9.1666667 - s).abs < 1.0e-05) }
200
+
201
+ # ---
202
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
203
+ s = a.variance
204
+ is_asserted_by { s.class == Float }
205
+ is_asserted_by { true == ((9.1666667 - s).abs < 1.0e-05) }
206
+
207
+ # ---
208
+ if CArray::HAVE_COMPLEX
209
+ a = CArray.cmplx128(10)
210
+ a.real.seq!(1) ### [1, 2, ..., 10]
211
+ a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
212
+ expect { s = a.variance }.to raise_error(CArray::DataTypeError)
213
+ end
214
+ end
215
+
216
+ example "wsum" do
217
+ # ---
218
+ a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
219
+ w = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
220
+ s = a.wsum(w)
221
+ is_asserted_by { s.class == Float }
222
+ is_asserted_by { 385 == s }
223
+
224
+ # ---
225
+ a = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
226
+ w = CArray.float64(10).seq!(1) ### [1, 2, ..., 10]
227
+ s = a.wsum(w)
228
+ is_asserted_by { s.class == Float }
229
+ is_asserted_by { 385 == s }
230
+
231
+ # ---
232
+ if CArray::HAVE_COMPLEX
233
+ a = CArray.cmplx128(10)
234
+ w = CArray.cmplx128(10)
235
+ a.real.seq!(1) ### [1, 2, ..., 10]
236
+ a.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
237
+ w.real.seq!(1) ### [1, 2, ..., 10]
238
+ w.imag.seq!(-1,-1) ### [-1, -2, ..., -10]
239
+ s = a.wsum(w.conj)
240
+ is_asserted_by { s.class == CComplex }
241
+ is_asserted_by { 770 == s.real }
242
+ is_asserted_by { 0 == s.imag }
243
+ end
244
+ end
245
+
246
+ example "accumulate" do
247
+ # ---
248
+ a = CArray.uint8(255) {1}
249
+ b = CArray.uint8(256) {1}
250
+ is_asserted_by { 255 == a.accumulate }
251
+ is_asserted_by { 0 == b.accumulate }
252
+
253
+ # ---
254
+ a = CArray.int8(127) {1}
255
+ b = CArray.int8(128) {1}
256
+ c = CArray.int8(255) {1}
257
+ d = CArray.int8(256) {1}
258
+ is_asserted_by { 127 == a.accumulate }
259
+ is_asserted_by { (-128) == b.accumulate }
260
+ is_asserted_by { (-1) == c.accumulate }
261
+ is_asserted_by { 0 == d.accumulate }
262
+ end
263
+
264
+ example "count_equal" do
265
+ # ---
266
+ a = CArray.int32(10)
267
+ a[0] = 0
268
+ a[1..3] = 1
269
+ a[4..7] = 2
270
+ a[8..9] = 3 ### [0,1,1,1,2,2,2,2,3,3]
271
+
272
+ is_asserted_by { 1 == a.count_equal(0) }
273
+ is_asserted_by { 3 == a.count_equal(1) }
274
+ is_asserted_by { 4 == a.count_equal(2) }
275
+ is_asserted_by { 2 == a.count_equal(3) }
276
+ is_asserted_by { 0 == a.count_equal(4) }
277
+ end
278
+
279
+ example "count_equiv" do
280
+ # ---
281
+ a = CArray.int32(10)
282
+ a[0] = 0
283
+ a[1..3] = 1
284
+ a[4..7] = 2
285
+ a[8..9] = 3 ### [0,1,1,1,2,2,2,2,3,3]
286
+
287
+ is_asserted_by { 1 == a.count_equiv(0, 0.001) }
288
+ is_asserted_by { 3 == a.count_equiv(1, 0.001) }
289
+ is_asserted_by { 4 == a.count_equiv(2, 0.001) }
290
+ is_asserted_by { 2 == a.count_equiv(3, 0.001) }
291
+ is_asserted_by { 0 == a.count_equiv(4, 0.001) }
292
+
293
+ # ---
294
+ a = CArray.float64(10)
295
+ a[0] = 0.00005
296
+ a[1..3] = 1.00005
297
+ a[4..7] = 2.0001
298
+ a[8..9] = 3.00015 ### [0,1,1,1,2,2,2,2,3,3]
299
+
300
+ is_asserted_by { 0 == a.count_equiv(0, 1.0e-05) }
301
+ is_asserted_by { 0 == a.count_equiv(1, 1.0e-05) }
302
+ is_asserted_by { 0 == a.count_equiv(2, 1.0e-05) }
303
+ is_asserted_by { 0 == a.count_equiv(3, 1.0e-05) }
304
+ is_asserted_by { 0 == a.count_equiv(4, 1.0e-05) }
305
+
306
+ is_asserted_by { 0 == a.count_equiv(0, 0.0001) }
307
+
308
+ is_asserted_by { 3 == a.count_equiv(1, 0.0001) }
309
+ is_asserted_by { 4 == a.count_equiv(2, 0.0001) }
310
+ is_asserted_by { 2 == a.count_equiv(3, 0.0001) }
311
+ is_asserted_by { 0 == a.count_equiv(4, 0.0001) }
312
+ end
313
+
314
+ example "count_close" do
315
+ # ---
316
+ a = CArray.int32(10)
317
+ a[0] = 0
318
+ a[1..3] = 1
319
+ a[4..7] = 2
320
+ a[8..9] = 3 ### [0,1,1,1,2,2,2,2,3,3]
321
+
322
+ is_asserted_by { 1 == a.count_close(0, 0.001) }
323
+ is_asserted_by { 3 == a.count_close(1, 0.001) }
324
+ is_asserted_by { 4 == a.count_close(2, 0.001) }
325
+ is_asserted_by { 2 == a.count_close(3, 0.001) }
326
+ is_asserted_by { 0 == a.count_close(4, 0.001) }
327
+
328
+ # ---
329
+ a = CArray.float64(10)
330
+ a[0] = 0.00005
331
+ a[1..3] = 1.00005
332
+ a[4..7] = 2.00005
333
+ a[8..9] = 3.00005 ### [0,1,1,1,2,2,2,2,3,3]
334
+
335
+ is_asserted_by { 0 == a.count_close(0, 1.0e-05) }
336
+ is_asserted_by { 0 == a.count_close(1, 1.0e-05) }
337
+ is_asserted_by { 0 == a.count_close(2, 1.0e-05) }
338
+ is_asserted_by { 0 == a.count_close(3, 1.0e-05) }
339
+ is_asserted_by { 0 == a.count_close(4, 1.0e-05) }
340
+
341
+ is_asserted_by { 1 == a.count_close(0, 0.0001) }
342
+ is_asserted_by { 3 == a.count_close(1, 0.0001) }
343
+ is_asserted_by { 4 == a.count_close(2, 0.0001) }
344
+ is_asserted_by { 2 == a.count_close(3, 0.0001) }
345
+ is_asserted_by { 0 == a.count_close(4, 0.0001) }
346
+ end
347
+
348
+ example "all_equal" do
349
+ #---
350
+ a = CArray.int32(3,3) {1}
351
+ is_asserted_by { false == a.all_equal?(0) }
352
+ is_asserted_by { true == a.all_equal?(1) }
353
+ end
354
+
355
+ example "all_equiv" do
356
+ #---
357
+ a = CArray.float64(3,3) {1.00005}
358
+ is_asserted_by { false == a.all_equiv?(0, 0.0001) }
359
+ is_asserted_by { true == a.all_equiv?(1, 0.0001) }
360
+ is_asserted_by { false == a.all_equiv?(1, 1.0e-05) }
361
+ end
362
+
363
+ example "all_close" do
364
+ #---
365
+ a = CArray.float64(3,3) {1.0001}
366
+ is_asserted_by { false == a.all_close?(0, 0.0001) }
367
+ is_asserted_by { true == a.all_close?(1, 0.0001) }
368
+ is_asserted_by { false == a.all_close?(1, 1.0e-05) }
369
+ end
370
+
371
+ example "any_equal" do
372
+ #---
373
+ a = CArray.int32(3,3)
374
+ a[1,1] = 1
375
+ is_asserted_by { true == a.any_equal?(0) }
376
+ is_asserted_by { true == a.any_equal?(1) }
377
+ end
378
+
379
+ example "any_equiv" do
380
+ #---
381
+ a = CArray.float64(3,3) { 0.00005 }
382
+ a[1,1] = 1.00005
383
+ is_asserted_by { false == a.any_equiv?(0, 0.0001) }
384
+ is_asserted_by { true == a.any_equiv?(1, 0.0001) }
385
+ is_asserted_by { false == a.any_equiv?(1, 1.0e-05) }
386
+ end
387
+
388
+ example "any_close" do
389
+ #---
390
+ a = CArray.float64(3,3) { 0.0001 }
391
+ a[1,1] = 1.0001
392
+ is_asserted_by { true == a.any_close?(0, 0.0001) }
393
+ is_asserted_by { true == a.any_close?(1, 0.0001) }
394
+ is_asserted_by { false == a.any_close?(1, 1.0e-05) }
395
+ end
396
+
397
+ end