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.
@@ -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
-