ruby-dcl 1.6.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (319) hide show
  1. data/ChangeLog +321 -0
  2. data/GenWrapper/Makefile +14 -0
  3. data/GenWrapper/cproto2init.rb +41 -0
  4. data/GenWrapper/dcl_narrayed_funcs.rb +181 -0
  5. data/GenWrapper/dcl_rb_footing +29 -0
  6. data/GenWrapper/dcl_rb_heading +129 -0
  7. data/GenWrapper/def.rb +20 -0
  8. data/GenWrapper/elim_ary_size.rb +398 -0
  9. data/GenWrapper/etc/Mk_proto +375 -0
  10. data/GenWrapper/etc/dcl_gen +13 -0
  11. data/GenWrapper/etc/p_header +63 -0
  12. data/GenWrapper/etc/p_init +410 -0
  13. data/GenWrapper/pparse.rb +137 -0
  14. data/GenWrapper/proto/Makefile +35 -0
  15. data/GenWrapper/proto/grph1_00 +17 -0
  16. data/GenWrapper/proto/grph1_csgi.fp +27 -0
  17. data/GenWrapper/proto/grph1_scpack.fp +442 -0
  18. data/GenWrapper/proto/grph1_sgpack.fp +1406 -0
  19. data/GenWrapper/proto/grph1_slpack.fp +316 -0
  20. data/GenWrapper/proto/grph1_stpack.fp +411 -0
  21. data/GenWrapper/proto/grph1_swpack.fp +746 -0
  22. data/GenWrapper/proto/grph2_00 +16 -0
  23. data/GenWrapper/proto/grph2_grpack.fp +139 -0
  24. data/GenWrapper/proto/grph2_ucpack.fp +522 -0
  25. data/GenWrapper/proto/grph2_udpack.fp +692 -0
  26. data/GenWrapper/proto/grph2_uepack.fp +601 -0
  27. data/GenWrapper/proto/grph2_ugpack.fp +436 -0
  28. data/GenWrapper/proto/grph2_uhpack.fp +246 -0
  29. data/GenWrapper/proto/grph2_ulpack.fp +516 -0
  30. data/GenWrapper/proto/grph2_umpack.fp +676 -0
  31. data/GenWrapper/proto/grph2_uspack.fp +989 -0
  32. data/GenWrapper/proto/grph2_uupack.fp +663 -0
  33. data/GenWrapper/proto/grph2_uvpack.fp +246 -0
  34. data/GenWrapper/proto/grph2_uwpack.fp +203 -0
  35. data/GenWrapper/proto/grph2_uxpack.fp +295 -0
  36. data/GenWrapper/proto/grph2_uypack.fp +295 -0
  37. data/GenWrapper/proto/grph2_uzpack.fp +574 -0
  38. data/GenWrapper/proto/math1_00 +12 -0
  39. data/GenWrapper/proto/math1_blklib.fp +58 -0
  40. data/GenWrapper/proto/math1_chrlib.fp +83 -0
  41. data/GenWrapper/proto/math1_fnclib.fp +79 -0
  42. data/GenWrapper/proto/math1_gnmlib.fp +118 -0
  43. data/GenWrapper/proto/math1_gt2dlib.fp +144 -0
  44. data/GenWrapper/proto/math1_ifalib.fp +123 -0
  45. data/GenWrapper/proto/math1_indxlib.fp +222 -0
  46. data/GenWrapper/proto/math1_intlib.fp +46 -0
  47. data/GenWrapper/proto/math1_lrllib.fp +276 -0
  48. data/GenWrapper/proto/math1_maplib.fp +24 -0
  49. data/GenWrapper/proto/math1_oslib.fp +28 -0
  50. data/GenWrapper/proto/math1_rfalib.fp +420 -0
  51. data/GenWrapper/proto/math1_rfblib.fp +51 -0
  52. data/GenWrapper/proto/math1_sublib.fp +60 -0
  53. data/GenWrapper/proto/math1_syslib.fp +767 -0
  54. data/GenWrapper/proto/math1_vialib.fp +339 -0
  55. data/GenWrapper/proto/math1_viblib.fp +264 -0
  56. data/GenWrapper/proto/math1_vralib.fp +339 -0
  57. data/GenWrapper/proto/math1_vrblib.fp +264 -0
  58. data/GenWrapper/proto/math1_xfclib.fp +68 -0
  59. data/GenWrapper/proto/math2_00 +12 -0
  60. data/GenWrapper/proto/math2_fftlib.fp +803 -0
  61. data/GenWrapper/proto/math2_intrlib.fp +32 -0
  62. data/GenWrapper/proto/math2_odelib.fp_notused +568 -0
  63. data/GenWrapper/proto/math2_rnmlib.fp +54 -0
  64. data/GenWrapper/proto/math2_shtlib.fp +1292 -0
  65. data/GenWrapper/proto/math2_shtlib.fp_old +1294 -0
  66. data/GenWrapper/proto/math2_vstlib.fp +84 -0
  67. data/GenWrapper/proto/misc1_00 +17 -0
  68. data/GenWrapper/proto/misc1_chnlib.fp +47 -0
  69. data/GenWrapper/proto/misc1_datelib.fp +352 -0
  70. data/GenWrapper/proto/misc1_fmtlib.fp +19 -0
  71. data/GenWrapper/proto/misc1_misclib.fp +26 -0
  72. data/GenWrapper/proto/misc1_randlib.fp +36 -0
  73. data/GenWrapper/proto/misc1_timelib.fp +151 -0
  74. data/GenWrapper/proto2c.rb +13 -0
  75. data/GenWrapper/prototype.rb +209 -0
  76. data/GenWrapper/util.rb +9 -0
  77. data/GenWrapper/variable.rb +671 -0
  78. data/README +29 -0
  79. data/Rakefile +44 -0
  80. data/ToDo +2 -0
  81. data/dcl_cary2obj.c +354 -0
  82. data/dcl_narrayed_funcs.c +518 -0
  83. data/dcl_obj2cary.c +512 -0
  84. data/dcl_rubydcloriginal.c +129 -0
  85. data/demo/gokuraku/hop/hop.rb +26 -0
  86. data/demo/gokuraku/jump/jump1.rb +64 -0
  87. data/demo/gokuraku/jump/jump2.rb +50 -0
  88. data/demo/gokuraku/layout/lay1.rb +46 -0
  89. data/demo/gokuraku/layout/lay2.rb +44 -0
  90. data/demo/gokuraku/step/step0.rb +34 -0
  91. data/demo/gokuraku/step/step1.rb +84 -0
  92. data/demo/gokuraku/step/step2.rb +62 -0
  93. data/demo/gokuraku/u1d/u1d1.rb +37 -0
  94. data/demo/gokuraku/u1d/u1d2.rb +50 -0
  95. data/demo/gokuraku/u2d/u2d1.rb +46 -0
  96. data/demo/gokuraku/u2d/u2d2.rb +49 -0
  97. data/demo/gokuraku/u2d/u2d3.rb +45 -0
  98. data/demo/gokuraku/u2d/u2d4.rb +61 -0
  99. data/demo/grph1/scpack/scpkt2.rb +40 -0
  100. data/demo/grph1/scpack/scpkt3.rb +102 -0
  101. data/demo/grph1/scpack/scpkt6.rb +84 -0
  102. data/demo/grph1/scpack/scpkt7.rb +59 -0
  103. data/demo/grph1/scpack/scpkt8.rb +60 -0
  104. data/demo/grph1/scpack/t810630.dat +296 -0
  105. data/demo/grph1/scpack/t811231.dat +296 -0
  106. data/demo/grph1/sgpack/sgfont.rb +52 -0
  107. data/demo/grph1/sgpack/sgfonz.rb +53 -0
  108. data/demo/grph1/sgpack/sgksx1.rb +44 -0
  109. data/demo/grph1/sgpack/sgksx2.rb +90 -0
  110. data/demo/grph1/sgpack/sgksx3.rb +75 -0
  111. data/demo/grph1/sgpack/sglidx.rb +40 -0
  112. data/demo/grph1/sgpack/sgltyp.rb +48 -0
  113. data/demo/grph1/sgpack/sgpk01.rb +39 -0
  114. data/demo/grph1/sgpack/sgpk02.rb +79 -0
  115. data/demo/grph1/sgpack/sgpk03.rb +107 -0
  116. data/demo/grph1/sgpack/sgpk04.rb +86 -0
  117. data/demo/grph1/sgpack/sgpk05.rb +66 -0
  118. data/demo/grph1/sgpack/sgpk06.rb +91 -0
  119. data/demo/grph1/sgpack/sgpk07.rb +66 -0
  120. data/demo/grph1/sgpack/sgpk08.rb +67 -0
  121. data/demo/grph1/sgpack/sgpk09.rb +61 -0
  122. data/demo/grph1/sgpack/sgpk10.rb +58 -0
  123. data/demo/grph1/sgpack/sgtclr.rb +59 -0
  124. data/demo/grph1/sgpack/sgtone.rb +70 -0
  125. data/demo/grph1/sgpack/sgtonz.rb +66 -0
  126. data/demo/grph1/slpack/slpk01.rb +37 -0
  127. data/demo/grph1/slpack/slpk02.rb +29 -0
  128. data/demo/grph1/slpack/slpk03.rb +36 -0
  129. data/demo/grph1/slpack/slpk04.rb +36 -0
  130. data/demo/grph2/g2pack/g2pk01.rb +78 -0
  131. data/demo/grph2/g2pack/g2pk02.rb +75 -0
  132. data/demo/grph2/grpack/grpk01.rb +65 -0
  133. data/demo/grph2/ucpack/ucpk01.rb +95 -0
  134. data/demo/grph2/ucpack/ucpk02.rb +95 -0
  135. data/demo/grph2/udegpk/u2df01.rb +50 -0
  136. data/demo/grph2/udegpk/u2df02.rb +76 -0
  137. data/demo/grph2/udegpk/u2df03.rb +80 -0
  138. data/demo/grph2/udegpk/u2df04.rb +52 -0
  139. data/demo/grph2/udegpk/u2df05.rb +70 -0
  140. data/demo/grph2/udegpk/u2df06.rb +59 -0
  141. data/demo/grph2/udegpk/u2df07.rb +68 -0
  142. data/demo/grph2/udegpk/u2df08.rb +91 -0
  143. data/demo/grph2/udegpk/u2df09b.rb +88 -0
  144. data/demo/grph2/udegpk/u2df09c.rb +89 -0
  145. data/demo/grph2/udegpk/u2df09e.rb +88 -0
  146. data/demo/grph2/udegpk/u2df09f.rb +89 -0
  147. data/demo/grph2/ulpack/ulpk01.rb +111 -0
  148. data/demo/grph2/ulpack/ulpk01n.rb +111 -0
  149. data/demo/grph2/ulpack/ulpk02.rb +111 -0
  150. data/demo/grph2/ulpack/ulpk02n.rb +111 -0
  151. data/demo/grph2/umpack/t811231.dat +296 -0
  152. data/demo/grph2/umpack/test01.rb +69 -0
  153. data/demo/grph2/umpack/test02.rb +64 -0
  154. data/demo/grph2/umpack/test03.rb +84 -0
  155. data/demo/grph2/umpack/test04.rb +54 -0
  156. data/demo/grph2/umpack/test05.rb +85 -0
  157. data/demo/grph2/umpack/test06.rb +44 -0
  158. data/demo/grph2/umpack/test07.rb +43 -0
  159. data/demo/grph2/umpack/test08.rb +37 -0
  160. data/demo/grph2/umpack/test09.rb +86 -0
  161. data/demo/grph2/umpack/test10.rb +52 -0
  162. data/demo/grph2/umpack/umpk01.rb +65 -0
  163. data/demo/grph2/umpack/umpk02.rb +54 -0
  164. data/demo/grph2/umpack/umpk03.rb +46 -0
  165. data/demo/grph2/umpack/umpk04.rb +65 -0
  166. data/demo/grph2/umpack/umpk05.rb +84 -0
  167. data/demo/grph2/uspack/uspk01.rb +39 -0
  168. data/demo/grph2/uspack/uspk02.rb +47 -0
  169. data/demo/grph2/uspack/uspk03.rb +38 -0
  170. data/demo/grph2/uspack/uspk04.rb +64 -0
  171. data/demo/grph2/uspack/uspk05.rb +68 -0
  172. data/demo/grph2/uspack/uspk06.rb +43 -0
  173. data/demo/grph2/uspack/uspk07.rb +62 -0
  174. data/demo/grph2/uspack/uspk08.rb +46 -0
  175. data/demo/grph2/uspack/uspk09.rb +89 -0
  176. data/demo/grph2/uspack/uspk10.rb +48 -0
  177. data/demo/grph2/uspack/uspk11.rb +71 -0
  178. data/demo/grph2/uspack/uspk12.rb +51 -0
  179. data/demo/grph2/uupack/uupk01.rb +50 -0
  180. data/demo/grph2/uupack/uupk02.rb +74 -0
  181. data/demo/grph2/uupack/uupk03.rb +56 -0
  182. data/demo/grph2/uupack/uupk04.rb +84 -0
  183. data/demo/grph2/uupack/uupk05.rb +74 -0
  184. data/demo/grph2/uupack/uupk06.rb +77 -0
  185. data/demo/grph2/uupack/uupk07.rb +88 -0
  186. data/demo/grph2/uxyzpk/uxyz01.rb +38 -0
  187. data/demo/grph2/uxyzpk/uxyz02.rb +36 -0
  188. data/demo/grph2/uxyzpk/uxyz03.rb +37 -0
  189. data/demo/grph2/uxyzpk/uxyz04.rb +54 -0
  190. data/demo/grph2/uxyzpk/uxyz05.rb +41 -0
  191. data/demo/grph2/uxyzpk/uxyz06.rb +42 -0
  192. data/demo/grph2/uxyzpk/uxyz07.rb +38 -0
  193. data/demo/grph2/uxyzpk/uxyz08.rb +60 -0
  194. data/demo/grph2/uxyzpk/uxyz09.rb +52 -0
  195. data/demo/grph2/uxyzpk/uxyz10.rb +57 -0
  196. data/demo/grph2/ximage/MEMO +23 -0
  197. data/demo/grph2/ximage/tomsclm.dat +1344 -0
  198. data/demo/grph2/ximage/ximg01.rb +145 -0
  199. data/demo/math1/gt2dlib/gt2d01.rb +33 -0
  200. data/demo/math2/fftlib/fftl01.rb +28 -0
  201. data/demo/math2/fftlib/fftl02.rb +224 -0
  202. data/demo/math2/intrlib/intr01.rb +22 -0
  203. data/demo/math2/rnmlib/rnml01.rb +21 -0
  204. data/demo/math2/vstlib/vstl01.rb +29 -0
  205. data/demo/math2/vstlib/vstl02.rb +38 -0
  206. data/demo/rakuraku/color/color1.rb +61 -0
  207. data/demo/rakuraku/color/color2.rb +57 -0
  208. data/demo/rakuraku/color/color3.rb +81 -0
  209. data/demo/rakuraku/kihon/kihon1.rb +39 -0
  210. data/demo/rakuraku/kihon/kihon2.rb +56 -0
  211. data/demo/rakuraku/kihon/kihon3.rb +57 -0
  212. data/demo/rakuraku/kihon/kihon4.rb +70 -0
  213. data/demo/rakuraku/kihon/kihon5.rb +56 -0
  214. data/demo/rakuraku/kihon/kihon6.rb +40 -0
  215. data/demo/rakuraku/kihon/kihon7.rb +98 -0
  216. data/demo/rakuraku/kihon/kihon8.rb +93 -0
  217. data/demo/rakuraku/kihon/kihon9.rb +77 -0
  218. data/demo/rakuraku/kihon/kihona.rb +123 -0
  219. data/demo/rakuraku/kihon/kihonb.rb +73 -0
  220. data/demo/rakuraku/kihon/kihonc.rb +110 -0
  221. data/demo/rakuraku/layout/lay1.rb +29 -0
  222. data/demo/rakuraku/layout/lay2.rb +28 -0
  223. data/demo/rakuraku/layout/lay3.rb +29 -0
  224. data/demo/rakuraku/map3d/map3d1.rb +68 -0
  225. data/demo/rakuraku/map3d/map3d2.rb +67 -0
  226. data/demo/rakuraku/map3d/map3d3.rb +66 -0
  227. data/demo/rakuraku/map3d/map3d4.rb +107 -0
  228. data/demo/rakuraku/map3d/map3d5.rb +118 -0
  229. data/demo/rakuraku/map3d/map3d6.rb +96 -0
  230. data/demo/rakuraku/map3d/map3d7.rb +98 -0
  231. data/demo/rakuraku/miss/miss1.rb +76 -0
  232. data/demo/rakuraku/miss/miss2.rb +60 -0
  233. data/demo/rakuraku/quick/quick1.rb +35 -0
  234. data/demo/rakuraku/quick/quick2.rb +33 -0
  235. data/demo/rakuraku/quick/quick3.rb +46 -0
  236. data/demo/rakuraku/quick/quick4.rb +47 -0
  237. data/demo/rakuraku/quick/quick5.rb +46 -0
  238. data/demo/rakuraku/u2d/u2d1.rb +50 -0
  239. data/demo/rakuraku/u2d/u2d2.rb +51 -0
  240. data/demo/rakuraku/u2d/u2d3.rb +70 -0
  241. data/demo/rakuraku/u2d/u2d4.rb +59 -0
  242. data/demo/rakuraku/u2d/u2d5.rb +61 -0
  243. data/demo/rakuraku/u2d/u2d6.rb +57 -0
  244. data/demo/rakuraku/u2d/u2d7.rb +103 -0
  245. data/demo/rakuraku/uspack/uspac1.rb +35 -0
  246. data/demo/rakuraku/uspack/uspac2.rb +66 -0
  247. data/demo/rakuraku/uspack/uspac3.rb +37 -0
  248. data/demo/rakuraku/uspack/uspac4.rb +49 -0
  249. data/demo/rakuraku/uxyz/uxyz1.rb +37 -0
  250. data/demo/rakuraku/uxyz/uxyz2.rb +35 -0
  251. data/demo/rakuraku/uxyz/uxyz3.rb +38 -0
  252. data/demo/rakuraku/uxyz/uxyz4.rb +36 -0
  253. data/demo/rakuraku/uxyz/uxyz5.rb +57 -0
  254. data/demo/rakuraku/uxyz/uxyz6.rb +43 -0
  255. data/demo/rakuraku/uxyz/uxyz7.rb +39 -0
  256. data/demo/rakuraku/uxyz/uxyz8.rb +56 -0
  257. data/demo/rubydcloriginal/uemrkz1.rb +37 -0
  258. data/dummy.c +9 -0
  259. data/extconf.rb +133 -0
  260. data/grph1_csgi.c +119 -0
  261. data/grph1_scpack.c +1388 -0
  262. data/grph1_sgpack.c +3876 -0
  263. data/grph1_slpack.c +432 -0
  264. data/grph1_stpack.c +1044 -0
  265. data/grph1_swpack.c +1922 -0
  266. data/grph1_zgpack.c.org +141 -0
  267. data/grph2_grpack.c +368 -0
  268. data/grph2_ucpack.c +1232 -0
  269. data/grph2_udpack.c +1404 -0
  270. data/grph2_uepack.c +1443 -0
  271. data/grph2_ugpack.c +1073 -0
  272. data/grph2_uhpack.c +943 -0
  273. data/grph2_ulpack.c +1304 -0
  274. data/grph2_umpack.c +1540 -0
  275. data/grph2_uspack.c +2481 -0
  276. data/grph2_uupack.c +1689 -0
  277. data/grph2_uvpack.c +943 -0
  278. data/grph2_uwpack.c +552 -0
  279. data/grph2_uxpack.c +555 -0
  280. data/grph2_uypack.c +555 -0
  281. data/grph2_uzpack.c +1455 -0
  282. data/init.c.default +187 -0
  283. data/init.c.gtk +189 -0
  284. data/lib/dcl.rb +5084 -0
  285. data/math1_blklib.c +227 -0
  286. data/math1_chrlib.c +239 -0
  287. data/math1_fnclib.c +254 -0
  288. data/math1_gnmlib.c +356 -0
  289. data/math1_gt2dlib.c +514 -0
  290. data/math1_ifalib.c +427 -0
  291. data/math1_indxlib.c +709 -0
  292. data/math1_intlib.c +167 -0
  293. data/math1_lrllib.c +817 -0
  294. data/math1_maplib.c +172 -0
  295. data/math1_oslib.c +111 -0
  296. data/math1_rfalib.c +1417 -0
  297. data/math1_rfblib.c +232 -0
  298. data/math1_sublib.c +221 -0
  299. data/math1_syslib.c +2025 -0
  300. data/math1_vialib.c +1156 -0
  301. data/math1_viblib.c +1027 -0
  302. data/math1_vralib.c +1156 -0
  303. data/math1_vrblib.c +1027 -0
  304. data/math1_xfclib.c +142 -0
  305. data/math2_fftlib.c +833 -0
  306. data/math2_intrlib.c +138 -0
  307. data/math2_rnmlib.c +226 -0
  308. data/math2_shtlib.c +1568 -0
  309. data/math2_vstlib.c +366 -0
  310. data/misc1_chnlib.c +179 -0
  311. data/misc1_datelib.c +953 -0
  312. data/misc1_fmtlib.c +99 -0
  313. data/misc1_misclib.c +118 -0
  314. data/misc1_randlib.c +142 -0
  315. data/misc1_timelib.c +380 -0
  316. data/obsolete/README +2 -0
  317. data/obsolete/dcl_ext.html +112 -0
  318. data/obsolete/dcl_ext.rb +313 -0
  319. metadata +433 -0
