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
@@ -0,0 +1,115 @@
1
+
2
+
3
+
4
+ class CArray
5
+
6
+ #
7
+ # Returns the 8-bit integer CArray object filled with 0 which
8
+ # dimension size is same as +self+. The resulted array represents
9
+ # the logical array which has +false+ for its all elements.
10
+ #
11
+ def false ()
12
+ return template(:boolean)
13
+ end
14
+
15
+ #
16
+ # Returns the 8-bit integer CArray object filled with 1 which
17
+ # dimension size is same as +self+. The resulted array represents
18
+ # the logical array which has +true+ for its all elements.
19
+ #
20
+ def true ()
21
+ return template(:boolean) { 1 }
22
+ end
23
+
24
+ # Returns map
25
+ def map (&block)
26
+ return self.convert(CA_OBJECT, &block).to_a
27
+ end
28
+
29
+ #
30
+ # ref = CA_INT([[0,1,2],[1,2,0],[2,0,1]])
31
+ # a = CArray.int(3,3).seq(1)
32
+ # b = CArray.int(3,3).seq(11)
33
+ # c = CArray.int(3,3).seq(21)
34
+ #
35
+ # CArray.pickup(CA_OBJECT, ref, [a,b,c])
36
+ # => <CArray.object(3,3): elem=9 mem=72b
37
+ # [ [ 1, 12, 23 ],
38
+ # [ 14, 25, 6 ],
39
+ # [ 27, 8, 19 ] ]>
40
+ #
41
+ # CArray.pickup(CA_OBJECT, ref, ["a","b","c"])
42
+ # => <CArray.object(3,3): elem=9 mem=36b
43
+ # [ [ "a", "b", "c" ],
44
+ # [ "b", "c", "a" ],
45
+ # [ "c", "a", "b" ] ]>
46
+ #
47
+ def self.pickup (data_type, ref, args)
48
+ out = ref.template(data_type)
49
+ args.each_with_index do |v, i|
50
+ s = ref.eq(i)
51
+ case v
52
+ when CArray
53
+ out[s] = v[s]
54
+ else
55
+ out[s] = v
56
+ end
57
+ end
58
+ return out
59
+ end
60
+
61
+ # Array#join like method
62
+ #
63
+ # > a = CArray.object(3,3).seq("a",:succ)
64
+ # => <CArray.object(3,3): elem=9 mem=72b
65
+ # [ [ "a", "b", "c" ],
66
+ # [ "d", "e", "f" ],
67
+ # [ "g", "h", "i" ] ]>
68
+ #
69
+ # > a.join("\n",",")
70
+ # => "a,b,c\nd,e,f\ng,h,i"
71
+ #
72
+
73
+ def join (*argv)
74
+ case argv.size
75
+ when 0
76
+ return to_a.join()
77
+ when 1
78
+ sep = argv.shift
79
+ return to_a.join(sep)
80
+ else
81
+ sep = argv.shift
82
+ return self[:i, false].map { |s|
83
+ s[0, false].join(*argv)
84
+ }.join(sep)
85
+ end
86
+ end
87
+
88
+ #
89
+ def to_bit_string (nb)
90
+ hex = CArray.uint8(((nb*elements)/8.0).ceil)
91
+ hex.bits[nil].paste([0], self.bits[false,[(nb-1)..0]].flatten)
92
+ hex.bits[] = hex.bits[nil,[-1..0]]
93
+ return hex.to_s
94
+ end
95
+
96
+ def from_bit_string (bstr, nb)
97
+ hex = CArray.uint8(bstr.length).load_binary(bstr)
98
+ hex.bits[] = hex.bits[nil,[-1..0]]
99
+ bits = hex.bits.flatten
100
+ self.bits[false,[(nb-1)..0]][nil].paste([0], bits)
101
+ return self
102
+ end
103
+
104
+ def self.from_bit_string (bstr, nb, data_type=CA_INT32, dim=nil)
105
+ if dim
106
+ obj = CArray.new(data_type, dim)
107
+ else
108
+ dim0 = ((bstr.length*8)/nb.to_f).floor
109
+ obj = CArray.new(data_type, [dim0])
110
+ end
111
+ obj.from_bit_string(bstr, nb)
112
+ return obj
113
+ end
114
+
115
+ end
@@ -3,10 +3,8 @@
3
3
  # carray/info.rb
