int row = nb_points + (k - 1) * nb_classes + c;
scalar_t tau = nb_samples_per_class[c-1] / scalar_t(_nb_clusters);
glp_set_row_bnds(lp, row, GLP_FX, tau, tau);
int row = nb_points + (k - 1) * nb_classes + c;
scalar_t tau = nb_samples_per_class[c-1] / scalar_t(_nb_clusters);
glp_set_row_bnds(lp, row, GLP_FX, tau, tau);