testcorr.cpp

Go to the documentation of this file.

#include <iostream>
#include <UniformRNG.h>
#include "miscfunc.h"
#include <numeric>      
#include <algorithm>
#include <boost/bind.hpp>
using namespace std;
int main()
{
        complex<double> C(1.,1.);
        UniformRNG Random;
        const int length = 100; 
        const int shift = 10;
        vector<double> TimeSeries1(3*length,0);
        vector<double> TimeSeries2(3*length,0);
        vector<double> Output;
        
        for (int i = length; i < 2*length; ++i)
        {
                TimeSeries1.at(i) = 200*(0.5- Random.GetNumber());
        }
        for (int i = length+shift; i < 2*length; ++i)
        {
                TimeSeries2.at(i) = TimeSeries1.at(i-shift);
        }
        double mean1=accumulate(TimeSeries1.begin(),TimeSeries1.end(),0.)/length;
        double mean2=accumulate(TimeSeries2.begin(),TimeSeries2.end(),0.)/length;
        //cout << "Mean1: " << mean1 << " Mean2: " << mean2 << endl;
        //transform(TimeSeries2.begin(),TimeSeries2.end(),TimeSeries2.begin(),bind2nd(minus<double>(),mean2));
        //transform(TimeSeries1.begin(),TimeSeries1.end(),TimeSeries1.begin(),bind2nd(minus<double>(),mean1));
        Correl(TimeSeries2,TimeSeries1,Output);
        //cout << "C: " << C << " C*: " << boost::bind<complex<double> >(&conj<double>,_1)(C) << endl;
        for (int i = 0; i < Output.size(); ++i)
            cout << i << " "<<  TimeSeries1.at(i) << " " << TimeSeries2.at(i) << " " <<Output.at(i)/double(length) << " " << endl;
        cout << endl;
}

---