werteunit.pas: bugfix fft

1 files changed, 17 insertions, 17 deletions

diff --git a/werteunit.pas b/werteunit.pas
index f81e864..9b0d0ad 100644
--- a/werteunit.pas
+++ b/werteunit.pas
@@ -1054,7 +1054,7 @@ begin
         if not senkrecht and (algo.inOrdnung = doAlleResIms) then
           algo.laden(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.inOrdnung = doAlleResSmi) then      // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.laden(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.laden(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.laden(invers,pSingle(@(werte[i*sStep])),pStep);
         algo.summen(true,vorher,in0);
@@ -1063,7 +1063,7 @@ begin
         if not senkrecht and (algo.outOrdnung = doAlleResIms) then
           algo.speichern(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.outOrdnung = doAlleResSmi) then     // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.speichern(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.speichern(invers,pSingle(@(werte[i*sStep])),pSingle(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.speichern(invers,pSingle(@(werte[i*sStep])),pStep);
       end;
@@ -1072,7 +1072,7 @@ begin
         if not senkrecht and (algo.inOrdnung = doAlleResIms) then
           algo.laden(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.inOrdnung = doAlleResSmi) then      // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.laden(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.laden(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.laden(invers,pDouble(@(werte[i*sStep])),pStep);
         algo.summen(true,vorher,in0);
@@ -1081,7 +1081,7 @@ begin
         if not senkrecht and (algo.outOrdnung = doAlleResIms) then
           algo.speichern(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.outOrdnung = doAlleResSmi) then     // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.speichern(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.speichern(invers,pDouble(@(werte[i*sStep])),pDouble(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.speichern(invers,pDouble(@(werte[i*sStep])),pStep);
       end;
@@ -1090,7 +1090,7 @@ begin
         if not senkrecht and (algo.inOrdnung = doAlleResIms) then
           algo.laden(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.inOrdnung = doAlleResSmi) then      // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.laden(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.laden(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.laden(invers,pExtended(@(werte[i*sStep])),pStep);
         algo.summen(true,vorher,in0);
@@ -1099,7 +1099,7 @@ begin
         if not senkrecht and (algo.outOrdnung = doAlleResIms) then
           algo.speichern(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[imShift   +   i*sStep])),pStep)
         else if not senkrecht and (algo.outOrdnung = doAlleResSmi) then     // "Smi" bezieht sich auf die Dimension, in der nicht transformiert wird!
-          algo.speichern(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[2*imShift-1 - i*sStep])),pStep)
+          algo.speichern(invers,pExtended(@(werte[i*sStep])),pExtended(@(werte[2*imShift-(1+i)*sStep])),pStep)
         else
           algo.speichern(invers,pExtended(@(werte[i*sStep])),pStep);
       end;
@@ -1329,18 +1329,18 @@ begin
           else
             imIm:=0;            // erste und mittlere Zeile und Spalte sind x-t-reell
 
-          werte[i     +        j *2*xS2]:=reRe-imIm;
-          werte[i     + (4*tS4-j)*2*xS2]:=reIm+imRe;
+          werte[i     +          j *2*xS2]:=reRe-imIm;
+          werte[i     + (4*tS4-1-j)*2*xS2]:=reIm+imRe;
           if (i<>0) and (i<>xS2) then begin
-            werte[2*xS2-i +        j *2*xS2]:= reRe+imIm;
-            werte[2*xS2-i + (4*tS4-j)*2*xS2]:=-imRe+reIm;
+            werte[2*xS2-i +          j *2*xS2]:= reRe+imIm;
+            werte[2*xS2-i + (4*tS4-1-j)*2*xS2]:=-imRe+reIm;
           end;
-          if (j<>0) and (j<>2*tS4) then begin
-            werte[i     + (2*tS4-j)*2*xS2]:=reRe+imIm;
-            werte[i     + (2*tS4+j)*2*xS2]:=imRe-reIm;
+          if (j<>0) and (j<>tS4) then begin
+            werte[i     + (2*tS4-1-j)*2*xS2]:=reRe+imIm;
+            werte[i     +   (2*tS4+j)*2*xS2]:=imRe-reIm;
             if (i<>0) and (i<>xS2) then begin
-              werte[2*xS2-i + (2*tS4-j)*2*xS2]:= reRe-imIm;
-              werte[2*xS2-i + (2*tS4+j)*2*xS2]:=-imRe-reIm;
+              werte[2*xS2-i + (2*tS4-1-j)*2*xS2]:= reRe-imIm;
+              werte[2*xS2-i +   (2*tS4+j)*2*xS2]:=-imRe-reIm;
             end;
           end;
         end;
@@ -1357,7 +1357,7 @@ begin
           else
             reIm:=0;            // erste und mittlere Zeile ist t-reell
           if (i<>0) and (i<>xS2) and (j<>0) and (j<>tS4) then
-            imIm:=werte[-i + (2*tS4-j+2)*2*xS2]
+            imIm:=werte[-i + (2*tS4-j+1)*2*xS2]
           else
             imIm:=0;            // erste und mittlere Zeile und Spalte sind x-t-reell
 
@@ -1367,7 +1367,7 @@ begin
             werte[2*xS2-i +        j *2*xS2]:= reRe+imIm;
             werte[2*xS2-i + (2*tS4+j)*2*xS2]:=-imRe+reIm;
           end;
-          if (j<>0) and (j<>2*tS4) then begin
+          if (j<>0) and (j<>tS4) then begin
             werte[i     + (2*tS4-j)*2*xS2]:=reRe+imIm;
             werte[i     + (4*tS4-j)*2*xS2]:=imRe-reIm;
             if (i<>0) and (i<>xS2) then begin