source: trunk/test/test_svm.cc @ 73

// $Id: test_svm.cc 73 2004-05-03 10:13:24Z peter $

// C++ tools include
////////////////////
#include "matrix.h"
#include "vector.h"
#include "SVM.h"
#include "Kernel.h"
#include "PolynomialKernelFunction.h"

// Standard includes
////////////////////
#include <fstream>
#include <iostream>
#include <cstdlib>

using namespace std;

int main()

{  
  cout << "reading data/nm_kernel.txt" << endl;
  ifstream is("data/nm_kernel.txt");
  theplu::gslapi::matrix kernel(is);
  is.close();

  theplu::gslapi::matrix k(kernel);
 
  cout << "reading data/nm_target_bin.txt" << endl;
  is.open("data/nm_target_bin.txt");
  theplu::gslapi::vector target(is);
  is.close();

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


  cout << "create svm" << endl;
  theplu::cpptools::SVM svm(kernel, target);
  cout << "train" << endl;
  svm.train();
  cout << "get alpha" << endl;
  theplu::gslapi::vector alpha = svm.get_alpha();
      
  // Comparing alpha to alpha_matlab
  theplu::gslapi::vector diff_alpha = alpha - alpha_matlab;
  cout << diff_alpha*diff_alpha << endl;
  
  // 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];
    }
  }
  cout << "slack:" << slack << endl;

  return 0;
}