testmodeldiff.cpp

Go to the documentation of this file.

#include <iostream>
#include "SeismicModelDiff.cpp"
#include "gentypes.h"

using namespace std;
int main()
{
        ResPkModel VRefModel,CompMod;
        string refmodname,compmodname;
        cout << "Reference Model: ";
        cin>> refmodname;
        VRefModel.GetData(refmodname);
        SeismicModelDiff Diff(VRefModel);
        cout << "Compare Model: ";
        cin >> compmodname;
        CompMod.GetData(compmodname);
        const unsigned int nlayers = CompMod.GetSVelocity().size();
        ttranscribed member(nlayers*3);
        
        for     (unsigned int i =0 ; i < nlayers; ++i)
        {
                member(i+nlayers) = CompMod.GetThickness().at(i);
                member(i+2*nlayers) = CompMod.GetSVelocity().at(i);
        }
        cout << "Roughness: " << Diff.CalcPerformance(member) << endl;
        
}

---