diff --git a/Filters/Points/Testing/Python/TestExtractPoints.py b/Filters/Points/Testing/Python/TestExtractPoints.py
index 37e8ebf153b742425b0d57b37dbb2e61f7086106..21300dde802c03525be6a4cf9ca9c3fab988e5f2 100755
--- a/Filters/Points/Testing/Python/TestExtractPoints.py
+++ b/Filters/Points/Testing/Python/TestExtractPoints.py
@@ -36,8 +36,8 @@ extract.Update()
 timer.StopTimer()
 time = timer.GetElapsedTime()
 print("Time to remove points: {0}".format(time))
-print("   Number removed: {0}".format(extract.GetNumberOfPointsRemoved()),
-      " (out of: {}".format(NPts))
+print("   Number removed: {0}".format(extract.GetNumberOfPointsRemoved()))
+print("   Original number of points: {0}".format(NPts))
 
 # First output are the non-outliers
 extMapper = vtk.vtkPointGaussianMapper()
diff --git a/Filters/Points/Testing/Python/TestFitImplicitFunction.py b/Filters/Points/Testing/Python/TestFitImplicitFunction.py
index 06cd9d7d96475147da067b68d73c2934a646ef7b..97e0e400f84e8fa044eefdeec00e23ed0518df72 100755
--- a/Filters/Points/Testing/Python/TestFitImplicitFunction.py
+++ b/Filters/Points/Testing/Python/TestFitImplicitFunction.py
@@ -37,8 +37,8 @@ extract.Update()
 timer.StopTimer()
 time = timer.GetElapsedTime()
 print("Time to extract points: {0}".format(time))
-print("   Number removed: {0}".format(extract.GetNumberOfPointsRemoved()),
-      " (out of: {}".format(NPts))
+print("   Number removed: {0}".format(extract.GetNumberOfPointsRemoved()))
+print("   Original number of points: {0}".format(NPts))
 
 # First output are the non-outliers
 extMapper = vtk.vtkPointGaussianMapper()
diff --git a/Filters/Points/Testing/Python/TestRadiusOutlierRemoval.py b/Filters/Points/Testing/Python/TestRadiusOutlierRemoval.py
index 100af2c3838cd20ac06e441358b8711bdfd628e8..2cb69b6d06d4f51b50da31a1a82c43f830fa27f4 100755
--- a/Filters/Points/Testing/Python/TestRadiusOutlierRemoval.py
+++ b/Filters/Points/Testing/Python/TestRadiusOutlierRemoval.py
@@ -39,8 +39,8 @@ removal.Update()
 timer.StopTimer()
 time = timer.GetElapsedTime()
 print("Time to remove points: {0}".format(time))
-print("   Number removed: {0}".format(removal.GetNumberOfPointsRemoved()),
-      " (out of: {}".format(NPts))
+print("   Number removed: {0}".format(removal.GetNumberOfPointsRemoved()))
+print("   Original number of points: {0}".format(NPts))
 
 # First output are the non-outliers
 remMapper = vtk.vtkPointGaussianMapper()
diff --git a/Filters/Points/vtkVoxelGrid.cxx b/Filters/Points/vtkVoxelGrid.cxx
index 5951bc7ea74e56573a33e6c080942cf58449511f..7b913c7f336cbc47609cd4b6a993ae5c001e5a5e 100644
--- a/Filters/Points/vtkVoxelGrid.cxx
+++ b/Filters/Points/vtkVoxelGrid.cxx
@@ -237,7 +237,7 @@ int vtkVoxelGrid::RequestData(
   // Grab the point data for interpolation
   vtkPointData *inPD = input->GetPointData();
   vtkPointData *outPD = output->GetPointData();
-  outPD->InterpolateAllocate(inPD,numBins);
+  outPD->InterpolateAllocate(inPD,numOutPts);
 
   // Finally run over all of the bins, and those that are not emoty are
   // processed. The processing consists of averaging all of the points found