testkhl.cpp

Go to the documentation of this file.

#include "KHL.h"
#include "statutils.h"
#include "UniformRNG.h"
#include "VecMat.h"
#include <iostream>
#include <fstream>
#include <algorithm>
#include <iterator>
#include <string>

using namespace std;
int main()
{
        UniformRNG Random;
        
        
        string name;
        cout << "Infile: ";
        cin >> name;
        ifstream file(name.c_str());
        vector<double> in((istream_iterator<double>(file)), istream_iterator<double>());
        
        
        const size_t nobs = in.size();
        rvec input(nobs), eval(nobs);
        rmat evec(nobs,nobs);
        copy(in.begin(),in.end(),input.begin());
        
   
        KHLExpand(input,evec,eval);
        int index = 0;
        double biggest = eval(0);
        for (int i = 0; i < eval.size(); ++i)
                if (eval(i) > biggest)
                {
                        index = i;
                        biggest = eval(i);
                }
        
        std::ofstream khloutfile("khl.out");
        copy(column(evec,index).begin(),column(evec,index).end(),std::ostream_iterator<double>(khloutfile,"\n"));
        std::ofstream khlinfile("khl.in");
        copy(input.begin(),input.end(),std::ostream_iterator<double>(khlinfile,"\n"));
        sort(eval.begin(),eval.end());
        std::ofstream khlevalfile("khl.eval");
        copy(eval.begin(),eval.end(),std::ostream_iterator<double>(khlevalfile,"\n"));
        //std::cout << "Input: " << input << std::endl;
        //std::cout << "Evec: " << evec << std::endl;
        //std::cout << "Eval: " << eval << std::endl;
        //std::cout << "Rec: " << prod(evec,eval) << std::endl;
}

---