#include"ranksvm.h"
#include"ranksvmtn.h"
#include"rankmisc.h"
#include"../tools/matrixIO.h"
#include<iostream>
#include<fstream>
#include<string>

using namespace Eigen;
using namespace std;

double C=1e-4;// Compensating & scaling

// Main terminating criteria
int maxiter = 1; // max iteration count
double prec=1e-15; // precision
// conjugate gradient
double cg_prec=1e-20; // precision
int cg_maxiter = 1000;
int ls_maxiter = 20;
// line search
double ls_prec=1e-20; // precision
double ls_turb=1e-30; // purturbation


int RSVM::saveModel(const string fname){

    std::ofstream fout(fname.c_str());
    fout<<this->getName()<<endl;
    fout<<this->fsize<<endl;
    Eigen::write_stream(fout, this->model.weight);
    fout<<this->model.beta<<endl;
    return 0;
}

RSVM* RSVM::loadModel(const string fname){
    std::ifstream fin(fname.c_str());
    std::string type;
    int fsize;
    fin>>type;
    fin>>fsize;

    RSVM* rsvm;

    if (type=="TN")
        rsvm = new RSVMTN();
    if (type=="L1")
        rsvm = new RSVML1();
    if (type=="BH")
        rsvm = new RSVMBH();

    rsvm->fsize=fsize;
    SVMModel model;
    Eigen::read_stream(fin, model.weight);
    fin>>model.beta;
    rsvm->setModel(model);

    return rsvm;
}

int RSVM::setModel(const SVMModel &model) {
    if (model.weight.rows()!=fsize)
        LOG(FATAL) << "Feature size mismatch: "<<fsize<<" "<<model.weight.rows();
    this->model.weight=model.weight;
    this->model.beta=model.beta;
    return 0;
}