@@ -0,0 +1,21 @@
1
+ #
2
+ # rnml01.rb
3
+ #
4
+ # $Id: rnml01.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ N = 12
15
+
16
+ rx = [0, 10, 20, 30, 40, 30, 20, 10, 0, -10, -20, -10].collect{|i| i.to_f}
17
+
18
+ printf " LIST OF RX : %s\n", "%5.0f"*N % rx.to_a
19
+ ry = DCL::vrrnm(rx, N, 1, 1, 3)
20
+ print " AFTER CALLING VRRNM(RX,RY,12,1,1,3)\n"
21
+ printf " LIST OF RY : %s\n", "%5.0f"*N % ry.to_a
@@ -0,0 +1,29 @@
1
+ #
2
+ # vstl01.rb
3
+ #
4
+ # $Id: vstl01.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ IX = 3
15
+ NN = 4
16
+
17
+ x = NArray[1, 3,-1, 2,-4, 0, 2, 1, 3, 2, 1, 0]
18
+ x.reshape!(IX, NN)
19
+
20
+ print " *** TEST FOR VS1INT / VS1DIN / VS1OUT\n"
21
+ wz1, nw1 = DCL::vs1int(IX)
22
+ for j in 0..NN-1
23
+ printf " INPUT = %s\n", sprintf("%6.1f"*IX, x[0,j], x[1,j], x[2,j])
24
+ wz1, nw1 = DCL::vs1din(wz1, nw1, x[true,j])
25
+ end
26
+ wz1 = DCL::vs1out(wz1, nw1)
27
+ printf " AVE. = %s\n", "%6.1f"*IX % wz1[0..IX-1].to_a
28
+ printf " VAR. = %s\n", "%6.1f"*IX % wz1[IX..2*IX-1].to_a
29
+
@@ -0,0 +1,38 @@
1
+ #
2
+ # vstl02.rb
3
+ #
4
+ # $Id: vstl02.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ IX = 3
15
+ IY = 3
16
+ NN = 4
17
+
18
+ x = NArray[1, 3,-1, 2,-4, 0, 2, 1, 3, 2, 1, 0]
19
+ y = NArray[0, 1, 0, 1,-2, 1, 3, 0, 2, 1,-1, 0]
20
+ x.reshape!(IX, NN)
21
+ y.reshape!(IY, NN)
22
+
23
+ fmt1 = " INPUT(%s) = " + "%6.1f"*3 + "\n"
24
+ fmt2 = "%6.1f"*3 + ("\n" + " "*12 + "%6.1f"*3)*2
25
+
26
+ print " *** TEST FOR VS2INT / VS2DIN / VS2OUT\n"
27
+ wz2, nw2 = DCL::vs2int(IX, IY)
28
+ for j in 0..NN-1
29
+ printf fmt1, "X", x[0,j], x[1,j], x[2,j]
30
+ printf fmt1, "Y", y[0,j], y[1,j], y[2,j]
31
+ wz2, nw2 = DCL::vs2din(wz2, nw2, x[true,j], y[true,j])
32
+ end
33
+ wz2 = DCL::vs2out(wz2, nw2)
34
+ printf " AVE. (X) = %s\n", fmt2 % wz2.to_a.flatten[0,IX*IY]
35
+ printf " AVE. (Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY,IX*IY]
36
+ printf " VAR. (X) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*2,IX*IY]
37
+ printf " VAR. (Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*3,IX*IY]
38
+ printf " COV(X,Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*4,IX*IY]
@@ -0,0 +1,61 @@
1
+ #
2
+ # color1.rb
3
+ # $Id: color1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+ include NumRu
10
+ include Math
11
+
12
+
13
+ np = 4
14
+ nt = 10
15
+ dd = 0.07
16
+ ds = 0.003
17
+ dx = dd+ds*2
18
+ dy = 1.0/nt
19
+
20
+ #-- data ---
21
+ xbox = NArray[0.0, 1.0, 1.0, 0.0]
22
+ ybox = NArray[0.0, 0.0, 1.0, 1.0]
23
+
24
+ #-- graph ---
25
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
26
+ DCL::gliset('MAXMSG', 300)
27
+
28
+ DCL::sgopn(-(iws).abs)
29
+
30
+ DCL::sglset('LFULL', true)
31
+ DCL::slrat(0.85, 1.0)
32
+ DCL::slmgn(0.0, 0.0, 0.05, 0.1)
33
+ DCL::slsttl('TEST OF COLORMAP', 'T', 0.0, -0.5, 0.025, 1)
34
+
35
+ DCL::sgfrm
36
+ DCL::sgswnd(0.0, 1.0, 0.0, 1.0)
37
+ DCL::sgstrn(1)
38
+
39
+ DCL::sgstxs(0.01)
40
+ DCL::sgstxc(1)
41
+ for j in 0..9
42
+ for i in 0..9
43
+ x1 = dx*i + ds + (1-dx*nt)/2
44
+ y1 = dy*(9-j) + ds
45
+ x2 = x1 + dd
46
+ y2 = y1 + dd
47
+ itpat = (i+j*10)*1000 + 999
48
+
49
+ DCL::sgsvpt(x1, x2, y1, y2)
50
+ DCL::sgstrf
51
+
52
+ DCL::slpvpr(1)
53
+ DCL::sgstnp(itpat)
54
+ DCL::sgtnu(xbox, ybox)
55
+ chr = format("%5d", itpat)
56
+ DCL::sgtxv(x2, y2+(dy-dd)/2, chr)
57
+ end
58
+ end
59
+
60
+ DCL::sgcls
61
+
@@ -0,0 +1,57 @@
1
+ #
2
+ # color2.rb
3
+ # $Id: color2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+ include NumRu
10
+ include Math
11
+
12
+
13
+ nmax = 50
14
+ x = NArray.sfloat(nmax+1)
15
+ y = NArray.sfloat(nmax+1)
16
+
17
+ #-- data ---
18
+ x.indgen(1950.0, 1.0)
19
+ r = 3.7
20
+ y[0] = 0.5
21
+ for n in 0..nmax-1
22
+ y[n+1] = r*y[n]*(1.0-y[n])
23
+ end
24
+ yave = DCL::rave(y, nmax+1, 1)
25
+ y = -4 * (y - yave)
26
+
27
+ #-- graph ---
28
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
29
+ DCL::gropn iws
30
+
31
+ DCL::grfrm
32
+
33
+ DCL::grswnd(x[0], x[nmax], -1.5, 2.0)
34
+ DCL::grsvpt(0.2, 0.9, 0.2, 0.9)
35
+ DCL::grstrn(1)
36
+ DCL::grstrf
37
+
38
+ DCL::sgiset('IFONT', 2)
39
+ DCL::uziset('INDEXT1', 402)
40
+ DCL::uziset('INDEXT2', 404)
41
+ DCL::uziset('INDEXL1', 303)
42
+ DCL::uziset('INDEXL2', 123)
43
+ DCL::uxaxdv('B', 2.0, 10.0)
44
+ DCL::uxaxdv('T', 2.0, 10.0)
45
+ DCL::uxsttl('B', 'YEAR', 0.0)
46
+
47
+ DCL::uyaxdv('L', 0.1, 0.5)
48
+ DCL::uyaxdv('R', 0.1, 0.5)
49
+ cttl = DCL::csgi(131) + 'T [K]'
50
+ DCL::uysttl('L', cttl, 0.0)
51
+
52
+ DCL::uxmttl('T', 'INTERANNUAL VARIATION', 0.0)
53
+
54
+ DCL::sgspli(24)
55
+ DCL::sgplu(x, y)
56
+ DCL::grcls
57
+
@@ -0,0 +1,81 @@
1
+ #
2
+ # color3.rb
3
+ # $Id: color3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+ include NumRu
10
+ include Math
11
+
12
+
13
+ nx = 21
14
+ ny = 21
15
+ xmin = -10
16
+ xmax = 10
17
+ ymin = -10
18
+ ymax = 10
19
+ dx1 = 1
20
+ dx2 = 5
21
+ dy1 = 1
22
+ dy2 = 5
23
+ kmax = 56
24
+ pmin = 0
25
+ pmax = 1
26
+ p = NArray.sfloat(nx, ny)
27
+ pi = NArray.sfloat(2, kmax+1)
28
+
29
+ #-- data ---
30
+ for j in 0..ny-1
31
+ for i in 0..nx-1
32
+ x = xmin + (xmax-xmin)*i/(nx-1).to_f
33
+ y = ymin + (ymax-ymin)*j/(ny-1).to_f
34
+ p[i,j] = exp(-x**2/64 -y**2/25)
35
+ end
36
+ end
37
+
38
+ #-- graph ---
39
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
40
+ DCL::gropn iws
41
+
42
+ DCL::grfrm
43
+
44
+ DCL::grswnd(xmin, xmax, ymin, ymax)
45
+ DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
46
+ DCL::grstrn(1)
47
+ DCL::grstrf
48
+
49
+ dp = (pmax-pmin)/kmax.to_f
50
+ for k in 1..kmax
51
+ tlev1 = (k-1)*dp
52
+ tlev2 = tlev1 + dp
53
+ ipat = (29+k)*1000 + 999
54
+ DCL::uestlv(tlev1, tlev2, ipat)
55
+ end
56
+
57
+ DCL::uetone(p)
58
+ DCL::usdaxs
59
+
60
+ #-- �ȡ���С� ----
61
+ DCL::grswnd(0.0, 1.0, pmin, pmax)
62
+ DCL::grsvpt(0.85, 0.9, 0.2, 0.8)
63
+ DCL::grstrn(1)
64
+ DCL::grstrf
65
+
66
+ for k in 0..kmax
67
+ pi[0,k] = pmin + k*dp
68
+ pi[1,k] = pmin + k*dp
69
+ end
70
+
71
+ DCL::uetone(pi)
72
+
73
+ DCL::slpvpr(3)
74
+ DCL::uzlset('LABELYR', true)
75
+
76
+ DCL::uzfact(0.8)
77
+ DCL::uysfmt('(F4.1)')
78
+ DCL::uyaxdv('R', 0.1, 0.2)
79
+
80
+ DCL::grcls
81
+
@@ -0,0 +1,39 @@
1
+ #
2
+ # kihon1.rb
3
+ # $Id: kihon1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include NMath
11
+
12
+ nmax = 50
13
+
14
+ x = NArray.sfloat(nmax+1)
15
+ y1 = NArray.sfloat(nmax+1)
16
+ y2 = NArray.sfloat(nmax+1)
17
+ y3 = NArray.sfloat(nmax+1)
18
+
19
+ #-- data ----
20
+ x = x.indgen * 1.0/nmax
21
+ y1 = x**3
22
+ y2 = x**2
23
+ y3 = sqrt(x)
24
+
25
+ #-- graph ----
26
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
27
+ DCL::sgopn iws
28
+ DCL::sgfrm
29
+
30
+ DCL::slpvpr(1)
31
+
32
+ DCL::sgplv(x, y1)
33
+ DCL::sgpmv(x, y2)
34
+ DCL::sgtxv(0.5, 0.5, 'SGTXV')
35
+ DCL::sglset('LSOFTF', true)
36
+ DCL::sgtnv(x, y3)
37
+
38
+ DCL::sgcls
39
+
@@ -0,0 +1,56 @@
1
+ #
2
+ # kihon2.rb
3
+ # $Id: kihon2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nmax = 40
13
+ imax = 4
14
+ dt = 4.0*PI/nmax
15
+
16
+ x = NArray.sfloat(nmax+1)
17
+ y = NArray.sfloat(nmax+1, imax)
18
+
19
+ #-- data ----
20
+ x = x.indgen * 1.0/nmax
21
+ for n in 0..nmax
22
+ for i in 0..imax-1
23
+ y[n,i] = 0.2*sin(n*dt) + 0.8 - 0.2*i
24
+ end
25
+ end
26
+
27
+ #-- graph ----
28
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
29
+ DCL::sgopn iws
30
+
31
+ #-- line index : frame 1 --
32
+ DCL::sgfrm
33
+ DCL::slpvpr(1)
34
+
35
+ DCL::sgplv(x, y[true,0])
36
+
37
+ for i in 1..imax-1
38
+ DCL::sgspli(2*(i+1))
39
+ DCL::sgplv(x, y[true,i])
40
+ end
41
+
42
+ DCL::sgspli(1)
43
+
44
+ #-- line type : frame 2 --
45
+ DCL::sgfrm
46
+ DCL::slpvpr(1)
47
+
48
+ DCL::sgplv(x, y[true,0])
49
+
50
+ for i in 1..imax-1
51
+ DCL::sgsplt(i+1)
52
+ DCL::sgplv(x, y[true,i])
53
+ end
54
+
55
+ DCL::sgcls
56
+
@@ -0,0 +1,57 @@
1
+ #
2
+ # kihon3.rb
3
+ # $Id: kihon3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
4
+ #
5
+
6
+ require "narray"
7
+ require "numru/dcl"
8
+
9
+ include NumRu
10
+ include Math
11
+
12
+ nmax = 40
13
+ imax = 7
14
+ dt = 4.0*PI/nmax
15
+
16
+ x = NArray.sfloat(nmax+1)
17
+ y = NArray.sfloat(nmax+1, imax)
18
+
19
+ #-- data ----
20
+ x = x.indgen * 1.0/nmax
21
+ for n in 0..nmax
22
+ for i in 0..imax-1
23
+ y[n,i] = 0.2*sin(n*dt) + 0.8 - 0.1*i
24
+ end
25
+ end
26
+
27
+ #-- graph ----
28
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
29
+ DCL::sgopn iws
30
+
31
+ #-- marker type : frame 1 --
32
+ DCL::sgfrm
33
+ DCL::slpvpr(1)
34
+
35
+ DCL::sgpmv(x, y[true,0])
36
+
37
+ for i in 1..imax-1
38
+ itype = i + 1
39
+ DCL::sgspmt(itype)
40
+ DCL::sgpmv(x, y[true,i])
41
+ end
42
+
43
+ #-- marker size : frame 2 --
44
+ DCL::sgfrm
45
+ DCL::slpvpr(1)
46
+
47
+ DCL::sgspmt(10)
48
+ DCL::sgpmv(x, y[true,0])
49
+
50
+ for i in 1..imax-1
51
+ rsize = 0.005*i
52
+ DCL::sgspms(rsize)
53
+ DCL::sgpmv(x, y[true,i])
54
+ end
55
+
56
+ DCL::sgcls
57
+