Corrected approach to restarting the optimizer.

Author: richardbeare

include/itkArrivalFunctionToPathFilter.h

@@ -238,6 +238,10 @@ class ITK_EXPORT ArrivalFunctionToPathFilter :
   typename OptimizerType::MeasureType m_TerminationValue;
   std::vector<PointsContainerType> m_PointList;
   unsigned int m_CurrentOutput;
+  // variable for communicating state of a multisegment path
+  // Sometimes an optimizer may terminate during one segment of the
+  // path, and needs to be restarted
+  bool m_TrackingPath;
 };
 
 } // end namespace itk

include/itkArrivalFunctionToPathFilter.hxx

@@ -207,7 +207,10 @@ ArrivalFunctionToPathFilter<TInputImage,TOutputPath>
     m_Optimizer->SetInitialPosition( end );
 
     // Use optimizer to back propagate from end point
-    m_Optimizer->StartOptimization( );
+    m_TrackingPath = false;
+    do {
+        m_Optimizer->StartOptimization( );
+    } while ( m_TrackingPath );
     }
 
   // Clean up by removing observer