ruby-dcl 1.6.3.1 → 1.7.0

Sign up to get free protection for your applications and to get access to all the features.
@@ -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
-