STYLE: Replace T var; var = x with T var = x

Using Notepad++, Replace in Files, doing:

    Find what: ^( [ ]+)([^ ][^=\r\n]*)([ ]+)(\w+);[\r\n]+\1\4\ =
    Replace with: $1$2$3$4 =
    Filters: itk*.* !+\test
    Directory: D:\src\ITK\Modules
    [v] Match case
    (*) Regular expression

Excluded ConceptChecking, ImageKmeansModelEstimator, MahalanobisDistanceThresholdImageFunction::SetCovariance, QuasiNewtonOptimizerv4Template, VariableSizeMatrix::operator*

Modules/Core/Common/include/itkAnnulusOperator.hxx

@@ -30,9 +30,7 @@ template <typename TPixel, unsigned int TDimension, typename TAllocator>
 void
 AnnulusOperator<TPixel, TDimension, TAllocator>::CreateOperator()
 {
-  CoefficientVector coefficients;
-
-  coefficients = this->GenerateCoefficients();
+  CoefficientVector coefficients = this->GenerateCoefficients();
 
   this->Fill(coefficients);
 }

Modules/Core/Common/include/itkConstNeighborhoodIteratorWithOnlyIndex.h

@@ -409,8 +409,7 @@ inline ConstNeighborhoodIteratorWithOnlyIndex<TImage>
 operator+(const ConstNeighborhoodIteratorWithOnlyIndex<TImage> &                      it,
           const typename ConstNeighborhoodIteratorWithOnlyIndex<TImage>::OffsetType & ind)
 {
-  ConstNeighborhoodIteratorWithOnlyIndex<TImage> ret;
-  ret = it;
+  ConstNeighborhoodIteratorWithOnlyIndex<TImage> ret = it;
   ret += ind;
   return ret;
 }
@@ -428,8 +427,7 @@ inline ConstNeighborhoodIteratorWithOnlyIndex<TImage>
 operator-(const ConstNeighborhoodIteratorWithOnlyIndex<TImage> &                      it,
           const typename ConstNeighborhoodIteratorWithOnlyIndex<TImage>::OffsetType & ind)
 {
-  ConstNeighborhoodIteratorWithOnlyIndex<TImage> ret;
-  ret = it;
+  ConstNeighborhoodIteratorWithOnlyIndex<TImage> ret = it;
   ret -= ind;
   return ret;
 }

Modules/Core/Common/include/itkLaplacianOperator.hxx

@@ -35,9 +35,7 @@ template <typename TPixel, unsigned int VDimension, typename TAllocator>
 void
 LaplacianOperator<TPixel, VDimension, TAllocator>::CreateOperator()
 {
-  CoefficientVector coefficients;
-
-  coefficients = this->GenerateCoefficients();
+  CoefficientVector coefficients = this->GenerateCoefficients();
 
   this->Fill(coefficients);
 }

Modules/Core/GPUCommon/include/itkGPUInPlaceImageFilter.hxx

@@ -73,9 +73,7 @@ GPUInPlaceImageFilter<TInputImage, TOutputImage, TParentImageFilter>::AllocateOu
       {
         // if we cannot cast the input to an output type, then allocate
         // an output usual.
-        OutputImagePointer outputPtr;
-
-        outputPtr = this->GetOutput(0);
+        OutputImagePointer outputPtr = this->GetOutput(0);
         outputPtr->SetBufferedRegion(outputPtr->GetRequestedRegion());
         outputPtr->Allocate();
       }

Modules/Core/GPUCommon/include/itkGPUKernelManager.h

@@ -124,9 +124,7 @@ class GPUKernelManager : public LightObject
       return false;
     }
 
-    cl_int errid;
-
-    errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, sizeof(cl_mem), manager->GetGPUBufferPointer());
+    cl_int errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, sizeof(cl_mem), manager->GetGPUBufferPointer());
     OpenCLCheckError(errid, __FILE__, __LINE__, ITK_LOCATION);
 
     m_KernelArgumentReady[kernelIdx][argIdx].m_IsReady = true;

