Method update

JSAT/src/jsat/classifiers/linear/kernelized/DUOL.java:203–309

(DataPoint dataPoint, int targetClass)

    }

    @Override
    public synchronized void update(DataPoint dataPoint, int targetClass)
    {
        final Vec x_t = dataPoint.getNumericalValues();
        final double y_t = targetClass*2-1;
        final List<Double> qi = k.getQueryInfo(x_t);
        double score = score(x_t, qi, true);

        final double loss_t = max(0, 1-y_t*score);
        
        if(loss_t <= 0)
            return;

        //start of line 8:
        int b = -1;
        double w_min = Double.POSITIVE_INFINITY;
        for(int i = 0; i < S.size(); i++)
        {
            if(f_s.get(i) <= 1)
            {
                double tmp = signum(alphas.get(i))*y_t*kTmp.get(i);
                if(tmp <= w_min)
                {
                    w_min = tmp;
                    b = i;
                }
            }
        }
        
        final double k_t = k.eval(0, 0, Arrays.asList(x_t), qi);
        if(w_min <= -rho)
        {
            
            final double k_b = k.eval(b, b, S, accelCache);
            final double k_tb = kTmp.get(b);
            final double alpha_b = alphas.get(b);
            final double w_tb = y_t*signum(alpha_b)*k_tb;
            final double gamma_hat_b = abs(alpha_b);
            final double loss_b = 1-signum(alpha_b)*f_s.get(b);
            
            //(C-gamma_hat_b) is common expression bellow, so use this insted
            final double CmGhb = (C-gamma_hat_b);

            final double gamma_t;
            final double gamma_b;
            final double gamma_b_delta;
            
            if(k_t*C +w_tb*CmGhb-loss_t < 0 && k_b*CmGhb+w_tb*C-loss_b < 0)
            {
                gamma_t = C;
                gamma_b_delta = CmGhb;
            }
            else if(  (w_tb*w_tb*C-w_tb*loss_b-k_t*k_b*C+k_b*loss_t)/k_b > 0 && isIn((loss_b-w_tb*C)/k_b, -gamma_hat_b, CmGhb) )
            {
                gamma_t = C;
                gamma_b_delta = (loss_b-w_tb*C)/k_b;
            }
            else if(  isIn((loss_t-w_tb*CmGhb)/k_t, 0, C) && loss_b-k_b*CmGhb-w_tb*(loss_t-w_tb*CmGhb)/k_t > 0  )