summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorErich Eckner <git@eckner.net>2015-09-30 15:20:56 +0200
committerErich Eckner <git@eckner.net>2015-09-30 15:20:56 +0200
commit07e43942e2038fd8375d7835d5e21eab3040e734 (patch)
tree8521d1b388a0ca3046c6f8dc7642980eb49ea383
parent18d3285ee1611015c0aa2dea92a986c8a5cfdd36 (diff)
downloadPlasmapropagation-07e43942e2038fd8375d7835d5e21eab3040e734.tar.xz
FFT bei niedrigeren Frequenzen abgeschnitten
-rw-r--r--Physikunit.pas4
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]);