testsurfobjective.cpp

Go to the documentation of this file.

#include <iostream>
#include "SurfaceWaveObjective.h"
#include "SurfaceWaveModel.h"
#include "SurfaceWaveData.h"
#include <string>

using namespace std;
int main()
{
        string datafilename, modelfilename, backgroundname;
        cout << "Name of Datafile: ";
        cin >> datafilename;
        cout << "Name of Modelfile:";
        cin >> modelfilename;
        cout << "Name of Backgroundmodel: ";
        cin >> backgroundname;
        
        SurfaceWaveData Data;
        Data.ReadSurf96(datafilename);
        SurfaceWaveModel Model;

        SurfaceWaveObjective SeisObjective(Data);
        
        SeisObjective.SetFitExponent(2);
        SeisObjective.SetPoisson(sqrt(3));
        SeisObjective.SetErrorLevel(0.02);
        const unsigned int nlayers = Model.GetPvVelocities().size();
        const unsigned int nparams = nlayers * 2;
        ttranscribed member(nparams);
        for (unsigned i = 0; i < nlayers; ++i)
        {
                member(i) = Model.GetThicknesses().at(i);
                member(i+nlayers) = Model.GetPvVelocities().at(i);
        }
        cout << "Misfit: " << SeisObjective.CalcPerformance(member) << endl;
}

---