go home Home | Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Data Structures | File List | Namespace Members | Data Fields | Globals | Related Pages
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes
elastix::ConjugateGradient< TElastix > Class Template Reference

#include <elxConjugateGradient.h>

Detailed Description

template<class TElastix>
class elastix::ConjugateGradient< TElastix >

An optimizer based on the itk::GenericConjugateGradientOptimizer.

A ConjugateGradient optimizer, using the itk::MoreThuenteLineSearchOptimizer. Different conjugate gradient methods can be selected with this optimizer.

This optimizer support the NewSamplesEveryIteration option. It requests new samples for the computation of each search direction (not during the line search). Actually this makes no sense for a conjugate gradient optimizer. So, think twice before using the NewSamplesEveryIteration option.

The parameters used in this class are:

Parameters:

Optimizer: Select this optimizer as follows:
(Optimizer "ConjugateGradient")

GenerateLineSearchIterations: Whether line search iteration should be counted as elastix-iterations.
example: (GenerateLineSearchIterations "true")
Can only be specified for all resolutions at once.
Default value: "false".

MaximumNumberOfIterations: The maximum number of iterations in each resolution.
example: (MaximumNumberOfIterations 100 100 50)
Default value: 100.

MaximumNumberOfLineSearchIterations: The maximum number of iterations in each resolution.
example: (MaximumNumberOfIterations 10 10 5)
Default value: 10.

StepLength: Set the length of the initial step tried by the itk::MoreThuenteLineSearchOptimizer.
example: (StepLength 2.0 1.0 0.5)
Default value: 1.0.

LineSearchValueTolerance: Determine the Wolfe conditions that the itk::MoreThuenteLineSearchOptimizer tries to satisfy.
example: (LineSearchValueTolerance 0.0001 0.0001 0.0001)
Default value: 0.0001.

LineSearchGradientTolerance: Determine the Wolfe conditions that the itk::MoreThuenteLineSearchOptimizer tries to satisfy.
example: (LineSearchGradientTolerance 0.9 0.9 0.9)
Default value: 0.9.

ValueTolerance: Stopping criterion. See the documentation of the itk::GenericConjugateGradientOptimizer for more information.
example: (ValueTolerance 0.001 0.0001 0.000001)
Default value: 0.00001.

GradientMagnitudeTolerance: Stopping criterion. See the documentation of the itk::GenericConjugateGradientOptimizer for more information.
example: (GradientMagnitudeTolerance 0.001 0.0001 0.000001)
Default value: 0.000001.

ConjugateGradientType: a string that defines how 'beta' is computed in each resolution. The following methods are implemented: "SteepestDescent", "FletcherReeves", "PolakRibiere", "DaiYuan", "HestenesStiefel", and "DaiYuanHestenesStiefel". "SteepestDescent" simply sets beta=0. See the source code of the GenericConjugateGradientOptimizer for more information.
example: (ConjugateGradientType "FletcherReeves" "PolakRibiere")
Default value: "DaiYuanHestenesStiefel".

StopIfWolfeNotSatisfied: Whether to stop the optimisation if in one iteration the Wolfe conditions can not be satisfied by the itk::MoreThuenteLineSearchOptimizer.
In general it is wise to do so.
example: (StopIfWolfeNotSatisfied "true" "false")
Default value: "true".

Definition at line 91 of file elxConjugateGradient.h.

Inheritance diagram for elastix::ConjugateGradient< TElastix >:
Inheritance graph
[legend]

Public Types

using ConstPointer = itk::SmartPointer< const Self >
 
using EventPassThroughPointer = typename EventPassThroughType::Pointer
 
using EventPassThroughType = itk::ReceptorMemberCommand< Self >
 
using ITKBaseType = typename Superclass2::ITKBaseType
 
using LineOptimizerPointer = LineOptimizerType::Pointer
 
using LineOptimizerType = itk::MoreThuenteLineSearchOptimizer
 
using Pointer = itk::SmartPointer< Self >
 
