diff options
| author | Erich Eckner <git@eckner.net> | 2017-11-07 15:28:27 +0100 |
|---|---|---|
| committer | Erich Eckner <git@eckner.net> | 2017-11-07 15:28:27 +0100 |
| commit | 1caef83fc5136f70fce3ce38c68cb1ed5ff579c4 (patch) | |
| tree | 96ca9be7b0a6034229a05883a0915ae9393f174b /werteunit.inc | |
| parent | 6f2f524dc5aab53084ce0b9774e89995584ca3fc (diff) | |
| download | epost-1caef83fc5136f70fce3ce38c68cb1ed5ff579c4.tar.xz | |
werteunit.inc: fix fft ordering
Diffstat (limited to 'werteunit.inc')
| -rw-r--r-- | werteunit.inc | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/werteunit.inc b/werteunit.inc index 8e768ca..306be7d 100644 --- a/werteunit.inc +++ b/werteunit.inc @@ -342,12 +342,12 @@ begin // | // reRe+imIm + | reRe-imIm + - // i*(imRe-reIm) | i*(-imRe-reIm) + // i*(reIm-imRe) | i*(-reIm-imRe) // | // ------------------------------- // | // reRe-imIm + | reRe+imIm + - // i*(reIm+imRe) | i*(-imRe+reIm) + // i*(reIm+imRe) | i*(-reIm+imRe) // | // Mathematica (komplexe-Matrix-Multiplikation.nb) behauptet: @@ -408,12 +408,12 @@ begin // | // reRe+imIm + | reRe-imIm + - // i*(imRe-reIm) | i*(-imRe-reIm) + // i*(reIm-imRe) | i*(-reIm-imRe) // | // ------------------------------- // | // reRe-imIm + | reRe+imIm + - // i*(reIm+imRe) | i*(-imRe+reIm) + // i*(reIm+imRe) | i*(-reIm+imRe) // | // Mathematica (komplexe-Matrix-Multiplikation.nb) behauptet: @@ -424,11 +424,11 @@ begin werte[i+j*params.xSteps]:=reRe1*reRe2 + imIm1*imIm2 - imRe1*imRe2 - reIm1*reIm2; if (i<>0) and (i<>x2) then begin - werte[-i+2*x2+j*params.xSteps]:=-reIm1*imIm2 - imIm1*reIm2 + reRe1*imRe2 + imRe1*reRe2; + werte[-i+2*x2+j*params.xSteps]:=-imRe1*imIm2 - imIm1*imRe2 + reRe1*reIm2 + reIm1*reRe2; o0:=o0 and (werte[-i+2*x2+j*params.xSteps]=0); end; if (j<>0) and (j<>t2) then begin - werte[i+(-j+2*t2)*params.xSteps]:=-imRe1*imIm2 - imIm1*imRe2 + reRe1*reIm2 + reIm1*reRe2; + werte[i+(-j+2*t2)*params.xSteps]:=-reIm1*imIm2 - imIm1*reIm2 + reRe1*imRe2 + imRe1*reRe2; o0:=o0 and (werte[i+(-j+2*t2)*params.xSteps]=0); if (i<>0) and (i<>x2) then begin werte[-i+2*x2+(-j+2*t2)*params.xSteps]:=reIm1*imRe2 + imRe1*reIm2 + reRe1*imIm2 + imIm1*reRe2; |