carray 1.3.6 → 1.5.3

Sign up to get free protection for your applications and to get access to all the features.
Files changed (315) hide show
  1. checksums.yaml +4 -4
  2. data/.yardopts +32 -0
  3. data/LICENSE +22 -0
  4. data/NEWS.md +65 -0
  5. data/README.md +23 -26
  6. data/Rakefile +31 -0
  7. data/TODO.md +17 -0
  8. data/carray.gemspec +24 -31
  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} +38 -40
  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} +90 -59
  31. data/{carray_access.c → ext/carray_access.c} +244 -109
  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} +350 -149
  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} +58 -70
  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} +137 -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} +95 -114
  46. data/{carray_math.rb → ext/carray_math.rb} +47 -31
  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} +72 -36
  50. data/{carray_order.c → ext/carray_order.c} +73 -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} +59 -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} +1 -6
  58. data/{mkmath.rb → ext/mkmath.rb} +13 -3
  59. data/{ruby_carray.c → ext/ruby_carray.c} +20 -8
  60. data/{ruby_ccomplex.c → ext/ruby_ccomplex.c} +1 -3
  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 +52 -10
  65. data/lib/carray/{base/autoload.rb → autoload.rb} +48 -8
  66. data/lib/carray/autoload/autoload_base.rb +1 -5
  67. data/lib/carray/autoload/autoload_gem_cairo.rb +9 -0
  68. data/lib/carray/autoload/autoload_gem_ffi.rb +9 -0
  69. data/lib/carray/autoload/autoload_gem_gnuplot.rb +2 -0
  70. data/lib/carray/autoload/autoload_gem_io_csv.rb +14 -0
  71. data/lib/carray/autoload/autoload_gem_io_pg.rb +6 -0
  72. data/lib/carray/autoload/autoload_gem_io_sqlite3.rb +12 -0
  73. data/lib/carray/autoload/autoload_gem_narray.rb +10 -0
  74. data/lib/carray/autoload/autoload_gem_numo_narray.rb +15 -0
  75. data/lib/carray/autoload/autoload_gem_opencv.rb +16 -0
  76. data/lib/carray/autoload/autoload_gem_random.rb +8 -0
  77. data/lib/carray/autoload/autoload_gem_rmagick.rb +23 -0
  78. data/lib/carray/autoload/{autoload_graphics_zimg.rb → autoload_gem_zimg.rb} +0 -0
  79. data/lib/carray/basic.rb +191 -0
  80. data/lib/carray/broadcast.rb +82 -0
  81. data/lib/carray/compose.rb +315 -0
  82. data/lib/carray/construct.rb +460 -0
  83. data/lib/carray/convert.rb +115 -0
  84. data/lib/carray/info.rb +1 -3
  85. data/lib/carray/{base/inspect.rb → inspect.rb} +9 -11
  86. data/lib/carray/io/imagemagick.rb +2 -4
  87. data/lib/carray/{base/iterator.rb → iterator.rb} +6 -6
  88. data/lib/carray/mask.rb +102 -0
  89. data/lib/carray/{base/math.rb → math.rb} +21 -53
  90. data/lib/carray/math/histogram.rb +8 -10
  91. data/lib/carray/math/recurrence.rb +1 -3
  92. data/lib/carray/mkmf.rb +9 -3
  93. data/lib/carray/object/ca_obj_iterator.rb +1 -3
  94. data/lib/carray/object/ca_obj_link.rb +1 -3
  95. data/lib/carray/object/ca_obj_pack.rb +9 -11
  96. data/lib/carray/obsolete.rb +270 -0
  97. data/lib/carray/ordering.rb +181 -0
  98. data/lib/carray/{base/serialize.rb → serialize.rb} +50 -66
  99. data/lib/carray/string.rb +188 -0
  100. data/lib/carray/{base/struct.rb → struct.rb} +19 -21
  101. data/lib/carray/{io/table.rb → table.rb} +1 -9
  102. data/lib/carray/testing.rb +51 -0
  103. data/lib/carray/time.rb +76 -0
  104. data/lib/carray/transform.rb +98 -0
  105. data/misc/Methods.ja.md +182 -0
  106. data/{NOTE → misc/NOTE} +0 -0
  107. data/spec/Classes/CABitfield_spec.rb +58 -0
  108. data/spec/Classes/CABlockIterator_spec.rb +114 -0
  109. data/spec/Classes/CABlock_spec.rb +205 -0
  110. data/spec/Classes/CAField_spec.rb +39 -0
  111. data/spec/Classes/CAGrid_spec.rb +75 -0
  112. data/spec/Classes/CAMap_spec.rb +0 -0
  113. data/{test/test_CAMapping.rb → spec/Classes/CAMapping_spec.rb} +35 -36
  114. data/spec/Classes/CAObject_attribute_spec.rb +33 -0
  115. data/spec/Classes/CAObject_spec.rb +33 -0
  116. data/spec/Classes/CARefer_spec.rb +93 -0
  117. data/spec/Classes/CARepeat_spec.rb +65 -0
  118. data/spec/Classes/CASelect_spec.rb +22 -0
  119. data/spec/Classes/CAShift_spec.rb +16 -0
  120. data/spec/Classes/CAStruct_spec.rb +71 -0
  121. data/{test/test_CATranspose.rb → spec/Classes/CATranspose_spec.rb} +20 -21
  122. data/spec/Classes/CAUnboudRepeat_spec.rb +102 -0
  123. data/spec/Classes/CAWindow_spec.rb +54 -0
  124. data/spec/Classes/CAWrap_spec.rb +8 -0
  125. data/{test/test_CArray.rb → spec/Classes/CArray_spec.rb} +48 -92
  126. data/spec/Classes/CScalar_spec.rb +55 -0
  127. data/spec/Features/feature_130_spec.rb +19 -0
  128. data/spec/Features/feature_attributes_spec.rb +280 -0
  129. data/spec/Features/feature_boolean_spec.rb +98 -0
  130. data/spec/Features/feature_broadcast.rb +116 -0
  131. data/spec/Features/feature_cast_function.rb +19 -0
  132. data/spec/Features/feature_cast_spec.rb +33 -0
  133. data/spec/Features/feature_class_spec.rb +84 -0
  134. data/spec/Features/feature_complex_spec.rb +42 -0
  135. data/{test/test_composite.rb → spec/Features/feature_composite_spec.rb} +17 -18
  136. data/spec/Features/feature_convert_spec.rb +46 -0
  137. data/spec/Features/feature_copy_spec.rb +123 -0
  138. data/spec/Features/feature_creation_spec.rb +84 -0
  139. data/spec/Features/feature_element_spec.rb +144 -0
  140. data/spec/Features/feature_extream_spec.rb +54 -0
  141. data/spec/Features/feature_generate_spec.rb +74 -0
  142. data/spec/Features/feature_index_spec.rb +69 -0
  143. data/spec/Features/feature_mask_spec.rb +580 -0
  144. data/spec/Features/feature_math_spec.rb +97 -0
  145. data/spec/Features/feature_order_spec.rb +146 -0
  146. data/spec/Features/feature_ref_store_spec.rb +209 -0
  147. data/spec/Features/feature_serialization_spec.rb +125 -0
  148. data/spec/Features/feature_stat_spec.rb +397 -0
  149. data/spec/Features/feature_virtual_spec.rb +48 -0
  150. data/spec/Features/method_eq_spec.rb +81 -0
  151. data/spec/Features/method_is_nan_spec.rb +12 -0
  152. data/spec/Features/method_map_spec.rb +54 -0
  153. data/spec/Features/method_max_with.rb +20 -0
  154. data/spec/Features/method_min_with.rb +19 -0
  155. data/spec/Features/method_ne_spec.rb +18 -0
  156. data/spec/Features/method_project_spec.rb +188 -0
  157. data/spec/Features/method_ref_spec.rb +27 -0
  158. data/spec/Features/method_round_spec.rb +11 -0
  159. data/spec/Features/method_s_linspace_spec.rb +48 -0
  160. data/spec/Features/method_s_span_spec.rb +14 -0
  161. data/spec/Features/method_seq_spec.rb +47 -0
  162. data/spec/Features/method_sort_with.rb +43 -0
  163. data/spec/Features/method_sorted_with.rb +29 -0
  164. data/spec/Features/method_span_spec.rb +42 -0
  165. data/spec/Features/method_wrap_readonly_spec.rb +43 -0
  166. data/{test → spec/UnitTest}/test_CAVirtual.rb +0 -0
  167. data/spec/spec_all.rb +0 -1
  168. data/utils/convert_test.rb +73 -0
  169. data/utils/{extract_rdoc.rb → extract_yard.rb} +7 -12
  170. data/{devel → utils}/guess_shape.rb +0 -0
  171. data/utils/{diff_method.rb → monkey_patch_methods.rb} +17 -7
  172. metadata +173 -316
  173. data/COPYING +0 -56
  174. data/GPL +0 -340
  175. data/Gemfile +0 -8
  176. data/Gemfile.lock +0 -33
  177. data/LEGAL +0 -50
  178. data/Makefile +0 -451
  179. data/TODO +0 -5
  180. data/carray_cast_func.c +0 -8670
  181. data/carray_config.h +0 -26
  182. data/carray_math.c +0 -24805
  183. data/carray_random.c +0 -531
  184. data/carray_stat_proc.c +0 -12170
  185. data/ext/calculus/carray_calculus.c +0 -931
  186. data/ext/calculus/carray_interp.c +0 -358
  187. data/ext/calculus/extconf.rb +0 -12
  188. data/ext/calculus/lib/autoload/autoload_math_calculus.rb +0 -2
  189. data/ext/calculus/lib/math/calculus.rb +0 -119
  190. data/ext/calculus/lib/math/interp/adapter_interp1d.rb +0 -31
  191. data/ext/fortio/extconf.rb +0 -3
  192. data/ext/fortio/lib/carray/autoload/autoload_fortran_format.rb +0 -5
  193. data/ext/fortio/lib/carray/io/fortran_format.rb +0 -43
  194. data/ext/fortio/lib/fortio.rb +0 -3
  195. data/ext/fortio/lib/fortio/fortran_format.rb +0 -605
  196. data/ext/fortio/lib/fortio/fortran_format.tab.rb +0 -536
  197. data/ext/fortio/lib/fortio/fortran_format.y +0 -215
  198. data/ext/fortio/lib/fortio/fortran_namelist.rb +0 -151
  199. data/ext/fortio/lib/fortio/fortran_namelist.tab.rb +0 -470
  200. data/ext/fortio/lib/fortio/fortran_namelist.y +0 -213
  201. data/ext/fortio/lib/fortio/fortran_sequential.rb +0 -345
  202. data/ext/fortio/ruby_fortio.c +0 -182
  203. data/ext/fortio/test/test_H.rb +0 -5
  204. data/ext/fortio/test/test_T.rb +0 -7
  205. data/ext/fortio/test/test_fortran_format.rb +0 -86
  206. data/ext/fortio/test/test_namelist.rb +0 -25
  207. data/ext/fortio/test/test_namelist_write.rb +0 -10
  208. data/ext/fortio/test/test_sequential.rb +0 -13
  209. data/ext/fortio/test/test_sequential2.rb +0 -13
  210. data/ext/fortio/work/test.rb +0 -10
  211. data/ext/fortio/work/test_e.rb +0 -19
  212. data/ext/fortio/work/test_ep.rb +0 -10
  213. data/ext/fortio/work/test_parse.rb +0 -12
  214. data/ext/imagemap/carray_imagemap.c +0 -495
  215. data/ext/imagemap/doc/call_graph.dot +0 -28
  216. data/ext/imagemap/draw.c +0 -567
  217. data/ext/imagemap/extconf.rb +0 -13
  218. data/ext/imagemap/lib/autoload/autoload_graphics_imagemap.rb +0 -1
  219. data/ext/imagemap/lib/graphics/imagemap.rb +0 -273
  220. data/ext/imagemap/lib/image_map.rb +0 -4
  221. data/ext/imagemap/test/swath_index.rb +0 -83
  222. data/ext/imagemap/test/swath_warp.rb +0 -99
  223. data/ext/imagemap/test/test.rb +0 -23
  224. data/ext/imagemap/test/test_image.rb +0 -42
  225. data/ext/imagemap/test/test_line.rb +0 -14
  226. data/ext/imagemap/test/test_rotate.rb +0 -17
  227. data/ext/imagemap/test/test_triangle.rb +0 -20
  228. data/ext/imagemap/test/test_warp.rb +0 -26
  229. data/ext/mathfunc/extconf.rb +0 -18
  230. data/ext/mathfunc/lib/autoload/autoload_math_mathfunc.rb +0 -1
  231. data/ext/mathfunc/lib/math/mathfunc.rb +0 -15
  232. data/ext/mathfunc/test/test_hypot.rb +0 -5
  233. data/ext/mathfunc/test/test_j0.rb +0 -22
  234. data/ext/mathfunc/test/test_jn.rb +0 -8
  235. data/ext/mathfunc/test/test_sph.rb +0 -9
  236. data/ext/narray/README +0 -22
  237. data/ext/narray/ca_wrap_narray.c +0 -500
  238. data/ext/narray/carray_narray.c +0 -21
  239. data/ext/narray/extconf.rb +0 -57
  240. data/ext/narray/lib/autoload/autoload_math_narray.rb +0 -1
  241. data/ext/narray/lib/autoload/autoload_math_narray_miss.rb +0 -11
  242. data/ext/narray/lib/math/narray.rb +0 -17
  243. data/ext/narray/lib/math/narray_miss.rb +0 -45
  244. data/lib/carray/autoload/autoload_graphics_gnuplot.rb +0 -2
  245. data/lib/carray/autoload/autoload_io_csv.rb +0 -14
  246. data/lib/carray/autoload/autoload_io_pg.rb +0 -6
  247. data/lib/carray/autoload/autoload_io_sqlite3.rb +0 -12
  248. data/lib/carray/autoload/autoload_io_table.rb +0 -1
  249. data/lib/carray/autoload/autoload_math_interp.rb +0 -4
  250. data/lib/carray/base/basic.rb +0 -1090
  251. data/lib/carray/base/obsolete.rb +0 -131
  252. data/lib/carray/graphics/gnuplot.rb +0 -2131
  253. data/lib/carray/graphics/zimg.rb +0 -296
  254. data/lib/carray/io/csv.rb +0 -572
  255. data/lib/carray/io/pg.rb +0 -101
  256. data/lib/carray/io/sqlite3.rb +0 -215
  257. data/lib/carray/math/interp.rb +0 -57
  258. data/lib/carray/math/interp/adapter_gsl_spline.rb +0 -47
  259. data/mkmf.log +0 -375
  260. data/mt19937ar.c +0 -182
  261. data/mt19937ar.h +0 -86
  262. data/rdoc_ext.rb +0 -1115
  263. data/rdoc_main.rb +0 -27
  264. data/rdoc_math.rb +0 -5
  265. data/rdoc_stat.rb +0 -31
  266. data/spec/CABlockIterator/CABlockIterator_spec.rb +0 -113
  267. data/spec/CArray/bug/store_spec.rb +0 -27
  268. data/spec/CArray/index/repeat_spec.rb +0 -10
  269. data/spec/CArray/method/eq_spec.rb +0 -80
  270. data/spec/CArray/method/is_nan_spec.rb +0 -12
  271. data/spec/CArray/method/ne_spec.rb +0 -18
  272. data/spec/CArray/method/round_spec.rb +0 -11
  273. data/spec/CArray/object/_attribute_spec.rb +0 -32
  274. data/spec/CArray/object/s_new_spec.rb +0 -31
  275. data/spec/CArray/serialize/Serialization_spec.rb +0 -89
  276. data/test/test_130.rb +0 -23
  277. data/test/test_ALL.rb +0 -51
  278. data/test/test_CABitfield.rb +0 -59
  279. data/test/test_CABlock.rb +0 -208
  280. data/test/test_CAField.rb +0 -40
  281. data/test/test_CAGrid.rb +0 -76
  282. data/test/test_CAMmap.rb +0 -11
  283. data/test/test_CARefer.rb +0 -94
  284. data/test/test_CARepeat.rb +0 -66
  285. data/test/test_CASelect.rb +0 -23
  286. data/test/test_CAShift.rb +0 -17
  287. data/test/test_CAWindow.rb +0 -55
  288. data/test/test_CAWrap.rb +0 -9
  289. data/test/test_CComplex.rb +0 -83
  290. data/test/test_CScalar.rb +0 -91
  291. data/test/test_attribute.rb +0 -281
  292. data/test/test_block_iterator.rb +0 -17
  293. data/test/test_boolean.rb +0 -99
  294. data/test/test_cast.rb +0 -33
  295. data/test/test_class.rb +0 -85
  296. data/test/test_complex.rb +0 -43
  297. data/test/test_convert.rb +0 -79
  298. data/test/test_copy.rb +0 -141
  299. data/test/test_creation.rb +0 -85
  300. data/test/test_element.rb +0 -146
  301. data/test/test_extream.rb +0 -55
  302. data/test/test_generate.rb +0 -75
  303. data/test/test_index.rb +0 -71
  304. data/test/test_mask.rb +0 -578
  305. data/test/test_math.rb +0 -98
  306. data/test/test_narray.rb +0 -64
  307. data/test/test_order.rb +0 -147
  308. data/test/test_random.rb +0 -15
  309. data/test/test_ref_store.rb +0 -211
  310. data/test/test_stat.rb +0 -414
  311. data/test/test_struct.rb +0 -72
  312. data/test/test_virtual.rb +0 -49
  313. data/utils/create_rdoc.sh +0 -9
  314. data/utils/make_tgz.sh +0 -3
  315. data/version.h +0 -18