using Self = ConjugateGradient
 
enum  StopConditionType
 
using Superclass1 = GenericConjugateGradientOptimizer
 
using Superclass2 = OptimizerBase< TElastix >
 
- Public Types inherited from itk::GenericConjugateGradientOptimizer
using BetaDefinitionMapType = std::map< BetaDefinitionType, ComputeBetaFunctionType >
 
using BetaDefinitionType = std::string
 
using ComputeBetaFunctionType = double(Self::*)(const DerivativeType &, const DerivativeType &, const ParametersType &)
 
using ConstPointer = SmartPointer< const Self >
 
using LineSearchOptimizerPointer = LineSearchOptimizerType::Pointer
 
using LineSearchOptimizerType = LineSearchOptimizer
 
using Pointer = SmartPointer< Self >
 
using ScaledCostFunctionType = ScaledSingleValuedCostFunction
 
using ScalesType = NonLinearOptimizer::ScalesType
 
using Self = GenericConjugateGradientOptimizer
 
enum  StopConditionType {
  MetricError , LineSearchError , MaximumNumberOfIterations , GradientMagnitudeTolerance ,
  ValueTolerance , InfiniteBeta , Unknown
}
 
using Superclass = ScaledSingleValuedNonLinearOptimizer
 
- Public Types inherited from itk::ScaledSingleValuedNonLinearOptimizer
using ConstPointer = SmartPointer< const Self >
 
using Pointer = SmartPointer< Self >
 
using ScaledCostFunctionPointer = ScaledCostFunctionType::Pointer
 
using ScaledCostFunctionType = ScaledSingleValuedCostFunction
 
using ScalesType = NonLinearOptimizer::ScalesType
 
using Self = ScaledSingleValuedNonLinearOptimizer
 
using Superclass = SingleValuedNonLinearOptimizer
 
- Public Types inherited from elastix::OptimizerBase< TElastix >
using ElastixType = TElastix
 
using ITKBaseType = itk::Optimizer
 
using ParametersType = typename ITKBaseType::ParametersType
 
using RegistrationType = typename ElastixType::RegistrationBaseType
 
using Self = OptimizerBase
 
using Superclass = BaseComponentSE< TElastix >
 
- Public Types inherited from elastix::BaseComponentSE< TElastix >
using ConfigurationPointer = Configuration::Pointer
 
using ElastixType = TElastix
 
using RegistrationType = typename ElastixType::RegistrationBaseType
 
using Self = BaseComponentSE
 
using Superclass = BaseComponent
 

Public Member Functions

void AfterEachIteration () override
 
void AfterEachResolution () override
 
void AfterRegistration () override
 
void BeforeEachResolution () override
 
void BeforeRegistration () override
 
 elxClassNameMacro ("ConjugateGradient")
 
virtual const char * GetClassName () const
 
virtual bool GetStartLineSearch () const
 
 ITK_DISALLOW_COPY_AND_MOVE (ConjugateGradient)
 
void StartOptimization () override
 
- Public Member Functions inherited from itk::GenericConjugateGradientOptimizer
virtual const BetaDefinitionTypeGetBetaDefinition ()
 
virtual const char * GetClassName () const
 
virtual const DerivativeType & GetCurrentGradient ()
 
virtual unsigned long GetCurrentIteration () const
 
virtual double GetCurrentStepLength () const
 
virtual MeasureType GetCurrentValue () const
 
virtual double GetGradientMagnitudeTolerance () const
 
virtual bool GetInLineSearch () const
 
virtual unsigned long GetMaximumNumberOfIterations () const
 
virtual unsigned long GetMaxNrOfItWithoutImprovement () const
 
virtual const StopConditionTypeGetStopCondition ()
 
virtual double GetValueTolerance () const
 
 ITK_DISALLOW_COPY_AND_MOVE (GenericConjugateGradientOptimizer)
 
 itkGetModifiableObjectMacro (LineSearchOptimizer, LineSearchOptimizerType)
 