4
4
  #
5
5
  # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
6
  #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
10
8
  #
11
9
  # ----------------------------------------------------------------------------
12
10
 
@@ -1,12 +1,10 @@
1
1
  # ----------------------------------------------------------------------------
2
2
  #
3
- # carray/base/inspect.rb
3
+ # carray/inspect.rb
4
4
  #
5
5
  # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
6
  #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
10
8
  #
11
9
  # ----------------------------------------------------------------------------
12
10
 
@@ -17,14 +15,14 @@ class CArray::Inspector # :nodoc:
17
15
  end
18
16
 
19
17
  def inspect_string
20
- if @carray.rank == 0
21
- raise "can't inspect CArray of rank == 0"
18
+ if @carray.ndim == 0
19
+ raise "can't inspect CArray of ndim == 0"
22
20
  end
23
21
  formatter = get_formatter()
24
22
  class_name = get_class_name()
25
23
  type_name = get_type_name()
26
24
  shape = get_shape()
27
- data_spec = get_data_spec(0, Array.new(@carray.rank){0}, formatter)
25
+ data_spec = get_data_spec(0, Array.new(@carray.ndim){0}, formatter)
28
26
  info_list = get_info_list()
29
27
  output = ["<",
30
28
  format("%s.%s(%s)", class_name, type_name, shape.join(",")),
@@ -152,14 +150,14 @@ class CArray::Inspector # :nodoc:
152
150
 
153
151
  def get_data_spec (level, idx, formatter)
154
152
  io = "[ "
155
- rank = @carray.rank
153
+ ndim = @carray.ndim
156
154
  dim = @carray.dim
157
- if level == rank - 1
155
+ if level == ndim - 1
158
156
  over = false
159
157
  dim[level].times do |i|
160
158
  idx[level] = i
161
159
  v = @carray[*idx]
162
- if v.undef?
160
+ if v == UNDEF
163
161
  io << '_'
164
162
  else
165
163
  io << formatter[v]
@@ -178,7 +176,7 @@ class CArray::Inspector # :nodoc:
178
176
  if over
179
177
  idx[level] = dim[level] - 1
180
178
  v = @carray[*idx]
181
- if v.undef?
179
+ if v == UNDEF
182
180
  io << ", _"
183
181
  else
184
182
  io << ", " + formatter[v]
@@ -3,10 +3,8 @@
3
3
  # carray/io/imagemagick.rb
4
4
  #
5
5
  # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
6
  #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
10
8
  #
11
9
  # ----------------------------------------------------------------------------
12
10
 
@@ -121,7 +119,7 @@ class CArray
121
119
 
122
120
  def magick_guess_image_type
123
121
  image_type = nil
124
- case rank
122
+ case ndim
125
123
  when 2
126
124
  if fixlen?
127
125
  case bytes
@@ -1,12 +1,10 @@
1
1
  # ----------------------------------------------------------------------------
2
2
  #
3
- # carray/base/iterator.rb
3
+ # carray/iterator.rb
4
4
  #
5
5
  # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
6
  #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
10
8
  #
11
9
  # ----------------------------------------------------------------------------
12
10
 
@@ -280,6 +278,7 @@ class CAIterator
280
278
  end
281
279
 
282
280
  def sort_with (&block)
281
+ warn "CAIterator#sort_with will be obsolete"
283
282
  out = reference.template
284
283
  idx = CA.sort_addr(*yield(self))
285
284
  ca[idx].each_with_addr do |blk, i|
@@ -317,8 +316,9 @@ class CAClassIterator < CAIterator # :nodoc:
317
316
  @reference = reference
318
317
  @classifier = classifier || @reference.uniq.sort
319
318
  @null = CArray.new(@reference.data_type,[0])
319
+ @null.data_class = @reference.data_class if @reference.has_data_class?
320
320
  @table = {}
321
- @rank = 1
321
+ @ndim = 1
322
322
  @dim = [0]
323
323
  if @classifier.all_masked? or @classifier.size == 0
324
324
  @dim = [0]
@@ -343,7 +343,7 @@ class CAClassIterator < CAIterator # :nodoc:
343
343
  return self
344
344
  end
345
345
 
346
- def ranking (&block)
346
+ def ndiming (&block)
347
347
  block ||= lambda {|a| a.size }
348
348
  values = self.to_a.map{|v| block[v] }.to_ca
349
349
  addrs = values.sort_addr.reverse
@@ -0,0 +1,102 @@
1
+ # ----------------------------------------------------------------------------
2
+ #
3
+ # carray/mask.rb
4
+ #
5
+ # This file is part of Ruby/CArray extension library.
6
+ #
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
8
+ #
9
+ # ----------------------------------------------------------------------------
10
+
11
+ class CArray
12
+
13
+ # Guard methods for handling variables that can be UNDEF values.
14
+ # Provide different processing depending on whether
15
+ # the given value is UNDEF or not.
16
+ #
17
+ # @param value [Object] target object
18
+ # @param fill_value [Object] alternative value if the given value is UNDEF
19
+ #
20
+ # @return fill_value if the given value is UNDEF, block return value if block
21
+ # is given, or value itself
22
+ #
23
+ def self.guard_undef (*values, fill_value: UNDEF, &block)
24
+ return fill_value if values.any?{|v| v == UNDEF }
25
+ return block.(*values)
26
+ end
27
+
28
+ # mask
29
+ #
30
+ # Returns the number of masked elements.
31
+ #
32
+ def count_masked (*axis)
33
+ if has_mask?
34
+ return mask.int64.accumulate(*axis)
35
+ else
36
+ if axis.empty?
37
+ return 0
38
+ else
39
+ spec = shape.map{:i}
40
+ axis.each do |k|
41
+ spec[k] = nil
42
+ end
43
+ return self[*spec].ca.template(:int64) { 0 }
44
+ end
45
+ end
46
+ end
47
+
48
+ #
49
+ # Returns the number of not-masked elements.
50
+ #
51
+ def count_not_masked (*axis)
52
+ if has_mask?
53
+ return mask.not.int64.accumulate(*axis)
54
+ else
55
+ if axis.empty?
56
+ return elements
57
+ else
58
+ spec = shape.map {:i}
59
+ axis.each do |k|
60
+ spec[k] = nil
61
+ end
62
+ it = self[*spec].ca
63
+ count = self.elements/it.elements
64
+ return it.template(:int64) { count }
65
+ end
66
+ end
67
+ end
68
+
69
+ def maskout! (*argv)
70
+ case argv.size
71
+ when 1
72
+ val = argv.first
73
+ case val
74
+ when CArray, Symbol
75
+ self[val] = UNDEF
76
+ else
77
+ self[:eq, val] = UNDEF
78
+ end
79
+ else
80
+ self[*argv] = UNDEF
81
+ end
82
+ return self
83
+ end
84
+
85
+ def maskout (*argv)
86
+ obj = self.to_ca
87
+ case argv.size
88
+ when 1
89
+ val = argv.first
90
+ case val
91
+ when CArray, Symbol
92
+ obj[val] = UNDEF
93
+ else
94
+ obj[:eq, val] = UNDEF
95
+ end
96
+ else
97
+ obj[*argv] = UNDEF
98
+ end
99
+ return obj
100
+ end
101
+
102
+ end
@@ -1,12 +1,10 @@
1
1
  # ----------------------------------------------------------------------------
2
2
  #
3
- # carray/base/math.rb
3
+ # carray/math.rb
4
4
  #
5
5
  # This file is part of Ruby/CArray extension library.
6
- # You can redistribute it and/or modify it under the terms of
7
- # the Ruby Licence.
8
6
  #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
7
+ # Copyright (C) 2005-2020 Hiroki Motoyoshi
10
8
  #
11
9
  # ----------------------------------------------------------------------------
12
10
 
@@ -148,7 +146,7 @@ class CArray
148
146
  def quo (other)
149
147
  case
150
148
  when integer?
151
- return double/other
149
+ return object.quo_i(other)
152
150
  when object?
153
151
  return quo_i(other)
154
152
  else
@@ -226,7 +224,7 @@ class CArray
226
224
 
227
225
  def anomaly (*argv)
228
226
  opt = argv.last.is_a?(Hash) ? argv.pop : {}
229
- idxs = Array.new(self.rank) { |i| argv.include?(i) ? :* : nil }
227
+ idxs = Array.new(self.ndim) { |i| argv.include?(i) ? :* : nil }
230
228
  if mn = opt[:mean]
231
229
  return self - mn[*idxs]
232
230
  else
@@ -260,7 +258,7 @@ class CArray
260
258
  return (c[(n-1)/2] + c[n/2])/2.0
261
259
  end
262
260
  else
263
- raise "CArray#median is not implemented for multiple ranks"
261
+ raise "CArray#median is not implemented for multiple ndims"
264
262
  end
265
263
 
266
264
  end
@@ -285,27 +283,23 @@ class CArray
285
283
  end
286
284
  out = []
287
285
  begin
288
- pers.each do |per|
289
- if per == 100
290
- out << ca[n-1]
291
- elsif per >= 0 and per < 100
292
- if n > 1
293
- f = (n-1)*per/100.0
294
- k = f.floor
295
- r = f - k
296
- out << (1-r)*ca[k] + r*ca[k+1]
286
+ pers.each do |per|
287
+ if per == 100
288
+ out << ca[n-1]
289
+ elsif per >= 0 and per < 100
290
+ if n > 1
291
+ f = (n-1)*per/100.0
292
+ k = f.floor
293
+ r = f - k
294
+ out << (1-r)*ca[k] + r*ca[k+1]
295
+ else
296
+ out << ca[0]
297
+ end
297
298
  else
298
- out << ca[0]
299
+ out << CA_NAN
299
300
  end
300
- else
301
- out << CA_NAN
302
301
  end
303
302
  end
304
- rescue
305
- p self[:is_not_masked]
306
- p n
307
- raise
308
- end
309
303
  return out
310
304
  end
311
305
 
@@ -322,7 +316,7 @@ class CArray
322
316
  count = x.count_not_masked
323
317
  xm = x.mean(:min_count => min_count)
324
318
  ym = y.mean(:min_count => min_count)
325
- if ( xm.undef? or ym.undef? )
319
+ if ( xm == UNDEF or ym == UNDEF )
326
320
  return fill_value || UNDEF
327
321
  else
328
322
  return (x-xm).wsum(y-ym)/count
@@ -341,7 +335,7 @@ class CArray
341
335
  count = x.count_not_masked
342
336
  xm = x.mean(:min_count=>min_count)
343
337
  ym = y.mean(:min_count=>min_count)
344
- if ( xm.undef? or ym.undef? )
338
+ if ( xm == UNDEF or ym == UNDEF )
345
339
  return fill_value || UNDEF
346
340
  else
347
341
  return (x-xm).wsum(y-ym)/(count-1)
@@ -359,7 +353,7 @@ class CArray
359
353
  y.inherit_mask(x)
360
354
  xm = x.mean(:min_count=>min_count)
361
355
  ym = y.mean(:min_count=>min_count)
362
- if ( xm.undef? or ym.undef? )
356
+ if ( xm == UNDEF or ym == UNDEF )
363
357
  return fill_value || UNDEF
364
358
  else
365
359
  xd, yd = x-xm, y-ym
@@ -373,31 +367,5 @@ class CArray
373
367
 
374
368
  end
375
369
 
376
- class CArray
377
-
378
- def self.summation (*dim)
379
- out = nil
380
- first = true
381
- CArray.each_index(*dim) { |*idx|
382
- if first
383
- out = yield(*idx)
384
- first = false
385
- else
386
- out += yield(*idx)
387
- end
388
- }
389
- return out
390
- end
391
-
392
- def by (other)
393
- case other
394
- when CArray
395
- return (self[nil][nil,:*]*other[nil][:*,nil]).reshape(*(dim+other.dim))
396
- else
397
- return self * other
398
- end
399
- end
400
-
401
- end
402
370
 
403
371