spectest.cpp

Go to the documentation of this file.

#include "TsSpectrum.h"
#include "UniformRNG.h"
#include "types.h"
#include <iostream>
#include <gsl/gsl_math.h>

using namespace std;
int main()
{
        const double tolerance = 1e-6;
        const bool MultiCalc = false;
        const int length = 20000;
        const int runs = 100;
        ttsdata Input(length), Recovered(length);
        tcompdata Spectrum(length/2+1);
        UniformRNG Random;
        for (int r = 0; r < runs; ++r)
        {
                cout << "Run: " << r << "       ";
                bool runok = true;
                for (int i = 0; i < length; ++i)
                        Input.at(i) = Random.GetNumber();
                TsSpectrum SpecCalc(MultiCalc);
                SpecCalc.CalcSpectrum(Input.begin(),Input.end(),Spectrum.begin(),Spectrum.end());
                SpecCalc.CalcTimeSeries(Spectrum.begin(),Spectrum.end(),Recovered.begin(),Recovered.end());
                for (int i = 0; i < length; ++i)
                {
                        if (gsl_fcmp(Input.at(i),Recovered.at(i),tolerance) != 0)
                        {
                                cout << endl;
                                cout << i << " " << Input.at(i) << " " << Recovered.at(i) << endl;
                                runok = false;
                        }
                }
                if (runok)
                        cout << "Run OK !" << endl;
                else
                        cout << "Run failed !" << endl;
        }
        cout << "All finished !" << endl;
}

---