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,23 +0,0 @@
1
- require 'carray'
2
-
3
- #sc = ImageMap::ColorScale.pm3d(13,10,14,256,:interpolation => false)
4
- #sc = ImageMap::ColorScale.GMT("GMT_globe.cpt")
5
- #sc.set_range(0..Math::PI)
6
- #sc.set_upper_color(127,127,127)
7
- #sc.set_lower_color(127,127,127)
8
- #sc.set_mask_color( 127,127,127)
9
-
10
- ca = CArray.float(640, 480)
11
- ca[:i,nil].scale!(-0.1, Math::PI+0.1)
12
- ca.mask = 0
13
- ca.mask[90..110,90..110] = 1
14
-
15
- if true
16
- im = ImageMap.create(ca, sc)
17
- require 'carray/RMagick'
18
- il = Magick::Image.new_from_ca(im)
19
- il.display
20
- else
21
- im = ImageMap.display(ca, sc)
22
- end
23
-
@@ -1,42 +0,0 @@
1
- require 'carray'
2
- require 'carray/image_map'
3
- require 'benchmark'
4
-
5
- lon = CArray.float(37).span!(-90..180)
6
- lat = CArray.float(12).span!(0..90)
7
-
8
- u = lon[:%,lat.dim0]
9
- v = lat[lon.dim0,:%]
10
- val = u.rad.sin! + v.rad
11
-
12
- image = ImageMap.new(CA_DOUBLE, [600, 600]) { -9999 }
13
-
14
- image.set_xrange(-103..103)
15
- image.set_yrange(-103..103)
16
-
17
- Benchmark.bm do |o|
18
- o.report {
19
- image.warp(val, u, v, :grid=>"point", :gradation => true) {
20
- |_u,_v,_x,_y|
21
- _x[] = (100 - _v/3) * _u.rad.cos!
22
- _y[] = (100 - _v/3) * _u.rad.sin!
23
- }
24
- }
25
- end
26
-
27
- #image.draw_line(0, 0, 0, 90, 0)
28
- image.line(0, 0, 0, 90, 0)
29
-
30
- poly = CArray.float(4,2) {[[-90,90],[0,0],[45,45],[0,90]]}
31
-
32
- image.polyline(poly[nil,1],poly[nil,0], 0)
33
-
34
- image[:lt, -9998] = UNDEF
35
-
36
- #sc = nil
37
- #sc = ImageMap::ColorScale.pm3d(11,7,22,12)
38
- sc = ColorScale.GMT("GMT_polar.cpt")
39
- #sc.set_mask_color(64,64,64)
40
- sc.set_range(nil, nil)
41
-
42
- ImageMap.display(image, sc)
@@ -1,14 +0,0 @@
1
- require 'carray'
2
- require 'carray/image_map'
3
-
4
- p img = ImageMap.new(CA_INT8, [200,200]) {0}
5
- p img.class
6
-
7
- gx = CArray.float(2,2) { [[10,30],
8
- [10,30]] }
9
- gy = CArray.float(2,2) { [[10,10],
10
- [80,50]] }
11
-
12
- img.fill_rectangle(gy, gx, 255)
13
-
14
- img.display_by_magick
@@ -1,17 +0,0 @@
1
- require "carray"
2
-
3
- img = ImageMap.new(:int16, [300, 300]) {0}
4
- img.set_xrange 0...300
5
- img.set_yrange 0...300
6
-
7
- data = CArray.int16(16,16).seq!
8
-
9
- x = CArray.float(16).seq![16,:%]
10
- y = CArray.float(16).seq!.reverse![:%,16]
11
-
12
- img.warp(data, x, y, :grid=>"area0", :gradation=>false) {|gx,gy,ix,iy|
13
- ix[] = 8*(gy+gx) + 24
14
- iy[] = 8*(gy-gx) + 150
15
- }
16
-
17
- img.as_uint8.display_by_magick
@@ -1,20 +0,0 @@
1
- require 'carray'
2
-
3
- image = CArray.float(200, 200) { nan }
4
-
5
- y = CArray.float(2,2){[50,50, 65,100]}
6
- x = CArray.float(2,2){[70,100,50,90]}
7
- z = CArray.float(2,2){[30,40, 50,20]}
8
- ImageMap::Image.draw_rectangle_gradation(image, y, x, z)
9
-
10
- #image.draw_line(y[0], x[0], y[1], x[1])
11
- #image.draw_line(y[1], x[1], y[2], x[2])
12
- #image.draw_line(y[2], x[2], y[0], x[0])
13
- #image.draw_line(y[2], x[2], y[0], x[0])
14
-
15
- image.mask = image.is_nan
16
-
17
- sc = ImageMap::ColorScale.pm3d(13,6,10)
18
-
19
- ImageMap.display(image, sc)
20
-
@@ -1,26 +0,0 @@
1
- require 'carray'
2
-
3
- lon = CArray.float(121).span!(-180..180)
4
- lat = CArray.float(12).span!(0..90)
5
- p val = CArray.uint16(121,12).seq!.quantize(256)
6
-
7
- img = ImageMap.uint16(600,600) { 9999 }
8
-
9
- img.warp(val, lon[:%,lat.dim0], lat[lon.dim0,:%],
10
- :grid=>"area0", :gradation=>false) { |u,v,x,y|
11
- x[] = (100 - v/2) * u.rad.cos!
12
- y[] = (100 - v/2) * u.rad.sin!
13
- x.mul!(2).add!(300)
14
- y.mul!(2).add!(300)
15
- }
16
-
17
- img[:eq,9999] = UNDEF
18
-
19
- RGB = CA.struct{ uint8 :r,:g,:b }
20
- rgb = CArray.new(RGB, [255,3])
21
- rgb["r"].seq!
22
-
23
- out = rgb.project(img)
24
-
25
- out.display_by_magick("rgb")
26
-
@@ -1,18 +0,0 @@
1
- require "mkmf"
2
-
3
- dir_config("carray", "../..", "../..")
4
-
5
- have_header("carray.h")
6
-
7
- if /cygwin|mingw/ =~ RUBY_PLATFORM
8
- have_library("carray")
9
- end
10
-
11
- have_header("tgmath.h")
12
-
13
- have_func("atan2", "math.h")
14
- have_func("hypot", "math.h")
15
- have_func("lgamma", "math.h")
16
- have_func("expm1", "math.h")
17
-
18
- create_makefile("carray/carray_mathfunc")
@@ -1 +0,0 @@
1
- require "carray/math/mathfunc"
@@ -1,15 +0,0 @@
1
- # ----------------------------------------------------------------------------
2
- #
3
- # lib/carray/graphics/imagemap.rb
4
- #
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
- #
9
- # Copyright (C) 2005 Hiroki Motoyoshi
10
- #
11
- # ----------------------------------------------------------------------------
12
-
13
- require "carray/carray_mathfunc"
14
-
15
-
@@ -1,5 +0,0 @@
1
- require 'carray'
2
-
3
- p x = CArray.object(101).span(0..1)
4
-
5
- p CAMath.hypot(x[:*,nil], x[nil,:*])
@@ -1,22 +0,0 @@
1
- require 'carray'
2
-
3
- x = CArray.float(101).span(0..10)
4
-
5
- CA.gnuplot { |g|
6
- g.plot2d(
7
- [x, x.erf],
8
- [x, x.erfc]
9
- )
10
- g.plot2d(
11
- [x, x.j0],
12
- [x, x.j1],
13
- [x, x.jn(2)],
14
- [x, x.jn(3)],
15
- [x, x.jn(4)],
16
- [x, x.jn(5)]
17
- )
18
- g.plot2d(
19
- [x, x.y0],
20
- [x, x.y1]
21
- )
22
- }
@@ -1,8 +0,0 @@
1
- require "carray"
2
-
3
- p x = CAMath.yn(1, CA_FLOAT(0..100))
4
-
5
- CA.gnuplot { |g|
6
- g.plot2d([x], :with=>"lines")
7
- }
8
-
@@ -1,9 +0,0 @@
1
- require "carray"
2
-
3
- r = CA_DOUBLE(1)
4
- theta = CA_DOUBLE(30).rad
5
- phi = CA_DOUBLE(45).rad
6
-
7
- x,y,z = CAMath.spherical_to_xyz(r, theta, phi)
8
- p CAMath.xyz_to_spherical(x,y,z)
9
- p [r, theta, phi]
@@ -1,22 +0,0 @@
1
-
2
- This directory contains the plugin library used as the interface to NArray
3
- and NArrayMiss library from Ruby/CArray.
4
-
5
- ---
6
- NArray Library (http://narray.rubyforge.org/) by Masahiro Tanaka.
7
-
8
- NArray is an Numerical N-dimensional Array class. Supported element types
9
- are 1/2/4-byte Integer, single/double-precision Real/Complex, and Ruby
10
- Object. This extension library incorporates fast calculation and easy
11
- manipulation of large numerical arrays into the Ruby language. NArray
12
- has features similar to NumPy, but NArray has vector and matrix subclasses.
13
-
14
- (License: Ruby License)
15
-
16
- ---
17
- NArrayMiss Library (http://ruby.gfd-dennou.org/index.htm) by Seiya Nishizawa.
18
-
19
- NArrayMiss is a additional class processing of missing value with to NArray
20
- for Ruby.
21
-
22
- (License: GNU Lesser General Public License (LGPL) version 2.)
@@ -1,500 +0,0 @@
1
- /* ---------------------------------------------------------------------------
2
-
3
- carray/ca_wrap_narray.c
4
-
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
-
9
- Copyright (C) 2005 Hiroki Motoyoshi
10
-
11
- ---------------------------------------------------------------------------- */
12
-
13
- #include "ruby.h"
14
- #include "carray.h"
15
-
16
- #include "narray_config.h"
17
- #ifndef HAVE_U_INT8_T
18
- #define HAVE_U_INT8_T
19
- #endif
20
- #ifndef HAVE_INT16_T
21
- #define HAVE_INT16_T
22
- #endif
23
- #ifndef HAVE_INT16_T
24
- #define HAVE_INT32_T
25
- #endif
26
- #ifndef HAVE_U_INT32_T
27
- #define HAVE_U_INT32_T
28
- #endif
29
- #include "narray.h"
30
-
31
- #include <math.h>
32
-
33
- extern VALUE cNMatrix, cNVector, cNMatrixLU;
34
-
35
- static char EMPTY_ARRAY_PTR;
36
-
37
- /* -------------------------------------------------------------------- */
38
-
39
- static int8_t
40
- na_typecode_to_ca_data_type (int typecode)
41
- {
42
- int8_t data_type;
43
-
44
- switch ( typecode ) {
45
- case NA_BYTE:
46
- data_type = CA_UINT8; break;
47
- case NA_SINT:
48
- data_type = CA_INT16; break;
49
- case NA_LINT:
50
- data_type = CA_INT32; break;
51
- case NA_SFLOAT:
52
- data_type = CA_FLOAT32; break;
53
- case NA_DFLOAT:
54
- data_type = CA_FLOAT64; break;
55
- #ifdef HAVE_COMPLEX_H
56
- case NA_SCOMPLEX:
57
- data_type = CA_CMPLX64; break;
58
- case NA_DCOMPLEX:
59
- data_type = CA_CMPLX128; break;
60
- #endif
61
- case NA_ROBJ:
62
- data_type = CA_OBJECT; break;
63
- default:
64
- rb_raise(rb_eCADataTypeError,
65
- "unknown NArray typecode <%i>", typecode);
66
- }
67
- return data_type;
68
- }
69
-
70
- static int
71
- ca_data_type_to_na_typecode(int8_t data_type)
72
- {
73
- int typecode;
74
-
75
- switch ( data_type ) {
76
- case CA_BOOLEAN:
77
- case CA_UINT8:
78
- typecode = NA_BYTE; break;
79
- case CA_INT16:
80
- typecode = NA_SINT; break;
81
- case CA_INT32:
82
- typecode = NA_LINT; break;
83
- case CA_FLOAT32:
84
- typecode = NA_SFLOAT; break;
85
- case CA_FLOAT64:
86
- typecode = NA_DFLOAT; break;
87
- #ifdef HAVE_COMPLEX_H
88
- case CA_CMPLX64:
89
- typecode = NA_SCOMPLEX; break;
90
- case CA_CMPLX128:
91
- typecode = NA_DCOMPLEX; break;
92
- #endif
93
- case CA_OBJECT:
94
- typecode = NA_ROBJ; break;
95
- default:
96
- rb_raise(rb_eRuntimeError,
97
- "no corresponding NArray typecode for CArray data type <%s>",
98
- ca_type_name[data_type]);
99
- }
100
-
101
- return typecode;
102
- }
103
-
104
- /* -------------------------------------------------------------------- */
105
-
106
- static void
107
- cary_na_ref_free (struct NARRAY *ary)
108
- {
109
- xfree(ary->shape);
110
- xfree(ary);
111
- }
112
-
113
- static void
114
- cary_na_ref_mark (struct NARRAY *ary)
115
- {
116
- rb_gc_mark( ary->ref );
117
- }
118
-
119
- /* -------------------------------------------------------------------- */
120
-
121
- static VALUE
122
- rb_cary_na_ref_new_i (int argc, VALUE *argv, VALUE self, VALUE klass)
123
- {
124
- CArray *orig;
125
- struct NARRAY *ary;
126
- int8_t i;
127
-
128
- Data_Get_Struct(self, CArray, orig);
129
-
130
- if ( ! ca_is_attached(orig) ) {
131
- rb_raise(rb_eRuntimeError,
132
- "cannot create NArray reference for not attached CArray");
133
- }
134
-
135
- ary = ALLOC(struct NARRAY);
136
- if ( !ary ) {
137
- rb_raise(rb_eRuntimeError, "cannot allocate NArray");
138
- }
139
-
140
- if ( orig->elements == 0 ) {
141
- ary->rank = 0;
142
- ary->shape = NULL;
143
- ary->total = 0;
144
- ary->ptr = NULL;
145
- }
146
- else if ( argc == 0 ) {
147
- ary->rank = orig->rank;
148
- ary->shape = ALLOC_N(int, ary->rank);
149
- for (i=0; i<ary->rank; i++) {
150
- ary->shape[i] = (int) orig->dim[ary->rank-1-i];
151
- }
152
- ary->total = (int) orig->elements;
153
- ary->ptr = orig->ptr;
154
- }
155
- else {
156
- ca_size_t elements = 1;
157
- ary->rank = argc;
158
- ary->shape = ALLOC_N(int, ary->rank);
159
- for (i=0; i<ary->rank; i++) {
160
- ary->shape[i] = NUM2INT(argv[i]);
161
- elements *= ary->shape[i];
162
- }
163
- ary->total = (int) elements;
164
- ary->ptr = orig->ptr;
165
- if ( elements != orig->elements ) {
166
- free(ary->shape);
167
- rb_raise(rb_eRuntimeError, "elements mismatch");
168
- }
169
- }
170
-
171
- ary->type = ca_data_type_to_na_typecode(orig->data_type);
172
- ary->ref = self;
173
-
174
- return Data_Wrap_Struct(klass, cary_na_ref_mark, cary_na_ref_free, ary);
175
- }
176
-
177
- /* rdoc:
178
- class CArray
179
- def na
180
- end
181
- def nv
182
- end
183
- def nm
184
- end
185
- end
186
- */
187
-
188
- static VALUE
189
- rb_cary_na_ref_new (int argc, VALUE *argv, VALUE self)
190
- {
191
- return rb_cary_na_ref_new_i(argc, argv, self, cNArray);
192
- }
193
-
194
- static VALUE
195
- rb_cary_nv_ref_new (int argc, VALUE *argv, VALUE self)
196
- {
197
- return rb_cary_na_ref_new_i(argc, argv, self, cNVector);
198
- }
199
-
200
- static VALUE
201
- rb_cary_nm_ref_new (int argc, VALUE *argv, VALUE self)
202
- {
203
- return rb_cary_na_ref_new_i(argc, argv, self, cNMatrix);
204
- }
205
-
206
- /* -------------------------------------------------------------------- */
207
-
208
- static VALUE
209
- rb_cary_to_na_bang_i (VALUE self, VALUE klass)
210
- {
211
- volatile VALUE obj, fary;
212
- CArray *ca;
213
- struct NARRAY *na;
214
- int type;
215
-
216
- Data_Get_Struct(self, CArray, ca);
217
-
218
- type = ca_data_type_to_na_typecode(ca->data_type);
219
-
220
- if ( ca->elements == 0 ) {
221
- obj = na_make_empty(type, klass);
222
- }
223
- else {
224
- int shape[127];
225
- int i;
226
-
227
- for (i=0; i<ca->rank; i++) {
228
- shape[i] = (int) ca->dim[ca->rank-i-1];
229
- }
230
-
231
- obj = na_make_object(type, ca->rank, shape, klass);
232
- GetNArray(obj, na);
233
-
234
- fary = rb_ca_farray(self);
235
- Data_Get_Struct(fary, CArray, ca);
236
-
237
- ca_copy_data(ca, na->ptr);
238
- }
239
-
240
- return obj;
241
- }
242
-
243
- /* rdoc:
244
- class CArray
245
- def to_na!
246
- end
247
- def to_nv!
248
- end
249
- def to_nm!
250
- end
251
- end
252
- */
253
-
254
- static VALUE
255
- rb_cary_to_na_bang (VALUE self)
256
- {
257
- return rb_cary_to_na_bang_i(self, cNArray);
258
- }
259
-
260
- static VALUE
261
- rb_cary_to_nv_bang (VALUE self)
262
- {
263
- return rb_cary_to_na_bang_i(self, cNVector);
264
- }
265
-
266
- static VALUE
267
- rb_cary_to_nm_bang (VALUE self)
268
- {
269
- return rb_cary_to_na_bang_i(self, cNMatrix);
270
- }
271
-
272
- /* -------------------------------------------------------------------- */
273
-
274
- static VALUE
275
- rb_cary_to_na_i (VALUE self, VALUE klass)
276
- {
277
- volatile VALUE obj;
278
- CArray *ca;
279
- struct NARRAY *na;
280
- int type;
281
- int i;
282
-
283
- Data_Get_Struct(self, CArray, ca);
284
-
285
- type = ca_data_type_to_na_typecode(ca->data_type);
286
-
287
- if ( ca->elements == 0 ) {
288
- obj = na_make_empty(type, klass);
289
- }
290
- else {
291
- int shape[127];
292
-
293
- for (i=0; i<ca->rank; i++) {
294
- shape[i] = (int) ca->dim[ca->rank-i-1];
295
- }
296
-
297
- obj = na_make_object(type, ca->rank, shape, klass);
298
- GetNArray(obj, na);
299
-
300
- ca_copy_data(ca, na->ptr);
301
- }
302
-
303
- return obj;
304
- }
305
-
306
- /* rdoc:
307
- class CArray
308
- def to_na
309
- end
310
- def to_nv
311
- end
312
- def to_nm
313
- end
314
- end
315
- */
316
-
317
- static VALUE
318
- rb_cary_to_na (VALUE self, VALUE klass)
319
- {
320
- return rb_cary_to_na_i(self, cNArray);
321
- }
322
-
323
- static VALUE
324
- rb_cary_to_nv (VALUE self, VALUE klass)
325
- {
326
- return rb_cary_to_na_i(self, cNVector);
327
- }
328
-
329
- static VALUE
330
- rb_cary_to_nm (VALUE self, VALUE klass)
331
- {
332
- return rb_cary_to_na_i(self, cNMatrix);
333
- }
334
-
335
- /* -------------------------------------------------------------------- */
336
-
337
- /* rdoc:
338
- class NArray
339
- def ca
340
- end
341
- end
342
- */
343
-
344
- static VALUE
345
- rb_na_ca_ref_new (int argc, VALUE *argv, VALUE self)
346
- {
347
- volatile VALUE obj;
348
- struct NARRAY *na;
349
- ca_size_t dim[CA_RANK_MAX];
350
- int8_t data_type;
351
- int i;
352
-
353
- GetNArray(self, na);
354
-
355
- data_type = na_typecode_to_ca_data_type(na->type);
356
-
357
- if ( na->total == 0 ) {
358
- ca_size_t zero = 0;
359
- obj = rb_carray_wrap_ptr(data_type, 1, &zero, 0, NULL,
360
- &EMPTY_ARRAY_PTR, self); /* avoid ca->ptr == NULL */
361
- }
362
- else {
363
- if ( argc == 0 ) {
364
- CA_CHECK_RANK(na->rank);
365
- for (i=0; i<na->rank; i++) {
366
- dim[i] = na->shape[na->rank-i-1];
367
- }
368
- obj = rb_carray_wrap_ptr(data_type,
369
- na->rank, dim, 0, NULL, na->ptr, self);
370
- }
371
- else {
372
- CA_CHECK_RANK(argc);
373
- for (i=0; i<argc; i++) {
374
- dim[i] = NUM2INT(argv[i]);
375
- }
376
- obj = rb_carray_wrap_ptr(data_type, argc, dim, 0, NULL, na->ptr, self);
377
- }
378
- }
379
-
380
- return obj;
381
- }
382
-
383
- /* -------------------------------------------------------------------- */
384
-
385
- /* rdoc:
386
- class NArray
387
- def to_ca
388
- end
389
- end
390
- */
391
-
392
-
393
- static VALUE
394
- rb_na_to_ca (VALUE self)
395
- {
396
- volatile VALUE obj;
397
- CArray *ca;
398
- struct NARRAY *na;
399
- ca_size_t dim[CA_RANK_MAX];
400
- ca_size_t data_type;
401
- ca_size_t i;
402
-
403
- GetNArray(self, na);
404
-
405
- data_type = na_typecode_to_ca_data_type(na->type);
406
-
407
- if ( na->total == 0 ) {
408
- ca_size_t zero = 0;
409
- obj = rb_carray_new(data_type, 1, &zero, 0, NULL);
410
- }
411
- else {
412
- CA_CHECK_RANK(na->rank);
413
-
414
- for (i=0; i<na->rank; i++) {
415
- dim[i] = na->shape[na->rank-i-1];
416
- }
417
-
418
- obj = rb_carray_new(data_type, na->rank, dim, 0, NULL);
419
- Data_Get_Struct(obj, CArray, ca);
420
-
421
- ca_sync_data(ca, na->ptr);
422
- }
423
-
424
- return obj;
425
- }
426
-
427
- /* -------------------------------------------------------------------- */
428
-
429
- /* rdoc:
430
- class NArray
431
- def to_ca!
432
- end
433
- end
434
- */
435
-
436
- static VALUE
437
- rb_na_to_ca_bang (VALUE self)
438
- {
439
- volatile VALUE obj, fary;
440
- CArray *ca;
441
- struct NARRAY *na;
442
- ca_size_t dim[CA_RANK_MAX];
443
- ca_size_t data_type;
444
- ca_size_t i;
445
-
446
- GetNArray(self, na);
447
-
448
- data_type = na_typecode_to_ca_data_type(na->type);
449
-
450
- if ( na->total == 0 ) {
451
- ca_size_t zero = 0;
452
- obj = rb_carray_new(data_type, 1, &zero, 0, NULL);
453
- }
454
- else {
455
- CA_CHECK_RANK(na->rank);
456
-
457
- for (i=0; i<na->rank; i++) {
458
- dim[i] = na->shape[na->rank-i-1];
459
- }
460
-
461
- obj = rb_carray_new(data_type, na->rank, dim, 0, NULL);
462
- fary = rb_ca_farray(obj);
463
- Data_Get_Struct(fary, CArray, ca);
464
-
465
- ca_sync_data(ca, na->ptr);
466
- }
467
-
468
- return obj;
469
- }
470
-
471
- /* -------------------------------------------------------------------- */
472
-
473
- void
474
- Init_ca_wrap_narray ()
475
- {
476
- /* rb_require("narray"); */ /* "narray" should be loaded in config.rb */
477
-
478
- rb_define_const(rb_cCArray, "HAVE_NARRAY", Qtrue);
479
-
480
- /* CArray -> NArray */
481
- rb_define_method(rb_cCArray, "na", rb_cary_na_ref_new, -1);
482
- rb_define_method(rb_cCArray, "to_na", rb_cary_to_na, 0);
483
- rb_define_method(rb_cCArray, "to_na!", rb_cary_to_na_bang, 0);
484
-
485
- /* CArray -> NVector */
486
- rb_define_method(rb_cCArray, "nv", rb_cary_nv_ref_new, -1);
487
- rb_define_method(rb_cCArray, "to_nv", rb_cary_to_nv, 0);
488
- rb_define_method(rb_cCArray, "to_nv!", rb_cary_to_nv_bang, 0);
489
-
490
- /* CArray -> NMatrix */
491
- rb_define_method(rb_cCArray, "nm", rb_cary_nm_ref_new, -1);
492
- rb_define_method(rb_cCArray, "to_nm", rb_cary_to_nm, 0);
493
- rb_define_method(rb_cCArray, "to_nm!", rb_cary_to_nm_bang, 0);
494
-
495
- /* NArray -> CArray */
496
- rb_define_method(cNArray, "ca", rb_na_ca_ref_new, -1);
497
- rb_define_method(cNArray, "to_ca", rb_na_to_ca, 0);
498
- rb_define_method(cNArray, "to_ca!", rb_na_to_ca_bang, 0);
499
- }
500
-