diff --git a/Imaging/Core/vtkImageDifference.cxx b/Imaging/Core/vtkImageDifference.cxx
index 368fe37ada..3a2558a6ea 100644
--- a/Imaging/Core/vtkImageDifference.cxx
+++ b/Imaging/Core/vtkImageDifference.cxx
@@ -25,6 +25,7 @@
 #include <numeric>
 
 constexpr int MAX_NCOMPS = 4;
+constexpr double DEFAULT_ERROR = 1000.;
 using max_ncomps_array_t = std::array<int, MAX_NCOMPS>;
 
 vtkStandardNewMacro(vtkImageDifference);
@@ -317,7 +318,6 @@ void vtkImageDifference::ThreadedRequestData(vtkInformation* vtkNotUsed(request)
         max_ncomps_array_t rgbaMax;
         rgbaMax.fill(0);
         max_ncomps_array_t rgbaTresh;
-        rgbaTresh.fill(1000);
 
         // ignore the boundary within two pixels as we cannot
         // do a good average calc on the boundary
@@ -326,6 +326,7 @@ void vtkImageDifference::ThreadedRequestData(vtkInformation* vtkNotUsed(request)
         {
           for (int direction = 0; direction <= 1; ++direction)
           {
+            rgbaTresh.fill(DEFAULT_ERROR);
             unsigned char* dir1Ptr0 = direction == 0 ? in1Ptr0 : in2Ptr0;
             unsigned char* dir2Ptr0 = direction == 0 ? in2Ptr0 : in1Ptr0;
             bool haveValues = false;
@@ -401,7 +402,7 @@ void vtkImageDifference::ThreadedRequestData(vtkInformation* vtkNotUsed(request)
 
         std::copy(rgbaMax.begin(), rgbaMax.end(), rgbaTresh.begin());
 
-        error += std::accumulate(rgbaTresh.begin(), rgbaTresh.end(), 0) / (nComp * 255);
+        error += std::accumulate(rgbaTresh.begin(), rgbaTresh.end(), 0.) / (nComp * 255.);
 
         for (int i = 0; i < nComp; i++)
         {
@@ -410,7 +411,7 @@ void vtkImageDifference::ThreadedRequestData(vtkInformation* vtkNotUsed(request)
           *outPtr0++ = static_cast<unsigned char>(rgbaTresh[i]);
         }
         thresholdedError +=
-          std::accumulate(rgbaTresh.begin(), rgbaTresh.end(), 0) / (nComp * 255.0);
+          std::accumulate(rgbaTresh.begin(), rgbaTresh.end(), 0.) / (nComp * 255.);
 
         in1Ptr0 += nComp;
         in2Ptr0 += nComp;
@@ -489,8 +490,8 @@ int vtkImageDifference::RequestData(
     // Report errors here, do not report errors while multithreading!
     vtkErrorMacro("RequestData: " << this->ErrorMessage);
     this->ErrorMessage = nullptr;
-    this->Error = 1000.0;
-    this->ThresholdedError = 1000.0;
+    this->Error = DEFAULT_ERROR;
+    this->ThresholdedError = DEFAULT_ERROR;
     ret = 0;
   }
 
@@ -515,8 +516,8 @@ int vtkImageDifference::RequestInformation(vtkInformation* vtkNotUsed(request),
   if (in1Ext[0] != in2Ext[0] || in1Ext[1] != in2Ext[1] || in1Ext[2] != in2Ext[2] ||
     in1Ext[3] != in2Ext[3] || in1Ext[4] != in2Ext[4] || in1Ext[5] != in2Ext[5])
   {
-    this->Error = 1000.0;
-    this->ThresholdedError = 1000.0;
+    this->Error = DEFAULT_ERROR;
+    this->ThresholdedError = DEFAULT_ERROR;
 
     vtkErrorMacro("ExecuteInformation: Input are not the same size.\n"
       << " Input1 is: " << in1Ext[0] << "," << in1Ext[1] << "," << in1Ext[2] << "," << in1Ext[3]