Seismic Tools/testrecobjective.cpp

Go to the documentation of this file.

#include <iostream>
#include "AbsVelRecObjective.h"
#include "SurfaceWaveData.h"
#include "ResPkModel.h"
#include "SeismicDataComp.h"
#include "CCalcRecConf.h"
#include <string>

using namespace std;
int main()
{
        string datafilename, modelfilename;
        cout << "Name of Datafile: ";
        cin >> datafilename;
        cout << "Name of Modelfile:";
        cin >> modelfilename;
        
        CCalcRecConf Config;
        Config.GetData("calcrec.conf");
        SeismicDataComp RecData;
        SurfaceWaveData SurfData;
        RecData.GetData(datafilename+".rec");
        SurfData.ReadAscii(datafilename+".vel");
        ResPkModel Model;
        Model.GetData(modelfilename);
        
        AbsVelRecObjective SeisObjective(RecData,SurfData,Config.shift,Config.omega,Config.sigma,Config.cc,0.07,RecCalc::iterdecon);
        
        SeisObjective.SetFitExponent(2);
        SeisObjective.SetPoisson(sqrt(3));
        SeisObjective.SetErrorLevel(0.02);
        SeisObjective.SetAbsVelWeight(1);
        SeisObjective.SetRecWeight(1);
        const unsigned int nlayers = Model.Thickness.size();
        const unsigned int nparams = nlayers * 2;
        ttranscribed member(nparams);
        for (unsigned i = 0; i < nlayers; ++i)
        {
                member(i) = Model.Thickness.at(i);
                member(i+nlayers) = Model.SVelocity.at(i);
        }
        cout << "Misfit: " << SeisObjective.CalcPerformance(member) << endl;
}

---