Seismic Tools/testrecobjective.cpp File Reference

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

Include dependency graph for Seismic Tools/testrecobjective.cpp:

Go to the source code of this file.

Functions
int	main ()

---

Function Documentation

int main

(

void

)

Definition at line 10 of file Seismic Tools/testrecobjective.cpp.

References GeneralObjective::CalcPerformance(), CCalcRecConf::cc, ResPkModel::GetData(), SeismicDataComp::GetData(), CCalcRecConf::GetData(), RecCalc::iterdecon, CCalcRecConf::omega, SurfaceWaveData::ReadAscii(), RecData, AbsVelRecObjective::SetAbsVelWeight(), C1DRecObjective::SetErrorLevel(), GeneralObjective::SetFitExponent(), C1DRecObjective::SetPoisson(), AbsVelRecObjective::SetRecWeight(), CCalcRecConf::shift, CCalcRecConf::sigma, SeismicModel::SVelocity, and SeismicModel::Thickness.

{
        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;
}

Here is the call graph for this function:

---