fftw_loop.cpp

Go to the documentation of this file.

 
#include "sigpack.h"
using namespace arma;
 
int main()
{
    int      N = 45191; // Large prime number
    sp::FFTW ss( N, FFTW_ESTIMATE );
    vec      x( N );
    cx_vec   Sxx( N );
 
    vec     P( N );
    clock_t tic;
    x.randn();
    tic = clock();
    Sxx = fft( x );
    cout << "Armadillo FFT Time = " << ( clock() - tic ) / double( CLOCKS_PER_SEC ) << endl;
 
    tic = clock();
    Sxx = fft( x );
    cout << "Armadillo FFT Time = " << ( clock() - tic ) / double( CLOCKS_PER_SEC ) << endl;
 
    for( int n = 1; n < 5; n++ )
    {
        x.randn();
        tic = clock();
        Sxx = ss.fft( x );
        cout << "FFTW Loop[" << n << "]. Time = " << ( clock() - tic ) / double( CLOCKS_PER_SEC ) << endl;
    }
    return 0;
}