#pragma omp parallel for private(x, j, tid) schedule(static) reduction(+:cluster0IndexCounter)

for (x=0; x<=698; x++)
{

   for (j=0; j<9; j++) //populate cluster0 with points
    {
       cluster0[cluster0IndexCounter][j]=matxOriginal[x][j];

    }

    cluster0IndexCounter=cluster0IndexCounter+1;

}


printf("cluster0: %f, %f, %f, %f, %fn", cluster0[9][0], cluster0[9][1], cluster0[9][2], cluster0[9][3], cluster0[9][4]);