head 1.5; access; symbols pkgsrc-2026Q1:1.5.0.148 pkgsrc-2026Q1-base:1.5 pkgsrc-2025Q4:1.5.0.146 pkgsrc-2025Q4-base:1.5 pkgsrc-2025Q3:1.5.0.144 pkgsrc-2025Q3-base:1.5 pkgsrc-2025Q2:1.5.0.142 pkgsrc-2025Q2-base:1.5 pkgsrc-2025Q1:1.5.0.140 pkgsrc-2025Q1-base:1.5 pkgsrc-2024Q4:1.5.0.138 pkgsrc-2024Q4-base:1.5 pkgsrc-2024Q3:1.5.0.136 pkgsrc-2024Q3-base:1.5 pkgsrc-2024Q2:1.5.0.134 pkgsrc-2024Q2-base:1.5 pkgsrc-2024Q1:1.5.0.132 pkgsrc-2024Q1-base:1.5 pkgsrc-2023Q4:1.5.0.130 pkgsrc-2023Q4-base:1.5 pkgsrc-2023Q3:1.5.0.128 pkgsrc-2023Q3-base:1.5 pkgsrc-2023Q2:1.5.0.126 pkgsrc-2023Q2-base:1.5 pkgsrc-2023Q1:1.5.0.124 pkgsrc-2023Q1-base:1.5 pkgsrc-2022Q4:1.5.0.122 pkgsrc-2022Q4-base:1.5 pkgsrc-2022Q3:1.5.0.120 pkgsrc-2022Q3-base:1.5 pkgsrc-2022Q2:1.5.0.118 pkgsrc-2022Q2-base:1.5 pkgsrc-2022Q1:1.5.0.116 pkgsrc-2022Q1-base:1.5 pkgsrc-2021Q4:1.5.0.114 pkgsrc-2021Q4-base:1.5 pkgsrc-2021Q3:1.5.0.112 pkgsrc-2021Q3-base:1.5 pkgsrc-2021Q2:1.5.0.110 pkgsrc-2021Q2-base:1.5 pkgsrc-2021Q1:1.5.0.108 pkgsrc-2021Q1-base:1.5 pkgsrc-2020Q4:1.5.0.106 pkgsrc-2020Q4-base:1.5 pkgsrc-2020Q3:1.5.0.104 pkgsrc-2020Q3-base:1.5 pkgsrc-2020Q2:1.5.0.100 pkgsrc-2020Q2-base:1.5 pkgsrc-2020Q1:1.5.0.80 pkgsrc-2020Q1-base:1.5 pkgsrc-2019Q4:1.5.0.102 pkgsrc-2019Q4-base:1.5 pkgsrc-2019Q3:1.5.0.98 pkgsrc-2019Q3-base:1.5 pkgsrc-2019Q2:1.5.0.96 pkgsrc-2019Q2-base:1.5 pkgsrc-2019Q1:1.5.0.94 pkgsrc-2019Q1-base:1.5 pkgsrc-2018Q4:1.5.0.92 pkgsrc-2018Q4-base:1.5 pkgsrc-2018Q3:1.5.0.90 pkgsrc-2018Q3-base:1.5 pkgsrc-2018Q2:1.5.0.88 pkgsrc-2018Q2-base:1.5 pkgsrc-2018Q1:1.5.0.86 pkgsrc-2018Q1-base:1.5 pkgsrc-2017Q4:1.5.0.84 pkgsrc-2017Q4-base:1.5 pkgsrc-2017Q3:1.5.0.82 pkgsrc-2017Q3-base:1.5 pkgsrc-2017Q2:1.5.0.78 pkgsrc-2017Q2-base:1.5 pkgsrc-2017Q1:1.5.0.76 pkgsrc-2017Q1-base:1.5 pkgsrc-2016Q4:1.5.0.74 pkgsrc-2016Q4-base:1.5 pkgsrc-2016Q3:1.5.0.72 pkgsrc-2016Q3-base:1.5 pkgsrc-2016Q2:1.5.0.70 pkgsrc-2016Q2-base:1.5 pkgsrc-2016Q1:1.5.0.68 pkgsrc-2016Q1-base:1.5 pkgsrc-2015Q4:1.5.0.66 pkgsrc-2015Q4-base:1.5 pkgsrc-2015Q3:1.5.0.64 pkgsrc-2015Q3-base:1.5 pkgsrc-2015Q2:1.5.0.62 pkgsrc-2015Q2-base:1.5 pkgsrc-2015Q1:1.5.0.60 pkgsrc-2015Q1-base:1.5 pkgsrc-2014Q4:1.5.0.58 pkgsrc-2014Q4-base:1.5 pkgsrc-2014Q3:1.5.0.56 pkgsrc-2014Q3-base:1.5 pkgsrc-2014Q2:1.5.0.54 pkgsrc-2014Q2-base:1.5 pkgsrc-2014Q1:1.5.0.52 pkgsrc-2014Q1-base:1.5 pkgsrc-2013Q4:1.5.0.50 pkgsrc-2013Q4-base:1.5 pkgsrc-2013Q3:1.5.0.48 pkgsrc-2013Q3-base:1.5 pkgsrc-2013Q2:1.5.0.46 pkgsrc-2013Q2-base:1.5 pkgsrc-2013Q1:1.5.0.44 pkgsrc-2013Q1-base:1.5 pkgsrc-2012Q4:1.5.0.42 pkgsrc-2012Q4-base:1.5 pkgsrc-2012Q3:1.5.0.40 pkgsrc-2012Q3-base:1.5 pkgsrc-2012Q2:1.5.0.38 pkgsrc-2012Q2-base:1.5 pkgsrc-2012Q1:1.5.0.36 pkgsrc-2012Q1-base:1.5 pkgsrc-2011Q4:1.5.0.34 pkgsrc-2011Q4-base:1.5 pkgsrc-2011Q3:1.5.0.32 pkgsrc-2011Q3-base:1.5 pkgsrc-2011Q2:1.5.0.30 pkgsrc-2011Q2-base:1.5 pkgsrc-2011Q1:1.5.0.28 pkgsrc-2011Q1-base:1.5 pkgsrc-2010Q4:1.5.0.26 pkgsrc-2010Q4-base:1.5 pkgsrc-2010Q3:1.5.0.24 pkgsrc-2010Q3-base:1.5 pkgsrc-2010Q2:1.5.0.22 pkgsrc-2010Q2-base:1.5 pkgsrc-2010Q1:1.5.0.20 pkgsrc-2010Q1-base:1.5 pkgsrc-2009Q4:1.5.0.18 pkgsrc-2009Q4-base:1.5 pkgsrc-2009Q3:1.5.0.16 pkgsrc-2009Q3-base:1.5 pkgsrc-2009Q2:1.5.0.14 pkgsrc-2009Q2-base:1.5 pkgsrc-2009Q1:1.5.0.12 pkgsrc-2009Q1-base:1.5 pkgsrc-2008Q4:1.5.0.10 pkgsrc-2008Q4-base:1.5 pkgsrc-2008Q3:1.5.0.8 pkgsrc-2008Q3-base:1.5 cube-native-xorg:1.5.0.6 cube-native-xorg-base:1.5 pkgsrc-2008Q2:1.5.0.4 pkgsrc-2008Q2-base:1.5 cwrapper:1.5.0.2 pkgsrc-2008Q1:1.4.0.20 pkgsrc-2008Q1-base:1.4 pkgsrc-2007Q4:1.4.0.18 pkgsrc-2007Q4-base:1.4 pkgsrc-2007Q3:1.4.0.16 pkgsrc-2007Q3-base:1.4 pkgsrc-2007Q2:1.4.0.14 pkgsrc-2007Q2-base:1.4 pkgsrc-2007Q1:1.4.0.12 pkgsrc-2007Q1-base:1.4 pkgsrc-2006Q4:1.4.0.10 pkgsrc-2006Q4-base:1.4 pkgsrc-2006Q3:1.4.0.8 pkgsrc-2006Q3-base:1.4 pkgsrc-2006Q2:1.4.0.6 pkgsrc-2006Q2-base:1.4 pkgsrc-2006Q1:1.4.0.4 pkgsrc-2006Q1-base:1.4 pkgsrc-2005Q4:1.4.0.2 pkgsrc-2005Q4-base:1.4 pkgsrc-2005Q3:1.1.0.24 pkgsrc-2005Q3-base:1.1 pkgsrc-2005Q2:1.1.0.22 pkgsrc-2005Q2-base:1.1 pkgsrc-2005Q1:1.1.0.20 pkgsrc-2005Q1-base:1.1 pkgsrc-2004Q4:1.1.0.18 pkgsrc-2004Q4-base:1.1 pkgsrc-2004Q3:1.1.0.16 pkgsrc-2004Q3-base:1.1 pkgsrc-2004Q2:1.1.0.14 pkgsrc-2004Q2-base:1.1 pkgsrc-2004Q1:1.1.0.12 pkgsrc-2004Q1-base:1.1 pkgsrc-2003Q4:1.1.0.10 pkgsrc-2003Q4-base:1.1 netbsd-1-6-1:1.1.0.6 netbsd-1-6-1-base:1.1 netbsd-1-6:1.1.0.8 netbsd-1-6-RELEASE-base:1.1 pkgviews:1.1.0.4 pkgviews-base:1.1 buildlink2:1.1.0.2 buildlink2-base:1.1 netbsd-1-5-PATCH003:1.1; locks; strict; comment @# @; 1.5 date 2008.06.21.20.00.22; author joerg; state Exp; branches; next 1.4; 1.4 date 2005.12.20.17.58.51; author joerg; state dead; branches; next 1.3; 1.3 date 2005.12.13.17.20.53; author joerg; state Exp; branches; next 1.2; 1.2 date 2005.11.05.21.04.37; author joerg; state dead; branches; next 1.1; 1.1 date 2001.12.07.19.46.07; author jmc; state Exp; branches; next ; desc @@ 1.5 log @Fix build with GCC 4 and f2c. @ text @$NetBSD$ Reorder functions in use order to fix compilation with f2c. --- f_source/sciport/fftd.f.orig 2008-06-21 13:43:11.000000000 +0000 +++ f_source/sciport/fftd.f @@@@ -1,3 +1,192 @@@@ +C----------------------------------------------- ************ +C CABLE2 +C ************ + SUBROUTINE SPDCABLE2(NN,WORK) +C + DOUBLE PRECISION WORK(2,NN),P2,TWOPI + DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ +C + N = 2 * NN + P2 = TWOPI/N + DO 10 I=1, NN + WORK(1,I) = DCOS(P2 * (I-1)) + WORK(2,I) = DSIN(P2 * (I-1)) + 10 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C RCONV2 +C ************ + SUBROUTINE SPDRCONV2(N,CY,C,CH) +C + DOUBLE COMPLEX CY(1) + DOUBLE PRECISION CH(N/2,2),P(2,1),C(2,N/2) + DOUBLE PRECISION X,Y,Z,Z1 +C + N2 = N/2 + P(1,1) = (C(1,1) + C(2,1)) * 2 + P(2,1) = (C(1,1) - C(2,1)) * 2 + CY(1) = DCMPLX(P(1,1),0.0D0) + CY(N2+1) = DCMPLX(P(2,1),0.0D0) + K = N2 + DO 10 I=2, N2 + X = C(1,I)+C(1,K) + Y = C(2,I)+C(2,K) + Z = C(1,I)-C(1,K) + Z1= C(2,I)-C(2,K) + P(1,1) = X + CH(I,1) * Y - CH(I,2) * Z + P(2,1) = Z1 - CH(I,2) * Y - CH(I,1) * Z + CY(I) = DCMPLX(P(1,1),P(2,1)) + K = K - 1 + 10 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C RCONV1 +C ************ + SUBROUTINE SPDRCONV1(N,CY,C,CH) +C + DOUBLE COMPLEX CY(1) + DOUBLE PRECISION CH(N/2,2),P(2,1),C(2,N/2) + DOUBLE PRECISION X,Y,Z,Z1 +C + N2 = N/2 + P(1,1) = (C(1,1) + C(2,1)) * 2 + P(2,1) = (C(1,1) - C(2,1)) * 2 + CY(1) = DCMPLX(P(1,1),0.0D0) + CY(N2+1) = DCMPLX(P(2,1),0.0D0) + K = N2 + DO 10 I=2, N2 + X = C(1,I)+C(1,K) + Y = C(2,I)+C(2,K) + Z = C(1,I)-C(1,K) + Z1= C(2,I)-C(2,K) + P(1,1) = X + CH(I,1) * Y + CH(I,2) * Z + P(2,1) = Z1 + CH(I,2) * Y - CH(I,1) * Z + CY(I) = DCMPLX(P(1,1),P(2,1)) + K = K - 1 + 10 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C RTOCK3 +C ************ + SUBROUTINE SPDRTOCK3(LS,NS,C,CH,CH2) +C + DOUBLE COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS) + DOUBLE PRECISION CH2(2,NS,LS,2) +C + IF (LS .GT. NS) GOTO 30 + DO 600 I=1, LS + DO 600 J=1, NS + WYK =CONJG(DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2))) * CH(J,2,I) + C(J,I,1) = CH(J,1,I) + WYK + C(J,I,2) = CH(J,1,I) - WYK + 600 CONTINUE + RETURN + 30 CONTINUE + DO 800 J=1, NS + DO 800 I=1, LS + WYK =CONJG(DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2))) * CH(J,2,I) + C(J,I,1) = CH(J,1,I) + WYK + C(J,I,2) = CH(J,1,I) - WYK + 800 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C RABLE1 +C ************ + SUBROUTINE SPDRABLE1(NN,WORK) +C + DOUBLE PRECISION WORK(NN,2),P2,TWOPI + DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ +C + N = 2 * NN + P2 = TWOPI/N + DO 10 I=1, NN + WORK(I,1) = DCOS(P2 * (I-1)) + WORK(I,2) = DSIN(P2 * (I-1)) + 10 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C RTOCK2 +C ************ + SUBROUTINE SPDRTOCK2(LS,NS,C,CH,CH2) +C + DOUBLE COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS) + DOUBLE PRECISION CH2(2,NS,LS,2) +C + IF (LS .GT. NS) GOTO 20 + DO 200 I=1, LS + DO 200 J=1, NS + WYK = DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2)) * CH(J,2,I) + C(J,I,1) = CH(J,1,I) + WYK + C(J,I,2) = CH(J,1,I) - WYK + 200 CONTINUE + RETURN + 20 CONTINUE + DO 400 J=1, NS + DO 400 I=1, LS + WYK = DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2)) * CH(J,2,I) + C(J,I,1) = CH(J,1,I) + WYK + C(J,I,2) = CH(J,1,I) - WYK + 400 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C CRBLE1 +C ************ + SUBROUTINE SPDCRBLE1(NN,WORK) +C + DOUBLE PRECISION WORK(NN,2),P2,TWOPI + DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ +C + N = 2 * NN + P2 = TWOPI/N + DO 10 I=1, NN + WORK(I,1) = DCOS(P2 * (I-1)) + WORK(I,2) = DSIN(P2 * (I-1)) + 10 CONTINUE + RETURN + END + +C----------------------------------------------- ************ +C CRFORM +C ************ + SUBROUTINE SPDCRFORM(IX,NS,NDIV2,CX,C,CH2) +C + DOUBLE COMPLEX CX(1),WYK1,C(NS,2),WYK + DOUBLE PRECISION CH2(NDIV2,2) +C + IF (IX .GT. 0) GOTO 50 + K = NS + 1 + DO 10 I=1, NS + WYK = CONJG(CX(NDIV2-I+2)) + C(I,1)= CX(I)+WYK + (CX(I) - WYK) * DCMPLX(CH2(I,2),CH2(I,1)) + WYK1 = CONJG(CX(NDIV2-K+2)) + C(I,2)= CX(K)+WYK1+ (CX(K) -WYK1) * DCMPLX(CH2(K,2),CH2(K,1)) + K = K + 1 + 10 CONTINUE + RETURN + 50 CONTINUE + K = NS + 1 + DO 20 I=1, NS + WYK = CONJG(CX(NDIV2-I+2)) + C(I,1)= CX(I)+WYK + (CX(I) - WYK) * DCMPLX(-CH2(I,2),CH2(I,1)) + WYK1 = CONJG(CX(NDIV2-K+2)) + C(I,2)= CX(K)+WYK1 +(CX(K) -WYK1) * DCMPLX(-CH2(K,2),CH2(K,1)) + K = K + 1 + 20 CONTINUE + RETURN + END + C------------------------------------------------------------- ************ C CRFFT2 C ************ @@@@ -62,36 +251,6 @@@@ C END C----------------------------------------------- ************ -C CRFORM -C ************ - SUBROUTINE SPDCRFORM(IX,NS,NDIV2,CX,C,CH2) -C - DOUBLE COMPLEX CX(1),WYK1,C(NS,2),WYK - DOUBLE PRECISION CH2(NDIV2,2) -C - IF (IX .GT. 0) GOTO 50 - K = NS + 1 - DO 10 I=1, NS - WYK = CONJG(CX(NDIV2-I+2)) - C(I,1)= CX(I)+WYK + (CX(I) - WYK) * DCMPLX(CH2(I,2),CH2(I,1)) - WYK1 = CONJG(CX(NDIV2-K+2)) - C(I,2)= CX(K)+WYK1+ (CX(K) -WYK1) * DCMPLX(CH2(K,2),CH2(K,1)) - K = K + 1 - 10 CONTINUE - RETURN - 50 CONTINUE - K = NS + 1 - DO 20 I=1, NS - WYK = CONJG(CX(NDIV2-I+2)) - C(I,1)= CX(I)+WYK + (CX(I) - WYK) * DCMPLX(-CH2(I,2),CH2(I,1)) - WYK1 = CONJG(CX(NDIV2-K+2)) - C(I,2)= CX(K)+WYK1 +(CX(K) -WYK1) * DCMPLX(-CH2(K,2),CH2(K,1)) - K = K + 1 - 20 CONTINUE - RETURN - END - -C----------------------------------------------- ************ C CROCK1 C ************ SUBROUTINE SPDCROCK1(NS,C,CH) @@@@ -157,23 +316,6 @@@@ C RETURN END -C----------------------------------------------- ************ -C CRBLE1 -C ************ - SUBROUTINE SPDCRBLE1(NN,WORK) -C - DOUBLE PRECISION WORK(NN,2),P2,TWOPI - DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ -C - N = 2 * NN - P2 = TWOPI/N - DO 10 I=1, NN - WORK(I,1) = DCOS(P2 * (I-1)) - WORK(I,2) = DSIN(P2 * (I-1)) - 10 CONTINUE - RETURN - END - C------------------------------------------------------------- ************ C RCFFT2 C ************ @@@@ -236,32 +378,6 @@@@ C END C----------------------------------------------- ************ -C RTOCK2 -C ************ - SUBROUTINE SPDRTOCK2(LS,NS,C,CH,CH2) -C - DOUBLE COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS) - DOUBLE PRECISION CH2(2,NS,LS,2) -C - IF (LS .GT. NS) GOTO 20 - DO 200 I=1, LS - DO 200 J=1, NS - WYK = DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2)) * CH(J,2,I) - C(J,I,1) = CH(J,1,I) + WYK - C(J,I,2) = CH(J,1,I) - WYK - 200 CONTINUE - RETURN - 20 CONTINUE - DO 400 J=1, NS - DO 400 I=1, LS - WYK = DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2)) * CH(J,2,I) - C(J,I,1) = CH(J,1,I) + WYK - C(J,I,2) = CH(J,1,I) - WYK - 400 CONTINUE - RETURN - END - -C----------------------------------------------- ************ C RTOCK1 C ************ SUBROUTINE SPDRTOCK1(NS,C,CH) @@@@ -275,105 +391,6 @@@@ C RETURN END -C----------------------------------------------- ************ -C RTOCK3 -C ************ - SUBROUTINE SPDRTOCK3(LS,NS,C,CH,CH2) -C - DOUBLE COMPLEX WYK,C(NS,LS,2),CH(NS,2,LS) - DOUBLE PRECISION CH2(2,NS,LS,2) -C - IF (LS .GT. NS) GOTO 30 - DO 600 I=1, LS - DO 600 J=1, NS - WYK =CONJG(DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2))) * CH(J,2,I) - C(J,I,1) = CH(J,1,I) + WYK - C(J,I,2) = CH(J,1,I) - WYK - 600 CONTINUE - RETURN - 30 CONTINUE - DO 800 J=1, NS - DO 800 I=1, LS - WYK =CONJG(DCMPLX(CH2(1,1,I,1),CH2(1,1,I,2))) * CH(J,2,I) - C(J,I,1) = CH(J,1,I) + WYK - C(J,I,2) = CH(J,1,I) - WYK - 800 CONTINUE - RETURN - END - -C----------------------------------------------- ************ -C RABLE1 -C ************ - SUBROUTINE SPDRABLE1(NN,WORK) -C - DOUBLE PRECISION WORK(NN,2),P2,TWOPI - DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ -C - N = 2 * NN - P2 = TWOPI/N - DO 10 I=1, NN - WORK(I,1) = DCOS(P2 * (I-1)) - WORK(I,2) = DSIN(P2 * (I-1)) - 10 CONTINUE - RETURN - END - -C----------------------------------------------- ************ -C RCONV1 -C ************ - SUBROUTINE SPDRCONV1(N,CY,C,CH) -C - DOUBLE COMPLEX CY(1) - DOUBLE PRECISION CH(N/2,2),P(2,1),C(2,N/2) - DOUBLE PRECISION X,Y,Z,Z1 -C - N2 = N/2 - P(1,1) = (C(1,1) + C(2,1)) * 2 - P(2,1) = (C(1,1) - C(2,1)) * 2 - CY(1) = DCMPLX(P(1,1),0.0D0) - CY(N2+1) = DCMPLX(P(2,1),0.0D0) - K = N2 - DO 10 I=2, N2 - X = C(1,I)+C(1,K) - Y = C(2,I)+C(2,K) - Z = C(1,I)-C(1,K) - Z1= C(2,I)-C(2,K) - P(1,1) = X + CH(I,1) * Y + CH(I,2) * Z - P(2,1) = Z1 + CH(I,2) * Y - CH(I,1) * Z - CY(I) = DCMPLX(P(1,1),P(2,1)) - K = K - 1 - 10 CONTINUE - RETURN - END - -C----------------------------------------------- ************ -C RCONV2 -C ************ - SUBROUTINE SPDRCONV2(N,CY,C,CH) -C - DOUBLE COMPLEX CY(1) - DOUBLE PRECISION CH(N/2,2),P(2,1),C(2,N/2) - DOUBLE PRECISION X,Y,Z,Z1 -C - N2 = N/2 - P(1,1) = (C(1,1) + C(2,1)) * 2 - P(2,1) = (C(1,1) - C(2,1)) * 2 - CY(1) = DCMPLX(P(1,1),0.0D0) - CY(N2+1) = DCMPLX(P(2,1),0.0D0) - K = N2 - DO 10 I=2, N2 - X = C(1,I)+C(1,K) - Y = C(2,I)+C(2,K) - Z = C(1,I)-C(1,K) - Z1= C(2,I)-C(2,K) - P(1,1) = X + CH(I,1) * Y - CH(I,2) * Z - P(2,1) = Z1 - CH(I,2) * Y - CH(I,1) * Z - CY(I) = DCMPLX(P(1,1),P(2,1)) - K = K - 1 - 10 CONTINUE - RETURN - END - C------------------------------------------------------------- ************ C CFFT2 C ************ @@@@ -516,23 +533,6 @@@@ C RETURN END -C----------------------------------------------- ************ -C CABLE2 -C ************ - SUBROUTINE SPDCABLE2(NN,WORK) -C - DOUBLE PRECISION WORK(2,NN),P2,TWOPI - DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ -C - N = 2 * NN - P2 = TWOPI/N - DO 10 I=1, NN - WORK(1,I) = DCOS(P2 * (I-1)) - WORK(2,I) = DSIN(P2 * (I-1)) - 10 CONTINUE - RETURN - END - C------------------------------------------------------------- ************ C ABORT C ************ @ 1.4 log @Backout last change. f2c is not correctly working here. @ text @d1 1 a1 1 $NetBSD: patch-ag,v 1.3 2005/12/13 17:20:53 joerg Exp $ d3 3 a5 1 --- f_source/sciport/fftd.f.orig 2005-12-13 13:09:28.000000000 +0000 d7 201 a207 4 @@@@ -66,8 +66,7 @@@@ C C ************ SUBROUTINE SPDCRFORM(IX,NS,NDIV2,CX,C,CH2) C d210 25 a234 5 + DOUBLE COMPLEX CX(1),WYK1,C(NS,2),WYK,CH2(NDIV2,2) C IF (IX .GT. 0) GOTO 50 K = NS + 1 @@@@ -162,7 +161,8 @@@@ C d236 10 a245 2 SUBROUTINE SPDCRBLE1(NN,WORK) C d247 22 a268 18 + DOUBLE COMPLEX WORK(NN,2) + DOUBLE PRECISION P2,TWOPI DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ C N = 2 * NN @@@@ -218,7 +218,7 @@@@ C LS = 2 NS = NS/2 60 CONTINUE - CALL SPDRTOCK3(LS,NS,CWORK(NN+1),CWORK(1),CWORK(N2)) + CALL SPDRTOCK3(LS,NS,CWORK(NN+1),CWORK(1),CWORK(N2)) IF (NS .EQ. 1) GOTO 70 LS = LS + LS NS = NS/2 @@@@ -240,8 +240,7 @@@@ C C ************ SUBROUTINE SPDRTOCK2(LS,NS,C,CH,CH2) C d271 21 a291 5 + DOUBLE 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 +279,7 @@@@ C d293 10 a302 2 SUBROUTINE SPDRTOCK3(LS,NS,C,CH,CH2) C d305 24 a328 8 + DOUBLE 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 +304,8 @@@@ C C ************ SUBROUTINE SPDRABLE1(NN,WORK) C d330 16 a345 9 + DOUBLE COMPLEX WORK(NN,2) + DOUBLE PRECISION P2,TWOPI DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ C N = 2 * NN @@@@ -323,8 +322,8 @@@@ C C ************ SUBROUTINE SPDRCONV1(N,CY,C,CH) C d348 26 a373 9 + DOUBLE COMPLEX CY(1),CH(N/2,2),C(2,N/2) + DOUBLE PRECISION P(2,1) DOUBLE PRECISION X,Y,Z,Z1 C N2 = N/2 @@@@ -351,8 +350,8 @@@@ C C ************ SUBROUTINE SPDRCONV2(N,CY,C,CH) C d376 33 a408 9 + DOUBLE COMPLEX CY(1),CH(N/2,2),C(2,N/2) + DOUBLE PRECISION P(2,1) DOUBLE PRECISION X,Y,Z,Z1 C N2 = N/2 @@@@ -521,7 +520,8 @@@@ C C ************ SUBROUTINE SPDCABLE2(NN,WORK) C d410 14 a423 5 + DOUBLE COMPLEX WORK(2,NN) + DOUBLE PRECISION P2,TWOPI DATA TWOPI /6.2831853071795864769252867665590057683943388D0/ C N = 2 * NN @ 1.3 log @Match function arguments with actual usage. Since neither me nor the maintainer is very fluent with Fortan, let's assume f2c knows what it was doing. @ text @d1 1 a1 1 $NetBSD$ @ 1.2 log @Do not intermingle build phases, CONFIGURE_DIRS can specify a list of directories to process. The removal of pre-configure fixes the config.* handling at the same time. Also specify the config.status overrides. Patch configure to not try to detect -lrfftw, just link the static version explicitly. Do the same for -lfftw. This ensures that we actually do get the static versions without having to worry about libtool being too clever. Retire manual config.sub hacks, this should be handled by config.guess already. just go ahead and commit that benchfft patch and let's never speak of it again :) @ text @d1 1 a1 1 $NetBSD: patch-ag,v 1.1 2001/12/07 19:46:07 jmc Exp $ d3 93 a95 11 --- config.sub.orig Fri Dec 7 02:05:59 2001 +++ config.sub Fri Dec 7 02:06:22 2001 @@@@ -182,7 +182,7 @@@@ | sparc64-* | mips64-* | mipsel-* \ | mips64el-* | mips64orion-* | mips64orionel-* \ | mipstx39-* | mipstx39el-* \ - | f301-*) + | f301-* | macppc-* ) ;; # Recognize the various machine names and aliases which stand # for a CPU type and a company and sometimes even an OS. @ 1.1 log @Add patches to config.sub so macppc is recognized @ text @d1 1 a1 1 $NetBSD$ @