         while( j == i ) 
           j = ceil( (float) trainingSize * 
                       float( rand() ) / ( float(RAND_MAX) + 1.0f ) );
-        printf("j is %d\n",j);
+//        j = 3918;
+ //       printf("j is %d\n",j);
 
         for( int k = 0; k < trainingSize; k++ )
         {  

          continue;
        } /* end if */

        float b1 = b - E[i]
                     - y[i] * (alphas[i] - alphaIOld) * 
                            K[INDX(i,j,trainingSize)]
                     - y[j] * (alphas[j] - alphaJOld) * 
                            K[INDX(i,j,trainingSize)];


      } /* end if */

      exit(911);
    } /* end for i */ 
  } /* end while */

  free(E); 
  free(alphas);


        float alphaIOld = alphas[i];
        float alphaJOld = alphas[j];

//        printf("alphaIOld %f alphaJOld %f\n",alphaIOld,alphaJOld);
 //       printf("yi %f yj %f\n",y[i],y[j]);

        if( y[i] == y[j] )
        {