Modules/Core/GPUCommon/src/itkGPUKernelManager.cxx

@@ -317,9 +317,7 @@ GPUKernelManager::SetKernelArg(int kernelIdx, cl_uint argIdx, size_t argSize, co
     return false;
   }
 
-  cl_int errid;
-
-  errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, argSize, argVal);
+  cl_int errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, argSize, argVal);
   OpenCLCheckError(errid, __FILE__, __LINE__, ITK_LOCATION);
 
   m_KernelArgumentReady[kernelIdx][argIdx].m_IsReady = true;
@@ -413,9 +411,7 @@ GPUKernelManager::SetKernelArgWithImage(int kernelIdx, cl_uint argIdx, GPUDataMa
     return false;
   }
 
-  cl_int errid;
-
-  errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, sizeof(cl_mem), manager->GetGPUBufferPointer());
+  cl_int errid = clSetKernelArg(m_KernelContainer[kernelIdx], argIdx, sizeof(cl_mem), manager->GetGPUBufferPointer());
   OpenCLCheckError(errid, __FILE__, __LINE__, ITK_LOCATION);
 
   m_KernelArgumentReady[kernelIdx][argIdx].m_IsReady = true;
@@ -624,16 +620,15 @@ GPUKernelManager::LaunchKernel(int kernelIdx, int dim, size_t * globalWorkSize,
   // localWorkSize[0] = localWorkSize[1] = localWorkSize[2] = 1;
   //
 
-  cl_int errid;
-  errid = clEnqueueNDRangeKernel(m_Manager->GetCommandQueue(m_CommandQueueId),
-                                 m_KernelContainer[kernelIdx],
-                                 (cl_uint)dim,
-                                 nullptr,
-                                 globalWorkSize,
-                                 localWorkSize,
-                                 0,
-                                 nullptr,
-                                 nullptr);
+  cl_int errid = clEnqueueNDRangeKernel(m_Manager->GetCommandQueue(m_CommandQueueId),
+                                        m_KernelContainer[kernelIdx],
+                                        (cl_uint)dim,
+                                        nullptr,
+                                        globalWorkSize,
+                                        localWorkSize,
+                                        0,
+                                        nullptr,
+                                        nullptr);
   OpenCLCheckError(errid, __FILE__, __LINE__, ITK_LOCATION);
 
   /*

Modules/Core/GPUCommon/src/itkOpenCLUtil.cxx

@@ -54,9 +54,7 @@ OpenCLGetAvailableDevices(cl_platform_id platform, cl_device_type devType, cl_ui
   cl_uint        totalNumDevices;
 
   // get total # of devices
-  cl_int errid;
-
-  errid = clGetDeviceIDs(platform, devType, 0, nullptr, &totalNumDevices);
+  cl_int errid = clGetDeviceIDs(platform, devType, 0, nullptr, &totalNumDevices);
   OpenCLCheckError(errid, __FILE__, __LINE__, ITK_LOCATION);
 
   auto * totalDevices = (cl_device_id *)malloc(totalNumDevices * sizeof(cl_device_id));

Modules/Core/ImageAdaptors/include/itkNthElementPixelAccessor.h

@@ -156,9 +156,7 @@ class NthElementPixelAccessor<TOutputPixelType, itk::VariableLengthVector<TPixel
   inline ExternalType
   Get(const ActualPixelType & input) const
   {
-    ExternalType output;
-
-    output = static_cast<ExternalType>(input[m_ElementNumber]);
+    ExternalType output = static_cast<ExternalType>(input[m_ElementNumber]);
     return output;
   }
 

Modules/Core/ImageAdaptors/include/itkVectorImageToImageAdaptor.h

@@ -74,9 +74,7 @@ class VectorImageToImagePixelAccessor : private DefaultVectorPixelAccessor<TType
   inline ExternalType
   Get(const ActualPixelType & input) const
   {
-    ExternalType output;
-
-    output = input[m_ComponentIdx];
+    ExternalType output = input[m_ComponentIdx];
     return output;
   }
 

Modules/Core/ImageFunction/include/itkMeanImageFunction.hxx

@@ -34,9 +34,7 @@ template <typename TInputImage, typename TCoordRep>
 auto
 MeanImageFunction<TInputImage, TCoordRep>::EvaluateAtIndex(const IndexType & index) const -> RealType
 {
-  RealType sum;
-
-  sum = RealType{};
+  RealType sum = RealType{};
 
   const InputImageType * const image = this->GetInputImage();
 

Modules/Core/ImageFunction/include/itkNeighborhoodOperatorImageFunction.h

@@ -90,8 +90,7 @@ class ITK_TEMPLATE_EXPORT NeighborhoodOperatorImageFunction : public ImageFuncti
   Evaluate(const PointType &) const override
   {
     std::cout << "NeighborhoodOperatorImageFunction::Evaluate(): Not implemented!" << std::endl;
-    TOutput out;
-    out = 0;
+    TOutput out = 0;
     return out;
   }
 
@@ -105,8 +104,7 @@ class ITK_TEMPLATE_EXPORT NeighborhoodOperatorImageFunction : public ImageFuncti
   EvaluateAtContinuousIndex(const ContinuousIndexType &) const override
   {
     std::cout << "NeighborhoodOperatorImageFunction::EvaluateAtContinuousIndex():Not implemented!" << std::endl;
-    TOutput out;
-    out = 0;
+    TOutput out = 0;
     return out;
   }
 

Modules/Core/ImageFunction/include/itkSumOfSquaresImageFunction.hxx

@@ -35,9 +35,7 @@ template <typename TInputImage, typename TCoordRep>
 auto
 SumOfSquaresImageFunction<TInputImage, TCoordRep>::EvaluateAtIndex(const IndexType & index) const -> RealType
 {
-  RealType sumOfSquares;
-
-  sumOfSquares = RealType{};
+  RealType sumOfSquares = RealType{};
 
   const InputImageType * const image = this->GetInputImage();
 

Modules/Core/Mesh/include/itkConnectedRegionsMeshFilter.hxx

@@ -253,8 +253,7 @@ ConnectedRegionsMeshFilter<TInputMesh, TOutputMesh>::GenerateData()
       }
 
       // get the cells using the closest point and use them as seeds
-      InputMeshCellLinksContainerConstPointer cellLinks;
-      cellLinks = input->GetCellLinks();
+      InputMeshCellLinksContainerConstPointer cellLinks = input->GetCellLinks();
 
       auto links = cellLinks->ElementAt(minId);
       for (auto citer = links.cbegin(); citer != links.cend(); ++citer)

Modules/Core/SpatialObjects/include/itkContourSpatialObject.hxx

@@ -105,8 +105,7 @@ ContourSpatialObject<TDimension>::SetControlPoints(const ContourPointListType &
 {
   m_ControlPoints.clear();
 
-  typename ContourPointListType::const_iterator it;
-  it = points.begin();
+  typename ContourPointListType::const_iterator it = points.begin();
   while (it != points.end())
   {
     m_ControlPoints.push_back(*it);
@@ -219,8 +218,7 @@ ContourSpatialObject<TDimension>::Update()
             newPoint[d] = pnt[d] + i * step[d];
           }
         }
-        typename Superclass::SpatialObjectPointType newSOPoint;
-        newSOPoint = (*it);
+        typename Superclass::SpatialObjectPointType newSOPoint = (*it);
         newSOPoint.SetSpatialObject(this);
         newSOPoint.SetPositionInObjectSpace(newPoint);
         this->m_Points.push_back(newSOPoint);

Modules/Core/Transform/include/itkQuaternionRigidTransform.hxx

@@ -198,9 +198,7 @@ QuaternionRigidTransform<TParametersValueType>::ComputeMatrix()
   VnlQuaternionType conjugateRotation = m_Rotation.conjugate();
   // this is done to compensate for the transposed representation
   // between VNL and ITK
-  MatrixType newMatrix;
-
-  newMatrix = conjugateRotation.rotation_matrix_transpose();
+  MatrixType newMatrix = conjugateRotation.rotation_matrix_transpose();
   this->SetVarMatrix(newMatrix);
 }
 

Modules/Filtering/BinaryMathematicalMorphology/include/itkBinaryPruningImageFilter.hxx

@@ -111,8 +111,7 @@ BinaryPruningImageFilter<TInputImage, TOutputImage>::ComputePruneImage()
     {
       if (ot.GetCenterPixel())
       {
-        PixelType genus;
-        genus = ot.GetPixel(offset1) + ot.GetPixel(offset2);
+        PixelType genus = ot.GetPixel(offset1) + ot.GetPixel(offset2);
         genus += ot.GetPixel(offset3) + ot.GetPixel(offset4);
         genus += ot.GetPixel(offset5) + ot.GetPixel(offset6);
         genus += ot.GetPixel(offset7) + ot.GetPixel(offset8);

Modules/Filtering/Convolution/include/itkMaskedFFTNormalizedCorrelationImageFilter.hxx

@@ -646,15 +646,13 @@ MaskedFFTNormalizedCorrelationImageFilter<TInputImage, TOutputImage, TMaskImage>
   // is bigger than our input.
 
   // Cast away the constness so we can set the requested region.
-  InputImagePointer inputPtr;
-  inputPtr = const_cast<InputImageType *>(this->GetFixedImage());
+  InputImagePointer inputPtr = const_cast<InputImageType *>(this->GetFixedImage());
   inputPtr->SetRequestedRegion(this->GetFixedImage()->GetLargestPossibleRegion());
 
   inputPtr = const_cast<InputImageType *>(this->GetMovingImage());
   inputPtr->SetRequestedRegion(this->GetMovingImage()->GetLargestPossibleRegion());
 
-  MaskImagePointer maskPtr;
-  maskPtr = const_cast<MaskImageType *>(this->GetFixedImageMask());
+  MaskImagePointer maskPtr = const_cast<MaskImageType *>(this->GetFixedImageMask());
   if (maskPtr)
   {
     maskPtr->SetRequestedRegion(this->GetFixedImageMask()->GetLargestPossibleRegion());

Modules/Filtering/CurvatureFlow/include/itkBinaryMinMaxCurvatureFlowImageFilter.hxx

@@ -28,8 +28,7 @@ BinaryMinMaxCurvatureFlowImageFilter<TInputImage, TOutputImage>::BinaryMinMaxCur
 {
   m_Threshold = 0.0;
 
-  typename BinaryMinMaxCurvatureFlowFunctionType::Pointer cffp;
-  cffp = BinaryMinMaxCurvatureFlowFunctionType::New();
+  typename BinaryMinMaxCurvatureFlowFunctionType::Pointer cffp = BinaryMinMaxCurvatureFlowFunctionType::New();
 
   this->SetDifferenceFunction(static_cast<FiniteDifferenceFunctionType *>(cffp.GetPointer()));
 }

Modules/Filtering/CurvatureFlow/include/itkCurvatureFlowImageFilter.hxx

@@ -29,8 +29,7 @@ CurvatureFlowImageFilter<TInputImage, TOutputImage>::CurvatureFlowImageFilter()
   this->SetNumberOfIterations(0);
   m_TimeStep = 0.05f;
 
-  typename CurvatureFlowFunctionType::Pointer cffp;
-  cffp = CurvatureFlowFunctionType::New();
+  typename CurvatureFlowFunctionType::Pointer cffp = CurvatureFlowFunctionType::New();
 
   this->SetDifferenceFunction(static_cast<FiniteDifferenceFunctionType *>(cffp.GetPointer()));
 }

Modules/Filtering/CurvatureFlow/include/itkMinMaxCurvatureFlowFunction.hxx

@@ -404,8 +404,7 @@ MinMaxCurvatureFlowFunction<TImage>::ComputeUpdate(const NeighborhoodType & it,
     return update;
   }
 
-  PixelType threshold;
-  threshold = this->ComputeThreshold(Dispatch<ImageDimension>(), it);
+  PixelType threshold = this->ComputeThreshold(Dispatch<ImageDimension>(), it);
 
   NeighborhoodInnerProduct<ImageType> innerProduct;
   PixelType                           avgValue = innerProduct(it, m_StencilOperator);

Modules/Filtering/CurvatureFlow/include/itkMinMaxCurvatureFlowImageFilter.hxx

@@ -29,8 +29,7 @@ MinMaxCurvatureFlowImageFilter<TInputImage, TOutputImage>::MinMaxCurvatureFlowIm
 {
   m_StencilRadius = 2;
 
-  typename MinMaxCurvatureFlowFunctionType::Pointer cffp;
-  cffp = MinMaxCurvatureFlowFunctionType::New();
+  typename MinMaxCurvatureFlowFunctionType::Pointer cffp = MinMaxCurvatureFlowFunctionType::New();
 
   this->SetDifferenceFunction(static_cast<FiniteDifferenceFunctionType *>(cffp.GetPointer()));
 }

Modules/Filtering/Denoising/include/itkPatchBasedDenoisingBaseImageFilter.hxx

@@ -211,8 +211,7 @@ PatchBasedDenoisingBaseImageFilter<TInputImage, TOutputImage>::GetPatchRadiusInV
     thisPtr->m_InputImage = this->GetInput();
   }
   const typename InputImageType::SpacingType & spacing = this->m_InputImage->GetSpacing();
-  typename InputImageType::SpacingValueType    maxSpacing;
-  maxSpacing = spacing[0];
+  typename InputImageType::SpacingValueType    maxSpacing = spacing[0];
   for (unsigned int dim = 1; dim < ImageDimension; ++dim)
   {
     if (spacing[dim] > maxSpacing)

Modules/Filtering/FastMarching/include/itkFastMarchingImageFilter.hxx

@@ -136,8 +136,7 @@ FastMarchingImageFilter<TLevelSet, TSpeedImage>::Initialize(LevelSetImageType *
   // set all output value to infinity
   using OutputIterator = ImageRegionIterator<LevelSetImageType>;
 
-  PixelType outputPixel;
-  outputPixel = m_LargeValue;
+  PixelType outputPixel = m_LargeValue;
 
   for (OutputIterator outIt(output, output->GetBufferedRegion()); !outIt.IsAtEnd(); ++outIt)
   {

Modules/Filtering/ImageFeature/include/itkHessian3DToVesselnessMeasureImageFilter.hxx

@@ -53,8 +53,8 @@ Hessian3DToVesselnessMeasureImageFilter<TPixel>::GenerateData()
 
   // walk the region of eigen values and get the vesselness measure
   EigenValueArrayType                                 eigenValue;
-  ImageRegionConstIterator<EigenValueOutputImageType> it;
-  it = ImageRegionConstIterator<EigenValueOutputImageType>(eigenImage, eigenImage->GetRequestedRegion());
+  ImageRegionConstIterator<EigenValueOutputImageType> it =
+    ImageRegionConstIterator<EigenValueOutputImageType>(eigenImage, eigenImage->GetRequestedRegion());
   ImageRegionIterator<OutputImageType> oit;
   this->AllocateOutputs();
   oit = ImageRegionIterator<OutputImageType>(output, output->GetRequestedRegion());

Modules/Filtering/ImageFilterBase/include/itkMovingHistogramImageFilterBase.hxx

@@ -52,8 +52,8 @@ MovingHistogramImageFilterBase<TInputImage, TOutputImage, TKernel>::SetKernel(co
   tmpSEImage->SetRegions(kernel.GetSize());
   tmpSEImage->Allocate();
   RegionType                                  tmpSEImageRegion = tmpSEImage->GetRequestedRegion();
-  ImageRegionIteratorWithIndex<BoolImageType> kernelImageIt;
-  kernelImageIt = ImageRegionIteratorWithIndex<BoolImageType>(tmpSEImage, tmpSEImageRegion);
+  ImageRegionIteratorWithIndex<BoolImageType> kernelImageIt =
+    ImageRegionIteratorWithIndex<BoolImageType>(tmpSEImage, tmpSEImageRegion);
   KernelIteratorType kernel_it = kernel.Begin();
   OffsetListType     kernelOffsets;
 

Modules/Filtering/ImageGradient/include/itkVectorGradientMagnitudeImageFilter.h

@@ -318,9 +318,7 @@ class ITK_TEMPLATE_EXPORT VectorGradientMagnitudeImageFilter : public ImageToIma
   NonPCEvaluateAtNeighborhood(const ConstNeighborhoodIteratorType & it) const
   {
     unsigned int i, j;
-    TRealType    dx, sum, accum;
-
-    accum = TRealType{};
+    TRealType    dx, sum, accum = TRealType{};
     for (i = 0; i < ImageDimension; ++i)
     {
       sum = TRealType{};

Modules/Filtering/ImageGrid/include/itkBSplineControlPointImageFilter.hxx

@@ -114,8 +114,7 @@ BSplineControlPointImageFilter<TInputImage, TOutputImage>::SetSplineOrder(ArrayT
 
     if (this->m_DoMultilevel)
     {
-      typename KernelType::MatrixType C;
-      C = this->m_Kernel[i]->GetShapeFunctionsInZeroToOneInterval();
+      typename KernelType::MatrixType C = this->m_Kernel[i]->GetShapeFunctionsInZeroToOneInterval();
 
       vnl_matrix<RealType> R;
       vnl_matrix<RealType> S;

Modules/Filtering/ImageGrid/include/itkBSplineResampleImageFilterBase.hxx

@@ -348,8 +348,7 @@ BSplineResampleImageFilterBase<TInputImage, TOutputImage>::ReduceNDImage(OutputI
 
   // Initialize scratchImage space and allocate memory
   InitializeScratch(startSize);
-  typename TOutputImage::Pointer scratchImage;
-  scratchImage = TOutputImage::New();
+  typename TOutputImage::Pointer scratchImage = TOutputImage::New();
   scratchImage->CopyInformation(inputPtr);
   RegionType scratchRegion = inputPtr->GetBufferedRegion();
   currentSize = startSize;
@@ -378,8 +377,7 @@ BSplineResampleImageFilterBase<TInputImage, TOutputImage>::ReduceNDImage(OutputI
    * OutputImage for direct writing into the final variable. */
 
   // The first time through the loop our input image is inputPtr
-  typename TInputImage::ConstPointer workingImage;
-  workingImage = inputPtr;
+  typename TInputImage::ConstPointer workingImage = inputPtr;
 
   unsigned int     count = scratchRegion.GetNumberOfPixels() * ImageDimension;
   ProgressReporter progress(this, 0, count, 10);
@@ -462,8 +460,7 @@ BSplineResampleImageFilterBase<TInputImage, TOutputImage>::ExpandNDImage(OutputI
 
   // Initialize scratchImage space and allocate memory
   InitializeScratch(startSize);
-  typename TOutputImage::Pointer scratchImage;
-  scratchImage = TOutputImage::New();
+  typename TOutputImage::Pointer scratchImage = TOutputImage::New();
   scratchImage->CopyInformation(inputPtr);
   RegionType scratchRegion = inputPtr->GetBufferedRegion();
   currentSize = startSize;
@@ -493,8 +490,7 @@ BSplineResampleImageFilterBase<TInputImage, TOutputImage>::ExpandNDImage(OutputI
   **/
 
   // The first time through the loop our input image is m_Image
-  typename TInputImage::ConstPointer workingImage;
-  workingImage = inputPtr;
+  typename TInputImage::ConstPointer workingImage = inputPtr;
 
   RegionType workingRegion = validRegion;