/*=auto=========================================================================
(c) Copyright 2003 Massachusetts Institute of Technology (MIT) All Rights Reserved.

This software ("3D Slicer") is provided by The Brigham and Women's 
Hospital, Inc. on behalf of the copyright holders and contributors.
Permission is hereby granted, without payment, to copy, modify, display 
and distribute this software and its documentation, if any, for  
research purposes only, provided that (1) the above copyright notice and 
the following four paragraphs appear on all copies of this software, and 
(2) that source code to any modifications to this software be made 
publicly available under terms no more restrictive than those in this 
License Agreement. Use of this software constitutes acceptance of these 
terms and conditions.

3D Slicer Software has not been reviewed or approved by the Food and 
Drug Administration, and is for non-clinical, IRB-approved Research Use 
Only.  In no event shall data or images generated through the use of 3D 
Slicer Software be used in the provision of patient care.

IN NO EVENT SHALL THE COPYRIGHT HOLDERS AND CONTRIBUTORS BE LIABLE TO 
ANY PARTY FOR DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL 
DAMAGES ARISING OUT OF THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, 
EVEN IF THE COPYRIGHT HOLDERS AND CONTRIBUTORS HAVE BEEN ADVISED OF THE 
POSSIBILITY OF SUCH DAMAGE.

THE COPYRIGHT HOLDERS AND CONTRIBUTORS SPECIFICALLY DISCLAIM ANY EXPRESS 
OR IMPLIED WARRANTIES INCLUDING, BUT NOT LIMITED TO, THE IMPLIED 
WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, AND 
NON-INFRINGEMENT.

THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS 
IS." THE COPYRIGHT HOLDERS AND CONTRIBUTORS HAVE NO OBLIGATION TO 
PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.

=========================================================================auto=*/
#ifndef __KLRegistration_txx
#define __KLRegistration_txx

#include "KLRegistration.h"

#include "itkCommand.h"

#include "vtkMatrix4x4.h"
// #include "NewStoppingCondition.h"

namespace itk
{

template <typename TFixedImage, typename TMovingImage>
KLRegistration<TFixedImage,TMovingImage>::KLRegistration()
{

  // Default number of bins
  unsigned int nBins = 32;
  HistogramSizeType histSize;
  histSize[0] = nBins;
  histSize[1] = nBins;
  this->m_Metric->UsePaddingValue(false);
  this->m_Metric->SetHistogramSize(histSize);
  this->m_Metric->SetDerivativeStepLength(0.1); // 0.1 mm

  // set the step length scales
  typedef typename MetricType::ScalesType ScalesType;
  ScalesType DerivativeStepLengthScales( this->GetNumberOfParameters() );
  DerivativeStepLengthScales.fill(1.0);
  this->m_Metric->SetDerivativeStepLengthScales(DerivativeStepLengthScales);

  // possible memory leak
  typedef typename InterpolatorType::Pointer  InterpolatorPointer;
  InterpolatorPointer  TrainingInterpolator = InterpolatorType::New();
  this->m_Metric->SetTrainingInterpolator(TrainingInterpolator);

  // Default parameters
  this->m_MovingImageStandardDeviation = 0.4;
  this->m_FixedImageStandardDeviation = 0.4;
  this->m_NumberOfSpatialSamples = 50;
}

//----------------------------------------------------------------------

// Clean Up: Remove the observer
template <typename TFixedImage, typename TMovingImage>
KLRegistration<TFixedImage,TMovingImage>::~KLRegistration()
{
}

//----------------------------------------------------------------------------

template < typename TFixedImage, typename TMovingImage  >
void KLRegistration<TFixedImage,TMovingImage>
::PrintSelf(std::ostream& os, Indent indent) const
{
  Superclass::PrintSelf(os, indent);

  os << indent << "MovingImageStandardDeviation: " 
        << m_MovingImageStandardDeviation   << endl;
  os << indent << "FixedImageStandardDeviation: " 
        <<   m_FixedImageStandardDeviation  << endl; 
  os << indent << "NumberOfSpatialSamples: " 
        <<   m_NumberOfSpatialSamples       << endl;
}

//----------------------------------------------------------------------

template < typename TFixedImage, typename TMovingImage  >
void KLRegistration<TFixedImage,TMovingImage>
::SetMetricParam()
{
  this->m_Optimizer->MinimizeOn();
}

} // namespace itk


#endif /* __KLRegistration_txx */