w4kek
/
WSJT-X
mirror of https://github.com/saitohirga/WSJT-X.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Further refinements to complex to real FFT initialization

This commit is contained in:
Bill Somerville 2021-10-20 19:28:58 +01:00
parent 7695b7d331
commit 30e11106c8
No known key found for this signature in database
GPG Key ID: D864B06D1E81618F
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
lib/four2a.f90
View File

@ -21,7 +21,7 @@ subroutine four2a(a,nfft,ndim,isign,iform)
use fftw3 use fftw3
parameter (NPMAX=2100) !Max numberf of stored plans parameter (NPMAX=2100) !Max numberf of stored plans
parameter (NSMALL=16384) !Max size of "small" FFTs parameter (NSMALL=16385) !Max half complex size of "small" FFTs
complex a(nfft) !Array to be transformed complex a(nfft) !Array to be transformed
complex aa(NSMALL) !Local copy of "small" a() complex aa(NSMALL) !Local copy of "small" a()
integer nn(NPMAX),ns(NPMAX),nf(NPMAX) !Params of stored plans integer nn(NPMAX),ns(NPMAX),nf(NPMAX) !Params of stored plans
@ -67,7 +67,7 @@ subroutine four2a(a,nfft,ndim,isign,iform)
if(nfft.le.NSMALL) then if(nfft.le.NSMALL) then
jz=nfft jz=nfft
if(iform.le.0) jz=nfft/2 if(iform.le.0) jz=nfft/2+1
aa(1:jz)=a(1:jz) aa(1:jz)=a(1:jz)
endif endif
@ -87,7 +87,7 @@ subroutine four2a(a,nfft,ndim,isign,iform)
if(nfft.le.NSMALL) then if(nfft.le.NSMALL) then
jz=nfft jz=nfft
if(iform.le.0) jz=nfft/2 if(iform.le.0) jz=nfft/2+1
a(1:jz)=aa(1:jz) a(1:jz)=aa(1:jz)
endif endif
end if end if