int l = model->l;
		int i;

                double *kvalue = Malloc(double,l);
                for(i=0;i<l;i++)
                        kvalue[i] = Kernel::k_function(x,model->SV[i],model->param);

                int *start = Malloc(int,nr_class);
                start[0] = 0;
                for(i=1;i<nr_class;i++)
                        start[i] = start[i-1]+model->nSV[i-1];


                i = 0;
                int j = nr_class - 1;
                while (i < j)
                {
                        double sum = 0;
                        int si = start[i];
                        int sj = start[j];
                        int ci = model->nSV[i];
                        int cj = model->nSV[j];

                        int k;
                        double *coef1 = model->sv_coef[j-1];
                        double *coef2 = model->sv_coef[i];
                        for(k=0;k<ci;k++)
                                sum += coef1[si+k] * kvalue[si+k];
                        for(k=0;k<cj;k++)
                                sum += coef2[sj+k] * kvalue[sj+k];
                        sum -= model->rho[i*(nr_class - 1) -
                                i*(i + 1)/2 + j - 1];
                        if(sum > 0)
                                j--;
                        else
                                i++;
                }

                free(kvalue);
                free(start);
                return model->label[i];