00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkLevenbergMarquardtOptimizer.h,v $ 00005 Language: C++ 00006 Date: $Date: 2009-06-24 12:02:52 $ 00007 Version: $Revision: 1.24 $ 00008 00009 Copyright (c) Insight Software Consortium. All rights reserved. 00010 See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details. 00011 00012 This software is distributed WITHOUT ANY WARRANTY; without even 00013 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 00014 PURPOSE. See the above copyright notices for more information. 00015 00016 =========================================================================*/ 00017 #ifndef __itkLevenbergMarquardtOptimizer_h 00018 #define __itkLevenbergMarquardtOptimizer_h 00019 00020 #include "itkMultipleValuedNonLinearVnlOptimizer.h" 00021 #include "vnl/algo/vnl_levenberg_marquardt.h" 00022 00023 namespace itk 00024 { 00025 00031 class ITK_EXPORT LevenbergMarquardtOptimizer : 00032 public MultipleValuedNonLinearVnlOptimizer 00033 { 00034 public: 00036 typedef LevenbergMarquardtOptimizer Self; 00037 typedef MultipleValuedNonLinearVnlOptimizer Superclass; 00038 typedef SmartPointer<Self> Pointer; 00039 typedef SmartPointer<const Self> ConstPointer; 00040 00042 itkNewMacro(Self); 00043 00045 itkTypeMacro( LevenbergMarquardtOptimizer, MultipleValuedNonLinearVnlOptimizer ); 00046 00048 typedef vnl_vector<double> InternalParametersType; 00049 00051 typedef vnl_levenberg_marquardt InternalOptimizerType; 00052 00054 vnl_levenberg_marquardt * GetOptimizer(void) const; 00055 00057 void StartOptimization( void ); 00058 00060 virtual void SetCostFunction( MultipleValuedCostFunction * costFunction ); 00061 00062 void SetNumberOfIterations(unsigned int iterations); 00063 void SetValueTolerance(double tol); 00064 void SetGradientTolerance(double tol); 00065 void SetEpsilonFunction(double epsilon); 00066 00068 MeasureType GetValue() const; 00069 00070 const std::string GetStopConditionDescription() const; 00071 00072 protected: 00073 LevenbergMarquardtOptimizer(); 00074 virtual ~LevenbergMarquardtOptimizer(); 00075 00076 typedef Superclass::CostFunctionAdaptorType CostFunctionAdaptorType; 00077 00078 private: 00079 LevenbergMarquardtOptimizer(const Self&); //purposely not implemented 00080 void operator=(const Self&); //purposely not implemented 00081 00082 bool m_OptimizerInitialized; 00083 InternalOptimizerType * m_VnlOptimizer; 00084 unsigned int m_NumberOfIterations; 00085 double m_ValueTolerance; 00086 double m_GradientTolerance; 00087 double m_EpsilonFunction; 00088 00089 }; 00090 00091 } // end namespace itk 00092 00093 00094 #endif 00095