ruby-dcl 1.6.3.1 → 1.7.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- data/.gitignore +17 -0
- data/ChangeLog +57 -1
- data/Gemfile +4 -0
- data/GenWrapper/Makefile +1 -0
- data/GenWrapper/dcl_rb_heading +2 -8
- data/GenWrapper/elim_ary_size.rb +1 -1
- data/GenWrapper/proto/Makefile +17 -16
- data/GenWrapper/proto/grph1_sgpack.fp +367 -1
- data/GenWrapper/proto/grph1_stpack.fp +45 -1
- data/GenWrapper/proto/grph1_swpack.fp +70 -1
- data/GenWrapper/proto/grph2_uipack.fp +508 -0
- data/GenWrapper/proto/math1_clsplib.fp +74 -0
- data/GenWrapper/proto/misc1_randlib.fp +15 -1
- data/README.md +29 -0
- data/Rakefile +18 -38
- data/dcl_narrayed_funcs.c +23 -22
- data/dcl_rubydcloriginal.c +21 -21
- data/demo/grph2/uipack/uipk01.rb +38 -0
- data/demo/grph2/uipack/uipk02.rb +48 -0
- data/demo/grph2/umpack/test09_2.rb +89 -0
- data/demo/math1/clsplib/cslb01.rb +55 -0
- data/extconf.rb +16 -48
- data/grph1_sgpack.c +1166 -0
- data/grph1_stpack.c +141 -0
- data/grph1_swpack.c +182 -0
- data/grph2_uipack.c +1456 -0
- data/init.c.default +8 -0
- data/lib/.keepme +0 -0
- data/lib/dcl/dcl.rb +5578 -0
- data/lib/dcl/version.rb +3 -0
- data/lib/dcl.rb +2 -5084
- data/math1_clsplib.c +329 -0
- data/misc1_randlib.c +37 -0
- data/ruby-dcl.gemspec +30 -0
- metadata +343 -352
- data/GenWrapper/etc/Mk_proto +0 -375
- data/GenWrapper/etc/dcl_gen +0 -13
- data/GenWrapper/etc/p_init +0 -410
- data/GenWrapper/proto/math2_shtlib.fp_old +0 -1294
- data/README_old +0 -79
- data/dummy.c +0 -9
- data/grph1_zgpack.c.org +0 -141
- data/init.c.gtk +0 -189
@@ -1,1294 +0,0 @@
|
|
1
|
-
/*
|
2
|
-
* math2_shtlib.fp
|
3
|
-
*
|
4
|
-
* $Id: math2_shtlib.fp_old,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 o W */
|
370
|
-
/* REAL o S */
|
371
|
-
REAL ((2*JM+1)*(2*MM+1)) o W
|
372
|
-
REAL ((MM+1)*(MM+1)) o S
|
373
|
-
REAL ((JM+1)*(4*JM+5*MM+14)+(MM+1)*(MM+1)+MM+2+6*IM+15) WORK
|
374
|
-
)
|
375
|
-
|
376
|
-
/*****************************************************
|
377
|
-
* math2/shtlib/shinit.f
|
378
|
-
****************************************************/
|
379
|
-
/*
|
380
|
-
* SUBROUTINE SHINIT(MM,JM,IM,Y,YS,YC,X,XS,XC,XW,Z,PY,PX,R,WFFT)
|
381
|
-
* REAL Y(0:JM),YS(0:JM),YC(0:JM)
|
382
|
-
* REAL X(0:JM),XS(0:JM),XC(0:JM),XW(0:JM)
|
383
|
-
* REAL Z(JM,0:JM,4)
|
384
|
-
* REAL PY(2,0:JM,0:MM)
|
385
|
-
* REAL PX(2,0:JM,0:MM)
|
386
|
-
* REAL R((MM+1)*(MM+1))
|
387
|
-
* REAL WFFT(*)
|
388
|
-
*/
|
389
|
-
/*
|
390
|
-
SUBROUTINE SHINIT(
|
391
|
-
INTEGER MM
|
392
|
-
INTEGER JM
|
393
|
-
INTEGER IM
|
394
|
-
REAL o Y
|
395
|
-
REAL (0:JM) o YS
|
396
|
-
REAL (0:JM) o YC
|
397
|
-
REAL o X
|
398
|
-
REAL (0:JM) o XS
|
399
|
-
REAL (0:JM) o XC
|
400
|
-
REAL (0:JM) o XW
|
401
|
-
REAL (JM,0:JM,4) o Z
|
402
|
-
REAL (2,0:JM,0:MM) o PY
|
403
|
-
REAL (2,0:JM,0:MM) o PX
|
404
|
-
REAL ((MM+1)*(MM+1)) o R
|
405
|
-
REAL (*) o WFFT
|
406
|
-
)
|
407
|
-
*/
|
408
|
-
|
409
|
-
/*****************************************************
|
410
|
-
* math2/shtlib/shinix.f
|
411
|
-
****************************************************/
|
412
|
-
/*
|
413
|
-
* SUBROUTINE SHINIX(JM,X,XS,XC,XW)
|
414
|
-
* IMPLICIT REAL*8(A-H,O-Z)
|
415
|
-
* REAL X(0:JM),XS(0:JM),XC(0:JM),XW(0:JM)
|
416
|
-
*/
|
417
|
-
/*
|
418
|
-
SUBROUTINE SHINIX(
|
419
|
-
INTEGER JM
|
420
|
-
REAL o X
|
421
|
-
REAL (0:JM) o XS
|
422
|
-
REAL (0:JM) o XC
|
423
|
-
REAL (0:JM) o XW
|
424
|
-
)
|
425
|
-
*/
|
426
|
-
|
427
|
-
/*****************************************************
|
428
|
-
* math2/shtlib/shiniy.f
|
429
|
-
****************************************************/
|
430
|
-
/*
|
431
|
-
* SUBROUTINE SHINIY(JM,Y,YS,YC)
|
432
|
-
* REAL Y(0:JM),YS(0:JM),YC(0:JM)
|
433
|
-
*/
|
434
|
-
/*
|
435
|
-
SUBROUTINE SHINIY(
|
436
|
-
INTEGER JM
|
437
|
-
REAL o Y
|
438
|
-
REAL (0:JM) o YS
|
439
|
-
REAL (0:JM) o YC
|
440
|
-
)
|
441
|
-
*/
|
442
|
-
|
443
|
-
/*****************************************************
|
444
|
-
* math2/shtlib/shiniz.f
|
445
|
-
****************************************************/
|
446
|
-
/*
|
447
|
-
* SUBROUTINE SHINIZ(JM,X,Y,Z)
|
448
|
-
* REAL X(0:JM),Y(0:JM)
|
449
|
-
* REAL Z(JM,0:JM,4)
|
450
|
-
*/
|
451
|
-
/*
|
452
|
-
SUBROUTINE SHINIZ(
|
453
|
-
INTEGER JM
|
454
|
-
REAL (0:JM) X
|
455
|
-
REAL (0:JM) Y
|
456
|
-
REAL (JM,0:JM,4) o Z
|
457
|
-
)
|
458
|
-
*/
|
459
|
-
|
460
|
-
/*****************************************************
|
461
|
-
* math2/shtlib/shinip.f
|
462
|
-
****************************************************/
|
463
|
-
/*
|
464
|
-
* SUBROUTINE SHINIP(MM,JM,YS,YC,PY)
|
465
|
-
* REAL YS(0:JM),YC(0:JM)
|
466
|
-
* REAL PY(2,0:JM,0:MM)
|
467
|
-
*/
|
468
|
-
/*
|
469
|
-
SUBROUTINE SHINIP(
|
470
|
-
INTEGER MM
|
471
|
-
INTEGER JM
|
472
|
-
REAL (0:JM) YS
|
473
|
-
REAL (0:JM) YC
|
474
|
-
REAL (2,0:JM,0:MM) o PY
|
475
|
-
)
|
476
|
-
*/
|
477
|
-
|
478
|
-
/*****************************************************
|
479
|
-
* math2/shtlib/shinic.f
|
480
|
-
****************************************************/
|
481
|
-
/*
|
482
|
-
* SUBROUTINE SHINIC(MM,JM,XW,PX)
|
483
|
-
* REAL XW(0:JM),PX(2,0:JM,0:MM)
|
484
|
-
*/
|
485
|
-
/*
|
486
|
-
SUBROUTINE SHINIC(
|
487
|
-
INTEGER MM
|
488
|
-
INTEGER JM
|
489
|
-
REAL (0:JM) XW
|
490
|
-
REAL (2,0:JM,0:MM) o PX
|
491
|
-
)
|
492
|
-
*/
|
493
|
-
|
494
|
-
/*****************************************************
|
495
|
-
* math2/shtlib/shinir.f
|
496
|
-
****************************************************/
|
497
|
-
/*
|
498
|
-
* SUBROUTINE SHINIR(MM,R)
|
499
|
-
* REAL R((MM+1)*(MM+1))
|
500
|
-
*/
|
501
|
-
/*
|
502
|
-
SUBROUTINE SHINIR(
|
503
|
-
INTEGER MM
|
504
|
-
REAL ((MM+1)*(MM+1)) o R
|
505
|
-
)
|
506
|
-
*/
|
507
|
-
|
508
|
-
/*****************************************************
|
509
|
-
* math2/shtlib/shinif.f
|
510
|
-
****************************************************/
|
511
|
-
/*
|
512
|
-
* SUBROUTINE SHINIF(NFFT,WFFT)
|
513
|
-
* REAL WFFT(*)
|
514
|
-
*/
|
515
|
-
/*
|
516
|
-
SUBROUTINE SHINIF(
|
517
|
-
INTEGER NFFT
|
518
|
-
REAL (*) o WFFT
|
519
|
-
)
|
520
|
-
*/
|
521
|
-
|
522
|
-
/*****************************************************
|
523
|
-
* math2/shtlib/shnm2l.f
|
524
|
-
****************************************************/
|
525
|
-
/*
|
526
|
-
* SUBROUTINE SHNM2L(MM,N,M,LR,LI)
|
527
|
-
*/
|
528
|
-
/*
|
529
|
-
SUBROUTINE SHNM2L(
|
530
|
-
INTEGER MM
|
531
|
-
INTEGER N
|
532
|
-
INTEGER M
|
533
|
-
INTEGER o LR
|
534
|
-
INTEGER o LI
|
535
|
-
)
|
536
|
-
*/
|
537
|
-
|
538
|
-
/*****************************************************
|
539
|
-
* math2/shtlib/sholap.f
|
540
|
-
****************************************************/
|
541
|
-
/*
|
542
|
-
* SUBROUTINE SHOLAP(MM,IND,A,B)
|
543
|
-
* REAL A((MM+1)*(MM+1)),B((MM+1)*(MM+1))
|
544
|
-
*/
|
545
|
-
/*
|
546
|
-
SUBROUTINE SHOLAP(
|
547
|
-
INTEGER MM
|
548
|
-
INTEGER IND
|
549
|
-
REAL ((MM+1)*(MM+1)) A
|
550
|
-
REAL ((MM+1)*(MM+1)) o B
|
551
|
-
)
|
552
|
-
*/
|
553
|
-
|
554
|
-
/*****************************************************
|
555
|
-
* math2/shtlib/shfw2g.f
|
556
|
-
****************************************************/
|
557
|
-
/*
|
558
|
-
* SUBROUTINE SHFW2G(MM,JM,IM,W,G,H,WFFT)
|
559
|
-
* REAL W(-JM:JM,-MM:MM),G(-IM:IM,-JM:JM)
|
560
|
-
* REAL H(0:2*IM-1)
|
561
|
-
* REAL WFFT(*)
|
562
|
-
*/
|
563
|
-
/*
|
564
|
-
SUBROUTINE SHFW2G(
|
565
|
-
INTEGER MM
|
566
|
-
INTEGER JM
|
567
|
-
INTEGER IM
|
568
|
-
REAL (-JM:JM,-MM:MM) W
|
569
|
-
REAL (-IM:IM,-JM:JM) o G
|
570
|
-
REAL (0:2*IM-1) t H
|
571
|
-
REAL (*) t WFFT
|
572
|
-
)
|
573
|
-
*/
|
574
|
-
|
575
|
-
/*****************************************************
|
576
|
-
* math2/shtlib/shfwga.f
|
577
|
-
****************************************************/
|
578
|
-
/*
|
579
|
-
* SUBROUTINE SHFWGA(MM,JM,IM,M1,M2,W,G,H,WFFT)
|
580
|
-
* REAL W(-JM:JM,-MM:MM),G(-IM:IM,-JM:JM)
|
581
|
-
* REAL H(0:2*IM-1)
|
582
|
-
* REAL WFFT(*)
|
583
|
-
*/
|
584
|
-
/*
|
585
|
-
SUBROUTINE SHFWGA(
|
586
|
-
INTEGER MM
|
587
|
-
INTEGER JM
|
588
|
-
INTEGER IM
|
589
|
-
INTEGER M1
|
590
|
-
INTEGER M2
|
591
|
-
REAL (-JM:JM,-MM:MM) W
|
592
|
-
REAL (-IM:IM,-JM:JM) o G
|
593
|
-
REAL (0:2*IM-1) t H
|
594
|
-
REAL (*) t WFFT
|
595
|
-
)
|
596
|
-
*/
|
597
|
-
|
598
|
-
/*****************************************************
|
599
|
-
* math2/shtlib/shfwgm.f
|
600
|
-
****************************************************/
|
601
|
-
/*
|
602
|
-
* SUBROUTINE SHFWGM(JM,IM,M,WR,WI,G,H,WFFT)
|
603
|
-
* REAL WR(-JM:JM),WI(-JM:JM),G(-IM:IM,-JM:JM)
|
604
|
-
* REAL H(0:2*IM-1)
|
605
|
-
* REAL WFFT(*)
|
606
|
-
*/
|
607
|
-
/* SUBROUTINE SHFWGM(
|
608
|
-
INTEGER JM
|
609
|
-
INTEGER IM
|
610
|
-
INTEGER M
|
611
|
-
REAL (-JM:JM) WR
|
612
|
-
REAL (-JM:JM) WI
|
613
|
-
REAL (-IM:IM,-JM:JM) o G
|
614
|
-
REAL (0:2*IM-1) t H
|
615
|
-
REAL (*) t WFFT
|
616
|
-
)
|
617
|
-
*/
|
618
|
-
|
619
|
-
/*****************************************************
|
620
|
-
* math2/shtlib/shfwgz.f
|
621
|
-
****************************************************/
|
622
|
-
/*
|
623
|
-
* SUBROUTINE SHFWGZ(JM,IM,WZ,G)
|
624
|
-
* REAL WZ(-JM:JM),G(-IM:IM,-JM:JM)
|
625
|
-
*/
|
626
|
-
/*
|
627
|
-
SUBROUTINE SHFWGZ(
|
628
|
-
INTEGER JM
|
629
|
-
INTEGER IM
|
630
|
-
REAL (-JM:JM) WZ
|
631
|
-
REAL (-IM:IM,-JM:JM) o G
|
632
|
-
)
|
633
|
-
*/
|
634
|
-
|
635
|
-
/*****************************************************
|
636
|
-
* math2/shtlib/shfwgj.f
|
637
|
-
****************************************************/
|
638
|
-
/*
|
639
|
-
* SUBROUTINE SHFWGJ(MM,IM,M1,M2,WJ,GJ,H,WFFT)
|
640
|
-
* REAL WJ(-MM:MM),GJ(-IM:IM)
|
641
|
-
* REAL H(0:2*IM-1)
|
642
|
-
* REAL WFFT(*)
|
643
|
-
*/
|
644
|
-
/*
|
645
|
-
SUBROUTINE SHFWGJ(
|
646
|
-
INTEGER MM
|
647
|
-
INTEGER IM
|
648
|
-
INTEGER M1
|
649
|
-
INTEGER M2
|
650
|
-
REAL (-MM:MM) WJ
|
651
|
-
REAL (-IM:IM) o GJ
|
652
|
-
REAL (0:2*IM-1) t H
|
653
|
-
REAL (*) t WFFT
|
654
|
-
)
|
655
|
-
*/
|
656
|
-
|
657
|
-
/*****************************************************
|
658
|
-
* math2/shtlib/shfg2w.f
|
659
|
-
****************************************************/
|
660
|
-
/*
|
661
|
-
* SUBROUTINE SHFG2W(MM,JM,IM,G,W,H,WFFT)
|
662
|
-
* REAL G(-IM:IM,-JM:JM),W(-JM:JM,-MM:MM)
|
663
|
-
* REAL H(0:2*IM-1)
|
664
|
-
* REAL WFFT(*)
|
665
|
-
*/
|
666
|
-
/*
|
667
|
-
SUBROUTINE SHFG2W(
|
668
|
-
INTEGER MM
|
669
|
-
INTEGER JM
|
670
|
-
INTEGER IM
|
671
|
-
REAL (-IM:IM,-JM:JM) G
|
672
|
-
REAL (-JM:JM,-MM:MM) o W
|
673
|
-
REAL (0:2*IM-1) t H
|
674
|
-
REAL (*) t WFFT
|
675
|
-
)
|
676
|
-
*/
|
677
|
-
|
678
|
-
/*****************************************************
|
679
|
-
* math2/shtlib/shfftf.f
|
680
|
-
****************************************************/
|
681
|
-
/*
|
682
|
-
* SUBROUTINE SHFFTF(NFFT,XFFT,WFFT)
|
683
|
-
* REAL XFFT(NFFT)
|
684
|
-
* REAL WFFT(*)
|
685
|
-
*/
|
686
|
-
/*
|
687
|
-
SUBROUTINE SHFFTF(
|
688
|
-
INTEGER NFFT
|
689
|
-
REAL (NFFT) io XFFT
|
690
|
-
REAL (*) t WFFT
|
691
|
-
)
|
692
|
-
*/
|
693
|
-
|
694
|
-
/*****************************************************
|
695
|
-
* math2/shtlib/shfftb.f
|
696
|
-
****************************************************/
|
697
|
-
/*
|
698
|
-
* SUBROUTINE SHFFTB(NFFT,XFFT,WFFT)
|
699
|
-
* REAL XFFT(NFFT)
|
700
|
-
* REAL WFFT(*)
|
701
|
-
*/
|
702
|
-
/*
|
703
|
-
SUBROUTINE SHFFTB(
|
704
|
-
INTEGER NFFT
|
705
|
-
REAL (NFFT) io XFFT
|
706
|
-
REAL (*) t WFFT
|
707
|
-
)
|
708
|
-
*/
|
709
|
-
|
710
|
-
/*****************************************************
|
711
|
-
* math2/shtlib/shms2w.f
|
712
|
-
****************************************************/
|
713
|
-
/*
|
714
|
-
* SUBROUTINE SHMS2W(MM,JM,ISW,S,W,SD,PM,YS,YC,PY,R)
|
715
|
-
* REAL S((MM+1)*(MM+1)),W(-JM:JM,-MM:MM)
|
716
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
717
|
-
* REAL YS(0:JM),YC(0:JM)
|
718
|
-
* REAL PY(2,0:JM,0:MM)
|
719
|
-
* REAL R((MM+1)*(MM+1))
|
720
|
-
*/
|
721
|
-
/*
|
722
|
-
SUBROUTINE SHMS2W(
|
723
|
-
INTEGER MM
|
724
|
-
INTEGER JM
|
725
|
-
INTEGER ISW
|
726
|
-
REAL ((MM+1)*(MM+1)) S
|
727
|
-
REAL (-JM:JM,-MM:MM) o W
|
728
|
-
REAL (0:MM+1) t SD
|
729
|
-
REAL (0:MM+1,0:JM) t PM
|
730
|
-
REAL (0:JM) t YS
|
731
|
-
REAL (0:JM) t YC
|
732
|
-
REAL (2,0:JM,0:MM) t PY
|
733
|
-
REAL ((MM+1)*(MM+1)) t R
|
734
|
-
)
|
735
|
-
*/
|
736
|
-
|
737
|
-
/*****************************************************
|
738
|
-
* math2/shtlib/shmswa.f
|
739
|
-
****************************************************/
|
740
|
-
/*
|
741
|
-
* SUBROUTINE SHMSWA(MM,JM,ISW,M1,M2,S,W,SD,PM,YS,YC,PY,R)
|
742
|
-
* REAL S((MM+1)*(MM+1)),W(-JM:JM,-MM:MM)
|
743
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
744
|
-
* REAL YS(0:JM),YC(0:JM)
|
745
|
-
* REAL PY(2,0:JM,0:MM)
|
746
|
-
* REAL R((MM+1)*(MM+1))
|
747
|
-
*/
|
748
|
-
/*
|
749
|
-
SUBROUTINE SHMSWA(
|
750
|
-
INTEGER MM
|
751
|
-
INTEGER JM
|
752
|
-
INTEGER ISW
|
753
|
-
INTEGER M1
|
754
|
-
INTEGER M2
|
755
|
-
REAL ((MM+1)*(MM+1)) S
|
756
|
-
REAL (-JM:JM,-MM:MM) o W
|
757
|
-
REAL (0:MM+1) t SD
|
758
|
-
REAL (0:MM+1,0:JM) t PM
|
759
|
-
REAL (0:JM) t YS
|
760
|
-
REAL (0:JM) t YC
|
761
|
-
REAL (2,0:JM,0:MM) t PY
|
762
|
-
REAL ((MM+1)*(MM+1)) t R
|
763
|
-
)
|
764
|
-
*/
|
765
|
-
|
766
|
-
/*****************************************************
|
767
|
-
* math2/shtlib/shmswz.f
|
768
|
-
****************************************************/
|
769
|
-
/*
|
770
|
-
* SUBROUTINE SHMSWZ(MM,JM,ISW,S,WZ,SD,PM,YS,YC,PY,R)
|
771
|
-
* REAL S((MM+1)*(MM+1)),WZ(-JM:JM)
|
772
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
773
|
-
* REAL YS(0:JM),YC(0:JM)
|
774
|
-
* REAL PY(2,0:JM,0:MM)
|
775
|
-
* REAL R((MM+1)*(MM+1))
|
776
|
-
*/
|
777
|
-
/*
|
778
|
-
SUBROUTINE SHMSWZ(
|
779
|
-
INTEGER MM
|
780
|
-
INTEGER JM
|
781
|
-
INTEGER ISW
|
782
|
-
REAL ((MM+1)*(MM+1)) S
|
783
|
-
REAL (-JM:JM) o WZ
|
784
|
-
REAL (0:MM+1) t SD
|
785
|
-
REAL (0:MM+1,0:JM) t PM
|
786
|
-
REAL (0:JM) t YS
|
787
|
-
REAL (0:JM) t YC
|
788
|
-
REAL (2,0:JM,0:MM) t PY
|
789
|
-
REAL ((MM+1)*(MM+1)) t R
|
790
|
-
)
|
791
|
-
*/
|
792
|
-
|
793
|
-
/*****************************************************
|
794
|
-
* math2/shtlib/shmswm.f
|
795
|
-
****************************************************/
|
796
|
-
/*
|
797
|
-
* SUBROUTINE SHMSWM(MM,JM,M,ISW,S,WR,WI,SD,PM,YS,YC,PY,R)
|
798
|
-
* REAL S((MM+1)*(MM+1)),WR(-JM:JM),WI(-JM:JM)
|
799
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
800
|
-
* REAL YS(0:JM),YC(0:JM)
|
801
|
-
* REAL PY(2,0:JM,0:MM)
|
802
|
-
* REAL R((MM+1)*(MM+1))
|
803
|
-
*/
|
804
|
-
/*
|
805
|
-
SUBROUTINE SHMSWM(
|
806
|
-
INTEGER MM
|
807
|
-
INTEGER JM
|
808
|
-
INTEGER M
|
809
|
-
INTEGER ISW
|
810
|
-
REAL ((MM+1)*(MM+1)) S
|
811
|
-
REAL (-JM:JM) o WR
|
812
|
-
REAL (-JM:JM) o WI
|
813
|
-
REAL (0:MM+1) t SD
|
814
|
-
REAL (0:MM+1,0:JM) t PM
|
815
|
-
REAL (0:JM) t YS
|
816
|
-
REAL (0:JM) t YC
|
817
|
-
REAL (2,0:JM,0:MM) t PY
|
818
|
-
REAL ((MM+1)*(MM+1)) t R
|
819
|
-
)
|
820
|
-
*/
|
821
|
-
|
822
|
-
/*****************************************************
|
823
|
-
* math2/shtlib/shmdxm.f
|
824
|
-
****************************************************/
|
825
|
-
/*
|
826
|
-
* SUBROUTINE SHMDXM(JM,WR,WI)
|
827
|
-
* REAL WR(-JM:JM),WI(-JM:JM)
|
828
|
-
*/
|
829
|
-
/*
|
830
|
-
SUBROUTINE SHMDXM(
|
831
|
-
INTEGER JM
|
832
|
-
REAL (-JM:JM) io WR
|
833
|
-
REAL (-JM:JM) io WI
|
834
|
-
)
|
835
|
-
*/
|
836
|
-
|
837
|
-
/*****************************************************
|
838
|
-
* math2/shtlib/shmswj.f
|
839
|
-
****************************************************/
|
840
|
-
/*
|
841
|
-
* SUBROUTINE SHMSWJ(MM,JM,ISW,J,M1,M2,S,WJ,SD,PMJ,YS,YC,PY,R)
|
842
|
-
* REAL S((MM+1)*(MM+1)),WJ(-MM:MM)
|
843
|
-
* REAL SD(0:MM+1),PMJ(0:MM+1)
|
844
|
-
* REAL YS(0:JM),YC(0:JM)
|
845
|
-
* REAL PY(2,0:JM,0:MM)
|
846
|
-
* REAL R((MM+1)*(MM+1))
|
847
|
-
*/
|
848
|
-
/*
|
849
|
-
SUBROUTINE SHMSWJ(
|
850
|
-
INTEGER MM
|
851
|
-
INTEGER JM
|
852
|
-
INTEGER ISW
|
853
|
-
INTEGER J
|
854
|
-
INTEGER M1
|
855
|
-
INTEGER M2
|
856
|
-
REAL ((MM+1)*(MM+1)) S
|
857
|
-
REAL (-MM:MM) o WJ
|
858
|
-
REAL (0:MM+1) t SD
|
859
|
-
REAL (0:MM+1) t PMJ
|
860
|
-
REAL (0:JM) t YS
|
861
|
-
REAL (0:JM) t YC
|
862
|
-
REAL (2,0:JM,0:MM) t PY
|
863
|
-
REAL ((MM+1)*(MM+1)) t R
|
864
|
-
)
|
865
|
-
*/
|
866
|
-
|
867
|
-
/*****************************************************
|
868
|
-
* math2/shtlib/shmwjz.f
|
869
|
-
****************************************************/
|
870
|
-
/*
|
871
|
-
* SUBROUTINE SHMWJZ(MM,JM,ISW,J,S,WJD,SD,PMJ,YS,YC,PY,R)
|
872
|
-
* REAL S((MM+1)*(MM+1)),WJD
|
873
|
-
* REAL SD(0:MM+1),PMJ(0:MM+1)
|
874
|
-
* REAL YS(0:JM),YC(0:JM)
|
875
|
-
* REAL PY(2,0:JM,0:MM)
|
876
|
-
* REAL R((MM+1)*(MM+1))
|
877
|
-
*/
|
878
|
-
/*
|
879
|
-
SUBROUTINE SHMWJZ(
|
880
|
-
INTEGER MM
|
881
|
-
INTEGER JM
|
882
|
-
INTEGER ISW
|
883
|
-
INTEGER J
|
884
|
-
REAL ((MM+1)*(MM+1)) S
|
885
|
-
REAL o WJD
|
886
|
-
REAL (0:MM+1) t SD
|
887
|
-
REAL (0:MM+1) t PMJ
|
888
|
-
REAL (0:JM) t YS
|
889
|
-
REAL (0:JM) t YC
|
890
|
-
REAL (2,0:JM,0:MM) t PY
|
891
|
-
REAL ((MM+1)*(MM+1)) t R
|
892
|
-
)
|
893
|
-
*/
|
894
|
-
|
895
|
-
/*****************************************************
|
896
|
-
* math2/shtlib/shmwjm.f
|
897
|
-
****************************************************/
|
898
|
-
/*
|
899
|
-
* SUBROUTINE SHMWJM(MM,JM,ISW,J,M,S,WJR,WJI,SD,PMJ,YS,YC,PY,R)
|
900
|
-
* REAL S((MM+1)*(MM+1)),WJR,WJI
|
901
|
-
* REAL SD(0:MM+1),PMJ(0:MM+1)
|
902
|
-
* REAL YS(0:JM),YC(0:JM)
|
903
|
-
* REAL PY(2,0:JM,0:MM)
|
904
|
-
* REAL R((MM+1)*(MM+1))
|
905
|
-
*/
|
906
|
-
/*
|
907
|
-
SUBROUTINE SHMWJM(
|
908
|
-
INTEGER MM
|
909
|
-
INTEGER JM
|
910
|
-
INTEGER ISW
|
911
|
-
INTEGER J
|
912
|
-
INTEGER M
|
913
|
-
REAL ((MM+1)*(MM+1)) S
|
914
|
-
REAL o WJR
|
915
|
-
REAL o WJI
|
916
|
-
REAL (0:MM+1) t SD
|
917
|
-
REAL (0:MM+1) t PMJ
|
918
|
-
REAL (0:JM) t YS
|
919
|
-
REAL (0:JM) t YC
|
920
|
-
REAL (2,0:JM,0:MM) t PY
|
921
|
-
REAL ((MM+1)*(MM+1)) t R
|
922
|
-
)
|
923
|
-
*/
|
924
|
-
|
925
|
-
/*****************************************************
|
926
|
-
* math2/shtlib/shmdxj.f
|
927
|
-
****************************************************/
|
928
|
-
/*
|
929
|
-
* SUBROUTINE SHMDXJ(WRD,WID)
|
930
|
-
* REAL WRD,WID
|
931
|
-
*/
|
932
|
-
/*
|
933
|
-
SUBROUTINE SHMDXJ(
|
934
|
-
REAL io WRD
|
935
|
-
REAL io WID
|
936
|
-
)
|
937
|
-
*/
|
938
|
-
|
939
|
-
/*****************************************************
|
940
|
-
* math2/shtlib/shmw2s.f
|
941
|
-
****************************************************/
|
942
|
-
/*
|
943
|
-
* SUBROUTINE SHMW2S(MM,JM,ISW,W,S,SD,PM,XS,XC,PX,WY,WX,R,Z)
|
944
|
-
* REAL W(-JM:JM,-MM:MM),S((MM+1)*(MM+1))
|
945
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
946
|
-
* REAL XS(0:JM),XC(0:JM)
|
947
|
-
* REAL PX(2,0:JM,0:MM)
|
948
|
-
* REAL WY(0:JM,2),WX(0:JM,2)
|
949
|
-
* REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
|
950
|
-
*/
|
951
|
-
/*
|
952
|
-
SUBROUTINE SHMW2S(
|
953
|
-
INTEGER MM
|
954
|
-
INTEGER JM
|
955
|
-
INTEGER ISW
|
956
|
-
REAL (-JM:JM,-MM:MM) W
|
957
|
-
REAL ((MM+1)*(MM+1)) o S
|
958
|
-
REAL (0:MM+1) t SD
|
959
|
-
REAL (0:MM+1,0:JM) t PM
|
960
|
-
REAL (0:JM) t XS
|
961
|
-
REAL (0:JM) t XC
|
962
|
-
REAL (2,0:JM,0:MM) t PX
|
963
|
-
REAL (0:JM,2) t WY
|
964
|
-
REAL (0:JM,2) t WX
|
965
|
-
REAL ((MM+1)*(MM+1)) t R
|
966
|
-
REAL (JM,0:JM,4) t Z
|
967
|
-
)
|
968
|
-
*/
|
969
|
-
|
970
|
-
/*****************************************************
|
971
|
-
* math2/shtlib/shmsri.f
|
972
|
-
****************************************************/
|
973
|
-
/*
|
974
|
-
* SUBROUTINE SHMSRI(MM,M,SR,SI)
|
975
|
-
* REAL SR(M:MM),SI(M:MM)
|
976
|
-
*/
|
977
|
-
/*
|
978
|
-
SUBROUTINE SHMSRI(
|
979
|
-
INTEGER MM
|
980
|
-
INTEGER M
|
981
|
-
REAL (M:MM) io SR
|
982
|
-
REAL (M:MM) io SI
|
983
|
-
)
|
984
|
-
*/
|
985
|
-
|
986
|
-
/*****************************************************
|
987
|
-
* math2/shtlib/shlfwu.f
|
988
|
-
****************************************************/
|
989
|
-
/*
|
990
|
-
* SUBROUTINE SHLFWU(MM,JM,M,ISW,WM,SM,SD,XS,XC,PM,WY,WX,PX,R,Z)
|
991
|
-
* REAL SM(M:MM),WM(-JM:JM)
|
992
|
-
* REAL SD(0:MM+1)
|
993
|
-
* REAL XS(0:JM),XC(0:JM)
|
994
|
-
* REAL PM(0:MM+1,0:JM)
|
995
|
-
* REAL WY(0:JM,2),WX(0:JM,2)
|
996
|
-
* REAL PX(2,0:JM,0:MM)
|
997
|
-
* REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
|
998
|
-
*/
|
999
|
-
/*
|
1000
|
-
SUBROUTINE SHLFWU(
|
1001
|
-
INTEGER MM
|
1002
|
-
INTEGER JM
|
1003
|
-
INTEGER M
|
1004
|
-
INTEGER ISW
|
1005
|
-
REAL (-JM:JM) WM
|
1006
|
-
REAL (M:MM) o SM
|
1007
|
-
REAL (0:MM+1) t SD
|
1008
|
-
REAL (0:JM) t XS
|
1009
|
-
REAL (0:JM) t XC
|
1010
|
-
REAL (0:MM+1,0:JM) t PM
|
1011
|
-
REAL (0:JM,2) t WY
|
1012
|
-
REAL (0:JM,2) t WX
|
1013
|
-
REAL (2,0:JM,0:MM) t PX
|
1014
|
-
REAL ((MM+1)*(MM+1)) t R
|
1015
|
-
REAL (JM,0:JM,4) t Z
|
1016
|
-
)
|
1017
|
-
*/
|
1018
|
-
|
1019
|
-
/*****************************************************
|
1020
|
-
* math2/shtlib/shlfwm.f
|
1021
|
-
****************************************************/
|
1022
|
-
/*
|
1023
|
-
* SUBROUTINE SHLFWM(MM,JM,M,ISW,WM,SM,SD,XC,PM,WY,WX,R,Z)
|
1024
|
-
* REAL SM(M:MM),WM(-JM:JM)
|
1025
|
-
* REAL SD(0:MM+1)
|
1026
|
-
* REAL XC(0:JM)
|
1027
|
-
* REAL PM(0:MM+1,0:JM)
|
1028
|
-
* REAL WY(0:JM,2),WX(0:JM,2)
|
1029
|
-
* REAL R((MM+1)*(MM+1)),Z(JM,0:JM,4)
|
1030
|
-
*/
|
1031
|
-
/*
|
1032
|
-
SUBROUTINE SHLFWM(
|
1033
|
-
INTEGER MM
|
1034
|
-
INTEGER JM
|
1035
|
-
INTEGER M
|
1036
|
-
INTEGER ISW
|
1037
|
-
REAL (-JM:JM) WM
|
1038
|
-
REAL (M:MM) SM
|
1039
|
-
REAL (0:MM+1) o SD
|
1040
|
-
REAL (0:JM) t XC
|
1041
|
-
REAL (0:MM+1,0:JM) t PM
|
1042
|
-
REAL (0:JM,2) t WY
|
1043
|
-
REAL (0:JM,2) t WX
|
1044
|
-
REAL ((MM+1)*(MM+1)) t R
|
1045
|
-
REAL (JM,0:JM,4) t Z
|
1046
|
-
)
|
1047
|
-
*/
|
1048
|
-
|
1049
|
-
/*****************************************************
|
1050
|
-
* math2/shtlib/shlfwl.f
|
1051
|
-
****************************************************/
|
1052
|
-
/*
|
1053
|
-
* SUBROUTINE SHLFWL(MM,JM,N,M,WX,SD,PM)
|
1054
|
-
* REAL WX(0:JM,2)
|
1055
|
-
* REAL SD,PM(0:MM+1,0:JM)
|
1056
|
-
*/
|
1057
|
-
/*
|
1058
|
-
SUBROUTINE SHLFWL(
|
1059
|
-
INTEGER MM
|
1060
|
-
INTEGER JM
|
1061
|
-
INTEGER N
|
1062
|
-
INTEGER M
|
1063
|
-
REAL (0:JM,2) WX
|
1064
|
-
REAL o SD
|
1065
|
-
REAL (0:MM+1,0:JM) PM
|
1066
|
-
)
|
1067
|
-
*/
|
1068
|
-
|
1069
|
-
/*****************************************************
|
1070
|
-
* math2/shtlib/shly2x.f
|
1071
|
-
****************************************************/
|
1072
|
-
/*
|
1073
|
-
* SUBROUTINE SHLY2X(JM,M,ISW,WM,WX,WY,Z)
|
1074
|
-
* REAL WM(-JM:JM)
|
1075
|
-
* REAL WY(0:JM,2),WX(0:JM,2)
|
1076
|
-
* REAL Z(JM,0:JM,4)
|
1077
|
-
*/
|
1078
|
-
/*
|
1079
|
-
SUBROUTINE SHLY2X(
|
1080
|
-
INTEGER JM
|
1081
|
-
INTEGER M
|
1082
|
-
INTEGER ISW
|
1083
|
-
REAL (-JM:JM) WM
|
1084
|
-
REAL (0:JM,2) o WX
|
1085
|
-
REAL (0:JM,2) o WY
|
1086
|
-
REAL (JM,0:JM,4) Z
|
1087
|
-
)
|
1088
|
-
*/
|
1089
|
-
|
1090
|
-
/*****************************************************
|
1091
|
-
* math2/shtlib/shlbwu.f
|
1092
|
-
****************************************************/
|
1093
|
-
/*
|
1094
|
-
* SUBROUTINE SHLBWU(MM,JM,M,ISW,SM,WM,SD,PM,YS,YC,PY,R)
|
1095
|
-
* REAL SM(M:MM)
|
1096
|
-
* REAL WM(-JM:JM)
|
1097
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
1098
|
-
* REAL YS(0:JM),YC(0:JM)
|
1099
|
-
* REAL PY(2,0:JM,0:MM)
|
1100
|
-
* REAL R((MM+1)*(MM+1))
|
1101
|
-
*/
|
1102
|
-
/*
|
1103
|
-
SUBROUTINE SHLBWU(
|
1104
|
-
INTEGER MM
|
1105
|
-
INTEGER JM
|
1106
|
-
INTEGER M
|
1107
|
-
INTEGER ISW
|
1108
|
-
REAL (M:MM) SM
|
1109
|
-
REAL (-JM:JM) o WM
|
1110
|
-
REAL (0:MM+1) t SD
|
1111
|
-
REAL (0:MM+1,0:JM) t PM
|
1112
|
-
REAL (0:JM) t YS
|
1113
|
-
REAL (0:JM) t YC
|
1114
|
-
REAL (2,0:JM,0:MM) t PY
|
1115
|
-
REAL ((MM+1)*(MM+1)) t R
|
1116
|
-
)
|
1117
|
-
*/
|
1118
|
-
|
1119
|
-
/*****************************************************
|
1120
|
-
* math2/shtlib/shlbwm.f
|
1121
|
-
****************************************************/
|
1122
|
-
/*
|
1123
|
-
* SUBROUTINE SHLBWM(MM,JM,M,ISW,SM,WM,SD,PM,YC,R)
|
1124
|
-
* REAL SM(M:MM)
|
1125
|
-
* REAL WM(-JM:JM)
|
1126
|
-
* REAL SD(0:MM+1),PM(0:MM+1,0:JM)
|
1127
|
-
* REAL YC(0:JM)
|
1128
|
-
* REAL R((MM+1)*(MM+1))
|
1129
|
-
*/
|
1130
|
-
/*
|
1131
|
-
SUBROUTINE SHLBWM(
|
1132
|
-
INTEGER MM
|
1133
|
-
INTEGER JM
|
1134
|
-
INTEGER M
|
1135
|
-
INTEGER ISW
|
1136
|
-
REAL (M:MM) SM
|
1137
|
-
REAL (-JM:JM) o WM
|
1138
|
-
REAL (0:MM+1) t SD
|
1139
|
-
REAL (0:MM+1,0:JM) t PM
|
1140
|
-
REAL (0:JM) t YC
|
1141
|
-
REAL ((MM+1)*(MM+1)) t R
|
1142
|
-
)
|
1143
|
-
*/
|
1144
|
-
|
1145
|
-
/*****************************************************
|
1146
|
-
* math2/shtlib/shlbwj.f
|
1147
|
-
****************************************************/
|
1148
|
-
/*
|
1149
|
-
* SUBROUTINE SHLBWJ(MM,JM,M,J,ISW,SM,WMJ,SD,PMJ,YC,R)
|
1150
|
-
* REAL SM(M:MM),WMJ
|
1151
|
-
* REAL SD(0:MM+1),PMJ(0:MM+1)
|
1152
|
-
* REAL YC(0:JM)
|
1153
|
-
* REAL R((MM+1)*(MM+1))
|
1154
|
-
*/
|
1155
|
-
/*
|
1156
|
-
SUBROUTINE SHLBWJ(
|
1157
|
-
INTEGER MM
|
1158
|
-
INTEGER JM
|
1159
|
-
INTEGER M
|
1160
|
-
INTEGER J
|
1161
|
-
INTEGER ISW
|
1162
|
-
REAL (M:MM) SM
|
1163
|
-
REAL o WMJ
|
1164
|
-
REAL (0:MM+1) t SD
|
1165
|
-
REAL (0:MM+1) t PMJ
|
1166
|
-
REAL (0:JM) t YC
|
1167
|
-
REAL ((MM+1)*(MM+1)) t R
|
1168
|
-
)
|
1169
|
-
*/
|
1170
|
-
|
1171
|
-
/*****************************************************
|
1172
|
-
* math2/shtlib/shlbwl.f
|
1173
|
-
****************************************************/
|
1174
|
-
/*
|
1175
|
-
* SUBROUTINE SHLBWL(MM,M,SD,WS,WA,PMJ)
|
1176
|
-
* REAL SD(0:MM+1)
|
1177
|
-
* REAL WS,WA
|
1178
|
-
* REAL PMJ(0:MM+1)
|
1179
|
-
*/
|
1180
|
-
/*
|
1181
|
-
SUBROUTINE SHLBWL(
|
1182
|
-
INTEGER MM
|
1183
|
-
INTEGER M
|
1184
|
-
REAL (0:MM+1) SD
|
1185
|
-
REAL o WS
|
1186
|
-
REAL o WA
|
1187
|
-
REAL (0:MM+1) PMJ
|
1188
|
-
)
|
1189
|
-
*/
|
1190
|
-
|
1191
|
-
/*****************************************************
|
1192
|
-
* math2/shtlib/shlssd.f
|
1193
|
-
****************************************************/
|
1194
|
-
/*
|
1195
|
-
* SUBROUTINE SHLSSD(MM,M,ISW,SM,SD,R)
|
1196
|
-
* REAL SM(M:MM)
|
1197
|
-
* REAL SD(0:MM+1)
|
1198
|
-
* REAL R((MM+1)*(MM+1))
|
1199
|
-
*/
|
1200
|
-
/*
|
1201
|
-
SUBROUTINE SHLSSD(
|
1202
|
-
INTEGER MM
|
1203
|
-
INTEGER M
|
1204
|
-
INTEGER ISW
|
1205
|
-
REAL (M:MM) SM
|
1206
|
-
REAL (0:MM+1) o SD
|
1207
|
-
REAL ((MM+1)*(MM+1)) R
|
1208
|
-
)
|
1209
|
-
*/
|
1210
|
-
|
1211
|
-
/*****************************************************
|
1212
|
-
* math2/shtlib/shlsds.f
|
1213
|
-
****************************************************/
|
1214
|
-
/*
|
1215
|
-
* SUBROUTINE SHLSDS(MM,M,ISW,SD,SM,R)
|
1216
|
-
* REAL SD(0:MM+1)
|
1217
|
-
* REAL SM(M:MM)
|
1218
|
-
* REAL R((MM+1)*(MM+1))
|
1219
|
-
*/
|
1220
|
-
/*
|
1221
|
-
SUBROUTINE SHLSDS(
|
1222
|
-
INTEGER MM
|
1223
|
-
INTEGER M
|
1224
|
-
INTEGER ISW
|
1225
|
-
REAL (0:MM+1) SD
|
1226
|
-
REAL (M:MM) o SM
|
1227
|
-
REAL ((MM+1)*(MM+1)) R
|
1228
|
-
)
|
1229
|
-
*/
|
1230
|
-
|
1231
|
-
/*****************************************************
|
1232
|
-
* math2/shtlib/shpfun.f
|
1233
|
-
****************************************************/
|
1234
|
-
/*
|
1235
|
-
* SUBROUTINE SHPFUN(MM,JM,M,FUN,PM,YS,YC,PY,R)
|
1236
|
-
* REAL FUN(-JM:JM,M:MM)
|
1237
|
-
* REAL PM(0:MM+1,0:JM)
|
1238
|
-
* REAL YS(0:JM),YC(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
|
1239
|
-
*/
|
1240
|
-
/*
|
1241
|
-
SUBROUTINE SHPFUN(
|
1242
|
-
INTEGER MM
|
1243
|
-
INTEGER JM
|
1244
|
-
INTEGER M
|
1245
|
-
REAL (-JM:JM,M:MM) o FUN
|
1246
|
-
REAL (0:MM+1,0:JM) t PM
|
1247
|
-
REAL (0:JM) t YS
|
1248
|
-
REAL (0:JM) t YC
|
1249
|
-
REAL (2,0:JM,0:MM) t PY
|
1250
|
-
REAL ((MM+1)*(MM+1)) t R
|
1251
|
-
)
|
1252
|
-
*/
|
1253
|
-
|
1254
|
-
/*****************************************************
|
1255
|
-
* math2/shtlib/shppma.f
|
1256
|
-
****************************************************/
|
1257
|
-
/*
|
1258
|
-
* SUBROUTINE SHPPMA(MM,JM,M,PM,YS,PY,R)
|
1259
|
-
* REAL PM(0:MM+1,0:JM)
|
1260
|
-
* REAL YS(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
|
1261
|
-
*/
|
1262
|
-
/*
|
1263
|
-
SUBROUTINE SHPPMA(
|
1264
|
-
INTEGER MM
|
1265
|
-
INTEGER JM
|
1266
|
-
INTEGER M
|
1267
|
-
REAL (0:MM+1,0:JM) io PM
|
1268
|
-
REAL (0:JM) YS
|
1269
|
-
REAL (2,0:JM,0:MM) PY
|
1270
|
-
REAL ((MM+1)*(MM+1)) R
|
1271
|
-
)
|
1272
|
-
*/
|
1273
|
-
|
1274
|
-
/*****************************************************
|
1275
|
-
* math2/shtlib/shppmj.f
|
1276
|
-
****************************************************/
|
1277
|
-
/*
|
1278
|
-
* SUBROUTINE SHPPMJ(MM,JM,M,J,PMJ,YS,PY,R)
|
1279
|
-
* REAL PMJ(0:MM+1)
|
1280
|
-
* REAL YS(0:JM),PY(2,0:JM,0:MM),R((MM+1)*(MM+1))
|
1281
|
-
*/
|
1282
|
-
/*
|
1283
|
-
SUBROUTINE SHPPMJ(
|
1284
|
-
INTEGER MM
|
1285
|
-
INTEGER JM
|
1286
|
-
INTEGER M
|
1287
|
-
INTEGER J
|
1288
|
-
REAL (0:MM+1) io PMJ
|
1289
|
-
REAL (0:JM) YS
|
1290
|
-
REAL (2,0:JM,0:MM) PY
|
1291
|
-
REAL ((MM+1)*(MM+1)) R
|
1292
|
-
)
|
1293
|
-
*/
|
1294
|
-
|