virtual void ResumeOptimization ()
 
void SetBetaDefinition (const BetaDefinitionType &arg)
 
virtual void SetGradientMagnitudeTolerance (double _arg)
 
virtual void SetLineSearchOptimizer (LineSearchOptimizerType *_arg)
 
virtual void SetMaximumNumberOfIterations (unsigned long _arg)
 
virtual void SetMaxNrOfItWithoutImprovement (unsigned long arg)
 
virtual void SetValueTolerance (double _arg)
 
void StartOptimization () override
 
virtual void StopOptimization ()
 
- Public Member Functions inherited from itk::ScaledSingleValuedNonLinearOptimizer
virtual const char * GetClassName () const
 
const ParametersType & GetCurrentPosition () const override
 
virtual bool GetMaximize () const
 
virtual const ScaledCostFunctionTypeGetScaledCostFunction ()
 
virtual const ParametersType & GetScaledCurrentPosition ()
 
bool GetUseScales () const
 
virtual void InitializeScales ()
 
 ITK_DISALLOW_COPY_AND_MOVE (ScaledSingleValuedNonLinearOptimizer)
 
virtual void MaximizeOff ()
 
virtual void MaximizeOn ()
 
void SetCostFunction (CostFunctionType *costFunction) override
 
virtual void SetMaximize (bool _arg)
 
virtual void SetUseScales (bool arg)
 
- Public Member Functions inherited from elastix::OptimizerBase< TElastix >
void AfterRegistrationBase () override
 
void BeforeEachResolutionBase () override
 
ITKBaseTypeGetAsITKBaseType ()
 
const ITKBaseTypeGetAsITKBaseType () const
 
virtual const char * GetClassName () const
 
 ITK_DISALLOW_COPY_AND_MOVE (OptimizerBase)
 
virtual void SetCurrentPositionPublic (const ParametersType &param)
 
virtual void SetSinusScales (double amplitude, double frequency, unsigned long numberOfParameters)
 
- Public Member Functions inherited from elastix::BaseComponentSE< TElastix >
void AddTargetCellToIterationInfo (const char *const name)
 
ConfigurationGetConfiguration () const
 
ElastixTypeGetElastix () const
 
xl::xoutbaseGetIterationInfoAt (const char *const name)
 
RegistrationTypeGetRegistration () const
 
 ITK_DISALLOW_COPY_AND_MOVE (BaseComponentSE)
 
int RemoveTargetCellFromIterationInfo (const char *const name)
 
void SetConfiguration (Configuration *_arg)
 
void SetElastix (ElastixType *_arg)
 
- Public Member Functions inherited from elastix::BaseComponent
virtual void AfterEachIteration ()
 
virtual void AfterEachIterationBase ()
 
virtual void AfterEachResolution ()
 
virtual void AfterEachResolutionBase ()
 
virtual void AfterRegistration ()
 
virtual void AfterRegistrationBase ()
 
virtual int BeforeAll ()
 
virtual int BeforeAllBase ()
 
virtual void BeforeEachResolution ()
 
virtual void BeforeEachResolutionBase ()
 
virtual void BeforeRegistration ()
 
virtual void BeforeRegistrationBase ()
 
virtual const char * elxGetClassName () const
 
const char * GetComponentLabel () const
 
 ITK_DISALLOW_COPY_AND_MOVE (BaseComponent)
 
 itkTypeMacroNoParent (BaseComponent)
 
void SetComponentLabel (const char *label, unsigned int idx)
 

Static Public Member Functions

static Pointer New ()
 
- Static Public Member Functions inherited from itk::GenericConjugateGradientOptimizer
static Pointer New ()
 
- Static Public Member Functions inherited from itk::ScaledSingleValuedNonLinearOptimizer
static Pointer New ()
 
- Static Public Member Functions inherited from elastix::BaseComponent
template<typename TBaseComponent >
static auto AsITKBaseType (TBaseComponent *const baseComponent) -> decltype(baseComponent->GetAsITKBaseType())
 
