/* 
 * DSPIFFT_RT - DSP Blockset 1-D IFFT 
 * Algorithm: y = (conj(FFT(conj(y)))/n
 *
 *  Copyright (c) 1995-2000 The MathWorks, Inc.
 *  $Revision: 1.5 $  $Date: 2000/03/03 21:28:05 $
 */

#include "dspifft_rt.h"
#include "dspfft_rt.h"

void dspifft(const int_T n, creal_T *y)
{
    int_T i;

    /* Conjugate input */
    for (i=0; i<n; i++) {
#if 1
    /* Bug fix for SGI64 compiler */
      const real_T t = -(y->im);
      (y++)->im = t;
#else
      (y++)->im = -(y->im);
#endif
    }
    y -= n;

    dspfft(n, y);

    {
        const real_T nrecip = 1.0 / n;
        for (i=0; i<n; i++) {
            y->re     *=  nrecip;
            (y++)->im *= -nrecip;
        }
    }
}

/* [EOF] dspifft_rt.c */