Changeset 209

Timestamp:

Nov 3, 2004, 9:24:41 AM (19 years ago)

Author:

Peter

Message:

update of namespaces

File:

• trunk/test/test_svm.cc (modified) (2 diffs)

trunk/test/test_svm.cc

@@ -14 +14 @@
 #include <iostream>
 #include <cstdlib>
+#include <limits.h>
 
 using namespace std;
@@ -20 +21 @@
 
 {  
-  ifstream is("data/nm_kernel.txt");
-  theplu::gslapi::matrix kernel_matrix(is);
-  is.close();
+  std::cout << ULLONG_MAX << endl;
+
+//  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_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();
+//  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::cpptools::SVM svm(kernel_matrix, target);
+//    svm.train();
 
-  theplu::gslapi::vector alpha = svm.get_alpha();
+//  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 XOR problem
-  cout << "testing XOR\n";
-  theplu::gslapi::matrix data(4,2);
-  data(0,0)=1;
-  data(0,1)=1;
+//  // 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;
+//  }
   
-  data(1,0)=1;
-  data(1,1)=-1;
-  
-  data(2,0)=-1;
-  data(2,1)=1;
-  
-  data(3,0)=-1;
-  data(3,1)=-1;
-  
-  target = theplu::gslapi::vector(4,1.0);
-  target(1)=-1; 
+  // testing on non solvable problem
+  cout << "testing \n";
+  theplu::gslapi::matrix data(3,1);
+  data(0,0)=-1;
+  data(1,0)=10;
+  data(2,0)=0;
+    
+  theplu::gslapi::vector target(3,1.0);
   target(2)=-1;
-  
+  std::cout << "target:\n" << target << "\n";
   theplu::cpptools::KernelFunction* kf =
     new theplu::cpptools::PolynomialKernelFunction(); 
   theplu::cpptools::Kernel kernel(data,*kf); 
-  svm = theplu::cpptools::SVM(kernel.get(),target);
-  //  svm.train();
-  cout << kernel.get() << "\n";
+  theplu::gslapi::matrix m=kernel.get();
+  std::cout << "kernel:\n" << m << "\n";
+  theplu::cpptools::SVM svm(kernel.get(),target);
+  //  svm.set_c(1);
+  svm.train();
+  
   return 0;