1
SUBROUTINE DCOSQF (N,X,WSAVE)
2
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
3
DIMENSION X(*) ,WSAVE(*)
4
DATA SQRT2 /1.41421356237309504880D0/
12
103 CALL DCOSQF1 (N,X,WSAVE,WSAVE(N+1))
16
SUBROUTINE DCOSQF1 (N,X,W,XH)
17
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
18
DIMENSION X(1) ,W(1) ,XH(1)
27
IF (MODN .EQ. 0) XH(NS2+1) = X(NS2+1)+X(NS2+1)
30
X(K) = W(K-1)*XH(KC)+W(KC-1)*XH(K)
31
X(KC) = W(K-1)*XH(K)-W(KC-1)*XH(KC)
33
IF (MODN .EQ. 0) X(NS2+1) = W(NS2)*XH(NS2+1)