GPLIB++: testmtobjective.cpp Source File

00001 #include "C1DMTObjective.h"
00002 #include <iostream>
00003 #include <string>
00004 #include <boost/numeric/ublas/io.hpp>
00005 
00006 int main()
00007 {
00008         CMTStation Stat1, Stat2;
00009         
00010         Stat1.GetData("test1.mtt");
00011         
00012         Iso1DMTObjective::datafuncvector_t DataFunc, ErrFunc;
00013         std::vector<double> Errlevel;
00014         ttranscribed ModelVector(4);
00015         DataFunc.push_back(&MTTensor::GetRhoxy);
00016         ErrFunc.push_back(&MTTensor::GetdRhoxy);
00017         Errlevel.push_back(0.01);
00018         DataFunc.push_back(&MTTensor::GetPhixy);
00019         ErrFunc.push_back(&MTTensor::GetdPhixy);
00020         Errlevel.push_back(0.01);
00021         Iso1DMTObjective Objective(Stat1);
00022         //Objective.SetFitParameters(DataFunc,ErrFunc,Errlevel); 
00023         ModelVector(0) = 1.0;
00024         ModelVector(1) = 2.0;
00025         ModelVector(2) = 10.0;
00026         ModelVector(3) = 1.0;
00027         Objective.CalcPerformance(ModelVector);
00028         std::cout << Objective.GetSynthData() << std::endl;
00029         std::cout << Objective.GetMisfit() << std::endl;
00030 }