cbitrev.cpp

来自「Digital filter designer s handbook C++ c」· C++ 代码 · 共 36 行

 //
 //  File = cbitrev.cpp
 //

 #include "d_cmplx.h"
 #include "cbitrev.h"

 void ComplexBitReverse( double_complex *array, int size)
 {
  double_complex tt;
  int nv2, nm1, i, j, k;
  
  nv2 = size/2;
  nm1 = size - 1;
  
  j=0;
  for (i=0; i<nm1; i++)
    {
     if (i<j)
       {
        tt = array[j];
        array[j] = array[i];
        array[i] = tt;
       }
     k = nv2;
     while (k<=j)
       {
        j -= k;
        k /= 2;
       }
     j += k;
    }
 }
 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++