16 using namespace shogun;
17 #define INDEX(ROW,COL,DIM) (((COL)*(DIM))+(ROW))
42 SG_ERROR(
"Number of training vectors does not match number of labels\n")
47 SG_INFO(
"%d trainlabels\n", num_data)
50 for (int32_t i=0; i<num_data; i++)
58 SG_ERROR(
"label unknown (%f)\n", lab)
62 int32_t tmax=1000000000;
71 for(int32_t i=0; i<num_data; i++) {
77 memset(vector_c, 0, num_data*
sizeof(
float64_t));
87 npp.
gnpp_imdm(diagK, vector_c, vector_y, num_data,
88 tmax, tolabs, tolrel, thlb, alpha, &t, &aHa11, &aHa22,
95 for(int32_t i = 0; i < num_data; i++ )
97 if( alpha[i] != 0 ) num_sv++;
100 alpha[i] = alpha[i]*2/nconst;
101 if( alpha[i]/(2*C) >= 1 ) trnerr++;
105 alpha[i] = -alpha[i]*2/nconst;
106 if( alpha[i]/(2*C) <= -1 ) trnerr++;
110 float64_t b = 0.5*(aHa22 - aHa11)/nconst;;
117 for (int32_t i=0; i<num_data; i++)