00001 #include "PsDoublModel.h"
00002 #include <iostream>
00003 #include <fstream>
00004 #include <string>
00005
00006 using namespace std;
00007
00008 int main(void)
00009 {
00010 PsDoublModel Model;
00011 string inputfilename, outputfilename;
00012 ofstream outfile("psdoubl.conf");
00013 cout << "Inputfile: ";
00014 cin >> inputfilename;
00015
00016 Model.GetData(inputfilename);
00017 const int size = Model.Thickness.size();
00018 const double variation = 0.1;
00019 const int genesize = 5;
00020 double stepsize = 0;
00021 const int steps = pow(2.,genesize);
00022 vector<double> values;
00023
00024 for (int i = 0; i < size; ++i)
00025 {
00026 if (Model.Thickness.at(i) < 0)
00027 values.push_back(Model.Thickness.at(i));
00028 if (Model.PVelocity.at(i) < 0)
00029 values.push_back(Model.PVelocity.at(i));
00030 if (Model.SVelocity.at(i) < 0)
00031 values.push_back(Model.SVelocity.at(i));
00032 if (Model.Density.at(i) < 0)
00033 values.push_back(Model.Density.at(i));
00034 if (Model.Qp.at(i) < 0)
00035 values.push_back(Model.Qp.at(i));
00036 if (Model.Qs.at(i) < 0)
00037 values.push_back(Model.Qs.at(i));
00038 }
00039 for (int i = 0; i < values.size(); ++i)
00040 outfile << abs(values.at(i) - values.at(i)*variation )<< " ";
00041 outfile << endl;
00042 for (int i = 0; i < values.size(); ++i)
00043 outfile << genesize << " ";
00044 outfile << endl;
00045 for (int i = 0; i < values.size(); ++i)
00046 outfile << abs(2* values.at(i) * variation/steps) << " ";
00047
00048 }