@@ -1,55 +0,0 @@
1
- $:.unshift("../lib", "..")
2
-
3
- require 'test/unit'
4
- require 'carray'
5
-
6
- class TestCArrayExtream < Test::Unit::TestCase
7
-
8
- def test_zero_length_array
9
- # ---
10
- a = CArray.int(0)
11
- assert_equal(1, a.rank)
12
- assert_equal(0, a.elements)
13
- assert_equal([0], a.dim)
14
- assert_equal(0, a.dim0)
15
-
16
- # ---
17
- a = CArray.int(0, 3)
18
- assert_equal(2, a.rank)
19
- assert_equal(0, a.elements)
20
- assert_equal([0,3], a.dim)
21
- assert_equal(0, a.dim0)
22
- assert_equal(3, a.dim1)
23
-
24
- # ---
25
- a = CArray.int(3, 0)
26
- assert_equal(2, a.rank)
27
- assert_equal(0, a.elements)
28
- assert_equal([3,0], a.dim)
29
- assert_equal(3, a.dim0)
30
- assert_equal(0, a.dim1)
31
- end
32
-
33
- def test_zero_length_fixlen
34
- # ---
35
- assert_nothing_raised { CArray.new(CA_FIXLEN, [3,3], :bytes => 0) }
36
- assert_raise(RuntimeError){ CArray.new(CA_FIXLEN, [3,3], :bytes => -1) }
37
- end
38
-
39
- def test_large_rank_array
40
- # ---
41
- dim = [2] * (CA_RANK_MAX)
42
- assert_instance_of(CArray, CArray.int8(*dim))
43
-
44
- # ---
45
- dim = [2] * (CA_RANK_MAX+1)
46
- assert_raise(RuntimeError){ CArray.int8(*dim) }
47
- end
48
-
49
- def test_negative_dimension_size
50
- # ---
51
- assert_raise(RuntimeError){ CArray.int8(-1) }
52
- assert_raise(RuntimeError){ CArray.int8(-1, -1) }
53
- end
54
-
55
- end
@@ -1,75 +0,0 @@
1
- #$:.unshift("../lib", "..")
2
-
3
- require 'test/unit'
4
- require 'carray'
5
-
6
- class TestCArrayGenerate < Test::Unit::TestCase
7
-
8
- def test_swap_bytes
9
- #
10
- a = CArray.int16(3,3) { 0x1234 }
11
- b = CArray.int16(3,3) { 0x3412 }
12
- assert_equal(b, a.swap_bytes)
13
- assert_equal(a, b.swap_bytes)
14
- assert_equal(a, a.swap_bytes.swap_bytes)
15
-
16
- #
17
- a = CArray.int32(3,3) { 0x12345678 }
18
- b = CArray.int32(3,3) { 0x78563412 }
19
- assert_equal(b, a.swap_bytes)
20
- assert_equal(a, b.swap_bytes)
21
- assert_equal(a, a.swap_bytes.swap_bytes)
22
-
23
- #
24
- if CArray::HAVE_COMPLEX
25
- c = CArray.int32(1) { 0x12345678 }
26
- x = c.refer(CA_FLOAT32, [1])[0]
27
- y = c.swap_bytes.refer(CA_FLOAT32, [1])[0]
28
- a = CArray.complex(3,3) { x + y*CI }
29
- b = CArray.complex(3,3) { y + x*CI }
30
- assert_equal(b, a.swap_bytes)
31
- assert_equal(a, b.swap_bytes)
32
- assert_equal(a, a.swap_bytes.swap_bytes)
33
- end
34
- end
35
-
36
- def test_seq
37
- # ---
38
- a = CArray.object(3).seq!
39
- assert_equal(CA_OBJECT([0,1,2]), a)
40
-
41
- # ---
42
- a = CArray.object(3).seq!(1)
43
- assert_equal(CA_OBJECT([1,2,3]), a)
44
-
45
- # ---
46
- a = CArray.object(3).seq!(1,2)
47
- assert_equal(CA_OBJECT([1,3,5]), a)
48
-
49
- # ---
50
- a = CArray.object(3).seq!(3,-1)
51
- assert_equal(CA_OBJECT([3,2,1]), a)
52
-
53
- # ---
54
- a = CArray.object(3).seq!("a", "a")
55
- assert_equal(CA_OBJECT(["a", "aa", "aaa"]), a)
56
-
57
- # ---
58
- a = CArray.object(3).seq!("a", :succ)
59
- assert_equal(CA_OBJECT(["a", "b", "c"]), a)
60
-
61
- # ---
62
- a = CArray.object(3).seq!("a", :succ) { |x| "@" + x }
63
- assert_equal(CA_OBJECT(["@a", "@b", "@c"]), a)
64
-
65
- # ---
66
- a = CArray.object(3).seq!("a", :succ)
67
- assert_equal(CA_OBJECT(["@a", "@b", "@c"]), CA_OBJECT("@") + a)
68
-
69
- # ---
70
- a = CArray.object(3).seq!([], [nil])
71
- assert_equal(CA_OBJECT([[], [nil], [nil, nil]]), a)
72
- end
73
-
74
-
75
- end
@@ -1,71 +0,0 @@
1
- $:.unshift("../lib", "..")
2
-
3
- require 'test/unit'
4
- require 'carray'
5
-
6
- class TestCArrayIndex < Test::Unit::TestCase
7
-
8
- def test_index
9
- a = CArray.int(3,3).seq!
10
- assert_raise(IndexError) { a[0, 0, 0]} ### invalid rank
11
- assert_raise(IndexError) { a[[0,1,1,2], 0]}
12
- assert_raise(IndexError) { a[[0,4], 0]}
13
- assert_raise(IndexError) { a[[0,3,2], 0]}
14
- assert_raise(NoMethodError) { a[:abcdefg, 0]}
15
- assert_raise(IndexError) { a[0, :abcdefg]}
16
- end
17
-
18
- def test_rubber_dim
19
- a = CArray.int(3,3,3,3,3).seq!
20
-
21
- assert_equal(a[false], a[nil, nil, nil, nil, nil])
22
-
23
- assert_equal(a[0,false], a[0, nil, nil, nil, nil])
24
- assert_equal(a[false,0], a[nil, nil, nil, nil, 0])
25
- assert_equal(a[0,false,0], a[0, nil, nil, nil, 0])
26
-
27
- assert_equal(a[0..1,false], a[0..1, nil, nil, nil, nil])
28
- assert_equal(a[false,0..1], a[nil, nil, nil, nil, 0..1])
29
- assert_equal(a[0..1,false,0..1], a[0..1, nil, nil, nil, 0..1])
30
-
31
- assert_equal(a[nil,false], a[nil, nil, nil, nil, nil])
32
- assert_equal(a[false,nil], a[nil, nil, nil, nil, nil])
33
- assert_equal(a[nil,false,nil], a[nil, nil, nil, nil, nil])
34
-
35
- assert_raise(IndexError) { a[0,0,0,0,0,0,false] }
36
-
37
- end
38
-
39
-
40
- def test_scan_index
41
- info = CArray.scan_index([3,3],[])
42
- assert_equal([info.type, info.index], [CA_REG_ALL, []])
43
- info = CArray.scan_index([3,3],[1])
44
- assert_equal([info.type, info.index], [CA_REG_ADDRESS, [1]])
45
- info = CArray.scan_index([3,3],[1..2])
46
- assert_equal([info.type, info.index], [CA_REG_ADDRESS_COMPLEX, [[1,2,1]]])
47
- info = CArray.scan_index([3,3],[1,1])
48
- assert_equal([info.type, info.index], [CA_REG_POINT, [1,1]])
49
- info = CArray.scan_index([3,3],[1,nil])
50
- assert_equal([info.type, info.index], [CA_REG_BLOCK, [1,[0,3,1]]])
51
- a = CArray.int(3,3).seq!
52
- i = CArray.int(3).seq!
53
- info = CArray.scan_index([3,3],[a.eq(4)])
54
- assert_equal([info.type, info.index], [CA_REG_SELECT, []])
55
- info = CArray.scan_index([3,3],[a])
56
- assert_equal([info.type, info.index], [CA_REG_MAPPING, []])
57
- info = CArray.scan_index([3,3],[i,i])
58
- assert_equal([info.type, info.index], [CA_REG_GRID, []])
59
- info = CArray.scan_index([3,3],[:%, :%, 2])
60
- assert_equal([info.type, info.index], [CA_REG_REPEAT, []])
61
- info = CArray.scan_index([3,3],[nil, nil, :*])
62
- assert_equal([info.type, info.index], [CA_REG_UNBOUND_REPEAT, []])
63
- end
64
-
65
- def test_addr2index
66
- a = CArray.int(3,3,3)
67
- assert_equal([1,1,1], a.addr2index(13))
68
- assert_equal(13, a.index2addr(1,1,1))
69
- end
70
-
71
- end
@@ -1,578 +0,0 @@
1
- $:.unshift("../lib", "..")
2
-
3
- require 'test/unit'
4
- require 'carray'
5
-
6
- class TestCArrayMask < Test::Unit::TestCase
7
-
8
- def test_basic_features
9
-
10
- # ---
11
- a = CArray.int32(3,3).seq!
12
- assert_equal(false, a.has_mask?) ### mask array is not created
13
- assert_equal(CA_INT32([[0,1,2],
14
- [3,4,5],
15
- [6,7,8]]), a.value)
16
- assert_equal(0, a.mask) ### mask array is not created
17
- assert_equal(CA_BOOLEAN([[0,0,0],
18
- [0,0,0],
19
- [0,0,0]]), a.is_masked)
20
- assert_equal(CA_BOOLEAN([[1,1,1],
21
- [1,1,1],
22
- [1,1,1]]), a.is_not_masked)
23
- assert_equal(0, a.count_masked)
24
- assert_equal(9, a.count_not_masked)
25
-
26
- # ---
27
- a = CArray.int32(3,3).seq!
28
- a.mask = 0
29
- assert_equal(true, a.has_mask?) ### mask array is created
30
- assert_equal(CA_INT32([[0,1,2],
31
- [3,4,5],
32
- [6,7,8]]), a.value)
33
- assert_equal(CA_BOOLEAN([[0,0,0],
34
- [0,0,0],
35
- [0,0,0]]), a.mask) ### mask array is created
36
- assert_equal(CA_BOOLEAN([[0,0,0],
37
- [0,0,0],
38
- [0,0,0]]), a.is_masked)
39
- assert_equal(CA_BOOLEAN([[1,1,1],
40
- [1,1,1],
41
- [1,1,1]]), a.is_not_masked)
42
- assert_equal(0, a.count_masked)
43
- assert_equal(9, a.count_not_masked)
44
-
45
- # ---
46
- a = CArray.int32(3,3).seq!
47
- a.mask = CA_BOOLEAN([[0,1,0],
48
- [1,1,1],
49
- [0,1,0]])
50
- assert_equal(true, a.has_mask?)
51
- assert_equal(CA_INT32([[0,1,2],
52
- [3,4,5],
53
- [6,7,8]]), a.value)
54
- assert_equal(CA_BOOLEAN([[0,1,0],
55
- [1,1,1],
56
- [0,1,0]]), a.mask)
57
- assert_equal(CA_BOOLEAN([[0,1,0],
58
- [1,1,1],
59
- [0,1,0]]), a.is_masked)
60
- assert_equal(CA_BOOLEAN([[1,0,1],
61
- [0,0,0],
62
- [1,0,1]]), a.is_not_masked)
63
- assert_equal(5, a.count_masked)
64
- assert_equal(4, a.count_not_masked)
65
- end
66
-
67
- def test_equality_when_mask_is_added
68
- # ---
69
- a = CArray.int32(3,3).seq!
70
- a.mask = CA_BOOLEAN([[0,1,0],
71
- [1,1,1],
72
- [0,1,0]])
73
- b = CArray.int32(3,3).seq!
74
- b.mask = CA_BOOLEAN([[0,1,0],
75
- [1,1,1],
76
- [0,1,0]])
77
- assert_equal(a.mask, b.mask)
78
- assert_equal(a.is_masked, b.is_masked)
79
- assert_equal(a, b)
80
-
81
- # ---
82
- a = CArray.int32(3,3).seq!
83
- a.mask = CA_BOOLEAN([[0,1,0],
84
- [0,1,0],
85
- [0,1,0]])
86
- b = CArray.int32(3,3).seq!
87
- b.mask = CA_BOOLEAN([[0,0,0],
88
- [1,1,1],
89
- [0,0,0]])
90
- assert_not_equal(a.mask, b.mask)
91
- assert_not_equal(a.is_masked, b.is_masked)
92
- assert_not_equal(a, b)
93
-
94
- # ---
95
- a = CArray.int32(3,3).seq!
96
- b = CArray.int32(3,3).seq!
97
- b.mask = CA_BOOLEAN([[0,0,0],
98
- [0,0,0],
99
- [0,0,0]])
100
- assert_not_equal(a.mask, b.mask)
101
- assert_equal(a.is_masked, b.is_masked)
102
- assert_equal(a, b)
103
-
104
- # ---
105
- a = CArray.int32(3,3).seq!
106
- a.mask = CA_BOOLEAN([[0,0,0],
107
- [1,1,1],
108
- [0,0,0]])
109
- b = CArray.int32(3,3).seq!
110
- b.mask = CA_BOOLEAN([[0,0,0],
111
- [1,1,1],
112
- [0,0,0]])
113
- a.value[1,nil] = 1
114
- b.value[1,nil] = -1
115
- assert_not_equal(a.value, b.value)
116
- assert_equal(a.mask, b.mask)
117
- assert_equal(a.is_masked, b.is_masked)
118
- assert_equal(a, b)
119
- end
120
-
121
- def test_has_mask?
122
- # test_basic_feature
123
- end
124
-
125
- def test_any_masked?
126
- #
127
- a = CA_INT(0..2)
128
- a.mask = 0
129
- assert_equal(false, a.any_masked?)
130
- a[1] = UNDEF
131
- assert_equal(true, a.any_masked?)
132
- a.unmask
133
- assert_equal(false, a.any_masked?)
134
- end
135
-
136
- def test_value
137
- #
138
- a = CArray.int32(3,3)
139
- a.mask = 0
140
- assert_raise(RuntimeError) {
141
- a.value.mask = 1 ### can't set mask to value array
142
- }
143
- end
144
-
145
- def test_mask
146
- # test_basic_feature
147
- a = CArray.int32(3,3)
148
- a.mask = 0
149
- assert_raise(RuntimeError) {
150
- a.mask.mask = 1 ### can't set mask to mask array
151
- }
152
- end
153
-
154
- def test_mask=
155
- # ---
156
- a = CArray.int32(3,3).seq!
157
- a.mask = CA_BOOLEAN([[0,1,0],
158
- [1,1,1],
159
- [0,1,0]]) ### Though the mask array can hold the value
160
- ### other than 0 or 1, the value 0 means
161
- ### masked and otherwise means not-masked.
162
- assert_equal(CA_BOOLEAN([[0,1,0],
163
- [1,1,1],
164
- [0,1,0]]), a.mask)
165
- assert_equal(CA_BOOLEAN([[0,1,0],
166
- [1,1,1],
167
- [0,1,0]]), a.is_masked)
168
- assert_equal(CA_BOOLEAN([[1,0,1],
169
- [0,0,0],
170
- [1,0,1]]), a.is_not_masked)
171
- assert_equal(5, a.count_masked)
172
- assert_equal(4, a.count_not_masked)
173
-
174
- # ---
175
- a = CArray.int32(3,3).seq!
176
- m = CA_BOOLEAN([[0,1,0],
177
- [0,1,0],
178
- [0,1,0]])
179
- m.mask = CA_BOOLEAN([[0,0,0],
180
- [1,1,1],
181
- [0,0,0]])
182
- a.mask = m ### When masked array is set as mask,
183
- ### masked position is also set as
184
- ### masked elements.
185
- assert_not_equal(m, a.mask)
186
- assert_equal(CA_BOOLEAN([[0,1,0],
187
- [1,1,1],
188
- [0,1,0]]), a.mask)
189
- end
190
-
191
- def test_is_masked
192
- # test_basic_features
193
- end
194
-
195
- def test_is_not_masked
196
- # test_basic_features
197
- end
198
-
199
- def test_count_masked
200
- # test_basic_featrues
201
- end
202
-
203
- def test_count_not_masked
204
- # test_basic_featrues
205
- end
206
-
207
- def test_unmask
208
- # ---
209
- x = CArray.int32(3,3).seq!
210
- x.mask = CA_BOOLEAN([[0,1,0],
211
- [1,1,1],
212
- [0,1,0]])
213
-
214
- a = x.to_ca
215
- b = x.to_ca
216
-
217
- a.unmask()
218
- assert_equal(true, a.has_mask?) ### never mask eliminated from array
219
- assert_equal(CA_INT32([[0,1,2],
220
- [3,4,5],
221
- [6,7,8]]), a)
222
-
223
- b.unmask(-1)
224
- assert_equal(true, b.has_mask?) ### never mask eliminated from array
225
- assert_equal(CA_INT32([[0,-1,2], ### filled with fill_value
226
- [-1,-1,-1],
227
- [6,-1,8]]), b)
228
- end
229
-
230
- def test_unmask_copy
231
- # ---
232
- a = CArray.int32(3,3).seq!
233
- a.mask = CA_BOOLEAN([[0,1,0],
234
- [1,1,1],
235
- [0,1,0]])
236
- u = a.unmask_copy()
237
- assert_equal(false, u.has_mask?) ### never mask eliminated from array
238
- assert_equal(CA_INT32([[0,1,2],
239
- [3,4,5],
240
- [6,7,8]]), u)
241
-
242
- # ---
243
- a = CArray.int32(3,3).seq!
244
- a.mask = CA_BOOLEAN([[0,1,0],
245
- [1,1,1],
246
- [0,1,0]])
247
- u = a.unmask_copy(-1)
248
- assert_equal(false, u.has_mask?) ### never mask eliminated from array
249
- assert_equal(CA_INT32([[0,-1,2], ### filled with fill_value
250
- [-1,-1,-1],
251
- [6,-1,8]]), u)
252
- end
253
-
254
- def test_inherit_mask
255
- # ---
256
- a = CArray.int32(3,3).seq!
257
- b = CArray.int32(3,3).seq!
258
- c = CArray.int32(3,3).seq!
259
- a.mask = CA_BOOLEAN([[0,1,0],
260
- [0,0,0],
261
- [0,1,0]])
262
- b.mask = CA_BOOLEAN([[0,0,0],
263
- [1,0,1],
264
- [0,0,0]])
265
- c.mask = CA_BOOLEAN([[0,0,0],
266
- [0,1,0],
267
- [0,0,0]])
268
- c.inherit_mask(a,b)
269
- assert_equal(CA_BOOLEAN([[0,1,0],
270
- [1,1,1],
271
- [0,1,0]]), c.mask)
272
- end
273
-
274
- def test_inherit_mask_replace
275
- # ---
276
- a = CArray.int32(3,3).seq!
277
- b = CArray.int32(3,3).seq!
278
- c = CArray.int32(3,3).seq!
279
- a.mask = CA_BOOLEAN([[0,1,0],
280
- [0,0,0],
281
- [0,1,0]])
282
- b.mask = CA_BOOLEAN([[0,0,0],
283
- [1,0,1],
284
- [0,0,0]])
285
- c.mask = CA_BOOLEAN([[0,0,0],
286
- [0,1,0],
287
- [0,0,0]])
288
- c.inherit_mask_replace(a,b)
289
- assert_equal(CA_BOOLEAN([[0,1,0],
290
- [1,0,1],
291
- [0,1,0]]), c.mask)
292
-
293
- # ---
294
- a = CArray.int32(3,3).seq!
295
-
296
- end
297
-
298
- # ----------------------------------------------------------------------
299
-
300
- def test_monop
301
- # ---
302
- a = CArray.float64(3,3) {1}
303
- a.mask = CA_BOOLEAN([[0,1,0],
304
- [1,1,1],
305
- [0,1,0]])
306
- b = -a ### example of monop
307
- assert_equal(a.mask, b.mask)
308
- end
309
-
310
- def test_binop
311
- # ---
312
- a = CArray.float64(3,3) {1}
313
- b = CArray.float64(3,3) {1}
314
- a.mask = CA_BOOLEAN([[0,1,0],
315
- [0,1,0],
316
- [0,1,0]])
317
- b.mask = CA_BOOLEAN([[0,0,0],
318
- [1,1,1],
319
- [0,0,0]])
320
- c = a + b ### example of binop
321
- assert_equal(CA_BOOLEAN([[0,1,0],
322
- [1,1,1],
323
- [0,1,0]]), c.mask)
324
- end
325
-
326
- def test_cmpop
327
- # ---
328
- a = CArray.float64(3,3) {1}
329
- b = CArray.float64(3,3) {1}
330
- a.mask = CA_BOOLEAN([[0,1,0],
331
- [0,1,0],
332
- [0,1,0]])
333
- b.mask = CA_BOOLEAN([[0,0,0],
334
- [1,1,1],
335
- [0,0,0]])
336
- c = a.le(b) ### example of binop
337
- assert_equal(CA_BOOLEAN([[0,1,0],
338
- [1,1,1],
339
- [0,1,0]]), c.mask)
340
- end
341
-
342
- # ----------------------------------------------------------------------
343
-
344
- def test_basic_stat
345
- # ---
346
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
347
-
348
- assert_equal(1, a.min)
349
- assert_equal(1, a.min(:mask_limit => 0))
350
- assert_equal(1, a.min(:mask_limit => 0, :fill_value => -9999))
351
-
352
- assert_equal(10, a.max)
353
- assert_equal(10, a.max(:mask_limit => 0))
354
- assert_equal(10, a.max(:mask_limit => 0, :fill_value => -9999))
355
-
356
- assert_equal(55, a.sum)
357
- assert_equal(55, a.sum(:mask_limit => 0))
358
- assert_equal(55, a.sum(:mask_limit => 0, :fill_value => -9999))
359
-
360
- assert_equal(5.5, a.mean)
361
- assert_equal(5.5, a.mean(:mask_limit => 0))
362
- assert_equal(5.5, a.mean(:mask_limit => 0, :fill_value => -9999))
363
-
364
- # ---
365
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
366
- a.mask = 0
367
- a[0] = UNDEF
368
- a[9] = UNDEF ### [_, 2, ..., 9,_]
369
-
370
- assert_equal(2, a.min)
371
- assert_equal(2, a.min(:mask_limit => 0))
372
- assert_equal(UNDEF, a.min(:mask_limit => 1))
373
- assert_equal(UNDEF, a.min(:mask_limit => 2))
374
- assert_equal(2, a.min(:mask_limit => 3))
375
- assert_equal(2, a.min(:mask_limit => 0, :fill_value => -9999))
376
- assert_equal(-9999, a.min(:mask_limit => 1, :fill_value => -9999))
377
- assert_equal(-9999, a.min(:mask_limit => 2, :fill_value => -9999))
378
- assert_equal(2, a.min(:mask_limit => 3, :fill_value => -9999))
379
-
380
- assert_equal(9, a.max)
381
- assert_equal(9, a.max(:mask_limit => 0))
382
- assert_equal(UNDEF, a.max(:mask_limit => 1))
383
- assert_equal(UNDEF, a.max(:mask_limit => 2))
384
- assert_equal(9, a.max(:mask_limit => 3))
385
- assert_equal(9, a.max(:mask_limit => 0, :fill_value => -9999))
386
- assert_equal(-9999, a.max(:mask_limit => 1, :fill_value => -9999))
387
- assert_equal(-9999, a.max(:mask_limit => 2, :fill_value => -9999))
388
- assert_equal(9, a.max(:mask_limit => 3, :fill_value => -9999))
389
-
390
- assert_equal(44, a.sum)
391
- assert_equal(44, a.sum(:mask_limit => 0))
392
- assert_equal(UNDEF, a.sum(:mask_limit => 1))
393
- assert_equal(UNDEF, a.sum(:mask_limit => 2))
394
- assert_equal(44, a.sum(:mask_limit => 3))
395
- assert_equal(44, a.sum(:mask_limit => 0, :fill_value => -9999))
396
- assert_equal(-9999, a.sum(:mask_limit => 1, :fill_value => -9999))
397
- assert_equal(-9999, a.sum(:mask_limit => 2, :fill_value => -9999))
398
- assert_equal(44, a.sum(:mask_limit => 3, :fill_value => -9999))
399
-
400
- assert_equal(5.5, a.mean)
401
- assert_equal(5.5, a.mean(:mask_limit => 0))
402
- assert_equal(UNDEF, a.mean(:mask_limit => 1))
403
- assert_equal(UNDEF, a.mean(:mask_limit => 2))
404
- assert_equal(5.5, a.mean(:mask_limit => 3))
405
- assert_equal(5.5, a.mean(:mask_limit => 0, :fill_value => -9999))
406
- assert_equal(-9999, a.mean(:mask_limit => 1, :fill_value => -9999))
407
- assert_equal(-9999, a.mean(:mask_limit => 2, :fill_value => -9999))
408
- assert_equal(5.5, a.mean(:mask_limit => 3, :fill_value => -9999))
409
-
410
- end
411
-
412
- def test_accumulate
413
- # ---
414
- a = CArray.uint8(256) { 1 }
415
-
416
- assert_equal(0, a.accumulate)
417
- assert_equal(0, a.accumulate(:mask_limit => 0))
418
- assert_equal(0, a.accumulate(:mask_limit => 0, :fill_value => -9999))
419
-
420
- # ---
421
- a = CArray.uint8(256) { 1 }
422
- a[0] = UNDEF
423
- a[255] = UNDEF
424
-
425
- assert_equal(254, a.accumulate())
426
- assert_equal(254, a.accumulate(:mask_limit => 0))
427
- assert_equal(UNDEF, a.accumulate(:mask_limit => 1))
428
- assert_equal(UNDEF, a.accumulate(:mask_limit => 2))
429
- assert_equal(254, a.accumulate(:mask_limit => 0, :fill_value => -9999))
430
- assert_equal(-9999, a.accumulate(:mask_limit => 1, :fill_value => -9999))
431
- assert_equal(-9999, a.accumulate(:mask_limit => 2, :fill_value => -9999))
432
- end
433
-
434
- def test_variance
435
- # ---
436
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
437
-
438
- assert_equal(true, (9.16666666666667 - a.variance()) < 0.000001 )
439
- assert_equal(true, (9.16666666666667 - a.variance(:mask_limit=>0)) < 0.000001 )
440
- assert_equal(true, (9.16666666666667 - a.variance(:mask_limit=>0, :fill_value=>-9999)) < 0.000001 )
441
-
442
- # ---
443
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
444
- a.mask = 0
445
- a[0] = UNDEF
446
- a[9] = UNDEF ### [_, 2, ..., 9,_]
447
-
448
- assert_equal(6, a.variance())
449
- assert_equal(6, a.variance(:mask_limit=>0))
450
- assert_equal(UNDEF, a.variance(:mask_limit=>1))
451
- assert_equal(UNDEF, a.variance(:mask_limit=>2))
452
- assert_equal(6, a.variance(:mask_limit=>3))
453
- assert_equal(6, a.variance(:mask_limit=>0, :fill_value=>-9999))
454
- assert_equal(-9999, a.variance(:mask_limit=>1, :fill_value=>-9999))
455
- assert_equal(-9999, a.variance(:mask_limit=>2, :fill_value=>-9999))
456
- assert_equal(6, a.variance(:mask_limit=>3, :fill_value=>-9999))
457
-
458
- end
459
-
460
- def test_variancep
461
- # ---
462
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
463
-
464
- assert_equal(true, (8.25 - a.variancep()) < 0.000001 )
465
- assert_equal(true, (8.25 - a.variancep(:mask_limit=>0)) < 0.000001 )
466
- assert_equal(true, (8.25 - a.variancep(:mask_limit=>0, :fill_value=>-9999)) < 0.000001 )
467
-
468
- # ---
469
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
470
- a.mask = 0
471
- a[0] = UNDEF
472
- a[9] = UNDEF ### [_, 2, ..., 9,_]
473
-
474
- assert_equal(true, (5.25 - a.variancep()) < 0.000001 )
475
- assert_equal(true, (5.25 - a.variancep(:mask_limit=>0)) < 0.000001 )
476
- assert_equal(UNDEF, a.variancep(:mask_limit=>1))
477
- assert_equal(UNDEF, a.variancep(:mask_limit=>2))
478
- assert_equal(true, (5.25 - a.variancep(:mask_limit=>3)) < 0.000001 )
479
- assert_equal(true, (5.25 - a.variancep(:mask_limit=>0, :fill_value=>-9999)) < 0.000001 )
480
- assert_equal(-9999, a.variancep(:mask_limit=>1, :fill_value=>-9999))
481
- assert_equal(-9999, a.variancep(:mask_limit=>2, :fill_value=>-9999))
482
- assert_equal(true, (5.25 - a.variancep(:mask_limit=>3, :fill_value=>-9999)) < 0.000001 )
483
- end
484
-
485
- def test_prod
486
- # ---
487
- a = CArray.int32(5).seq!(1) ### [1, 2, ..., 5]
488
-
489
- assert_equal(120, a.prod())
490
- assert_equal(120, a.prod(:mask_limit => 0))
491
- assert_equal(120, a.prod(:mask_limit => 0, :fill_value => -9999))
492
-
493
- # ---
494
- a = CArray.int32(5).seq!(1) ### [1, 2, ..., 5]
495
- a[0] = UNDEF
496
- a[2] = UNDEF ### [0, 2, _, 4, 5]
497
-
498
- assert_equal(40, a.prod())
499
- assert_equal(40, a.prod(:mask_limit => 0))
500
- assert_equal(UNDEF, a.prod(:mask_limit => 1))
501
- assert_equal(UNDEF, a.prod(:mask_limit => 2))
502
- assert_equal(40, a.prod(:mask_limit => 3))
503
- assert_equal(40, a.prod(:mask_limit => 0, :fill_value => -9999))
504
- assert_equal(-9999, a.prod(:mask_limit => 1, :fill_value => -9999))
505
- assert_equal(-9999, a.prod(:mask_limit => 2, :fill_value => -9999))
506
- assert_equal(40, a.prod(:mask_limit => 3, :fill_value => -9999))
507
- end
508
-
509
- def test_wsum
510
-
511
- # ---
512
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
513
- w = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
514
-
515
- assert_equal(385, a.wsum(w))
516
- assert_equal(385, a.wsum(w, 0))
517
- assert_equal(385, a.wsum(w, 0, -9999))
518
-
519
- # ---
520
- a = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
521
- w = CArray.int32(10).seq!(1) ### [1, 2, ..., 10]
522
- a.mask = 0
523
- w.mask = 0
524
- a[0] = UNDEF ### [_, 2, ..., 9,10]
525
- w[9] = UNDEF ### [1, 2, ..., 9,_]
526
-
527
- # assert_equal(284, a.wsum(w))
528
- # assert_equal(UNDEF, a.wsum(w, 0))
529
- # assert_equal(UNDEF, a.wsum(w, 1))
530
- # assert_equal(284, a.wsum(w, 2))
531
-
532
- # assert_equal(-9999, a.wsum(w, 0, -9999))
533
- # assert_equal(-9999, a.wsum(w, 1, -9999))
534
- # assert_equal(284, a.wsum(w, 2, -9999))
535
-
536
- end
537
-
538
- def test_cumsum
539
- # ---
540
- a = CArray.int32(10).seq!(1)
541
- a[2] = UNDEF
542
- assert_equal(CA_DOUBLE([1, 3, 3, 7, 12, 18, 25, 33, 42, 52]),
543
- a.cumsum())
544
- assert_equal(CA_DOUBLE([1, 3, *([UNDEF]*8)]),
545
- a.cumsum(0))
546
- assert_equal(CA_DOUBLE([1, 3, *([-9999]*8)]),
547
- a.cumsum(0, -9999))
548
- assert_equal(CA_DOUBLE([1, 3, 3, 7, 12, 18, 25, 33, 42, 52]),
549
- a.cumsum(1))
550
- end
551
-
552
- def test_cumprod
553
- # ---
554
- a = CArray.int32(5).seq!(1)
555
- a[2] = UNDEF
556
- assert_equal(CA_DOUBLE([1, 2, 2, 8, 40]),
557
- a.cumprod())
558
- assert_equal(CA_DOUBLE([1, 2, *([UNDEF]*3)]),
559
- a.cumprod(0))
560
- assert_equal(CA_DOUBLE([1, 2, *([-9999]*3)]),
561
- a.cumprod(0, -9999))
562
- assert_equal(CA_DOUBLE([1, 2, 2, 8, 40]),
563
- a.cumprod(1))
564
- end
565
-
566
- def test_count_xxx
567
- # TODO
568
- end
569
-
570
- def test_all_xxx?
571
- # TODO
572
- end
573
-
574
- def test_any_xxx?
575
- # TODO
576
- end
577
-
578
- end