testwal.cpp

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include "MTStation.h"
#include <boost/bind.hpp>
#include <gsl/gsl_math.h>
using namespace std;

int main(int argc, char *argv[])
{
        string infilename;
        MTStation Data;
        if (argc == 2)
        {
                infilename = argv[1];
        }
        else
        {
                cout << "Infilename: ";
                cin >> infilename;
        }
        
        Data.GetData(infilename);
        
        
        for (size_t i = 0; i < Data.GetMTData().size(); ++i)
        {
                double r = Data.GetMTData().at(i).GetI2()/Data.GetMTData().at(i).GetI1();
                double a = gsl_pow_2(Data.GetMTData().at(i).GetI5() - Data.GetMTData().at(i).GetI6());
                double b = 1.0 - Data.GetMTData().at(i).GetI5() * Data.GetMTData().at(i).GetI6() 
                                + sqrt(1.0 + gsl_pow_2(Data.GetMTData().at(i).GetI5())*gsl_pow_2(Data.GetMTData().at(i).GetI6())
                                                - gsl_pow_2(Data.GetMTData().at(i).GetI5())- gsl_pow_2(Data.GetMTData().at(i).GetI6()));
                //double walkappa = sqrt((1.0+r*r)/(1- a/(2*b) + r*r*b/2)-1.0);
                double walkappa = sqrt((1.0+gsl_pow_2(Data.GetMTData().at(i).Getr()))/
                                        (1.0- Data.GetMTData().at(i).Geta()/(2*Data.GetMTData().at(i).Getb()) 
                                                        + Data.GetMTData().at(i).Getr()*Data.GetMTData().at(i).Getr()*Data.GetMTData().at(i).Getb()/2)-1.0);
                double walsigma = (gsl_pow_2(Data.GetMTData().at(i).GetI3()) + gsl_pow_2(Data.GetMTData().at(i).GetI4()) * r*r)
                                                /(1- a/(2*b) + r*r*b/2);
                cout << 1./Data.GetMTData().at(i).GetFrequency() << " ";
                cout << Data.GetMTData().at(i).GetKappa() << " ";
                cout << walkappa << " ";
                cout << Data.GetMTData().at(i).GetSigma() << " ";
                cout <<walsigma<< " ";
                cout << endl;
        }
        
}

---