text: /*=========================================================================
   *
   *  Copyright UMC Utrecht and contributors
   *
   *  Licensed under the Apache License, Version 2.0 (the "License");
   *  you may not use this file except in compliance with the License.
   *  You may obtain a copy of the License at
   *
   *        http://www.apache.org/licenses/LICENSE-2.0.txt
   *
   *  Unless required by applicable law or agreed to in writing, software
   *  distributed under the License is distributed on an "AS IS" BASIS,
   *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
   *  See the License for the specific language governing permissions and
   *  limitations under the License.
   *
   *=========================================================================*/
  #ifndef __elxAdvancedMeanSquaresMetric_HXX__
  #define __elxAdvancedMeanSquaresMetric_HXX__
  
  #include "elxAdvancedMeanSquaresMetric.h"
  #include "itkTimeProbe.h"
  
  namespace elastix
  {
  
  /**
   * ******************* Initialize ***********************
   */
  
  template< class TElastix >
  void
  AdvancedMeanSquaresMetric< TElastix >
  ::Initialize( void ) throw ( itk::ExceptionObject )
  {
    itk::TimeProbe timer;
    timer.Start();
    this->Superclass1::Initialize();
    timer.Stop();
    elxout << "Initialization of AdvancedMeanSquares metric took: "
           << static_cast< long >( timer.GetMean() * 1000 ) << " ms." << std::endl;
  
  } // end Initialize()
  
  
  /**
   * ***************** BeforeEachResolution ***********************
   */
  
  template< class TElastix >
  void
  AdvancedMeanSquaresMetric< TElastix >
  ::BeforeEachResolution( void )
  {
    /** Get the current resolution level. */
    unsigned int level
      = ( this->m_Registration->GetAsITKBaseType() )->GetCurrentLevel();
  
    /** Get and set the normalization. */
    bool useNormalization = false;
    this->GetConfiguration()->ReadParameter( useNormalization,
      "UseNormalization", this->GetComponentLabel(), level, 0 );
    this->SetUseNormalization( useNormalization );
  
    /** Experimental options for SelfHessian */
  
    /** Set the number of samples used to compute the SelfHessian */
    unsigned int numberOfSamplesForSelfHessian = 100000;
    this->GetConfiguration()->ReadParameter( numberOfSamplesForSelfHessian,
      "NumberOfSamplesForSelfHessian", this->GetComponentLabel(), level, 0 );
    this->SetNumberOfSamplesForSelfHessian( numberOfSamplesForSelfHessian );
  
    /** Set the smoothing sigma used to compute the SelfHessian */
    double selfHessianSmoothingSigma = 1.0;
    this->GetConfiguration()->ReadParameter( selfHessianSmoothingSigma,
      "SelfHessianSmoothingSigma", this->GetComponentLabel(), level, 0 );
    this->SetSelfHessianSmoothingSigma( selfHessianSmoothingSigma );
  
    /** Set the smoothing sigma used to compute the SelfHessian */
    double selfHessianNoiseRange = 1.0;
    this->GetConfiguration()->ReadParameter( selfHessianNoiseRange,
      "SelfHessianNoiseRange", this->GetComponentLabel(), level, 0 );
    this->SetSelfHessianNoiseRange( selfHessianNoiseRange );
  
    /** Select the use of an OpenMP implementation for GetValueAndDerivative. */
    std::string useOpenMP = this->m_Configuration->GetCommandLineArgument( "-useOpenMP_SSD" );
    if( useOpenMP == "true" )
    {
      this->SetUseOpenMP( true );
    }
  
  } // end BeforeEachResolution()
  
  
  } // end namespace elastix
  
  #endif // end #ifndef __elxAdvancedMeanSquaresMetric_HXX__
  