Changeset 3384 in Sophya for trunk/SophyaLib/NTools/fftservintf.cc

Timestamp:

Nov 21, 2007, 11:46:19 AM (18 years ago)

Author:

ansari

Message:

Ajout test sur Imag(Nyquist) pour calcul taille tableau reel de sortie pour FFT-multidim complexe->reel, Reza 21/11/2007

File:

• trunk/SophyaLib/NTools/fftservintf.cc (modified) (2 diffs)

trunk/SophyaLib/NTools/fftservintf.cc

@@ -320 +320 @@
   }
   else {  // We have to resize the output array 
+    T thr = ZeroThreshold();   // Seuil pour tester Imag(Nyquist) == 0 
     sa_size_t sz[BASEARRAY_MAXNDIMS];
     if (ndg1 > 1) {
-      sz[0] = 2*in.Size(0)-1;
+      T imnyq = in(in.Size(0)-1,0,0).imag();
+      // Rz+cmc/Nov07 : 
+      // Calcul de la taille SizeX/Sz[0]  paire/impaire en fonction de Imag(Nyquist)
+      sz[0] = ((imnyq < -thr)||(imnyq > thr)) ? 2*in.Size(0)-1 : 2*in.Size(0)-2;
       for(k=1; k<in.NbDimensions(); k++) 
         sz[k] = in.Size(k);
@@ -329 +333 @@
     else {
       for(k=0; k<BASEARRAY_MAXNDIMS; k++)  sz[k] = 1; 
-      T thr = ZeroThreshold();
       sa_size_t n = in.Size(in.MaxSizeKA());
       sa_size_t ncs = ( (in[n-1].imag() < -thr) || (in[n-1].imag() > thr) )