36 for (c = 0; c < C->nclasses; c++) {
38 for (band = 0; band < C->nbands; band++)
39 C->sumdiff[band][c] = 0;
44 for (p = 0; p < C->npoints; p++) {
52 for (c = 0; c < C->nclasses; c++) {
57 for (band = 0; band < C->nbands; band++) {
58 z = C->points[band][p] * np - C->sum[band][c];
63 if (first || (d <
min)) {
70 if (C->class[p] !=
class) {
75 C->countdiff[
class]++;
78 for (band = 0; band < C->nbands; band++) {
79 q = C->points[band][p];
80 C->sumdiff[band][
class] += q;
81 C->sumdiff[band][old] -= q;
87 for (c = 0; c < C->nclasses; c++) {
88 C->count[c] += C->countdiff[c];
89 for (band = 0; band < C->nbands; band++)
90 C->sum[band][c] += C->sumdiff[band][c];