seisspectest.cpp

Go to the documentation of this file.

#include <iostream>
#include <gsl/gsl_math.h>
#include "SeismicDataComp.h"
#include "RecCalc.h"
using namespace std;

void Subroutine(const SeismicDataComp &SubTest)
{
        SeismicDataComp Local(SubTest);
}
int main()
{
        const double tolerance = 1e-6;
        const int runs = 20;
        for (int i = 0; i < runs; ++i)
        {
                RecCalc RecCalc(0,0,4,0.001);
                SeismicDataComp Test1;
                Test1.GetData("smp.rec");
                Subroutine(Test1);
                SeismicDataComp Test2(Test1);
                Subroutine(Test2);
                SeismicDataComp Test3;
                Test3 = Test1;
                Subroutine(Test3);
                
                for (int i = 0; i < Test1.Data.size(); ++i)
                {
                        if (gsl_fcmp(Test1.Data.at(i),Test2.Data.at(i),tolerance) != 0 
                                || gsl_fcmp(Test1.Data.at(i),Test2.Data.at(i),tolerance) != 0)
                                {
                                cout << Test1.Data.at(i) << " " << Test2.Data.at(i) << " "<< Test3.Data.at(i) << " " << endl;
                                }
                }
                RecCalc.CalcRec(Test1,Test2,Test3);
        }
        cout << "Done." << endl;
}

---