Initial import of new "fftpack" package:

Single precision Fortran FFT subroutines

13 files changed, 145 insertions, 0 deletions

diff --git a/math/fftpack/patches/patch-aa b/math/fftpack/patches/patch-aa
new file mode 100644
index 00000000000..8e9e8756d07
--- /dev/null
+++ b/math/fftpack/patches/patch-aa
@@ -0,0 +1,11 @@
+$NetBSD: patch-aa,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cfftb1.f.orig	Fri Nov 17 16:48:23 2000
++++ cfftb1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE CFFTB1 (N,C,CH,WA,IFAC)
+-      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(1)
++      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(*)
+       NF = IFAC(2)
+       NA = 0
+       L1 = 1
diff --git a/math/fftpack/patches/patch-ab b/math/fftpack/patches/patch-ab
new file mode 100644
index 00000000000..f617a4e5259
--- /dev/null
+++ b/math/fftpack/patches/patch-ab
@@ -0,0 +1,11 @@
+$NetBSD: patch-ab,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cfftf1.f.orig	Fri Nov 17 16:48:23 2000
++++ cfftf1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE CFFTF1 (N,C,CH,WA,IFAC)
+-      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(1)
++      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(*)
+       NF = IFAC(2)
+       NA = 0
+       L1 = 1
diff --git a/math/fftpack/patches/patch-ac b/math/fftpack/patches/patch-ac
new file mode 100644
index 00000000000..39f46278ec5
--- /dev/null
+++ b/math/fftpack/patches/patch-ac
@@ -0,0 +1,11 @@
+$NetBSD: patch-ac,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cffti1.f.orig	Fri Nov 17 16:48:23 2000
++++ cffti1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE CFFTI1 (N,WA,IFAC)
+-      DIMENSION       WA(1)      ,IFAC(1)    ,NTRYH(4)
++      DIMENSION       WA(1)      ,IFAC(*)    ,NTRYH(4)
+       DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/3,4,2,5/
+       NL = N
+       NF = 0
diff --git a/math/fftpack/patches/patch-ad b/math/fftpack/patches/patch-ad
new file mode 100644
index 00000000000..6028dc059c9
--- /dev/null
+++ b/math/fftpack/patches/patch-ad
@@ -0,0 +1,11 @@
+$NetBSD: patch-ad,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cosqb.f.orig	Fri Nov 17 16:48:23 2000
++++ cosqb.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE COSQB (N,X,WSAVE)
+-      DIMENSION       X(1)       ,WSAVE(1)
++      DIMENSION       X(*)       ,WSAVE(1)
+       DATA TSQRT2 /2.82842712474619/
+       IF (N-2) 101,102,103
+   101 X(1) = 4.*X(1)
diff --git a/math/fftpack/patches/patch-ae b/math/fftpack/patches/patch-ae
new file mode 100644
index 00000000000..59df70ffbb2
--- /dev/null
+++ b/math/fftpack/patches/patch-ae
@@ -0,0 +1,11 @@
+$NetBSD: patch-ae,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cosqf.f.orig	Fri Nov 17 16:48:23 2000
++++ cosqf.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE COSQF (N,X,WSAVE)
+-      DIMENSION       X(1)       ,WSAVE(1)
++      DIMENSION       X(*)       ,WSAVE(1)
+       DATA SQRT2 /1.4142135623731/
+       IF (N-2) 102,101,103
+   101 TSQX = SQRT2*X(2)
diff --git a/math/fftpack/patches/patch-af b/math/fftpack/patches/patch-af
new file mode 100644
index 00000000000..b2ead1bb56b
--- /dev/null
+++ b/math/fftpack/patches/patch-af
@@ -0,0 +1,11 @@
+$NetBSD: patch-af,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- cost.f.orig	Fri Nov 17 16:48:23 2000
++++ cost.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE COST (N,X,WSAVE)
+-      DIMENSION       X(1)       ,WSAVE(1)
++      DIMENSION       X(*)       ,WSAVE(1)
+       NM1 = N-1
+       NP1 = N+1
+       NS2 = N/2
diff --git a/math/fftpack/patches/patch-ag b/math/fftpack/patches/patch-ag
new file mode 100644
index 00000000000..94984780262
--- /dev/null
+++ b/math/fftpack/patches/patch-ag
@@ -0,0 +1,11 @@
+$NetBSD: patch-ag,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- ezfft1.f.orig	Fri Nov 17 16:48:23 2000
++++ ezfft1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE EZFFT1 (N,WA,IFAC)
+-      DIMENSION       WA(1)      ,IFAC(1)    ,NTRYH(4)
++      DIMENSION       WA(1)      ,IFAC(*)    ,NTRYH(4)
+       DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/4,2,3,5/
+      1    ,TPI/6.28318530717959/
+       NL = N
diff --git a/math/fftpack/patches/patch-ah b/math/fftpack/patches/patch-ah
new file mode 100644
index 00000000000..682ba905388
--- /dev/null
+++ b/math/fftpack/patches/patch-ah
@@ -0,0 +1,11 @@
+$NetBSD: patch-ah,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- ezfftb.f.orig	Fri Nov 17 16:48:23 2000
++++ ezfftb.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE EZFFTB (N,R,AZERO,A,B,WSAVE)
+-      DIMENSION       R(1)       ,A(1)       ,B(1)       ,WSAVE(1)
++      DIMENSION       R(*)       ,A(1)       ,B(1)       ,WSAVE(1)
+       IF (N-2) 101,102,103
+   101 R(1) = AZERO
+       RETURN
diff --git a/math/fftpack/patches/patch-ai b/math/fftpack/patches/patch-ai
new file mode 100644
index 00000000000..8de5ce8db6e
--- /dev/null
+++ b/math/fftpack/patches/patch-ai
@@ -0,0 +1,13 @@
+$NetBSD: patch-ai,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- ezfftf.f.orig	Fri Nov 17 16:48:24 2000
++++ ezfftf.f
+@@ -2,7 +2,7 @@
+ C
+ C                       VERSION 3  JUNE 1979
+ C
+-      DIMENSION       R(1)       ,A(1)       ,B(1)       ,WSAVE(1)
++      DIMENSION       R(*)       ,A(1)       ,B(1)       ,WSAVE(1)
+       IF (N-2) 101,102,103
+   101 AZERO = R(1)
+       RETURN
diff --git a/math/fftpack/patches/patch-aj b/math/fftpack/patches/patch-aj
new file mode 100644
index 00000000000..579c44604fc
--- /dev/null
+++ b/math/fftpack/patches/patch-aj
@@ -0,0 +1,11 @@
+$NetBSD: patch-aj,v 1.1.1.1 2000/11/28 02:17:02 jtb Exp $
+
+--- rfftb1.f.orig	Fri Nov 17 16:48:24 2000
++++ rfftb1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE RFFTB1 (N,C,CH,WA,IFAC)
+-      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(1)
++      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(*)
+       NF = IFAC(2)
+       NA = 0
+       L1 = 1
diff --git a/math/fftpack/patches/patch-ak b/math/fftpack/patches/patch-ak
new file mode 100644
index 00000000000..503446b7e6d
--- /dev/null
+++ b/math/fftpack/patches/patch-ak
@@ -0,0 +1,11 @@
+$NetBSD: patch-ak,v 1.1.1.1 2000/11/28 02:17:03 jtb Exp $
+
+--- rfftf1.f.orig	Fri Nov 17 16:48:24 2000
++++ rfftf1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE RFFTF1 (N,C,CH,WA,IFAC)
+-      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(1)
++      DIMENSION       CH(1)      ,C(1)       ,WA(1)      ,IFAC(*)
+       NF = IFAC(2)
+       NA = 1
+       L2 = N
diff --git a/math/fftpack/patches/patch-al b/math/fftpack/patches/patch-al
new file mode 100644
index 00000000000..996199f83c4
--- /dev/null
+++ b/math/fftpack/patches/patch-al
@@ -0,0 +1,11 @@
+$NetBSD: patch-al,v 1.1.1.1 2000/11/28 02:17:02 jtb Exp $
+
+--- rffti1.f.orig	Fri Nov 17 16:48:24 2000
++++ rffti1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE RFFTI1 (N,WA,IFAC)
+-      DIMENSION       WA(1)      ,IFAC(1)    ,NTRYH(4)
++      DIMENSION       WA(1)      ,IFAC(*)    ,NTRYH(4)
+       DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/4,2,3,5/
+       NL = N
+       NF = 0
diff --git a/math/fftpack/patches/patch-am b/math/fftpack/patches/patch-am
new file mode 100644
index 00000000000..b2e778b9799
--- /dev/null
+++ b/math/fftpack/patches/patch-am
@@ -0,0 +1,11 @@
+$NetBSD: patch-am,v 1.1.1.1 2000/11/28 02:17:02 jtb Exp $
+
+--- sint1.f.orig	Fri Nov 17 16:48:24 2000
++++ sint1.f
+@@ -1,5 +1,5 @@
+       SUBROUTINE SINT1(N,WAR,WAS,XH,X,IFAC)
+-      DIMENSION WAR(1),WAS(1),X(1),XH(1),IFAC(1)
++      DIMENSION WAR(1),WAS(1),X(1),XH(*),IFAC(1)
+       DATA SQRT3 /1.73205080756888/
+       DO 100 I=1,N
+       XH(I) = WAR(I)