source: trunk/test/svm_test.cc @ 307

// $Id: svm_test.cc 307 2005-05-03 13:28:29Z peter $

#include <c++_tools/gslapi/matrix.h>
#include <c++_tools/gslapi/vector.h>
#include <c++_tools/svm/SVM.h>
#include <c++_tools/svm/Kernel_SEV.h>
#include <c++_tools/svm/PolynomialKernelFunction.h>

#include <fstream>
#include <iostream>
#include <cstdlib>
#include <limits.h>

using namespace theplu;

int main()

{  

//  ifstream is("data/nm_kernel.txt");
//  theplu::gslapi::matrix kernel_matrix(is);
//  is.close();
 
//  is.open("data/nm_target_bin.txt");
//  theplu::gslapi::vector target(is);
//  is.close();

//  is.open("data/nm_alpha_linear_matlab.txt");
//  theplu::gslapi::vector alpha_matlab(is);
//  is.close();

//  theplu::cpptools::SVM svm(kernel_matrix, target);
//    svm.train();

//  theplu::gslapi::vector alpha = svm.get_alpha();
      
//  // Comparing alpha to alpha_matlab
//  theplu::gslapi::vector diff_alpha = alpha - alpha_matlab;
//  if (diff_alpha*diff_alpha> pow(10.0,-10)){
//    cerr << "Difference to matlab alphas too large/n";
//    return -1;
//  }
//  // Comparing output to target
//  theplu::gslapi::vector output = svm.get_output();
//  double slack = 0;
//  for (unsigned int i=0; i<target.size(); i++){
//    if (output[i]*target[i] < 1){
//      slack += 1 - output[i]*target[i];
//    }
//  }
//  if (slack > pow(10.0,-10)){
//    cerr << "Difference to matlab alphas too large/n";
//    return -1;
//  }
  
  // testing on non solvable problem
  gslapi::matrix data(3,1);
  data(0,0)=-1;
  data(1,0)=10;
  data(2,0)=0;
    
  gslapi::vector target(3,1.0);
  target(2)=-1;
  svm::KernelFunction* kf =
    new svm::PolynomialKernelFunction(); 
  svm::Kernel_SEV kernel(data,*kf); 
  //theplu::gslapi::matrix m=kernel.get();
  //std::cout << "kernel:\n" << m << "\n";
  svm::SVM svm(kernel,target);
  //  svm.set_c(1);
  //svm.train();
  
  return 0;
  
}