15 using namespace shogun;
16 #define INDEX(ROW,COL,DIM) (((COL)*(DIM))+(ROW))
40 SG_ERROR(
"Number of training vectors does not match number of labels\n");
45 SG_INFO(
"%d trainlabels\n", num_data);
48 for (int32_t i=0; i<num_data; i++)
59 int32_t tmax=1000000000;
68 for(int32_t i=0; i<num_data; i++) {
74 memset(vector_c, 0, num_data*
sizeof(
float64_t));
84 npp.
gnpp_imdm(diagK, vector_c, vector_y, num_data,
85 tmax, tolabs, tolrel, thlb, alpha, &t, &aHa11, &aHa22,
92 for(int32_t i = 0; i < num_data; i++ )
94 if( alpha[i] != 0 ) num_sv++;
97 alpha[i] = alpha[i]*2/nconst;
98 if( alpha[i]/(2*C) >= 1 ) trnerr++;
102 alpha[i] = -alpha[i]*2/nconst;
103 if( alpha[i]/(2*C) <= -1 ) trnerr++;
107 float64_t b = 0.5*(aHa22 - aHa11)/nconst;;
114 for (int32_t i=0; i<num_data; i++)