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,28 +0,0 @@
1
- digraph ImageMapCallGraph {
2
- graph [size="7,7", ];
3
- line -> draw_line;
4
- poly_line -> draw_polyline;
5
- draw_polyline -> draw_line;
6
- map_points -> draw_points;
7
- draw_hline_gradation;
8
- fill_triangle -> draw_hline;
9
- draw_triangle_gradation -> draw_hline_gradation;
10
- draw_rectangle_gradation -> draw_triangle_gradation;
11
- draw_rectangle_gradation_grid -> draw_rectangle_gradation;
12
- draw_hline;
13
- fill_rectangle -> fill_triangle;
14
- fill_rectangle_image -> fill_rectangle;
15
- fill_rectangle_grid -> fill_rectangle;
16
- warp_image0 -> fill_rectangle_image;
17
- warp_image1 -> fill_rectangle_image;
18
- warp_grid -> fill_rectangle_grid;
19
- warp_image_gradation -> draw_rectangle_gradation_grid;
20
- warp_grid_gradation -> draw_rectangle_gradation_grid;
21
- warp_grid_interpolation -> warp_grid_gradation;
22
- warp -> warp_image_gradation [label="area\ngradation"];
23
- warp -> warp_image0 [label="area0"];
24
- warp -> warp_image1 [label="area1"];
25
- warp -> warp_grid_interpolation [label="point\ninterpolation"];
26
- warp -> warp_grid_gradation [label="point\ngradation"];
27
- warp -> warp_grid [label="point"];
28
- }
@@ -1,567 +0,0 @@
1
- /* ---------------------------------------------------------------------------
2
-
3
- carray/draw.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 "carray.h"
14
- #include <math.h>
15
-
16
- #define ca_ptr_1d(ap, i) \
17
- ca_ptr_at_addr((ap), (i))
18
- #define ca_ptr_2d(ap, i, j) \
19
- ca_ptr_at_addr((ap), ((CArray*)(ap))->dim[1]*(i) + (j))
20
- #define ca_ptr_3d(ap, i, j, k) \
21
- ca_ptr_at_addr(ap, \
22
- ((CArray*)(ap))->dim[2]*( ((CArray*)(ap))->dim[1]*(i) + (j) ) + (k))
23
-
24
- #define ca_ref_1d(type, ap, i) \
25
- (*(type*) ca_ptr_at_addr((ap), (i)))
26
- #define ca_ref_2d(type, ap, i, j) \
27
- (*(type*) ca_ptr_at_addr((ap), ((CArray*)(ap))->dim[1]*(i) + (j)))
28
- #define ca_ref_3d(type, ap, i, j, k) \
29
- (*(type*) ca_ptr_at_addr(ap, \
30
- ((CArray*)(ap))->dim[2]*( ((CArray*)(ap))->dim[1]*(i) + (j) ) + (k)))
31
-
32
-
33
- #define ROUND(x) (int)(((x)>0) ? floor((x)+0.5) : ((x)<0) ? ceil((x)-0.5) : 0)
34
- #define FLOOR(x) (int)(floor((float)x)+0.5)
35
- #define CEIL(x) (int)(ceil((float)x)+0.5)
36
-
37
- enum { is=0, ib=1 };
38
-
39
- void
40
- draw_hline (CArray *image, int iy, float x0, float x1, char *ptr)
41
- {
42
- int ix;
43
-
44
- if ( iy < 0 || iy >= image->dim[0] )
45
- return;
46
-
47
- for (ix=FLOOR(x0); ix<=FLOOR(x1); ix++) {
48
- if ( ix >= 0 && ix < image->dim[1] ) {
49
- memcpy(ca_ptr_2d(image, iy, ix), ptr, image->bytes);
50
- }
51
- }
52
- }
53
-
54
- void
55
- fill_triangle (CArray *image, float y[3], float x[3], char *ptr)
56
- {
57
- int imin, imid, imax;
58
- float xmin, xmax;
59
- float ax[2], x0[2], xt[2], y0;
60
- int iy;
61
- int il, ir, pass1;
62
- int i;
63
-
64
- imin = imid = imax = 0;
65
- xmin = xmax = x[0];
66
-
67
- for (i=1; i<3; i++) {
68
- if ( y[i] <= y[imin] ) {
69
- imid = imin;
70
- imin = i;
71
- }
72
- else if ( y[i] >= y[imax] ) {
73
- imid = imax;
74
- imax = i;
75
- }
76
- else
77
- imid = i;
78
-
79
- if ( x[i] <= xmin )
80
- xmin = x[i];
81
- else if ( x[i] >= x[imax] )
82
- xmax = x[i];
83
- }
84
-
85
- if ( y[imid] != y[imin] ) {
86
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
87
- ax[ib] = (x[imid] - x[imin])/(y[imid] - y[imin]);
88
- if ( ax[is] < ax[ib] ) {
89
- il = is;
90
- ir = ib;
91
- }
92
- else {
93
- il = ib;
94
- ir = is;
95
- }
96
- }
97
- else {
98
- if ( x[imin] < x[imid] ) {
99
- il = is;
100
- ir = ib;
101
- }
102
- else {
103
- il = ib;
104
- ir = is;
105
- }
106
- }
107
-
108
- /* --------------------------------------------------------
109
- iy = FLOOR(y[imin]);
110
- draw_hline(image, iy, x[imin], x[imin], ptr);
111
- ---------------------------------------------------------- */
112
-
113
- pass1 = 0;
114
-
115
- if ( FLOOR(y[imid]) != FLOOR(y[imin]) ) {
116
-
117
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
118
- ax[ib] = (x[imid] - x[imin])/(y[imid] - y[imin]);
119
- x0[is] = x[imin];
120
- x0[ib] = x[imin];
121
- y0 = y[imin];
122
-
123
- iy = CEIL(y[imin]);
124
- while ( iy <= FLOOR(y[imid]) ) {
125
- xt[is] = ax[is]*((float)iy-y0)+x0[is];
126
- xt[ib] = ax[ib]*((float)iy-y0)+x0[ib];
127
- draw_hline(image, iy, xt[il], xt[ir], ptr);
128
- iy += 1;
129
- pass1 = 1;
130
- }
131
-
132
- }
133
-
134
- if ( FLOOR(y[imax]) != FLOOR(y[imid]) ) {
135
-
136
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
137
- ax[ib] = (x[imax] - x[imid])/(y[imax] - y[imid]);
138
- x0[is] = ax[is]*(y[imid]-y[imin])+x[imin];
139
- x0[ib] = x[imid];
140
- y0 = y[imid];
141
-
142
- /* --------------------------------------------------------
143
- if ( pass1 ) {
144
- iy = FLOOR(y[imid]);
145
- xt[is] = ax[is]*(iy-y0)+x0[is];
146
- xt[ib] = x0[ib];
147
- draw_hline(image, iy, xt[il], xt[ir], ptr);
148
- }
149
- ---------------------------------------------------------- */
150
-
151
- iy = FLOOR(y[imid]);
152
- xt[is] = x0[is];
153
- xt[ib] = x0[ib];
154
- draw_hline(image, iy, xt[il], xt[ir], ptr);
155
-
156
- iy = CEIL(y[imid]);
157
- while ( iy <= FLOOR(y[imax]) ) {
158
- xt[is] = ax[is]*(iy-y0)+x0[is];
159
- xt[ib] = ax[ib]*(iy-y0)+x0[ib];
160
- draw_hline(image, iy, xt[il], xt[ir], ptr);
161
- iy += 1;
162
- }
163
-
164
- }
165
-
166
- /* --------------------------------------------------------
167
- iy = FLOOR(y[imax]);
168
- draw_hline_gradation(image, iy, x[imax], x[imax], z[imax], z[imax]);
169
- ---------------------------------------------------------- */
170
- }
171
-
172
- void
173
- fill_rectangle (CArray *image, float y[4], float x[4], char *ptr)
174
- {
175
- float yt[3], xt[3], xmin, xmax, ymin, ymax;
176
- int i;
177
-
178
- xmin = x[0];
179
- xmax = x[0];
180
- ymin = y[0];
181
- ymax = y[0];
182
-
183
- for (i=0; i<4; i++) {
184
- if ( ! ( isfinite(y[i]) && isfinite(x[i]) ) )
185
- return;
186
- if ( x[i] < xmin )
187
- xmin = x[i];
188
- else if ( x[i] > xmax )
189
- xmax = x[i];
190
- if ( y[i] < ymin )
191
- ymin = y[i];
192
- else if ( y[i] > ymax )
193
- ymax = y[i];
194
- }
195
-
196
- if ( ymin >= image->dim[0] || ymax < 0 )
197
- return;
198
-
199
- if ( xmin >= image->dim[1] || xmax < 0 )
200
- return;
201
-
202
- yt[0] = y[0]; yt[1] = y[1]; yt[2] = y[3];
203
- xt[0] = x[0]; xt[1] = x[1]; xt[2] = x[3];
204
- fill_triangle(image, yt, xt, ptr);
205
- yt[0] = y[0]; yt[1] = y[2]; yt[2] = y[3];
206
- xt[0] = x[0]; xt[1] = x[2]; xt[2] = x[3];
207
- fill_triangle(image, yt, xt, ptr);
208
- }
209
-
210
- void
211
- fill_rectangle_image (CArray *image, CArray *cy, CArray *cx, CArray *cv)
212
- {
213
- float by[4], bx[4];
214
- int i,j;
215
-
216
- for (i=0; i<cv->dim[0]; i++)
217
- for (j=0; j<cv->dim[1]; j++) {
218
- by[0] = ca_ref_2d(float, cy, i, j);
219
- by[1] = ca_ref_2d(float, cy, i, j+1);
220
- by[2] = ca_ref_2d(float, cy, i+1, j);
221
- by[3] = ca_ref_2d(float, cy, i+1, j+1);
222
- bx[0] = ca_ref_2d(float, cx, i, j);
223
- bx[1] = ca_ref_2d(float, cx, i, j+1);
224
- bx[2] = ca_ref_2d(float, cx, i+1, j);
225
- bx[3] = ca_ref_2d(float, cx, i+1, j+1);
226
- fill_rectangle(image, by, bx, ca_ptr_2d(cv, i, j));
227
- }
228
- }
229
-
230
- void
231
- fill_rectangle_grid (CArray *image, CArray *cy, CArray *cx, CArray *cv)
232
- {
233
- float by[4], bx[4], bv[4], mean;
234
- int i,j;
235
- char *val;
236
-
237
- val = ALLOC_N(char, cv->bytes);
238
-
239
- for (i=0; i<cv->dim[0]-1; i++) {
240
- for (j=0; j<cv->dim[1]-1; j++) {
241
- by[0] = ca_ref_2d(float, cy, i, j);
242
- by[1] = ca_ref_2d(float, cy, i, j+1);
243
- by[2] = ca_ref_2d(float, cy, i+1, j);
244
- by[3] = ca_ref_2d(float, cy, i+1, j+1);
245
- bx[0] = ca_ref_2d(float, cx, i, j);
246
- bx[1] = ca_ref_2d(float, cx, i, j+1);
247
- bx[2] = ca_ref_2d(float, cx, i+1, j);
248
- bx[3] = ca_ref_2d(float, cx, i+1, j+1);
249
- ca_ptr2val(cv, ca_ptr_2d(cv, i, j), CA_FLOAT, &bv[0]);
250
- ca_ptr2val(cv, ca_ptr_2d(cv, i, j+1), CA_FLOAT, &bv[1]);
251
- ca_ptr2val(cv, ca_ptr_2d(cv, i+1, j), CA_FLOAT, &bv[2]);
252
- ca_ptr2val(cv, ca_ptr_2d(cv, i+1, j+1), CA_FLOAT, &bv[3]);
253
- mean = (bv[0] + bv[1] + bv[2] + bv[3])/4.0;
254
- ca_val2ptr(CA_FLOAT, &mean, cv, val);
255
- fill_rectangle(image, by, bx, val);
256
- }
257
- }
258
- free(val);
259
- }
260
-
261
- /* ----------------------------------------------------------------- */
262
-
263
- void
264
- draw_hline_gradation (CArray *image, int iy,
265
- float x0, float x1, float z0, float z1)
266
- {
267
- int ix;
268
- float az, gz;
269
-
270
- if ( iy < 0 || iy >= image->dim[0] ) {
271
- return;
272
- }
273
-
274
- ix = FLOOR(x0);
275
- if ( ix >= 0 && ix < image->dim[1] ) {
276
- ca_val2ptr(CA_FLOAT, &z0, image, ca_ptr_2d(image, iy, ix));
277
- }
278
-
279
- if ( x1 != x0 ) {
280
- az = (z1-z0)/(x1-x0);
281
- ix = CEIL(x0);
282
- while ( ix <= FLOOR(x1) ) {
283
- gz = az * (ix - x0) + z0;
284
- if ( ix >= 0 && ix < image->dim[1] ) {
285
- ca_val2ptr(CA_FLOAT, &gz, image, ca_ptr_2d(image, iy, ix));
286
- }
287
- ix += 1;
288
- }
289
- }
290
-
291
- ix = FLOOR(x1);
292
- if ( ix >= 0 && ix < image->dim[1] ) {
293
- ca_val2ptr(CA_FLOAT, &z1, image, ca_ptr_2d(image, iy, ix));
294
- }
295
- }
296
-
297
- void
298
- draw_triangle_gradation (CArray *image, float y[3], float x[3], float z[3])
299
- {
300
- int imin, imid, imax;
301
- float xmin, xmax;
302
- float ax[2], az[2], x0[2], z0[2], xt[2], zt[2], y0;
303
- int iy;
304
- int il, ir, pass1;
305
- int i;
306
-
307
- imin = imid = imax = 0;
308
- xmin = xmax = x[0];
309
-
310
- for (i=1; i<3; i++) {
311
- if ( y[i] <= y[imin] ) {
312
- imid = imin;
313
- imin = i;
314
- }
315
- else if ( y[i] >= y[imax] ) {
316
- imid = imax;
317
- imax = i;
318
- }
319
- else
320
- imid = i;
321
-
322
- if ( x[i] <= xmin )
323
- xmin = x[i];
324
- else if ( x[i] >= x[imax] )
325
- xmax = x[i];
326
- }
327
-
328
- if ( y[imid] != y[imin] ) {
329
-
330
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
331
- ax[ib] = (x[imid] - x[imin])/(y[imid] - y[imin]);
332
-
333
- if ( ax[is] < ax[ib] ) {
334
- il = is;
335
- ir = ib;
336
- }
337
- else {
338
- il = ib;
339
- ir = is;
340
- }
341
- }
342
- else {
343
- if ( x[imin] < x[imid] ) {
344
- il = is;
345
- ir = ib;
346
- }
347
- else {
348
- il = ib;
349
- ir = is;
350
- }
351
- }
352
-
353
- /* --------------------------------------------------------
354
- iy = FLOOR(y[imin]);
355
- draw_hline_gradation(image, iy, x[imin], x[imin], z[imin], z[imin]);
356
- ---------------------------------------------------------- */
357
-
358
- pass1 = 0;
359
-
360
- if ( FLOOR(y[imid]) != FLOOR(y[imin]) ) {
361
-
362
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
363
- ax[ib] = (x[imid] - x[imin])/(y[imid] - y[imin]);
364
- az[is] = (z[imax] - z[imin])/(y[imax] - y[imin]);
365
- az[ib] = (z[imid] - z[imin])/(y[imid] - y[imin]);
366
- x0[is] = x[imin];
367
- x0[ib] = x[imin];
368
- z0[is] = z[imin];
369
- z0[ib] = z[imin];
370
- y0 = y[imin];
371
-
372
- iy = CEIL(y[imin]);
373
- while ( iy <= FLOOR(y[imid]) ) {
374
- xt[is] = ax[is]*((float)iy-y0)+x0[is];
375
- xt[ib] = ax[ib]*((float)iy-y0)+x0[ib];
376
- zt[is] = az[is]*((float)iy-y0)+z0[is];
377
- zt[ib] = az[ib]*((float)iy-y0)+z0[ib];
378
- draw_hline_gradation(image, iy, xt[il], xt[ir], zt[il], zt[ir]);
379
- iy += 1;
380
- pass1 = 1;
381
- }
382
-
383
- }
384
-
385
- if ( FLOOR(y[imax]) != FLOOR(y[imid]) ) {
386
-
387
- ax[is] = (x[imax] - x[imin])/(y[imax] - y[imin]);
388
- ax[ib] = (x[imax] - x[imid])/(y[imax] - y[imid]);
389
- az[is] = (z[imax] - z[imin])/(y[imax] - y[imin]);
390
- az[ib] = (z[imax] - z[imid])/(y[imax] - y[imid]);
391
- x0[is] = ax[is]*(y[imid]-y[imin])+x[imin];
392
- x0[ib] = x[imid];
393
- z0[is] = az[is]*(y[imid]-y[imin])+z[imin];
394
- z0[ib] = z[imid];
395
- y0 = y[imid];
396
-
397
- /* --------------------------------------------------------
398
- if ( pass1 ) {
399
- iy = FLOOR(y[imid]);
400
- xt[is] = ax[is]*(iy-y0)+x0[is];
401
- xt[ib] = x0[ib];
402
- zt[is] = az[is]*(iy-y0)+z0[is];
403
- zt[ib] = z0[ib];
404
- draw_hline_gradation(image, iy, xt[il], xt[ir], zt[il], zt[ir]);
405
- }
406
- ---------------------------------------------------------- */
407
-
408
- iy = FLOOR(y[imid]);
409
- xt[is] = x0[is];
410
- xt[ib] = x0[ib];
411
- zt[is] = z0[is];
412
- zt[ib] = z0[ib];
413
- draw_hline_gradation(image, iy, xt[il], xt[ir], zt[il], zt[ir]);
414
-
415
- iy = CEIL(y[imid]);
416
- while ( iy <= FLOOR(y[imax]) ) {
417
- xt[is] = ax[is]*(iy-y0)+x0[is];
418
- xt[ib] = ax[ib]*(iy-y0)+x0[ib];
419
- zt[is] = az[is]*(iy-y0)+z0[is];
420
- zt[ib] = az[ib]*(iy-y0)+z0[ib];
421
- draw_hline_gradation(image, iy, xt[il], xt[ir], zt[il], zt[ir]);
422
- iy += 1;
423
- }
424
-
425
- }
426
-
427
- /* --------------------------------------------------------
428
- iy = FLOOR(y[imax]);
429
- draw_hline_gradation(image, iy, x[imax], x[imax], z[imax], z[imax]);
430
- ---------------------------------------------------------- */
431
- }
432
-
433
- #include <math.h>
434
-
435
- void
436
- draw_rectangle_gradation (CArray *image, float y[4], float x[4], float z[4])
437
- {
438
- float yt[3], xt[3], zt[3], xmin, xmax, ymin, ymax;
439
- float y0, x0, z0;
440
- int i;
441
-
442
- xmin = x[0];
443
- xmax = x[0];
444
- ymin = y[0];
445
- ymax = y[0];
446
-
447
- y0 = x0 = z0 = 0;
448
- for (i=0; i<4; i++) {
449
- if ( x[i] < xmin )
450
- xmin = x[i];
451
- else if ( x[i] > xmax )
452
- xmax = x[i];
453
- if ( y[i] < ymin )
454
- ymin = y[i];
455
- else if ( y[i] > ymax )
456
- ymax = y[i];
457
- y0 += y[i]; x0 += x[i]; z0 += z[i];
458
- }
459
-
460
- if ( ymin >= image->dim[0] || ymax < 0 )
461
- return;
462
-
463
- if ( xmin >= image->dim[1] || xmax < 0 )
464
- return;
465
-
466
- if ( ! ( isfinite(y0) && isfinite(x0) ) )
467
- return;
468
-
469
- y0 /= 4.0; x0 /= 4.0; z0 /= 4.0;
470
-
471
- yt[0] = y0; xt[0] = x0; zt[0] = z0;
472
-
473
- yt[1] = y[0]; yt[2] = y[1];
474
- xt[1] = x[0]; xt[2] = x[1];
475
- zt[1] = z[0]; zt[2] = z[1];
476
- draw_triangle_gradation(image, yt, xt, zt);
477
- yt[1] = y[3];
478
- xt[1] = x[3];
479
- zt[1] = z[3];
480
- draw_triangle_gradation(image, yt, xt, zt);
481
- yt[2] = y[2];
482
- xt[2] = x[2];
483
- zt[2] = z[2];
484
- draw_triangle_gradation(image, yt, xt, zt);
485
- yt[1] = y[0];
486
- xt[1] = x[0];
487
- zt[1] = z[0];
488
- draw_triangle_gradation(image, yt, xt, zt);
489
-
490
- }
491
-
492
- void
493
- draw_rectangle_gradation_grid (CArray *image,
494
- CArray *cy, CArray *cx, CArray *cz)
495
- {
496
- float by[4], bx[4], bz[4];
497
- int i,j;
498
-
499
- for (i=0; i<cz->dim[0]-1; i++) {
500
- for (j=0; j<cz->dim[1]-1; j++) {
501
- by[0] = ca_ref_2d(float, cy, i, j);
502
- by[1] = ca_ref_2d(float, cy, i, j+1);
503
- by[2] = ca_ref_2d(float, cy, i+1, j);
504
- by[3] = ca_ref_2d(float, cy, i+1, j+1);
505
- bx[0] = ca_ref_2d(float, cx, i, j);
506
- bx[1] = ca_ref_2d(float, cx, i, j+1);
507
- bx[2] = ca_ref_2d(float, cx, i+1, j);
508
- bx[3] = ca_ref_2d(float, cx, i+1, j+1);
509
- ca_ptr2val(cz, ca_ptr_2d(cz, i, j), CA_FLOAT, &bz[0]);
510
- ca_ptr2val(cz, ca_ptr_2d(cz, i, j+1), CA_FLOAT, &bz[1]);
511
- ca_ptr2val(cz, ca_ptr_2d(cz, i+1, j), CA_FLOAT, &bz[2]);
512
- ca_ptr2val(cz, ca_ptr_2d(cz, i+1, j+1), CA_FLOAT, &bz[3]);
513
- draw_rectangle_gradation(image, by, bx, bz);
514
- }
515
- }
516
- }
517
-
518
- void
519
- draw_points (CArray *image,
520
- CArray *cy, CArray *cx, CArray *cz)
521
- {
522
- float y, x;
523
- int32_t iy, ix;
524
- int32_t i,j;
525
-
526
- ca_check_type(cz, image->data_type);
527
- ca_check_same_shape (cz, cy);
528
- ca_check_same_shape (cz, cx);
529
-
530
- for (i=0; i<cz->dim[0]; i++)
531
- for (j=0; j<cz->dim[1]; j++) {
532
- y = ca_ref_2d(float, cy, i, j);
533
- x = ca_ref_2d(float, cx, i, j);
534
- iy = ROUND(y);
535
- ix = ROUND(x);
536
- if ( ix < 0 || ix >= image->dim[1] ||
537
- iy < 0 || iy >= image->dim[0] )
538
- continue;
539
- memcpy(ca_ptr_2d(image, iy, ix), ca_ptr_2d(cz, i, j), image->bytes);
540
- }
541
- }
542
-
543
- void
544
- draw_line (CArray *ca, float y0, float x0, float y1, float x1, char *ptr);
545
-
546
- void
547
- draw_polyline (CArray *image,
548
- CArray *cy, CArray *cx, char *ptr)
549
- {
550
- float y1, x1, y2, x2;
551
- int32_t i;
552
-
553
- ca_check_same_shape (cx, cy);
554
-
555
- y1 = *(float *) ca_ptr_at_addr(cy, 0);
556
- x1 = *(float *) ca_ptr_at_addr(cx, 0);
557
-
558
- for (i=1; i<cx->elements; i++) {
559
- y2 = *(float *) ca_ptr_at_addr(cy, i);
560
- x2 = *(float *) ca_ptr_at_addr(cx, i);
561
- draw_line(image, y1, x1, y2, x2, ptr);
562
- y1 = y2;
563
- x1 = x2;
564
- }
565
-
566
- }
567
-