ruby-dcl 1.7.1 → 1.8.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (327) hide show
  1. checksums.yaml +5 -5
  2. data/ChangeLog +1149 -365
  3. data/LICENSE.txt +1 -1
  4. data/Rakefile +7 -11
  5. data/demo/gokuraku/hop/hop.rb +1 -2
  6. data/demo/gokuraku/jump/jump1.rb +1 -2
  7. data/demo/gokuraku/jump/jump2.rb +2 -3
  8. data/demo/gokuraku/layout/lay1.rb +1 -2
  9. data/demo/gokuraku/layout/lay2.rb +1 -2
  10. data/demo/gokuraku/step/step0.rb +1 -2
  11. data/demo/gokuraku/step/step1.rb +1 -2
  12. data/demo/gokuraku/step/step2.rb +1 -2
  13. data/demo/gokuraku/u1d/u1d1.rb +1 -2
  14. data/demo/gokuraku/u1d/u1d2.rb +1 -2
  15. data/demo/gokuraku/u2d/u2d1.rb +1 -2
  16. data/demo/gokuraku/u2d/u2d2.rb +1 -2
  17. data/demo/gokuraku/u2d/u2d3.rb +1 -2
  18. data/demo/gokuraku/u2d/u2d4.rb +1 -2
  19. data/demo/grph1/scpack/scpkt2.rb +1 -2
  20. data/demo/grph1/scpack/scpkt3.rb +1 -2
  21. data/demo/grph1/scpack/scpkt6.rb +1 -2
  22. data/demo/grph1/scpack/scpkt7.rb +1 -2
  23. data/demo/grph1/scpack/scpkt8.rb +1 -2
  24. data/demo/grph1/sgpack/sgfont.rb +2 -3
  25. data/demo/grph1/sgpack/sgfonz.rb +2 -3
  26. data/demo/grph1/sgpack/sgksx1.rb +1 -2
  27. data/demo/grph1/sgpack/sgksx2.rb +1 -2
  28. data/demo/grph1/sgpack/sgksx3.rb +1 -2
  29. data/demo/grph1/sgpack/sglidx.rb +1 -2
  30. data/demo/grph1/sgpack/sgltyp.rb +1 -2
  31. data/demo/grph1/sgpack/sgpk01.rb +1 -2
  32. data/demo/grph1/sgpack/sgpk02.rb +1 -2
  33. data/demo/grph1/sgpack/sgpk03.rb +1 -2
  34. data/demo/grph1/sgpack/sgpk04.rb +1 -2
  35. data/demo/grph1/sgpack/sgpk05.rb +1 -2
  36. data/demo/grph1/sgpack/sgpk06.rb +1 -2
  37. data/demo/grph1/sgpack/sgpk07.rb +9 -8
  38. data/demo/grph1/sgpack/sgpk08.rb +1 -2
  39. data/demo/grph1/sgpack/sgpk09.rb +1 -2
  40. data/demo/grph1/sgpack/sgpk10.rb +1 -2
  41. data/demo/grph1/sgpack/sgtclr.rb +1 -2
  42. data/demo/grph1/sgpack/sgtone.rb +1 -2
  43. data/demo/grph1/sgpack/sgtonz.rb +1 -2
  44. data/demo/grph1/slpack/slpk01.rb +1 -2
  45. data/demo/grph1/slpack/slpk02.rb +1 -2
  46. data/demo/grph1/slpack/slpk03.rb +1 -2
  47. data/demo/grph1/slpack/slpk04.rb +2 -3
  48. data/demo/grph2/grpack/grpk01.rb +1 -2
  49. data/demo/grph2/ucpack/ucpk01.rb +1 -2
  50. data/demo/grph2/ucpack/ucpk02.rb +1 -2
  51. data/demo/grph2/udegpk/u2df01.rb +1 -2
  52. data/demo/grph2/udegpk/u2df02.rb +1 -2
  53. data/demo/grph2/udegpk/u2df03.rb +3 -4
  54. data/demo/grph2/udegpk/u2df04.rb +1 -2
  55. data/demo/grph2/udegpk/u2df05.rb +1 -2
  56. data/demo/grph2/udegpk/u2df06.rb +1 -2
  57. data/demo/grph2/udegpk/u2df07.rb +1 -2
  58. data/demo/grph2/udegpk/u2df08.rb +1 -2
  59. data/demo/grph2/udegpk/u2df09b.rb +0 -1
  60. data/demo/grph2/udegpk/u2df09c.rb +0 -1
  61. data/demo/grph2/udegpk/u2df09e.rb +0 -1
  62. data/demo/grph2/udegpk/u2df09f.rb +0 -1
  63. data/demo/grph2/uipack/uipk01.rb +1 -2
  64. data/demo/grph2/uipack/uipk02.rb +26 -18
  65. data/demo/grph2/uipack/uipk03.rb +50 -0
  66. data/demo/grph2/ulpack/ulpk01.rb +1 -2
  67. data/demo/grph2/ulpack/ulpk01n.rb +1 -2
  68. data/demo/grph2/ulpack/ulpk02.rb +1 -2
  69. data/demo/grph2/ulpack/ulpk02n.rb +1 -2
  70. data/demo/grph2/umpack/test01.rb +1 -2
  71. data/demo/grph2/umpack/test02.rb +1 -2
  72. data/demo/grph2/umpack/test03.rb +1 -2
  73. data/demo/grph2/umpack/test04.rb +2 -3
  74. data/demo/grph2/umpack/test05.rb +1 -2
  75. data/demo/grph2/umpack/test05r.rb +86 -0
  76. data/demo/grph2/umpack/test06.rb +1 -2
  77. data/demo/grph2/umpack/test07.rb +1 -2
  78. data/demo/grph2/umpack/test08.rb +1 -2
  79. data/demo/grph2/umpack/test10.rb +0 -1
  80. data/demo/grph2/umpack/umpk01.rb +1 -2
  81. data/demo/grph2/umpack/umpk02.rb +1 -2
  82. data/demo/grph2/umpack/umpk03.rb +1 -2
  83. data/demo/grph2/umpack/umpk04.rb +1 -2
  84. data/demo/grph2/umpack/umpk05.rb +1 -2
  85. data/demo/grph2/uspack/uspk01.rb +1 -2
  86. data/demo/grph2/uspack/uspk02.rb +1 -2
  87. data/demo/grph2/uspack/uspk03.rb +1 -2
  88. data/demo/grph2/uspack/uspk04.rb +1 -2
  89. data/demo/grph2/uspack/uspk05.rb +1 -2
  90. data/demo/grph2/uspack/uspk06.rb +2 -3
  91. data/demo/grph2/uspack/uspk07.rb +1 -2
  92. data/demo/grph2/uspack/uspk08.rb +2 -3
  93. data/demo/grph2/uspack/uspk09.rb +1 -2
  94. data/demo/grph2/uspack/uspk10.rb +1 -2
  95. data/demo/grph2/uspack/uspk11.rb +1 -2
  96. data/demo/grph2/uspack/uspk12.rb +1 -2
  97. data/demo/grph2/uupack/uupk01.rb +1 -2
  98. data/demo/grph2/uupack/uupk02.rb +1 -2
  99. data/demo/grph2/uupack/uupk03.rb +1 -2
  100. data/demo/grph2/uupack/uupk04.rb +1 -2
  101. data/demo/grph2/uupack/uupk05.rb +1 -2
  102. data/demo/grph2/uupack/uupk06.rb +1 -2
  103. data/demo/grph2/uupack/uupk07.rb +1 -2
  104. data/demo/grph2/uxyzpk/uxyz01.rb +1 -2
  105. data/demo/grph2/uxyzpk/uxyz02.rb +1 -2
  106. data/demo/grph2/uxyzpk/uxyz03.rb +1 -2
  107. data/demo/grph2/uxyzpk/uxyz04.rb +1 -2
  108. data/demo/grph2/uxyzpk/uxyz05.rb +1 -2
  109. data/demo/grph2/uxyzpk/uxyz06.rb +1 -2
  110. data/demo/grph2/uxyzpk/uxyz07.rb +1 -2
  111. data/demo/grph2/uxyzpk/uxyz08.rb +1 -2
  112. data/demo/grph2/uxyzpk/uxyz09.rb +1 -2
  113. data/demo/grph2/uxyzpk/uxyz10.rb +1 -2
  114. data/demo/grph2/ximage/ximg01.rb +2 -3
  115. data/demo/math1/clsplib/cslb01.rb +1 -2
  116. data/demo/math2/fftlib/fftl01.rb +1 -2
  117. data/demo/math2/fftlib/fftl02.rb +1 -2
  118. data/demo/math2/intrlib/intr01.rb +1 -2
  119. data/demo/math2/rnmlib/rnml01.rb +1 -2
  120. data/demo/math2/vstlib/vstl01.rb +1 -2
  121. data/demo/math2/vstlib/vstl02.rb +1 -2
  122. data/demo/misc2/tmpack/tmpk01.rb +46 -0
  123. data/demo/misc2/tmpack/tmpk02.rb +58 -0
  124. data/demo/misc2/tmpack/tmpk03.rb +56 -0
  125. data/demo/rakuraku/color/color1.rb +1 -2
  126. data/demo/rakuraku/color/color2.rb +5 -5
  127. data/demo/rakuraku/color/color3.rb +1 -2
  128. data/demo/rakuraku/kihon/kihon1.rb +1 -2
  129. data/demo/rakuraku/kihon/kihon2.rb +1 -2
  130. data/demo/rakuraku/kihon/kihon3.rb +1 -2
  131. data/demo/rakuraku/kihon/kihon4.rb +1 -2
  132. data/demo/rakuraku/kihon/kihon5.rb +1 -2
  133. data/demo/rakuraku/kihon/kihon6.rb +1 -2
  134. data/demo/rakuraku/kihon/kihon7.rb +10 -5
  135. data/demo/rakuraku/kihon/kihon8.rb +1 -2
  136. data/demo/rakuraku/kihon/kihon9.rb +1 -2
  137. data/demo/rakuraku/kihon/kihona.rb +1 -2
  138. data/demo/rakuraku/kihon/kihonb.rb +12 -11
  139. data/demo/rakuraku/kihon/kihonc.rb +1 -2
  140. data/demo/rakuraku/layout/lay1.rb +1 -2
  141. data/demo/rakuraku/layout/lay2.rb +1 -2
  142. data/demo/rakuraku/layout/lay3.rb +1 -2
  143. data/demo/rakuraku/map3d/map3d1.rb +1 -2
  144. data/demo/rakuraku/map3d/map3d2.rb +1 -2
  145. data/demo/rakuraku/map3d/map3d3.rb +1 -2
  146. data/demo/rakuraku/map3d/map3d4.rb +1 -2
  147. data/demo/rakuraku/map3d/map3d5.rb +1 -2
  148. data/demo/rakuraku/map3d/map3d6.rb +0 -1
  149. data/demo/rakuraku/map3d/map3d7.rb +0 -1
  150. data/demo/rakuraku/miss/miss1.rb +1 -2
  151. data/demo/rakuraku/miss/miss2.rb +1 -2
  152. data/demo/rakuraku/quick/quick1.rb +1 -2
  153. data/demo/rakuraku/quick/quick2.rb +1 -2
  154. data/demo/rakuraku/quick/quick3.rb +1 -2
  155. data/demo/rakuraku/quick/quick4.rb +1 -2
  156. data/demo/rakuraku/quick/quick5.rb +1 -2
  157. data/demo/rakuraku/u2d/u2d1.rb +1 -2
  158. data/demo/rakuraku/u2d/u2d2.rb +1 -2
  159. data/demo/rakuraku/u2d/u2d3.rb +1 -2
  160. data/demo/rakuraku/u2d/u2d4.rb +1 -2
  161. data/demo/rakuraku/u2d/u2d5.rb +1 -2
  162. data/demo/rakuraku/u2d/u2d6.rb +1 -2
  163. data/demo/rakuraku/u2d/u2d7.rb +1 -2
  164. data/demo/rakuraku/uspack/uspac1.rb +2 -3
  165. data/demo/rakuraku/uspack/uspac2.rb +1 -2
  166. data/demo/rakuraku/uspack/uspac3.rb +2 -3
  167. data/demo/rakuraku/uspack/uspac4.rb +1 -2
  168. data/demo/rakuraku/uxyz/uxyz1.rb +1 -2
  169. data/demo/rakuraku/uxyz/uxyz2.rb +1 -2
  170. data/demo/rakuraku/uxyz/uxyz3.rb +1 -2
  171. data/demo/rakuraku/uxyz/uxyz4.rb +1 -2
  172. data/demo/rakuraku/uxyz/uxyz5.rb +1 -2
  173. data/demo/rakuraku/uxyz/uxyz6.rb +1 -2
  174. data/demo/rakuraku/uxyz/uxyz7.rb +1 -2
  175. data/demo/rakuraku/uxyz/uxyz8.rb +1 -2
  176. data/demo/rubydcloriginal/uemrkz1.rb +0 -1
  177. data/ext/numru/dcl_cary2obj.c +29 -24
  178. data/ext/numru/dcl_narrayed_funcs.c +21 -21
  179. data/ext/numru/dcl_obj2cary.c +9 -5
  180. data/ext/numru/dcl_rubydcloriginal.c +14 -5
  181. data/ext/numru/extconf.rb +73 -61
  182. data/ext/numru/grph1_csgi.c +38 -8
  183. data/ext/numru/grph1_scpack.c +49 -46
  184. data/ext/numru/grph1_sgpack.c +167 -164
  185. data/ext/numru/grph1_slpack.c +23 -20
  186. data/ext/numru/grph1_stpack.c +44 -41
  187. data/ext/numru/grph1_swpack.c +155 -92
  188. data/ext/numru/grph2_grpack.c +22 -19
  189. data/ext/numru/grph2_ucpack.c +55 -52
  190. data/ext/numru/grph2_udpack.c +61 -58
  191. data/ext/numru/grph2_uepack.c +60 -57
  192. data/ext/numru/grph2_ugpack.c +49 -46
  193. data/ext/numru/grph2_uhpack.c +25 -22
  194. data/ext/numru/grph2_uipack.c +57 -54
  195. data/ext/numru/grph2_ulpack.c +59 -56
  196. data/ext/numru/grph2_umpack.c +75 -72
  197. data/ext/numru/grph2_uspack.c +97 -94
  198. data/ext/numru/grph2_uupack.c +78 -75
  199. data/ext/numru/grph2_uvpack.c +25 -22
  200. data/ext/numru/grph2_uwpack.c +30 -27
  201. data/ext/numru/grph2_uxpack.c +23 -20
  202. data/ext/numru/grph2_uypack.c +23 -20
  203. data/ext/numru/grph2_uzpack.c +69 -66
  204. data/ext/numru/init.c.default +6 -2
  205. data/ext/numru/math1_blklib.c +13 -10
  206. data/ext/numru/math1_chrlib.c +16 -13
  207. data/ext/numru/math1_clsplib.c +13 -10
  208. data/ext/numru/math1_fnclib.c +16 -13
  209. data/ext/numru/math1_gnmlib.c +21 -18
  210. data/ext/numru/math1_gt2dlib.c +18 -15
  211. data/ext/numru/math1_ifalib.c +18 -15
  212. data/ext/numru/math1_indxlib.c +24 -21
  213. data/ext/numru/math1_intlib.c +13 -10
  214. data/ext/numru/math1_lrllib.c +33 -30
  215. data/ext/numru/math1_maplib.c +14 -11
  216. data/ext/numru/math1_oslib.c +11 -8
  217. data/ext/numru/math1_rfalib.c +39 -36
  218. data/ext/numru/math1_rfblib.c +12 -9
  219. data/ext/numru/math1_sublib.c +13 -10
  220. data/ext/numru/math1_syslib.c +82 -79
  221. data/ext/numru/math1_vialib.c +31 -28
  222. data/ext/numru/math1_viblib.c +24 -21
  223. data/ext/numru/math1_vralib.c +31 -28
  224. data/ext/numru/math1_vrblib.c +24 -21
  225. data/ext/numru/math1_xfclib.c +12 -9
  226. data/ext/numru/math2_fftlib.c +28 -25
  227. data/ext/numru/math2_intrlib.c +11 -8
  228. data/ext/numru/math2_rnmlib.c +12 -9
  229. data/ext/numru/math2_shtlib.c +34 -31
  230. data/ext/numru/math2_vstlib.c +15 -12
  231. data/ext/numru/misc1_chnlib.c +12 -9
  232. data/ext/numru/misc1_datelib.c +36 -33
  233. data/ext/numru/misc1_fmtlib.c +10 -7
  234. data/ext/numru/misc1_misclib.c +11 -8
  235. data/ext/numru/misc1_randlib.c +13 -10
  236. data/ext/numru/misc1_timelib.c +20 -17
  237. data/ext/numru/misc2_tmpack.c +1200 -0
  238. data/lib/numru/dcl/dcl.rb +5676 -0
  239. data/lib/numru/dcl/version.rb +5 -0
  240. data/lib/numru/dcl.rb +6 -5580
  241. data/ruby-dcl.gemspec +13 -9
  242. metadata +25 -99
  243. data/.gitignore +0 -13
  244. data/.rspec +0 -2
  245. data/.travis.yml +0 -3
  246. data/GenWrapper/Makefile +0 -15
  247. data/GenWrapper/cproto2init.rb +0 -41
  248. data/GenWrapper/dcl_narrayed_funcs.rb +0 -181
  249. data/GenWrapper/dcl_rb_footing +0 -29
  250. data/GenWrapper/dcl_rb_heading +0 -123
  251. data/GenWrapper/def.rb +0 -20
  252. data/GenWrapper/elim_ary_size.rb +0 -398
  253. data/GenWrapper/etc/p_header +0 -64
  254. data/GenWrapper/pparse.rb +0 -137
  255. data/GenWrapper/proto/Makefile +0 -36
  256. data/GenWrapper/proto/grph1_00 +0 -17
  257. data/GenWrapper/proto/grph1_csgi.fp +0 -27
  258. data/GenWrapper/proto/grph1_scpack.fp +0 -442
  259. data/GenWrapper/proto/grph1_sgpack.fp +0 -1772
  260. data/GenWrapper/proto/grph1_slpack.fp +0 -316
  261. data/GenWrapper/proto/grph1_stpack.fp +0 -455
  262. data/GenWrapper/proto/grph1_swpack.fp +0 -826
  263. data/GenWrapper/proto/grph2_00 +0 -16
  264. data/GenWrapper/proto/grph2_grpack.fp +0 -139
  265. data/GenWrapper/proto/grph2_ucpack.fp +0 -522
  266. data/GenWrapper/proto/grph2_udpack.fp +0 -692
  267. data/GenWrapper/proto/grph2_uepack.fp +0 -601
  268. data/GenWrapper/proto/grph2_ugpack.fp +0 -436
  269. data/GenWrapper/proto/grph2_uhpack.fp +0 -246
  270. data/GenWrapper/proto/grph2_uipack.fp +0 -508
  271. data/GenWrapper/proto/grph2_ulpack.fp +0 -516
  272. data/GenWrapper/proto/grph2_umpack.fp +0 -676
  273. data/GenWrapper/proto/grph2_uspack.fp +0 -989
  274. data/GenWrapper/proto/grph2_uupack.fp +0 -663
  275. data/GenWrapper/proto/grph2_uvpack.fp +0 -246
  276. data/GenWrapper/proto/grph2_uwpack.fp +0 -203
  277. data/GenWrapper/proto/grph2_uxpack.fp +0 -295
  278. data/GenWrapper/proto/grph2_uypack.fp +0 -295
  279. data/GenWrapper/proto/grph2_uzpack.fp +0 -574
  280. data/GenWrapper/proto/math1_00 +0 -12
  281. data/GenWrapper/proto/math1_blklib.fp +0 -58
  282. data/GenWrapper/proto/math1_chrlib.fp +0 -83
  283. data/GenWrapper/proto/math1_clsplib.fp +0 -74
  284. data/GenWrapper/proto/math1_fnclib.fp +0 -79
  285. data/GenWrapper/proto/math1_gnmlib.fp +0 -118
  286. data/GenWrapper/proto/math1_gt2dlib.fp +0 -144
  287. data/GenWrapper/proto/math1_ifalib.fp +0 -123
  288. data/GenWrapper/proto/math1_indxlib.fp +0 -222
  289. data/GenWrapper/proto/math1_intlib.fp +0 -46
  290. data/GenWrapper/proto/math1_lrllib.fp +0 -276
  291. data/GenWrapper/proto/math1_maplib.fp +0 -24
  292. data/GenWrapper/proto/math1_oslib.fp +0 -28
  293. data/GenWrapper/proto/math1_rfalib.fp +0 -420
  294. data/GenWrapper/proto/math1_rfblib.fp +0 -51
  295. data/GenWrapper/proto/math1_sublib.fp +0 -60
  296. data/GenWrapper/proto/math1_syslib.fp +0 -767
  297. data/GenWrapper/proto/math1_vialib.fp +0 -339
  298. data/GenWrapper/proto/math1_viblib.fp +0 -264
  299. data/GenWrapper/proto/math1_vralib.fp +0 -339
  300. data/GenWrapper/proto/math1_vrblib.fp +0 -264
  301. data/GenWrapper/proto/math1_xfclib.fp +0 -68
  302. data/GenWrapper/proto/math2_00 +0 -12
  303. data/GenWrapper/proto/math2_fftlib.fp +0 -803
  304. data/GenWrapper/proto/math2_intrlib.fp +0 -32
  305. data/GenWrapper/proto/math2_odelib.fp_notused +0 -568
  306. data/GenWrapper/proto/math2_rnmlib.fp +0 -54
  307. data/GenWrapper/proto/math2_shtlib.fp +0 -1292
  308. data/GenWrapper/proto/math2_vstlib.fp +0 -84
  309. data/GenWrapper/proto/misc1_00 +0 -17
  310. data/GenWrapper/proto/misc1_chnlib.fp +0 -47
  311. data/GenWrapper/proto/misc1_datelib.fp +0 -352
  312. data/GenWrapper/proto/misc1_fmtlib.fp +0 -19
  313. data/GenWrapper/proto/misc1_misclib.fp +0 -26
  314. data/GenWrapper/proto/misc1_randlib.fp +0 -50
  315. data/GenWrapper/proto/misc1_timelib.fp +0 -151
  316. data/GenWrapper/proto2c.rb +0 -13
  317. data/GenWrapper/prototype.rb +0 -209
  318. data/GenWrapper/util.rb +0 -9
  319. data/GenWrapper/variable.rb +0 -671
  320. data/README +0 -29
  321. data/ToDo +0 -2
  322. data/bin/console +0 -14
  323. data/bin/setup +0 -7
  324. data/demo/grph2/ximage/MEMO +0 -23
  325. data/obsolete/README +0 -2
  326. data/obsolete/dcl_ext.html +0 -112
  327. data/obsolete/dcl_ext.rb +0 -313
