fftw_test_basic.cpp

Go to the documentation of this file.

#include "CTsSpectrum.h"
#include "../Abstract Classes/types.h"
#include <fstream>
#include <iostream>
#include <miscfunc.h>
#include <fstream>
#include <boost/lambda/lambda.hpp>
#include <boost/numeric/ublas/vector.hpp>
using namespace std;

int main ()
{
        CTsSpectrum Spectrum1(true);
        const int points = 200;
        const double PI = 3.14159265359;
        //ttsdata TimeSeries1(points);
        //ttsdata TimeSeries2(points);
        //ttsdata CorrTs(points,0);
        //tcompdata FreqDomain1(TimeSeries1.size()/2+1,0);
        
        ublas::vector<double> TimeSeries1(points);
        ublas::vector<double> TimeSeries2(points);
        ublas::vector<double> CorrTs(points);
        ublas::vector<complex<double> > FreqDomain1(TimeSeries1.size()/2+1);
        const int calculations = 100;
        //TimeSeries1.assign(points,0);
        //Spectrum1.TsData->insert(Spectrum1.TsData->end(),points,1);
        for (int i = 0; i < points; ++i)
                TimeSeries1(i)=sin(4 * PI * i/points);
        
        Spectrum1.CalcSpectrum(TimeSeries1.begin(),TimeSeries1.end(),FreqDomain1.begin(),FreqDomain1.end());
        Convolve(TimeSeries1,TimeSeries1,CorrTs);
        ofstream tsfile("tsold");
        for (int i = 0; i < TimeSeries1.size(); ++i)
                tsfile << TimeSeries1(i) << CorrTs(i) << endl;
        Spectrum1.CalcTimeSeries(FreqDomain1.begin(),FreqDomain1.end(),TimeSeries2.begin(),TimeSeries2.end());
        ofstream newfile("tsnew");
        for (int i = 0; i < TimeSeries2.size(); ++i)
                newfile << TimeSeries2(i) << endl;      
}