freVnlSetupOptimizationOptimizer.cxx

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Version:   v.1.0.0
  Date:      $Date: 2006/09/01 12:00:00 $
  Module:    $RCSfile: freVnlSetupOptimizationOptimizer.cxx,v $
  Language:  C++


  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/

#include "freVnlSetupOptimizationOptimizer.h"
#include "itkCommand.h"
#include "itkEventObject.h"

namespace FREE
{

VnlSetupOptimizationOptimizer
::VnlSetupOptimizationOptimizer()
{
  m_CostFunctionAdaptor = 0;
  m_Maximize = false;
  m_Command = CommandType::New();
  m_Command->SetCallbackFunction( this, &VnlSetupOptimizationOptimizer::IterationReport );
  m_CachedValue = 0;
  m_CachedCurrentPosition.Fill(0);
  m_CachedDerivative.Fill(0);
}

VnlSetupOptimizationOptimizer
::~VnlSetupOptimizationOptimizer()
{
        if( m_CostFunctionAdaptor )
  {
    delete m_CostFunctionAdaptor;
    m_CostFunctionAdaptor = 0;
  }
}

void
VnlSetupOptimizationOptimizer
::PrintSelf(std::ostream& os, itk::Indent indent) const
{
  Superclass::PrintSelf( os, indent );
  os << indent << "Maximize flag: "
     << (m_Maximize ? "On" : "Off") << std::endl;
  os << indent << "Cached Value: " << m_CachedValue << std::endl;
  os << indent << "Cached Derivative: " << m_CachedDerivative << std::endl;
  os << indent << "Cached current positiion: "
     << m_CachedCurrentPosition << std::endl;
  os << "Command observer " << m_Command.GetPointer() << std::endl;
  os << "Cost Function adaptor" << m_CostFunctionAdaptor << std::endl;
}

void 
VnlSetupOptimizationOptimizer
::SetCostFunctionAdaptor( CostFunctionAdaptorType * adaptor )
{

  if( m_CostFunctionAdaptor == adaptor ) 
    {
    return;
    }

  if( m_CostFunctionAdaptor )
    {
    delete m_CostFunctionAdaptor;
    }

  m_CostFunctionAdaptor = adaptor; 

        m_CostFunctionAdaptor->AddObserver( itk::IterationEvent(), m_Command );
}

const VnlSetupOptimizationOptimizer::CostFunctionAdaptorType * 
VnlSetupOptimizationOptimizer
::GetCostFunctionAdaptor( void ) const
{
  return m_CostFunctionAdaptor;
}


VnlSetupOptimizationOptimizer::CostFunctionAdaptorType * 
VnlSetupOptimizationOptimizer
::GetCostFunctionAdaptor( void )
{
  return m_CostFunctionAdaptor;
}


VnlSetupOptimizationOptimizer::CostFunctionAdaptorType * 
VnlSetupOptimizationOptimizer
::GetNonConstCostFunctionAdaptor( void ) const
{
  return m_CostFunctionAdaptor;
}

void
VnlSetupOptimizationOptimizer
::IterationReport( const itk::EventObject & event ) 
{
  const CostFunctionAdaptorType * adaptor = this->GetCostFunctionAdaptor();
  m_CachedValue = adaptor->GetCachedValue();
  m_CachedDerivative = adaptor->GetCachedDerivative();
  m_CachedCurrentPosition = adaptor->GetCachedCurrentParameters();
  m_CurrentValue = adaptor->GetCachedValue();
  m_CurrentDecomposedValue = adaptor->GetCachedDecomposedValue();
  SetCurrentPosition(adaptor->GetCachedCurrentParameters());
  this->InvokeEvent( event );
}

} // end namespace FREE

---