GPLIB++: testanisomtobjective.cpp Source File

00001 #include "Aniso1DMTObjective.h"
00002 #include <iostream>
00003 #include <string>
00004 #include <boost/numeric/ublas/io.hpp>
00005 
00006 int main()
00007 {
00008         MTStation Stat1, Stat2;
00009         
00010         Stat1.GetData("test1.mtt");
00011         
00012         Aniso1DMTObjective::datafuncvector_t DataFunc, ErrFunc;
00013         std::vector<double> Errlevel;
00014         ttranscribed ModelVector(8);
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         Aniso1DMTObjective 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         ModelVector(4) = 0.0;
00028         ModelVector(5) = 0.0;
00029         ModelVector(6) = 0.0;
00030         ModelVector(7) = 0.0;
00031         Objective.CalcPerformance(ModelVector);
00032         std::cout << Objective.GetSynthData() << std::endl;
00033         std::cout << Objective.GetMisfit() << std::endl;
00034 }