% The multiplicative constants for the 17 point FFT

I = sqrt(-1);
u = [
       -1.062500000000000
       -0.257694101601104
        0.723407977286057
       -0.089055591620606
       -0.317176192832725
        0.924380996081242
        0.676798496730885
       -0.440889073481754
       -1.517002366671939
       -0.797601020823318
        1.281092943422807
        0.296310685295348
        0.060401262046216
       -0.420101934970527
        1.462686052158509*I
        2.709842506062867*I
       -1.124438635937869*I
       -1.808356521480244*I
        2.958485673330231*I
        0.222952651355246*I
       -0.906077574510765*I
       -2.491481444635763*I
        0.634492510107882*I
        2.681907643666417*I
       -1.890642422994411*I
        0.499530681019060*I
        0.524082025323152*I
       -1.205277132872840*I
        0.867029716652200*I
        0.032526452324592*I
        1.423638194299944*I
       -1.356975842482187*I
       -2.072296847912463*I
       -0.409600041534228*I
        0.312453977459405*I
        0.642137248078546*I
       -0.876604270228695*I
       -0.544991184003723*I
        0.436775561093086*I
        0.533921625167910*I
        0.361241666187153*I
];


% The input permutation for the 17 point FFT

ip = [
   1
   2
   4
   10
   11
   14
   6
   16
   12
   17
   15
   9
   8
   5
   13
   3
   7
];


% The output permutation for the 17 point FFT

op = [
   1
   17
   15
   2
   13
   6
   16
   12
   11
   3
   4
   8
   14
   5
   10
   7
   9
];