Also use between-cluster-distance prune for points.

2015-03-12 15:42:45 -04:00
parent a0a6bc1479
commit a09c3e3314
1 changed files with 10 additions and 0 deletions
@@ -508,6 +508,16 @@ prunedPoints[index] << ", lastOwner " << lastOwners[index] << ": invalid "
        prunedCentroids.col(owner) += dataset.col(index);
        prunedCounts(owner)++;
      }
+      else if (upperBounds[index] + clusterDistances[owner] < 0.5 *
+               interclusterDistances[newFromOldCentroids[owner]])
+      {
+        prunedPoints[index] = true;
+        upperBounds[index] += clusterDistances[owner];
+        lastOwners[index] = owner;
+        lowerSecondBounds[index] -= clusterDistances[centroids.n_cols];
+        prunedCentroids.col(owner) += dataset.col(index);
+        prunedCounts(owner)++;
+      }
      else
      {
        prunedPoints[index] = false;