#include "kissfft.hh"

#include <iostream></iostream>

#include <cstdlib></cstdlib>

#include <typeinfo></typeinfo>

#include <sys time.h=""></sys>

static inline

double curtime(void)

{

    struct timeval tv;

    gettimeofday(&tv, NULL);

    return (double)tv.tv_sec + (double)tv.tv_usec*.000001;

}

using namespace std;

template <class t=""></class>

void dotest(int nfft)

{

    typedef kissfft<t> FFT;</t>

    typedef std::complex<t> cpx_type;</t>

    cout << "type:" << typeid(T).name() << " nfft:" << nfft;

    FFT fft(nfft,false);

    vector<cpx_type> inbuf(nfft);</cpx_type>

    vector<cpx_type> outbuf(nfft);</cpx_type>

    for (int k=0;k

        inbuf[k]= cpx_type( 

                (T)(rand()/(double)RAND_MAX - .5),

                (T)(rand()/(double)RAND_MAX - .5) );

    fft.transform( &inbuf[0] , &outbuf[0] );

    long double totalpower=0;

    long double difpower=0;

    for (int k0=0;k0

        complex<long double=""> acc = 0;</long>

        long double phinc = 2*k0* M_PIl / nfft;

        for (int k1=0;k1

            complex<long double=""> x(inbuf[k1].real(),inbuf[k1].imag()); </long>

            acc += x * exp( complex<long double="">(0,-k1*phinc) );</long>

        }

        totalpower += norm(acc);

        complex<long double=""> x(outbuf[k0].real(),outbuf[k0].imag()); </long>

        complex<long double=""> dif = acc - x;</long>

        difpower += norm(dif);

    }

    cout << " RMSE:" << sqrt(difpower/totalpower) << "\t";

    double t0 = curtime();

    int nits=20e6/nfft;

    for (int k=0;k

        fft.transform( &inbuf[0] , &outbuf[0] );

    }

    double t1 = curtime();

    cout << " MSPS:" << ( (nits*nfft)*1e-6/ (t1-t0) ) << endl;

}

int main(int argc,char ** argv)

{

    if (argc>1) {

        for (int k=1;k

            int nfft = atoi(argv[k]);

            dotest<float>(nfft); dotest<double>(nfft); dotest<long double="">(nfft);</long></double></float>

        }

    }else{

        dotest<float>(32); dotest<double>(32); dotest<long double="">(32);</long></double></float>

        dotest<float>(1024); dotest<double>(1024); dotest<long double="">(1024);</long></double></float>

        dotest<float>(840); dotest<double>(840); dotest<long double="">(840);</long></double></float>

    }

    return 0;

}