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,88 @@
1
+ require "narray"
2
+ require "numru/dcl"
3
+ include NumRu
4
+ include Math
5
+
6
+
7
+ nx = 21
8
+ ny = 21
9
+ xmin = -10
10
+ xmax = 10
11
+ ymin = -10
12
+ ymax = 10
13
+ dx1 = 1
14
+ dx2 = 5
15
+ dy1 = 1
16
+ dy2 = 5
17
+ kmax = 5
18
+ pmin = 0
19
+ pmax = 1
20
+ u = NArray.sfloat(nx, ny)
21
+ v = NArray.sfloat(nx, ny)
22
+ p = NArray.sfloat(nx, ny)
23
+ pi = NArray.sfloat(2, kmax+1)
24
+
25
+ #-- data ---
26
+ for j in 0..ny-1
27
+ for i in 0..nx-1
28
+ x = xmin + (xmax-xmin)*i/(nx-1).to_f
29
+ y = ymin + (ymax-ymin)*j/(ny-1).to_f
30
+ u[i,j] = x
31
+ v[i,j] = -y
32
+ p[i,j] = exp(-x**2/64 -y**2/25)
33
+ end
34
+ end
35
+
36
+ #-- graph ---
37
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
38
+ DCL::gropn iws
39
+
40
+ DCL::sglset('LSOFTF', true)
41
+ DCL::grfrm
42
+
43
+ DCL::grswnd(xmin, xmax, ymin, ymax)
44
+ DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
45
+ DCL::grstrn(1)
46
+ DCL::grstrf
47
+
48
+ DCL::uxaxdv('B', dx1, dx2)
49
+ DCL::uxaxdv('T', dx1, dx2)
50
+ DCL::uxsttl('B', 'X-axis', 0.0)
51
+ DCL::uxmttl('T', 'FIGURE TITLE', 0.0)
52
+
53
+ DCL::uyaxdv('L', dy1, dy2)
54
+ DCL::uyaxdv('R', dy1, dy2)
55
+ DCL::uysttl('L', 'Y-axis', 0.0)
56
+
57
+ dp = (pmax-pmin)/kmax.to_f
58
+ for k in 1..kmax
59
+ tlev1 = (k-1)*dp
60
+ tlev2 = tlev1 + dp
61
+ ipat = 600 + k - 1
62
+ DCL::uestlv(tlev1, tlev2, ipat)
63
+ end
64
+ DCL::uetone(p)
65
+
66
+ #-- �ȡ���С� ----
67
+ DCL::grswnd(0.0, 1.0, pmin, pmax)
68
+ DCL::grsvpt(0.85, 0.9, 0.45, 0.75)
69
+ DCL::grstrn(1)
70
+ DCL::grstrf
71
+
72
+ for k in 1..kmax+1
73
+ pi[0,k-1] = pmin + (k-1)*dp
74
+ pi[1,k-1] = pmin + (k-1)*dp
75
+ end
76
+
77
+
78
+ DCL::udlset('LMSG', false)
79
+ DCL::udcntr(pi)
80
+ DCL::uetone(pi)
81
+
82
+ DCL::slpvpr(3)
83
+ DCL::uzlset('LABELYR', true)
84
+ DCL::uzfact(0.8)
85
+ DCL::uysfmt('(F4.1)')
86
+ DCL::uyaxdv('R', dp, dp)
87
+
88
+ DCL::grcls
@@ -0,0 +1,89 @@
1
+ require "narray"
2
+ require "numru/dcl"
3
+ include NumRu
4
+ include Math
5
+
6
+
7
+ nx = 21
8
+ ny = 21
9
+ xmin = -10
10
+ xmax = 10
11
+ ymin = -10
12
+ ymax = 10
13
+ dx1 = 1
14
+ dx2 = 5
15
+ dy1 = 1
16
+ dy2 = 5
17
+ kmax = 5
18
+ pmin = 0
19
+ pmax = 1
20
+ u = NArray.sfloat(nx, ny)
21
+ v = NArray.sfloat(nx, ny)
22
+ p = NArray.sfloat(nx, ny)
23
+ pi = NArray.sfloat(2, kmax+1)
24
+
25
+ #-- data ---
26
+ for j in 0..ny-1
27
+ for i in 0..nx-1
28
+ x = xmin + (xmax-xmin)*i/(nx-1).to_f
29
+ y = ymin + (ymax-ymin)*j/(ny-1).to_f
30
+ u[i,j] = x
31
+ v[i,j] = -y
32
+ p[i,j] = exp(-x**2/64 -y**2/25)
33
+ end
34
+ end
35
+
36
+ #-- graph ---
37
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
38
+ DCL::sgscmn(5)
39
+ DCL::gropn iws
40
+
41
+ DCL::sglset('LSOFTF', false)
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
+ DCL::uxaxdv('B', dx1, dx2)
50
+ DCL::uxaxdv('T', dx1, dx2)
51
+ DCL::uxsttl('B', 'X-axis', 0.0)
52
+ DCL::uxmttl('T', 'FIGURE TITLE', 0.0)
53
+
54
+ DCL::uyaxdv('L', dy1, dy2)
55
+ DCL::uyaxdv('R', dy1, dy2)
56
+ DCL::uysttl('L', 'Y-axis', 0.0)
57
+
58
+ dp = (pmax-pmin)/kmax.to_f
59
+ for k in 1..kmax
60
+ tlev1 = (k-1)*dp
61
+ tlev2 = tlev1 + dp
62
+ ipat = 999 + (99 - k*15)*1000
63
+ DCL::uestlv(tlev1, tlev2, ipat)
64
+ end
65
+ DCL::uetonf(p)
66
+
67
+ #-- �ȡ���С� ----
68
+ DCL::grswnd(0.0, 1.0, pmin, pmax)
69
+ DCL::grsvpt(0.85, 0.9, 0.45, 0.75)
70
+ DCL::grstrn(1)
71
+ DCL::grstrf
72
+
73
+ for k in 1..kmax+1
74
+ pi[0,k-1] = pmin + (k-1)*dp
75
+ pi[1,k-1] = pmin + (k-1)*dp
76
+ end
77
+
78
+
79
+ DCL::udlset('LMSG', false)
80
+ DCL::udcntr(pi)
81
+ DCL::uetonf(pi)
82
+
83
+ DCL::slpvpr(3)
84
+ DCL::uzlset('LABELYR', true)
85
+ DCL::uzfact(0.8)
86
+ DCL::uysfmt('(F4.1)')
87
+ DCL::uyaxdv('R', dp, dp)
88
+
89
+ DCL::grcls
@@ -0,0 +1,111 @@
1
+ #
2
+ # ulpk01.rb
3
+ #
4
+ # $Id: ulpk01.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ nbl = 2
15
+
16
+ #-- data ---
17
+ bl = NArray[1.0, 3.0]
18
+
19
+ #-- graph ---
20
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
21
+ DCL::sgopn(-(iws.abs))
22
+
23
+ DCL::sglset('LFULL', true)
24
+ DCL::uzfact(0.7)
25
+ DCL::slrat(0.75, 1.0)
26
+ DCL::sldiv('T', 1, 7)
27
+
28
+ DCL::sgfrm
29
+ DCL::uzinit
30
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
31
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
32
+ DCL::sgstrn(4)
33
+ DCL::sgstrf
34
+ DCL::uliset('IXTYPE', 1)
35
+ DCL::ulxlog('B', 1, 9)
36
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=1,NTCK=9)'
37
+ DCL::uxsttl('B', ctl, 0.0)
38
+
39
+ DCL::sglset('LCNTL', true)
40
+
41
+ DCL::sgfrm
42
+ DCL::uzinit
43
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
44
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
45
+ DCL::sgstrn(4)
46
+ DCL::sgstrf
47
+ DCL::uliset('IXTYPE', 1)
48
+ DCL::ulxlog('B', 3, 9)
49
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=3,NTCK=9)'
50
+ DCL::uxsttl('B', ctl, 0.0)
51
+
52
+ DCL::sgfrm
53
+ DCL::uzinit
54
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
55
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
56
+ DCL::sgstrn(4)
57
+ DCL::sgstrf
58
+ DCL::uliset('IXTYPE', 2)
59
+ DCL::ulxlog('B', 3, 9)
60
+ ctl = 'ULXLOG (IXTYPE=2,NLBL=3,NTCK=9)'
61
+ DCL::uxsttl('B', ctl, 0.0)
62
+
63
+ DCL::sgfrm
64
+ DCL::uzinit
65
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
66
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
67
+ DCL::sgstrn(4)
68
+ DCL::sgstrf
69
+ DCL::uliset('IXTYPE', 3)
70
+ DCL::ulxlog('B', 3, 9)
71
+ ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
72
+ DCL::uxsttl('B', ctl, 0.0)
73
+
74
+ DCL::sgfrm
75
+ DCL::uzinit
76
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
77
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
78
+ DCL::sgstrn(4)
79
+ DCL::sgstrf
80
+ DCL::uliset('IXTYPE', 4)
81
+ DCL::ulxlog('B', 3, 9)
82
+ ctl = 'ULXLOG (IXTYPE=4,NLBL=3,NTCK=9)'
83
+ DCL::uxsttl('B', ctl, 0.0)
84
+
85
+ DCL::sgfrm
86
+ DCL::uzinit
87
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
88
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
89
+ DCL::sgstrn(4)
90
+ DCL::sgstrf
91
+ DCL::ulxsfm('(F5.1)')
92
+ DCL::uliset('IXTYPE', 3)
93
+ DCL::ulxlog('B', 3, 9)
94
+ ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
95
+ DCL::uxsttl('B', ctl, 0.0)
96
+
97
+ DCL::sgfrm
98
+ DCL::uzinit
99
+ DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
100
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
101
+ DCL::sgstrn(4)
102
+ DCL::sgstrf
103
+ DCL::uliset('IXTYPE', 1)
104
+ DCL::uliset('IXCHR', 194)
105
+ DCL::ulsxbl(bl)
106
+ DCL::ulxlog('B', 4, 5)
107
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=4,NTCK=5)'
108
+ DCL::uxsttl('B', ctl, 0.0)
109
+
110
+ DCL::sgcls
111
+
@@ -0,0 +1,111 @@
1
+ #
2
+ # ulpk01.rb
3
+ #
4
+ # $Id: ulpk01n.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ nbl = 2
15
+
16
+ #-- data ---
17
+ bl = NArray[1.0, 3.0]
18
+
19
+ #-- graph ---
20
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
21
+ DCL::sgopn(-(iws.abs))
22
+
23
+ DCL::sglset('LFULL', true)
24
+ DCL::uzfact(0.6)
25
+ DCL::slrat(0.75, 1.0)
26
+ DCL::sldiv('T', 1, 7)
27
+
28
+ DCL::sgfrm
29
+ DCL::uzinit
30
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
31
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
32
+ DCL::sgstrn(4)
33
+ DCL::sgstrf
34
+ DCL::uliset('IXTYPE', 1)
35
+ DCL::ulxlog('B', 1, 9)
36
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=1,NTCK=9)'
37
+ DCL::uxsttl('B', ctl, 0.0)
38
+
39
+ DCL::sglset('LCNTL', true)
40
+
41
+ DCL::sgfrm
42
+ DCL::uzinit
43
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
44
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
45
+ DCL::sgstrn(4)
46
+ DCL::sgstrf
47
+ DCL::uliset('IXTYPE', 1)
48
+ DCL::ulxlog('B', 3, 9)
49
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=3,NTCK=9)'
50
+ DCL::uxsttl('B', ctl, 0.0)
51
+
52
+ DCL::sgfrm
53
+ DCL::uzinit
54
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
55
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
56
+ DCL::sgstrn(4)
57
+ DCL::sgstrf
58
+ DCL::uliset('IXTYPE', 2)
59
+ DCL::ulxlog('B', 3, 9)
60
+ ctl = 'ULXLOG (IXTYPE=2,NLBL=3,NTCK=9)'
61
+ DCL::uxsttl('B', ctl, 0.0)
62
+
63
+ DCL::sgfrm
64
+ DCL::uzinit
65
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
66
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
67
+ DCL::sgstrn(4)
68
+ DCL::sgstrf
69
+ DCL::uliset('IXTYPE', 3)
70
+ DCL::ulxlog('B', 3, 9)
71
+ ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
72
+ DCL::uxsttl('B', ctl, 0.0)
73
+
74
+ DCL::sgfrm
75
+ DCL::uzinit
76
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
77
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
78
+ DCL::sgstrn(4)
79
+ DCL::sgstrf
80
+ DCL::uliset('IXTYPE', 4)
81
+ DCL::ulxlog('B', 3, 9)
82
+ ctl = 'ULXLOG (IXTYPE=4,NLBL=3,NTCK=9)'
83
+ DCL::uxsttl('B', ctl, 0.0)
84
+
85
+ DCL::sgfrm
86
+ DCL::uzinit
87
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
88
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
89
+ DCL::sgstrn(4)
90
+ DCL::sgstrf
91
+ DCL::ulxsfm('(F5.1)')
92
+ DCL::uliset('IXTYPE', 3)
93
+ DCL::ulxlog('B', 3, 9)
94
+ ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
95
+ DCL::uxsttl('B', ctl, 0.0)
96
+
97
+ DCL::sgfrm
98
+ DCL::uzinit
99
+ DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
100
+ DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
101
+ DCL::sgstrn(4)
102
+ DCL::sgstrf
103
+ DCL::uliset('IXTYPE', 1)
104
+ DCL::uliset('IXCHR', 194)
105
+ DCL::ulsxbl(bl)
106
+ DCL::ulxlog('B', 4, 5)
107
+ ctl = 'ULXLOG (IXTYPE=1,NLBL=4,NTCK=5)'
108
+ DCL::uxsttl('B', ctl, 0.0)
109
+
110
+ DCL::sgcls
111
+
@@ -0,0 +1,111 @@
1
+ #
2
+ # ulpk02.rb
3
+ #
4
+ # $Id: ulpk02.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
5
+ #
6
+
7
+ require "narray"
8
+ require "numru/dcl"
9
+
10
+ include NumRu
11
+ include Math
12
+
13
+
14
+ nbl = 2
15
+
16
+ #-- data ---
17
+ bl = [1.0, 3.0]
18
+
19
+ #-- graph ---
20
+ iws = (ARGV[0] || (puts ' WORKSTATION ID (I) ? ;'; DCL::sgpwsn; gets)).to_i
21
+ DCL::sgopn(+iws.abs)
22
+
23
+ DCL::sglset('LFULL', true)
24
+ DCL::uzfact(0.7)
25
+ DCL::slrat(1.0, 0.75)
26
+ DCL::sldiv('Y', 7, 1)
27
+
28
+ DCL::sgfrm
29
+ DCL::uzinit
30
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
31
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
32
+ DCL::sgstrn(4)
33
+ DCL::sgstrf
34
+ DCL::uliset('IYTYPE', 1)
35
+ DCL::ulylog('L', 1, 9)
36
+ ctl = 'ULYLOG (IYTYPE=1,NLBL=1,NTCK=9)'
37
+ DCL::uysttl('L', ctl, 0.0)
38
+
39
+ DCL::sglset('LCNTL', true)
40
+
41
+ DCL::sgfrm
42
+ DCL::uzinit
43
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
44
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
45
+ DCL::sgstrn(4)
46
+ DCL::sgstrf
47
+ DCL::uliset('IYTYPE', 1)
48
+ DCL::ulylog('L', 3, 9)
49
+ ctl = 'ULYLOG (IYTYPE=1,NLBL=3,NTCK=9)'
50
+ DCL::uysttl('L', ctl, 0.0)
51
+
52
+ DCL::sgfrm
53
+ DCL::uzinit
54
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
55
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
56
+ DCL::sgstrn(4)
57
+ DCL::sgstrf
58
+ DCL::uliset('IYTYPE', 2)
59
+ DCL::ulylog('L', 3, 9)
60
+ ctl = 'ULYLOG (IYTYPE=2,NLBL=3,NTCK=9)'
61
+ DCL::uysttl('L', ctl, 0.0)
62
+
63
+ DCL::sgfrm
64
+ DCL::uzinit
65
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
66
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
67
+ DCL::sgstrn(4)
68
+ DCL::sgstrf
69
+ DCL::uliset('IYTYPE', 3)
70
+ DCL::ulylog('L', 3, 9)
71
+ ctl = 'ULYLOG (IYTYPE=3,NLBL=3,NTCK=9)'
72
+ DCL::uysttl('L', ctl, 0.0)
73
+
74
+ DCL::sgfrm
75
+ DCL::uzinit
76
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
77
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
78
+ DCL::sgstrn(4)
79
+ DCL::sgstrf
80
+ DCL::uliset('IYTYPE', 4)
81
+ DCL::ulylog('L', 3, 9)
82
+ ctl = 'ULYLOG (IYTYPE=4,NLBL=3,NTCK=9)'
83
+ DCL::uysttl('L', ctl, 0.0)
84
+
85
+ DCL::sgfrm
86
+ DCL::uzinit
87
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
88
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
89
+ DCL::sgstrn(4)
90
+ DCL::sgstrf
91
+ DCL::ulysfm('(F5.1)')
92
+ DCL::uliset('IYTYPE', 3)
93
+ DCL::ulylog('L', 3, 9)
94
+ ctl = 'ULYLOG (IYTYPE=3,NLBL=3,NTCK=9)'
95
+ DCL::uysttl('L', ctl, 0.0)
96
+
97
+ DCL::sgfrm
98
+ DCL::uzinit
99
+ DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
100
+ DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
101
+ DCL::sgstrn(4)
102
+ DCL::sgstrf
103
+ DCL::uliset('IYTYPE', 1)
104
+ DCL::uliset('IYCHR', 194)
105
+ DCL::ulsybl(bl)
106
+ DCL::ulylog('L', 4, 5)
107
+ ctl = 'ULYLOG (IYTYPE=1,NLBL=4,NTCK=5)'
108
+ DCL::uysttl('L', ctl, 0.0)
109
+
110
+ DCL::sgcls
111
+