Correct bug in Voss et al distortion correction

Author: ocommowi

Anima/registration/anatomical-commands/distortion_correction/animaDistortionCorrectionImageFilter.hxx

@@ -42,22 +42,24 @@ template< typename TInputImage >
 void DistortionCorrectionImageFilter < TInputImage >
 ::GenerateData()
 {
-  // Call a method that can be overridden by a subclass to allocate
-  // memory for the filter's outputs
-  this->AllocateOutputs();
-
-  // Call a method that can be overridden by a subclass to perform
-  // some calculations prior to splitting the main computations into
-  // separate threads
-  this->BeforeThreadedGenerateData();
-
-  using RegionType = itk::ImageRegion <TInputImage::ImageDimension>;
-  typename OutputImageType::Pointer outputImage(this->GetOutput());
-  const RegionType region = outputImage->GetRequestedRegion();
-
-  this->GetMultiThreader()->SetNumberOfWorkUnits(this->GetNumberOfWorkUnits());
-  this->GetMultiThreader()->template ParallelizeImageRegionRestrictDirection<OutputImageType::ImageDimension>(
-    this->m_Direction, region, [this](const RegionType & lambdaRegion) { this->DynamicThreadedGenerateData(lambdaRegion); }, this);
+    // Call a method that can be overridden by a subclass to allocate
+    // memory for the filter's outputs
+    this->AllocateOutputs();
+
+    // Call a method that can be overridden by a subclass to perform
+    // some calculations prior to splitting the main computations into
+    // separate threads
+    this->BeforeThreadedGenerateData();
+
+    using RegionType = itk::ImageRegion <TInputImage::ImageDimension>;
+    typename OutputImageType::Pointer outputImage(this->GetOutput());
+    const RegionType region = outputImage->GetRequestedRegion();
+
+    this->GetMultiThreader()->SetNumberOfWorkUnits(this->GetNumberOfWorkUnits());
+    this->GetMultiThreader()->template ParallelizeImageRegionRestrictDirection<OutputImageType::ImageDimension>(
+                this->m_Direction, region, [this](const RegionType & lambdaRegion) { this->DynamicThreadedGenerateData(lambdaRegion); }, this);
+
+    this->AfterThreadedGenerateData();
 }
 
 template< typename TInputImage >