static void InitializeElastixExecutable ()
 
static bool IsElastixLibrary ()
 

Protected Member Functions

 ConjugateGradient ()
 
virtual std::string DeterminePhase () const
 
virtual std::string GetLineSearchStopCondition () const
 
void LineSearch (const ParametersType searchDir, double &step, ParametersType &x, MeasureType &f, DerivativeType &g) override
 
bool TestConvergence (bool firstLineSearchDone) override
 
 ~ConjugateGradient () override=default
 
- Protected Member Functions inherited from itk::GenericConjugateGradientOptimizer
void AddBetaDefinition (const BetaDefinitionType &name, ComputeBetaFunctionType function)
 
virtual double ComputeBeta (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaDY (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaDYHS (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaFR (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaHS (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaPR (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
double ComputeBetaSD (const DerivativeType &previousGradient, const DerivativeType &gradient, const ParametersType &previousSearchDir)
 
virtual void ComputeSearchDirection (const DerivativeType &previousGradient, const DerivativeType &gradient, ParametersType &searchDir)
 
 GenericConjugateGradientOptimizer ()
 
virtual void LineSearch (const ParametersType searchDir, double &step, ParametersType &x, MeasureType &f, DerivativeType &g)
 
void PrintSelf (std::ostream &os, Indent indent) const override
 
virtual void SetInLineSearch (bool _arg)
 
virtual bool TestConvergence (bool firstLineSearchDone)
 
 ~GenericConjugateGradientOptimizer () override=default
 
- Protected Member Functions inherited from itk::ScaledSingleValuedNonLinearOptimizer
virtual void GetScaledDerivative (const ParametersType &parameters, DerivativeType &derivative) const
 
virtual MeasureType GetScaledValue (const ParametersType &parameters) const
 
virtual void GetScaledValueAndDerivative (const ParametersType &parameters, MeasureType &value, DerivativeType &derivative) const
 
void PrintSelf (std::ostream &os, Indent indent) const override
 
 ScaledSingleValuedNonLinearOptimizer ()
 
void SetCurrentPosition (const ParametersType &param) override
 
virtual void SetScaledCurrentPosition (const ParametersType &parameters)
 
 ~ScaledSingleValuedNonLinearOptimizer () override=default
 
- Protected Member Functions inherited from elastix::OptimizerBase< TElastix >
virtual bool GetNewSamplesEveryIteration () const
 
 OptimizerBase ()=default
 
virtual void SelectNewSamples ()
 
 ~OptimizerBase () override=default
 
- Protected Member Functions inherited from elastix::BaseComponentSE< TElastix >
 BaseComponentSE ()=default
 
 ~BaseComponentSE () override=default
 
- Protected Member Functions inherited from elastix::BaseComponent
 BaseComponent ()=default
 
virtual ~BaseComponent ()=default
 

Protected Attributes

LineOptimizerPointer m_LineOptimizer
 
- Protected Attributes inherited from itk::GenericConjugateGradientOptimizer
BetaDefinitionType m_BetaDefinition
 
BetaDefinitionMapType m_BetaDefinitionMap
 
DerivativeType m_CurrentGradient
 
unsigned long m_CurrentIteration { 0 }
 
double m_CurrentStepLength { 0.0 }
 
MeasureType m_CurrentValue { 0.0 }
 
bool m_InLineSearch { false }
 
bool m_PreviousGradientAndSearchDirValid { false }
 
bool m_Stop { false }
 
StopConditionType m_StopCondition { Unknown }
 
bool m_UseDefaultMaxNrOfItWithoutImprovement { true }
 
- Protected Attributes inherited from itk::ScaledSingleValuedNonLinearOptimizer
ScaledCostFunctionPointer m_ScaledCostFunction
 
ParametersType m_ScaledCurrentPosition
 
- Protected Attributes inherited from elastix::BaseComponentSE< TElastix >
ConfigurationPointer m_Configuration {}
 
itk::WeakPointer< TElastix > m_Elastix {}
 
RegistrationTypem_Registration {}
 

Private Member Functions

void InvokeIterationEvent (const itk::EventObject &event)
 

Private Attributes

 elxOverrideGetSelfMacro
 
EventPassThroughPointer m_EventPasser
 
bool m_GenerateLineSearchIterations
 
double m_SearchDirectionMagnitude
 
bool m_StartLineSearch
 
bool m_StopIfWolfeNotSatisfied
 
bool m_WolfeIsStopCondition
 

Member Typedef Documentation

◆ ConstPointer

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::ConstPointer = itk::SmartPointer<const Self>

Definition at line 103 of file elxConjugateGradient.h.

◆ EventPassThroughPointer

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::EventPassThroughPointer = typename EventPassThroughType::Pointer

Definition at line 134 of file elxConjugateGradient.h.

◆ EventPassThroughType

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::EventPassThroughType = itk::ReceptorMemberCommand<Self>

Definition at line 133 of file elxConjugateGradient.h.

◆ ITKBaseType

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::ITKBaseType = typename Superclass2::ITKBaseType

Definition at line 128 of file elxConjugateGradient.h.

◆ LineOptimizerPointer

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::LineOptimizerPointer = LineOptimizerType::Pointer

Definition at line 132 of file elxConjugateGradient.h.

◆ LineOptimizerType

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::LineOptimizerType = itk::MoreThuenteLineSearchOptimizer

Extra typedefs

Definition at line 131 of file elxConjugateGradient.h.

◆ Pointer

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::Pointer = itk::SmartPointer<Self>

Definition at line 102 of file elxConjugateGradient.h.

◆ Self

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::Self = ConjugateGradient

Standard ITK.

Definition at line 99 of file elxConjugateGradient.h.

◆ Superclass1

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::Superclass1 = GenericConjugateGradientOptimizer

Definition at line 100 of file elxConjugateGradient.h.

◆ Superclass2

template<class TElastix >
using elastix::ConjugateGradient< TElastix >::Superclass2 = OptimizerBase<TElastix>

Definition at line 101 of file elxConjugateGradient.h.

Member Enumeration Documentation

◆ StopConditionType

Definition at line 71 of file itkGenericConjugateGradientOptimizer.h.

Constructor & Destructor Documentation

◆ ConjugateGradient()

template<class TElastix >
elastix::ConjugateGradient< TElastix >::ConjugateGradient ( )
protected

◆ ~ConjugateGradient()

template<class TElastix >
elastix::ConjugateGradient< TElastix >::~ConjugateGradient ( )
overrideprotecteddefault

Member Function Documentation

◆ AfterEachIteration()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::AfterEachIteration ( )
overridevirtual

Reimplemented from elastix::BaseComponent.

◆ AfterEachResolution()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::AfterEachResolution ( )
overridevirtual

Reimplemented from elastix::BaseComponent.

◆ AfterRegistration()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::AfterRegistration ( )
overridevirtual

Reimplemented from elastix::BaseComponent.

◆ BeforeEachResolution()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::BeforeEachResolution ( )
overridevirtual

Reimplemented from elastix::BaseComponent.

◆ BeforeRegistration()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::BeforeRegistration ( )
overridevirtual

Methods to set parameters and print output at different stages in the registration process.

Reimplemented from elastix::BaseComponent.

◆ DeterminePhase()

template<class TElastix >
virtual std::string elastix::ConjugateGradient< TElastix >::DeterminePhase ( ) const
protectedvirtual

Generate a string, representing the phase of optimisation (line search, main)

◆ elxClassNameMacro()

template<class TElastix >
elastix::ConjugateGradient< TElastix >::elxClassNameMacro ( "ConjugateGradient< TElastix >"  )

Name of this class. Use this name in the parameter file to select this specific optimizer.
example: (Optimizer "ConjugateGradient")

◆ GetClassName()

template<class TElastix >
virtual const char * elastix::ConjugateGradient< TElastix >::GetClassName ( ) const
virtual

Run-time type information (and related methods).

Reimplemented from itk::GenericConjugateGradientOptimizer.

◆ GetLineSearchStopCondition()

template<class TElastix >
virtual std::string elastix::ConjugateGradient< TElastix >::GetLineSearchStopCondition ( ) const
protectedvirtual

Convert the line search stop condition to a string

◆ GetStartLineSearch()

template<class TElastix >
virtual bool elastix::ConjugateGradient< TElastix >::GetStartLineSearch ( ) const
virtual

◆ InvokeIterationEvent()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::InvokeIterationEvent ( const itk::EventObject &  event)
private

◆ ITK_DISALLOW_COPY_AND_MOVE()

template<class TElastix >
elastix::ConjugateGradient< TElastix >::ITK_DISALLOW_COPY_AND_MOVE ( ConjugateGradient< TElastix >  )

◆ LineSearch()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::LineSearch ( const ParametersType  searchDir,
double step,
ParametersType x,
MeasureType &  f,
DerivativeType &  g 
)
overrideprotected

Call the superclass' implementation. If an itk::ExceptionObject is caught, because the line search optimizer tried a too big step, the exception is printed, but ignored further. The optimizer stops, but elastix just goes on to the next resolution.

◆ New()

template<class TElastix >
static Pointer elastix::ConjugateGradient< TElastix >::New ( )
static

Method for creation through the object factory.

◆ StartOptimization()

template<class TElastix >
void elastix::ConjugateGradient< TElastix >::StartOptimization ( )
override

Check if any scales are set, and set the UseScales flag on or off; after that call the superclass' implementation

◆ TestConvergence()

template<class TElastix >
bool elastix::ConjugateGradient< TElastix >::TestConvergence ( bool  firstLineSearchDone)
overrideprotectedvirtual

Reimplement the superclass. Calls the superclass' implementation and checks if the MoreThuente line search routine has stopped with Wolfe conditions satisfied.

Reimplemented from itk::GenericConjugateGradientOptimizer.

Field Documentation

◆ elxOverrideGetSelfMacro

template<class TElastix >
elastix::ConjugateGradient< TElastix >::elxOverrideGetSelfMacro
private

Definition at line 190 of file elxConjugateGradient.h.

◆ m_EventPasser

template<class TElastix >
EventPassThroughPointer elastix::ConjugateGradient< TElastix >::m_EventPasser
private

Definition at line 195 of file elxConjugateGradient.h.

◆ m_GenerateLineSearchIterations

template<class TElastix >
bool elastix::ConjugateGradient< TElastix >::m_GenerateLineSearchIterations
private

Definition at line 198 of file elxConjugateGradient.h.

◆ m_LineOptimizer

template<class TElastix >
LineOptimizerPointer elastix::ConjugateGradient< TElastix >::m_LineOptimizer
protected

Definition at line 164 of file elxConjugateGradient.h.

◆ m_SearchDirectionMagnitude

template<class TElastix >
double elastix::ConjugateGradient< TElastix >::m_SearchDirectionMagnitude
private

Definition at line 196 of file elxConjugateGradient.h.

◆ m_StartLineSearch

template<class TElastix >
bool elastix::ConjugateGradient< TElastix >::m_StartLineSearch
private

Definition at line 197 of file elxConjugateGradient.h.

◆ m_StopIfWolfeNotSatisfied

template<class TElastix >
bool elastix::ConjugateGradient< TElastix >::m_StopIfWolfeNotSatisfied
private

Definition at line 199 of file elxConjugateGradient.h.

◆ m_WolfeIsStopCondition

template<class TElastix >
bool elastix::ConjugateGradient< TElastix >::m_WolfeIsStopCondition
private

Definition at line 200 of file elxConjugateGradient.h.



Generated on 2023-01-13 for elastix by doxygen 1.9.6 elastix logo