1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
$NetBSD: patch-af,v 1.3 2005/12/13 17:20:53 joerg Exp $
--- f_source/sciport/ffts.f.orig 2005-12-13 12:44:18.000000000 +0000
+++ f_source/sciport/ffts.f
@@ -66,8 +66,7 @@ C
C ************
SUBROUTINE SPSCRFORM(IX,NS,NDIV2,CX,C,CH2)
C
- COMPLEX CX(1),WYK1,C(NS,2),WYK
- REAL CH2(NDIV2,2)
+ COMPLEX CX(1),WYK1,C(NS,2),WYK, CH2(NDIV2,2)
C
IF (IX .GT. 0) GOTO 50
K = NS + 1
@@ -110,8 +109,7 @@ C
C ************
SUBROUTINE SPSCROCK2(LS,NS,C,CH,CH2)
C
- COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS)
- REAL CH2(2,NS,LS,2)
+ COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS),CH2(2,NS,LS,2)
C
IF (LS .GT. NS) GOTO 20
DO 200 I=1, LS
@@ -136,8 +134,7 @@ C
C ************
SUBROUTINE SPSCROCK3(LS,NS,C,CH,CH2)
C
- COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS)
- REAL CH2(2,NS,LS,2)
+ COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS),CH2(2,NS,LS,2)
C
IF (LS .GT. NS) GOTO 30
DO 600 I=1, LS
@@ -162,7 +159,8 @@ C
C ************
SUBROUTINE SPSCRBLE1(NN,WORK)
C
- REAL WORK(NN,2),TWOPI
+ COMPLEX WORK(NN,2)
+ REAL TWOPI
DATA TWOPI /6.28318530717958647692/
C
N = 2 * NN
@@ -240,8 +238,7 @@ C
C ************
SUBROUTINE SPSRTOCK2(LS,NS,C,CH,CH2)
C
- COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS)
- REAL CH2(2,NS,LS,2)
+ COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS),CH2(2,NS,LS,2)
C
IF (LS .GT. NS) GOTO 20
DO 200 I=1, LS
@@ -280,8 +277,7 @@ C
C ************
SUBROUTINE SPSRTOCK3(LS,NS,C,CH,CH2)
C
- COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS)
- REAL CH2(2,NS,LS,2)
+ COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS),CH2(2,NS,LS,2)
C
IF (LS .GT. NS) GOTO 30
DO 600 I=1, LS
@@ -306,7 +302,8 @@ C
C ************
SUBROUTINE SPSRABLE1(NN,WORK)
C
- REAL WORK(NN,2),TWOPI
+ COMPLEX WORK(NN,2)
+ REAL TWOPI
DATA TWOPI /6.28318530717958647692/
C
N = 2 * NN
@@ -323,8 +320,8 @@ C
C ************
SUBROUTINE SPSRCONV1(N,CY,C,CH)
C
- COMPLEX CY(1)
- REAL CH(N/2,2),P(2,1),C(2,N/2)
+ COMPLEX CY(1),CH(N/2,2),C(2,N/2)
+ REAL P(2,1)
C
N2 = N/2
P(1,1) = (C(1,1) + C(2,1)) * 2
@@ -350,8 +347,8 @@ C
C ************
SUBROUTINE SPSRCONV2(N,CY,C,CH)
C
- COMPLEX CY(1)
- REAL CH(N/2,2),P(2,1),C(2,N/2)
+ COMPLEX CY(1),CH(N/2,2),C(2,N/2)
+ REAL P(2,1)
C
N2 = N/2
P(1,1) = (C(1,1) + C(2,1)) * 2
@@ -519,7 +516,8 @@ C
C ************
SUBROUTINE SPSCABLE2(NN,WORK)
C
- REAL WORK(2,NN),TWOPI
+ REAL TWOPI
+ COMPLEX WORK(2,NN)
DATA TWOPI /6.28318530717958647692/
C
N = 2 * NN
|