@@ -1,1292 +0,0 @@
1
- /*
2
- * math2_shtlib.fp
3
- *
4
- * $Id: math2_shtlib.fp,v 1.1.1.1 2011/02/23 07:21:28 horinout Exp $
5
- */
6
-
7
- /*****************************************************
8
- * math2/shtlib/shtlib.f
9
- ****************************************************/
10
- /*
11
- * SUBROUTINE SHTLIB
12
- */
13
- SUBROUTINE SHTLIB(
14
- )
15
-
16
- /*
17
- * ENTRY SHTINT(MM,JM,IM,Q)
18
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
19
- */
20
- SUBROUTINE SHTINT(
21
- INTEGER MM
22
- INTEGER JM
23
- INTEGER IM
24
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) o WORK
25
- )
26
-
27
- /*
28
- * ENTRY SHTLAP(MM,IND,A,B)
29
- * REAL WJ(*),GJ(*),A(*),B(*)
30
- */
31
- SUBROUTINE SHTLAP(
32
- INTEGER MM
33
- INTEGER IND
34
- REAL ((MM+1)*(MM+1)) A
35
- REAL ((MM+1)*(MM+1)) o B
36
- )
37
-
38
- /*
39
- * ENTRY SHTNML(MM,N,M,LR,LI)
40
- */
41
- SUBROUTINE SHTNML(
42
- INTEGER MM
43
- INTEGER N
44
- INTEGER M
45
- INTEGER o LR
46
- INTEGER o LI
47
- )
48
-
49
- /*
50
- * ENTRY SHTFUN(MM,JM,M,FUN,Q)
51
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
52
- */
53
- SUBROUTINE SHTFUN(
54
- INTEGER MM
55
- INTEGER JM
56
- INTEGER M
57
- REAL ((2*JM+1),(MM-M+1)) o FUN
58
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
59
- )
60
-
61
- /*
62
- * ENTRY SHTLFW(MM,JM,M,ISW,WM,SM,Q)
63
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
64
- */
65
- SUBROUTINE SHTLFW(
66
- INTEGER MM
67
- INTEGER JM
68
- INTEGER M
69
- INTEGER ISW
70
- REAL (2*JM+1) WM
71
- REAL (MM-M+1) o SM
72
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
73
- )
74
-
75
- /*
76
- * ENTRY SHTLBW(MM,JM,M,ISW,SM,WM,Q)
77
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
78
- */
79
- SUBROUTINE SHTLBW(
80
- INTEGER MM
81
- INTEGER JM
82
- INTEGER M
83
- INTEGER ISW
84
- REAL (MM-M+1) SM
85
- REAL (2*JM+1) o WM
86
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
87
- )
88
-
89
- /*
90
- * ENTRY SHTS2W(MM,JM,ISW,S,W,Q)
91
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
92
- */
93
- SUBROUTINE SHTS2W(
94
- INTEGER MM
95
- INTEGER JM
96
- INTEGER ISW
97
- /* REAL S */
98
- /* REAL o W */
99
- REAL ((MM+1)*(MM+1)) S
100
- REAL ((2*JM+1),(2*MM+1)) o W
101
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
102
- )
103
-
104
- /*
105
- * ENTRY SHTSWA(MM,JM,ISW,M1,M2,S,W,Q)
106
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
107
- */
108
- SUBROUTINE SHTSWA(
109
- INTEGER MM
110
- INTEGER JM
111
- INTEGER ISW
112
- INTEGER M1
113
- INTEGER M2
114
- /* REAL S */
115
- /* REAL o W */
116
- REAL ((MM+1)*(MM+1)) S
117
- REAL ((2*JM+1),(2*MM+1)) o W
118
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
119
- )
120
-
121
- /*
122
- * ENTRY SHTSWZ(MM,JM,ISW,S,WZ,Q)
123
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
124
- */
125
- SUBROUTINE SHTSWZ(
126
- INTEGER MM
127
- INTEGER JM
128
- INTEGER ISW
129
- /* REAL S */
130
- REAL ((MM+1)*(MM+1)) S
131
- REAL (2*JM+1) o WZ
132
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
133
- )
134
-
135
- /*
136
- * ENTRY SHTSWM(MM,JM,M,ISW,S,WR,WI,Q)
137
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
138
- */
139
- SUBROUTINE SHTSWM(
140
- INTEGER MM
141
- INTEGER JM
142
- INTEGER M
143
- INTEGER ISW
144
- /* REAL S */
145
- REAL ((MM+1)*(MM+1)) S
146
- REAL (2*JM+1) o WR
147
- REAL (2*JM+1) o WI
148
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
149
- )
150
-
151
- /*
152
- * ENTRY SHTSWJ(MM,JM,ISW,J,M1,M2,S,WJ,Q)
153
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
154
- * REAL WJ(*),GJ(*),A(*),B(*)
155
- */
156
- SUBROUTINE SHTSWJ(
157
- INTEGER MM
158
- INTEGER JM
159
- INTEGER ISW
160
- INTEGER J
161
- INTEGER M1
162
- INTEGER M2
163
- /* REAL S */
164
- REAL ((MM+1)*(MM+1)) S
165
- REAL (2*MM+1) o WJ
166
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
167
- )
168
-
169
- /*
170
- * ENTRY SHTW2S(MM,JM,ISW,S,W,Q)
171
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
172
- */
173
- SUBROUTINE SHTW2S(
174
- INTEGER MM
175
- INTEGER JM
176
- INTEGER ISW
177
- /* REAL S */
178
- /* REAL o W */
179
- REAL ((2*JM+1)*(2*MM+1)) S
180
- REAL ((MM+1),(MM+1)) o W
181
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
182
- )
183
-
184
- /*
185
- * ENTRY SHTW2G(MM,JM,IM,W,G,Q)
186
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
187
- */
188
- SUBROUTINE SHTW2G(
189
- INTEGER MM
190
- INTEGER JM
191
- INTEGER IM
192
- /* REAL W */
193
- REAL ((2*JM+1),(2*MM+1)) W
194
- REAL ((2*IM+1),(2*JM+1)) o G
195
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
196
- )
197
-
198
- /*
199
- * ENTRY SHTWGA(MM,JM,IM,M1,M2,W,G,Q)
200
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
201
- */
202
- SUBROUTINE SHTWGA(
203
- INTEGER MM
204
- INTEGER JM
205
- INTEGER IM
206
- INTEGER M1
207
- INTEGER M2
208
- /* REAL W */
209
- REAL ((2*JM+1),(2*MM+1)) W
210
- REAL ((2*IM+1),(2*JM+1)) o G
211
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
212
- )
213
-
214
- /*
215
- * ENTRY SHTWGM(MM,JM,IM,M,WR,WI,G,Q)
216
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
217
- */
218
- SUBROUTINE SHTWGM(
219
- INTEGER MM
220
- INTEGER JM
221
- INTEGER IM
222
- INTEGER M
223
- REAL (2*JM+1) WR
224
- REAL (2*JM+1) WI
225
- REAL ((2*IM+1),(2*JM+1)) o G
226
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
227
- )
228
-
229
- /*
230
- * ENTRY SHTWGZ(JM,IM,WZ,G)
231
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
232
- */
233
- SUBROUTINE SHTWGZ(
234
- INTEGER JM
235
- INTEGER IM
236
- REAL (2*JM+1) WZ
237
- REAL ((2*IM+1),(2*JM+1)) o G
238
- )
239
-
240
- /*
241
- * ENTRY SHTWGJ(MM,IM,M1,M2,WJ,GJ,Q)
242
- * REAL WJ(*),GJ(*),A(*),B(*)
243
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
244
- */
245
- SUBROUTINE SHTWGJ(
246
- INTEGER MM
247
- INTEGER IM
248
- INTEGER M1
249
- INTEGER M2
250
- REAL (2*MM+1) WJ
251
- REAL (2*IM+1) o GJ
252
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
253
- )
254
-
255
- /*
256
- * ENTRY SHTG2W(MM,JM,IM,G,W,Q)
257
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
258
- */
259
- SUBROUTINE SHTG2W(
260
- INTEGER MM
261
- INTEGER JM
262
- INTEGER IM
263
- REAL ((2*IM+1),(2*JM+1)) G
264
- /* REAL o W */
265
- REAL ((2*JM+1),(2*MM+1)) o W
266
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
267
- )
268
-
269
- /*
270
- * ENTRY SHTS2G(MM,JM,IM,ISW,S,W,G,Q)
271
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
272
- */
273
- SUBROUTINE SHTS2G(
274
- INTEGER MM
275
- INTEGER JM
276
- INTEGER IM
277
- INTEGER ISW
278
- /* REAL S */
279
- /* REAL o W */
280
- REAL ((MM+1)*(MM+1)) S
281
- REAL ((2*JM+1),(2*MM+1)) o W
282
- REAL ((2*IM+1),(2*JM+1)) o G
283
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
284
- )
285
-
286
- /*
287
- * ENTRY SHTSGA(MM,JM,IM,ISW,M1,M2,S,W,G,Q)
288
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
289
- */
290
- SUBROUTINE SHTSGA(
291
- INTEGER MM
292
- INTEGER JM
293
- INTEGER IM
294
- INTEGER ISW
295
- INTEGER M1
296
- INTEGER M2
297
- /* REAL S */
298
- /* REAL o W */
299
- REAL ((MM+1)*(MM+1)) S
300
- REAL ((2*JM+1),(2*MM+1)) o W
301
- REAL ((2*IM+1),(2*JM+1)) o G
302
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
303
- )
304
-
305
- /*
306
- * ENTRY SHTSGZ(MM,JM,IM,ISW,S,WZ,G,Q)
307
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
308
- */
309
- SUBROUTINE SHTSGZ(
310
- INTEGER MM
311
- INTEGER JM
312
- INTEGER IM
313
- INTEGER ISW
314
- /* REAL S */
315
- REAL ((MM+1)*(MM+1)) S
316
- REAL (2*JM+1) o WZ
317
- REAL ((2*IM+1),(2*JM+1)) o G
318
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
319
- )
320
-
321
- /*
322
- * ENTRY SHTSGM(MM,JM,IM,M,ISW,S,WR,WI,G,Q)
323
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
324
- */
325
- SUBROUTINE SHTSGM(
326
- INTEGER MM
327
- INTEGER JM
328
- INTEGER IM
329
- INTEGER M
330
- INTEGER ISW
331
- /* REAL S */
332
- REAL ((MM+1)*(MM+1)) S
333
- REAL (2*JM+1) o WR
334
- REAL (2*JM+1) o WI
335
- REAL ((2*IM+1),(2*JM+1)) o G
336
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
337
- )
338
-
339
- /*
340
- * ENTRY SHTSGJ(MM,JM,IM,ISW,J,M1,M2,S,WJ,GJ,Q)
341
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
342
- * REAL WJ(*),GJ(*),A(*),B(*)
343
- */
344
- SUBROUTINE SHTSGJ(
345
- INTEGER MM
346
- INTEGER JM
347
- INTEGER IM
348
- INTEGER ISW
349
- INTEGER J
350
- INTEGER M1
351
- INTEGER M2
352
- /* REAL S */
353
- REAL ((MM+1)*(MM+1)) S
354
- REAL (2*MM+1) o WJ
355
- REAL (2*IM+1) o GJ
356
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
357
- )
358
-
359
- /*
360
- * ENTRY SHTG2S(MM,JM,IM,ISW,G,W,S,Q)
361
- * REAL Q(*),FUN(*),S(*),W(*),G(*),SM(*),WM(*),WZ(*),WR(*),WI(*)
362
- */
363
- SUBROUTINE SHTG2S(
364
- INTEGER MM
365
- INTEGER JM
366
- INTEGER IM
367
- INTEGER ISW
368
- REAL ((2*IM+1),(2*JM+1)) G
369
- REAL ((2*JM+1),(2*MM+1)) o W
370
- REAL ((MM+1)*(MM+1)) o S
371
- REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
372
- )
373
-
374
- /*****************************************************
375
- * math2/shtlib/shinit.f
376
- ****************************************************/
377
- /*
378
- * SUBROUTINE SHINIT(MM,JM,IM,Y,YS,YC,X,XS,XC,XW,Z,PY,PX,R,WFFT)
379
- * REAL Y(0:JM),YS(0:JM),YC(0:JM)
380
- * REAL X(0:JM),XS(0:JM),XC(0:JM),XW(0:JM)
381
- * REAL Z(JM,0:JM,4)
382
- * REAL PY(2,0:JM,0:MM)
383
- * REAL PX(2,0:JM,0:MM)
384
- * REAL R((MM+1)*(MM+1))
385
- * REAL WFFT(*)
386
- */
387
- /*
388
- SUBROUTINE SHINIT(
389
- INTEGER MM
390
- INTEGER JM
391
- INTEGER IM
392
- REAL o Y
393
- REAL (0:JM) o YS
394
- REAL (0:JM) o YC
395
- REAL o X
396
- REAL (0:JM) o XS
397
- REAL (0:JM) o XC
398
- REAL (0:JM) o XW
399
- REAL (JM,0:JM,4) o Z
400
- REAL (2,0:JM,0:MM) o PY
401
- REAL (2,0:JM,0:MM) o PX
402
- REAL ((MM+1)*(MM+1)) o R
403
- REAL (*) o WFFT
404
- )
405
- */
406
-
407
- /*****************************************************
408
- * math2/shtlib/shinix.f
409
- ****************************************************/
410
- /*
411
- * SUBROUTINE SHINIX(JM,X,XS,XC,XW)
412
- * IMPLICIT REAL*8(A-H,O-Z)
413
- * REAL X(0:JM),XS(0:JM),XC(0:JM),XW(0:JM)
414
- */
415
- /*
416
- SUBROUTINE SHINIX(
417
- INTEGER JM
418
- REAL o X
419
- REAL (0:JM) o XS
420
- REAL (0:JM) o XC
421
- REAL (0:JM) o XW
422
- )
423
- */
424
-
425
- /*****************************************************
426
- * math2/shtlib/shiniy.f
427
- ****************************************************/
428
- /*
429
- * SUBROUTINE SHINIY(JM,Y,YS,YC)
430
- * REAL Y(0:JM),YS(0:JM),YC(0:JM)
431
- */
432
- /*
433
- SUBROUTINE SHINIY(
434
- INTEGER JM
435
- REAL o Y
436
- REAL (0:JM) o YS
437
- REAL (0:JM) o YC
438
- )
439
- */
440
-
441
- /*****************************************************
442
- * math2/shtlib/shiniz.f
443
- ****************************************************/
444
- /*
445
- * SUBROUTINE SHINIZ(JM,X,Y,Z)
446
- * REAL X(0:JM),Y(0:JM)
447
- * REAL Z(JM,0:JM,4)
448
- */
449
- /*
450
- SUBROUTINE SHINIZ(
451
- INTEGER JM
452
- REAL (0:JM) X
453
- REAL (0:JM) Y
454
- REAL (JM,0:JM,4) o Z
455
- )
456
- */
457
-
458
- /*****************************************************
459
- * math2/shtlib/shinip.f
460
- ****************************************************/
461
- /*
462
- * SUBROUTINE SHINIP(MM,JM,YS,YC,PY)
463
- * REAL YS(0:JM),YC(0:JM)
464
- * REAL PY(2,0:JM,0:MM)
465
- */
466
- /*
467
- SUBROUTINE SHINIP(
468
- INTEGER MM
469
- INTEGER JM
470
- REAL (0:JM) YS
471
- REAL (0:JM) YC
472
- REAL (2,0:JM,0:MM) o PY
473
- )
474
- */
475
-
476
- /*****************************************************
477
- * math2/shtlib/shinic.f
478
- ****************************************************/
479
- /*
480
- * SUBROUTINE SHINIC(MM,JM,XW,PX)
481
- * REAL XW(0:JM),PX(2,0:JM,0:MM)
482
- */
483
- /*
484
- SUBROUTINE SHINIC(
485
- INTEGER MM
486
- INTEGER JM
487
- REAL (0:JM) XW
488
- REAL (2,0:JM,0:MM) o PX
489
- )
490
- */
491
-
492
- /*****************************************************
493
- * math2/shtlib/shinir.f
494
- ****************************************************/
495
- /*
496
- * SUBROUTINE SHINIR(MM,R)
497
- * REAL R((MM+1)*(MM+1))
498
- */
499
- /*
500
- SUBROUTINE SHINIR(
501
- INTEGER MM
502
- REAL ((MM+1)*(MM+1)) o R
503
- )
504
- */
505
-
506
- /*****************************************************
507
- * math2/shtlib/shinif.f
508
- ****************************************************/
509
- /*
510
- * SUBROUTINE SHINIF(NFFT,WFFT)
511
- * REAL WFFT(*)
512
- */
513
- /*
514
- SUBROUTINE SHINIF(
515
- INTEGER NFFT
516
- REAL (*) o WFFT
517
- )
518
- */
519
-
520
- /*****************************************************
521
- * math2/shtlib/shnm2l.f
522
- ****************************************************/
523
- /*
524
- * SUBROUTINE SHNM2L(MM,N,M,LR,LI)
525
- */
526
- /*
527
- SUBROUTINE SHNM2L(
528
- INTEGER MM
529
- INTEGER N
530
- INTEGER M
531
- INTEGER o LR
532
- INTEGER o LI
533
- )
534
- */
535
-
536
- /*****************************************************
537
- * math2/shtlib/sholap.f
538
- ****************************************************/
539
- /*
540
- * SUBROUTINE SHOLAP(MM,IND,A,B)
541
- * REAL A((MM+1)*(MM+1)),B((MM+1)*(MM+1))
542
- */
543
- /*
544
- SUBROUTINE SHOLAP(
545
- INTEGER MM
546
- INTEGER IND
547
- REAL ((MM+1)*(MM+1)) A
548
- REAL ((MM+1)*(MM+1)) o B
549
- )
550
- */
551
-
552
- /*****************************************************
553
- * math2/shtlib/shfw2g.f
554
- ****************************************************/
555
- /*
556
- * SUBROUTINE SHFW2G(MM,JM,IM,W,G,H,WFFT)
557
- * REAL W(-JM:JM,-MM:MM),G(-IM:IM,-JM:JM)
558
- * REAL H(0:2*IM-1)
559
- * REAL WFFT(*)
560
- */
561
- /*
562
- SUBROUTINE SHFW2G(
563
- INTEGER MM
564
- INTEGER JM
565
- INTEGER IM
566
- REAL (-JM:JM,-MM:MM) W
567
- REAL (-IM:IM,-JM:JM) o G
568
- REAL (0:2*IM-1) t H
569
- REAL (*) t WFFT
570
- )
571
- */
572
-
573
- /*****************************************************
574
- * math2/shtlib/shfwga.f
575
- ****************************************************/
576
- /*
577
- * SUBROUTINE SHFWGA(MM,JM,IM,M1,M2,W,G,H,WFFT)
578
- * REAL W(-JM:JM,-MM:MM),G(-IM:IM,-JM:JM)
579
- * REAL H(0:2*IM-1)
580
- * REAL WFFT(*)
581
- */
582
- /*
583
- SUBROUTINE SHFWGA(
584
- INTEGER MM
585
- INTEGER JM
586
- INTEGER IM
587
- INTEGER M1
588
- INTEGER M2
589
- REAL (-JM:JM,-MM:MM) W
590
- REAL (-IM:IM,-JM:JM) o G
591
- REAL (0:2*IM-1) t H
592
- REAL (*) t WFFT
593
- )
594
- */
595
-
596
- /*****************************************************
597
- * math2/shtlib/shfwgm.f
598
- ****************************************************/
599
- /*
600
- * SUBROUTINE SHFWGM(JM,IM,M,WR,WI,G,H,WFFT)
601
- * REAL WR(-JM:JM),WI(-JM:JM),G(-IM:IM,-JM:JM)
602
- * REAL H(0:2*IM-1)
603
- * REAL WFFT(*)
604
- */
605
- /* SUBROUTINE SHFWGM(
606
- INTEGER JM
607
- INTEGER IM
608
- INTEGER M
609
- REAL (-JM:JM) WR
610
- REAL (-JM:JM) WI
611
- REAL (-IM:IM,-JM:JM) o G
612
- REAL (0:2*IM-1) t H
613
- REAL (*) t WFFT
614
- )
615
- */
616
-
617
- /*****************************************************
618
- * math2/shtlib/shfwgz.f
619
- ****************************************************/
620
- /*
621
- * SUBROUTINE SHFWGZ(JM,IM,WZ,G)
622
- * REAL WZ(-JM:JM),G(-IM:IM,-JM:JM)
623
- */
624
- /*
625
- SUBROUTINE SHFWGZ(
626
- INTEGER JM
627
- INTEGER IM
628
- REAL (-JM:JM) WZ
629
- REAL (-IM:IM,-JM:JM) o G
630
- )
631
- */
632
-
633
- /*****************************************************
634
- * math2/shtlib/shfwgj.f
635
- ****************************************************/
636
- /*
637
- * SUBROUTINE SHFWGJ(MM,IM,M1,M2,WJ,GJ,H,WFFT)
638
- * REAL WJ(-MM:MM),GJ(-IM:IM)
639
- * REAL H(0:2*IM-1)
640
- * REAL WFFT(*)
641
- */
642
- /*
643
- SUBROUTINE SHFWGJ(
644
- INTEGER MM
645
- INTEGER IM
646
- INTEGER M1
647
- INTEGER M2
648
- REAL (-MM:MM) WJ
649
- REAL (-IM:IM) o GJ
650
- REAL (0:2*IM-1) t H
651
- REAL (*) t WFFT
652
- )
653
- */
654
-
655
- /*****************************************************
656
- * math2/shtlib/shfg2w.f
657
- ****************************************************/
658
- /*
659
- * SUBROUTINE SHFG2W(MM,JM,IM,G,W,H,WFFT)
660
- * REAL G(-IM:IM,-JM:JM),W(-JM:JM,-MM:MM)
661
- * REAL H(0:2*IM-1)
662
- * REAL WFFT(*)
663
- */
664
- /*
665
- SUBROUTINE SHFG2W(
666
- INTEGER MM
667
- INTEGER JM
668
- INTEGER IM
669
- REAL (-IM:IM,-JM:JM) G
670
- REAL (-JM:JM,-MM:MM) o W
671
- REAL (0:2*IM-1) t H
672
- REAL (*) t WFFT
673
- )
674
- */
675
-
676
- /*****************************************************
677
- * math2/shtlib/shfftf.f
678
- ****************************************************/
679
- /*
680
- * SUBROUTINE SHFFTF(NFFT,XFFT,WFFT)
681
- * REAL XFFT(NFFT)
682
- * REAL WFFT(*)
683
- */
684
- /*
685
- SUBROUTINE SHFFTF(
686
- INTEGER NFFT
687
- REAL (NFFT) io XFFT
688
- REAL (*) t WFFT
689
- )
690
- */
691
-
692
- /*****************************************************
693
- * math2/shtlib/shfftb.f
694
- ****************************************************/
695
- /*
696
- * SUBROUTINE SHFFTB(NFFT,XFFT,WFFT)
697
- * REAL XFFT(NFFT)
698
- * REAL WFFT(*)
699
- */
700
- /*
701
- SUBROUTINE SHFFTB(
702
- INTEGER NFFT
703
- REAL (NFFT) io XFFT
704
- REAL (*) t WFFT
705
- )
706
- */
707
-
708
- /*****************************************************
709
- * math2/shtlib/shms2w.f
710
- ****************************************************/
711
- /*
712
- * SUBROUTINE SHMS2W(MM,JM,ISW,S,W,SD,PM,YS,YC,PY,R)
713
- * REAL S((MM+1)*(MM+1)),W(-JM:JM,-MM:MM)
714
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
715
- * REAL YS(0:JM),YC(0:JM)
716
- * REAL PY(2,0:JM,0:MM)
717
- * REAL R((MM+1)*(MM+1))
718
- */
719
- /*
720
- SUBROUTINE SHMS2W(
721
- INTEGER MM
722
- INTEGER JM
723
- INTEGER ISW
724
- REAL ((MM+1)*(MM+1)) S
725
- REAL (-JM:JM,-MM:MM) o W
726
- REAL (0:MM+1) t SD
727
- REAL (0:MM+1,0:JM) t PM
728
- REAL (0:JM) t YS
729
- REAL (0:JM) t YC
730
- REAL (2,0:JM,0:MM) t PY
731
- REAL ((MM+1)*(MM+1)) t R
732
- )
733
- */
734
-
735
- /*****************************************************
736
- * math2/shtlib/shmswa.f
737
- ****************************************************/
738
- /*
739
- * SUBROUTINE SHMSWA(MM,JM,ISW,M1,M2,S,W,SD,PM,YS,YC,PY,R)
740
- * REAL S((MM+1)*(MM+1)),W(-JM:JM,-MM:MM)
741
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
742
- * REAL YS(0:JM),YC(0:JM)
743
- * REAL PY(2,0:JM,0:MM)
744
- * REAL R((MM+1)*(MM+1))
745
- */
746
- /*
747
- SUBROUTINE SHMSWA(
748
- INTEGER MM
749
- INTEGER JM
750
- INTEGER ISW
751
- INTEGER M1
752
- INTEGER M2
753
- REAL ((MM+1)*(MM+1)) S
754
- REAL (-JM:JM,-MM:MM) o W
755
- REAL (0:MM+1) t SD
756
- REAL (0:MM+1,0:JM) t PM
757
- REAL (0:JM) t YS
758
- REAL (0:JM) t YC
759
- REAL (2,0:JM,0:MM) t PY
760
- REAL ((MM+1)*(MM+1)) t R
761
- )
762
- */
763
-
764
- /*****************************************************
765
- * math2/shtlib/shmswz.f
766
- ****************************************************/
767
- /*
768
- * SUBROUTINE SHMSWZ(MM,JM,ISW,S,WZ,SD,PM,YS,YC,PY,R)
769
- * REAL S((MM+1)*(MM+1)),WZ(-JM:JM)
770
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
771
- * REAL YS(0:JM),YC(0:JM)
772
- * REAL PY(2,0:JM,0:MM)
773
- * REAL R((MM+1)*(MM+1))
774
- */
775
- /*
776
- SUBROUTINE SHMSWZ(
777
- INTEGER MM
778
- INTEGER JM
779
- INTEGER ISW
780
- REAL ((MM+1)*(MM+1)) S
781
- REAL (-JM:JM) o WZ
782
- REAL (0:MM+1) t SD
783
- REAL (0:MM+1,0:JM) t PM
784
- REAL (0:JM) t YS
785
- REAL (0:JM) t YC
786
- REAL (2,0:JM,0:MM) t PY
787
- REAL ((MM+1)*(MM+1)) t R
788
- )
789
- */
790
-
791
- /*****************************************************
792
- * math2/shtlib/shmswm.f
793
- ****************************************************/
794
- /*
795
- * SUBROUTINE SHMSWM(MM,JM,M,ISW,S,WR,WI,SD,PM,YS,YC,PY,R)
796
- * REAL S((MM+1)*(MM+1)),WR(-JM:JM),WI(-JM:JM)
797
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
798
- * REAL YS(0:JM),YC(0:JM)
799
- * REAL PY(2,0:JM,0:MM)
800
- * REAL R((MM+1)*(MM+1))
801
- */
802
- /*
803
- SUBROUTINE SHMSWM(
804
- INTEGER MM
805
- INTEGER JM
806
- INTEGER M
807
- INTEGER ISW
808
- REAL ((MM+1)*(MM+1)) S
809
- REAL (-JM:JM) o WR
810
- REAL (-JM:JM) o WI
811
- REAL (0:MM+1) t SD
812
- REAL (0:MM+1,0:JM) t PM
813
- REAL (0:JM) t YS
814
- REAL (0:JM) t YC
815
- REAL (2,0:JM,0:MM) t PY
816
- REAL ((MM+1)*(MM+1)) t R
817
- )
818
- */
819
-
820
- /*****************************************************
821
- * math2/shtlib/shmdxm.f
822
- ****************************************************/
823
- /*
824
- * SUBROUTINE SHMDXM(JM,WR,WI)
825
- * REAL WR(-JM:JM),WI(-JM:JM)
826
- */
827
- /*
828
- SUBROUTINE SHMDXM(
829
- INTEGER JM
830
- REAL (-JM:JM) io WR
831
- REAL (-JM:JM) io WI
832
- )
833
- */
834
-
835
- /*****************************************************
836
- * math2/shtlib/shmswj.f
837
- ****************************************************/
838
- /*
839
- * SUBROUTINE SHMSWJ(MM,JM,ISW,J,M1,M2,S,WJ,SD,PMJ,YS,YC,PY,R)
840
- * REAL S((MM+1)*(MM+1)),WJ(-MM:MM)
841
- * REAL SD(0:MM+1),PMJ(0:MM+1)
842
- * REAL YS(0:JM),YC(0:JM)
843
- * REAL PY(2,0:JM,0:MM)
844
- * REAL R((MM+1)*(MM+1))
845
- */
846
- /*
847
- SUBROUTINE SHMSWJ(
848
- INTEGER MM
849
- INTEGER JM
850
- INTEGER ISW
851
- INTEGER J
852
- INTEGER M1
853
- INTEGER M2
854
- REAL ((MM+1)*(MM+1)) S
855
- REAL (-MM:MM) o WJ
856
- REAL (0:MM+1) t SD
857
- REAL (0:MM+1) t PMJ
858
- REAL (0:JM) t YS
859
- REAL (0:JM) t YC
860
- REAL (2,0:JM,0:MM) t PY
861
- REAL ((MM+1)*(MM+1)) t R
862
- )
863
- */
864
-
865
- /*****************************************************
866
- * math2/shtlib/shmwjz.f
867
- ****************************************************/
868
- /*
869
- * SUBROUTINE SHMWJZ(MM,JM,ISW,J,S,WJD,SD,PMJ,YS,YC,PY,R)
870
- * REAL S((MM+1)*(MM+1)),WJD
871
- * REAL SD(0:MM+1),PMJ(0:MM+1)
872
- * REAL YS(0:JM),YC(0:JM)
873
- * REAL PY(2,0:JM,0:MM)
874
- * REAL R((MM+1)*(MM+1))
875
- */
876
- /*
877
- SUBROUTINE SHMWJZ(
878
- INTEGER MM
879
- INTEGER JM
880
- INTEGER ISW
881
- INTEGER J
882
- REAL ((MM+1)*(MM+1)) S
883
- REAL o WJD
884
- REAL (0:MM+1) t SD
885
- REAL (0:MM+1) t PMJ
886
- REAL (0:JM) t YS
887
- REAL (0:JM) t YC
888
- REAL (2,0:JM,0:MM) t PY
889
- REAL ((MM+1)*(MM+1)) t R
890
- )
891
- */
892
-
893
- /*****************************************************
894
- * math2/shtlib/shmwjm.f
895
- ****************************************************/
896
- /*
897
- * SUBROUTINE SHMWJM(MM,JM,ISW,J,M,S,WJR,WJI,SD,PMJ,YS,YC,PY,R)
898
- * REAL S((MM+1)*(MM+1)),WJR,WJI
899
- * REAL SD(0:MM+1),PMJ(0:MM+1)
900
- * REAL YS(0:JM),YC(0:JM)
901
- * REAL PY(2,0:JM,0:MM)
902
- * REAL R((MM+1)*(MM+1))
903
- */
904
- /*
905
- SUBROUTINE SHMWJM(
906
- INTEGER MM
907
- INTEGER JM
908
- INTEGER ISW
909
- INTEGER J
910
- INTEGER M
911
- REAL ((MM+1)*(MM+1)) S
912
- REAL o WJR
913
- REAL o WJI
914
- REAL (0:MM+1) t SD
915
- REAL (0:MM+1) t PMJ
916
- REAL (0:JM) t YS
917
- REAL (0:JM) t YC
918
- REAL (2,0:JM,0:MM) t PY
919
- REAL ((MM+1)*(MM+1)) t R
920
- )
921
- */
922
-
923
- /*****************************************************
924
- * math2/shtlib/shmdxj.f
925
- ****************************************************/
926
- /*
927
- * SUBROUTINE SHMDXJ(WRD,WID)
928
- * REAL WRD,WID
929
- */
930
- /*
931
- SUBROUTINE SHMDXJ(
932
- REAL io WRD
933
- REAL io WID
934
- )
935
- */
936
-
937
- /*****************************************************
938
- * math2/shtlib/shmw2s.f
939
- ****************************************************/
940
- /*
941
- * SUBROUTINE SHMW2S(MM,JM,ISW,W,S,SD,PM,XS,XC,PX,WY,WX,R,Z)
942
- * REAL W(-JM:JM,-MM:MM),S((MM+1)*(MM+1))
943
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
944
- * REAL XS(0:JM),XC(0:JM)
945
- * REAL PX(2,0:JM,0:MM)
946
- * REAL WY(0:JM,2),WX(0:JM,2)
947
- * REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
948
- */
949
- /*
950
- SUBROUTINE SHMW2S(
951
- INTEGER MM
952
- INTEGER JM
953
- INTEGER ISW
954
- REAL (-JM:JM,-MM:MM) W
955
- REAL ((MM+1)*(MM+1)) o S
956
- REAL (0:MM+1) t SD
957
- REAL (0:MM+1,0:JM) t PM
958
- REAL (0:JM) t XS
959
- REAL (0:JM) t XC
960
- REAL (2,0:JM,0:MM) t PX
961
- REAL (0:JM,2) t WY
962
- REAL (0:JM,2) t WX
963
- REAL ((MM+1)*(MM+1)) t R
964
- REAL (JM,0:JM,4) t Z
965
- )
966
- */
967
-
968
- /*****************************************************
969
- * math2/shtlib/shmsri.f
970
- ****************************************************/
971
- /*
972
- * SUBROUTINE SHMSRI(MM,M,SR,SI)
973
- * REAL SR(M:MM),SI(M:MM)
974
- */
975
- /*
976
- SUBROUTINE SHMSRI(
977
- INTEGER MM
978
- INTEGER M
979
- REAL (M:MM) io SR
980
- REAL (M:MM) io SI
981
- )
982
- */
983
-
984
- /*****************************************************
985
- * math2/shtlib/shlfwu.f
986
- ****************************************************/
987
- /*
988
- * SUBROUTINE SHLFWU(MM,JM,M,ISW,WM,SM,SD,XS,XC,PM,WY,WX,PX,R,Z)
989
- * REAL SM(M:MM),WM(-JM:JM)
990
- * REAL SD(0:MM+1)
991
- * REAL XS(0:JM),XC(0:JM)
992
- * REAL PM(0:MM+1,0:JM)
993
- * REAL WY(0:JM,2),WX(0:JM,2)
994
- * REAL PX(2,0:JM,0:MM)
995
- * REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
996
- */
997
- /*
998
- SUBROUTINE SHLFWU(
999
- INTEGER MM
1000
- INTEGER JM
1001
- INTEGER M
1002
- INTEGER ISW
1003
- REAL (-JM:JM) WM
1004
- REAL (M:MM) o SM
1005
- REAL (0:MM+1) t SD
1006
- REAL (0:JM) t XS
1007
- REAL (0:JM) t XC
1008
- REAL (0:MM+1,0:JM) t PM
1009
- REAL (0:JM,2) t WY
1010
- REAL (0:JM,2) t WX
1011
- REAL (2,0:JM,0:MM) t PX
1012
- REAL ((MM+1)*(MM+1)) t R
1013
- REAL (JM,0:JM,4) t Z
1014
- )
1015
- */
1016
-
1017
- /*****************************************************
1018
- * math2/shtlib/shlfwm.f
1019
- ****************************************************/
1020
- /*
1021
- * SUBROUTINE SHLFWM(MM,JM,M,ISW,WM,SM,SD,XC,PM,WY,WX,R,Z)
1022
- * REAL SM(M:MM),WM(-JM:JM)
1023
- * REAL SD(0:MM+1)
1024
- * REAL XC(0:JM)
1025
- * REAL PM(0:MM+1,0:JM)
1026
- * REAL WY(0:JM,2),WX(0:JM,2)
1027
- * REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
1028
- */
1029
- /*
1030
- SUBROUTINE SHLFWM(
1031
- INTEGER MM
1032
- INTEGER JM
1033
- INTEGER M
1034
- INTEGER ISW
1035
- REAL (-JM:JM) WM
1036
- REAL (M:MM) SM
1037
- REAL (0:MM+1) o SD
1038
- REAL (0:JM) t XC
1039
- REAL (0:MM+1,0:JM) t PM
1040
- REAL (0:JM,2) t WY
1041
- REAL (0:JM,2) t WX
1042
- REAL ((MM+1)*(MM+1)) t R
1043
- REAL (JM,0:JM,4) t Z
1044
- )
1045
- */
1046
-
1047
- /*****************************************************
1048
- * math2/shtlib/shlfwl.f
1049
- ****************************************************/
1050
- /*
1051
- * SUBROUTINE SHLFWL(MM,JM,N,M,WX,SD,PM)
1052
- * REAL WX(0:JM,2)
1053
- * REAL SD,PM(0:MM+1,0:JM)
1054
- */
1055
- /*
1056
- SUBROUTINE SHLFWL(
1057
- INTEGER MM
1058
- INTEGER JM
1059
- INTEGER N
1060
- INTEGER M
1061
- REAL (0:JM,2) WX
1062
- REAL o SD
1063
- REAL (0:MM+1,0:JM) PM
1064
- )
1065
- */
1066
-
1067
- /*****************************************************
1068
- * math2/shtlib/shly2x.f
1069
- ****************************************************/
1070
- /*
1071
- * SUBROUTINE SHLY2X(JM,M,ISW,WM,WX,WY,Z)
1072
- * REAL WM(-JM:JM)
1073
- * REAL WY(0:JM,2),WX(0:JM,2)
1074
- * REAL Z(JM,0:JM,4)
1075
- */
1076
- /*
1077
- SUBROUTINE SHLY2X(
1078
- INTEGER JM
1079
- INTEGER M
1080
- INTEGER ISW
1081
- REAL (-JM:JM) WM
1082
- REAL (0:JM,2) o WX
1083
- REAL (0:JM,2) o WY
1084
- REAL (JM,0:JM,4) Z
1085
- )
1086
- */
1087
-
1088
- /*****************************************************
1089
- * math2/shtlib/shlbwu.f
1090
- ****************************************************/
1091
- /*
1092
- * SUBROUTINE SHLBWU(MM,JM,M,ISW,SM,WM,SD,PM,YS,YC,PY,R)
1093
- * REAL SM(M:MM)
1094
- * REAL WM(-JM:JM)
1095
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
1096
- * REAL YS(0:JM),YC(0:JM)
1097
- * REAL PY(2,0:JM,0:MM)
1098
- * REAL R((MM+1)*(MM+1))
1099
- */
1100
- /*
1101
- SUBROUTINE SHLBWU(
1102
- INTEGER MM
1103
- INTEGER JM
1104
- INTEGER M
1105
- INTEGER ISW
1106
- REAL (M:MM) SM
1107
- REAL (-JM:JM) o WM
1108
- REAL (0:MM+1) t SD
1109
- REAL (0:MM+1,0:JM) t PM
1110
- REAL (0:JM) t YS
1111
- REAL (0:JM) t YC
1112
- REAL (2,0:JM,0:MM) t PY
1113
- REAL ((MM+1)*(MM+1)) t R
1114
- )
1115
- */
1116
-
1117
- /*****************************************************
1118
- * math2/shtlib/shlbwm.f
1119
- ****************************************************/
1120
- /*
1121
- * SUBROUTINE SHLBWM(MM,JM,M,ISW,SM,WM,SD,PM,YC,R)
1122
- * REAL SM(M:MM)
1123
- * REAL WM(-JM:JM)
1124
- * REAL SD(0:MM+1),PM(0:MM+1,0:JM)
1125
- * REAL YC(0:JM)
1126
- * REAL R((MM+1)*(MM+1))
1127
- */
1128
- /*
1129
- SUBROUTINE SHLBWM(
1130
- INTEGER MM
1131
- INTEGER JM
1132
- INTEGER M
1133
- INTEGER ISW
1134
- REAL (M:MM) SM
1135
- REAL (-JM:JM) o WM
1136
- REAL (0:MM+1) t SD
1137
- REAL (0:MM+1,0:JM) t PM
1138
- REAL (0:JM) t YC
1139
- REAL ((MM+1)*(MM+1)) t R
1140
- )
1141
- */
1142
-
1143
- /*****************************************************
1144
- * math2/shtlib/shlbwj.f
1145
- ****************************************************/
1146
- /*
1147
- * SUBROUTINE SHLBWJ(MM,JM,M,J,ISW,SM,WMJ,SD,PMJ,YC,R)
1148
- * REAL SM(M:MM),WMJ
1149
- * REAL SD(0:MM+1),PMJ(0:MM+1)
1150
- * REAL YC(0:JM)
1151
- * REAL R((MM+1)*(MM+1))
1152
- */
1153
- /*
1154
- SUBROUTINE SHLBWJ(
1155
- INTEGER MM
1156
- INTEGER JM
1157
- INTEGER M
1158
- INTEGER J
1159
- INTEGER ISW
1160
- REAL (M:MM) SM
1161
- REAL o WMJ
1162
- REAL (0:MM+1) t SD
1163
- REAL (0:MM+1) t PMJ
1164
- REAL (0:JM) t YC
1165
- REAL ((MM+1)*(MM+1)) t R
1166
- )
1167
- */
1168
-
1169
- /*****************************************************
1170
- * math2/shtlib/shlbwl.f
1171
- ****************************************************/
1172
- /*
1173
- * SUBROUTINE SHLBWL(MM,M,SD,WS,WA,PMJ)
1174
- * REAL SD(0:MM+1)
1175
- * REAL WS,WA
1176
- * REAL PMJ(0:MM+1)
1177
- */
1178
- /*
1179
- SUBROUTINE SHLBWL(
1180
- INTEGER MM
1181
- INTEGER M
1182
- REAL (0:MM+1) SD
1183
- REAL o WS
1184
- REAL o WA
1185
- REAL (0:MM+1) PMJ
1186
- )
1187
- */
1188
-
1189
- /*****************************************************
1190
- * math2/shtlib/shlssd.f
1191
- ****************************************************/
1192
- /*
1193
- * SUBROUTINE SHLSSD(MM,M,ISW,SM,SD,R)
1194
- * REAL SM(M:MM)
1195
- * REAL SD(0:MM+1)
1196
- * REAL R((MM+1)*(MM+1))
1197
- */
1198
- /*
1199
- SUBROUTINE SHLSSD(
1200
- INTEGER MM
1201
- INTEGER M
1202
- INTEGER ISW
1203
- REAL (M:MM) SM
1204
- REAL (0:MM+1) o SD
1205
- REAL ((MM+1)*(MM+1)) R
1206
- )
1207
- */
1208
-
1209
- /*****************************************************
1210
- * math2/shtlib/shlsds.f
1211
- ****************************************************/
1212
- /*
1213
- * SUBROUTINE SHLSDS(MM,M,ISW,SD,SM,R)
1214
- * REAL SD(0:MM+1)
1215
- * REAL SM(M:MM)
1216
- * REAL R((MM+1)*(MM+1))
1217
- */
1218
- /*
1219
- SUBROUTINE SHLSDS(
1220
- INTEGER MM
1221
- INTEGER M
1222
- INTEGER ISW
1223
- REAL (0:MM+1) SD
1224
- REAL (M:MM) o SM
1225
- REAL ((MM+1)*(MM+1)) R
1226
- )
1227
- */
1228
-
1229
- /*****************************************************
1230
- * math2/shtlib/shpfun.f
1231
- ****************************************************/
1232
- /*
1233
- * SUBROUTINE SHPFUN(MM,JM,M,FUN,PM,YS,YC,PY,R)
1234
- * REAL FUN(-JM:JM,M:MM)
1235
- * REAL PM(0:MM+1,0:JM)
1236
- * REAL YS(0:JM),YC(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
1237
- */
1238
- /*
1239
- SUBROUTINE SHPFUN(
1240
- INTEGER MM
1241
- INTEGER JM
1242
- INTEGER M
1243
- REAL (-JM:JM,M:MM) o FUN
1244
- REAL (0:MM+1,0:JM) t PM
1245
- REAL (0:JM) t YS
1246
- REAL (0:JM) t YC
1247
- REAL (2,0:JM,0:MM) t PY
1248
- REAL ((MM+1)*(MM+1)) t R
1249
- )
1250
- */
1251
-
1252
- /*****************************************************
1253
- * math2/shtlib/shppma.f
1254
- ****************************************************/
1255
- /*
1256
- * SUBROUTINE SHPPMA(MM,JM,M,PM,YS,PY,R)
1257
- * REAL PM(0:MM+1,0:JM)
1258
- * REAL YS(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
1259
- */
1260
- /*
1261
- SUBROUTINE SHPPMA(
1262
- INTEGER MM
1263
- INTEGER JM
1264
- INTEGER M
1265
- REAL (0:MM+1,0:JM) io PM
1266
- REAL (0:JM) YS
1267
- REAL (2,0:JM,0:MM) PY
1268
- REAL ((MM+1)*(MM+1)) R
1269
- )
1270
- */
1271
-
1272
- /*****************************************************
1273
- * math2/shtlib/shppmj.f
1274
- ****************************************************/
1275
- /*
1276
- * SUBROUTINE SHPPMJ(MM,JM,M,J,PMJ,YS,PY,R)
1277
- * REAL PMJ(0:MM+1)
1278
- * REAL YS(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
1279
- */
1280
- /*
1281
- SUBROUTINE SHPPMJ(
1282
- INTEGER MM
1283
- INTEGER JM
1284
- INTEGER M
1285
- INTEGER J
1286
- REAL (0:MM+1) io PMJ
1287
- REAL (0:JM) YS
1288
- REAL (2,0:JM,0:MM) PY
1289
- REAL ((MM+1)*(MM+1)) R
1290
- )
1291
- */
1292
-