epostunit.pas: tWerte.berechneRadonTransformation(): accellerate fft by using 2^x as edge length

1 files changed, 3 insertions, 2 deletions

diff --git a/epostunit.pas b/epostunit.pas
index 0f28c94..68a0308 100644
--- a/epostunit.pas
+++ b/epostunit.pas
@@ -6125,8 +6125,9 @@ begin
   if not sT then begin
     gibAus('kopiere Werte für Radon-Transformation ...',3);
     hilfsWerte:=tWerte.create(quelle,gExtended); // quelle zu hilfsWerte kopieren ...
-    hilfsWerte._xSteps:=2*_tSiz+1;               // ... und quadratisch ergänzen: Kantenlänge = 2*Diagonale von quelle + 1
-    hilfsWerte._tSiz:=2*_tSiz+1;
+                                                 // ... und quadratisch ergänzen: Kantenlänge >= 2*Diagonale von quelle + 1
+    hilfsWerte._xSteps:=round(power(2,ceil(ln(_tSiz+1)/ln(2))));
+    hilfsWerte._tSiz:=hilfsWerte._xSteps;
     hilfsWerte.holeRAM(3);
     hilfsWerte.eWerte.nullenEinfuegen(quelle._xSteps,quelle._tSiz);
     gibAus('berechne t-FFT für Radon-Transformation ...',3);