diff options
author | Erich Eckner <git@eckner.net> | 2015-09-30 15:20:56 +0200 |
---|---|---|
committer | Erich Eckner <git@eckner.net> | 2015-09-30 15:20:56 +0200 |
commit | 07e43942e2038fd8375d7835d5e21eab3040e734 (patch) | |
tree | 8521d1b388a0ca3046c6f8dc7642980eb49ea383 | |
parent | 18d3285ee1611015c0aa2dea92a986c8a5cfdd36 (diff) | |
download | Plasmapropagation-07e43942e2038fd8375d7835d5e21eab3040e734.tar.xz |
FFT bei niedrigeren Frequenzen abgeschnitten
-rw-r--r-- | Physikunit.pas | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Physikunit.pas b/Physikunit.pas index 8aec754..68ad65d 100644 --- a/Physikunit.pas +++ b/Physikunit.pas @@ -899,8 +899,8 @@ begin for j:=0 to aX-1 do for k:=0 to aP div 2 do begin // *-i*k*2*pi; tmp:=(fftTmp+j+k*aX)^.re; - (fftTmp+j+k*aX)^.re:= 2*pi/dP/aP * k * (fftTmp+j+k*aX)^.im/aP * byte(3*k < aP); - (fftTmp+j+k*aX)^.im:=- 2*pi/dP/aP * k * tmp/aP * byte(3*k < aP); + (fftTmp+j+k*aX)^.re:= 2*pi/dP/aP * k * (fftTmp+j+k*aX)^.im/aP * byte(4*k < aP); + (fftTmp+j+k*aX)^.im:=- 2*pi/dP/aP * k * tmp/aP * byte(4*k < aP); end; fftw_execute(ffts[i,1,true]); |