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 elxTransformBendingEnergyPenaltyTerm_hxx
  #define elxTransformBendingEnergyPenaltyTerm_hxx
  
  #include "elxTransformBendingEnergyPenaltyTerm.h"
  #include "itkTimeProbe.h"
  
  namespace elastix
  {
  
  /**
   * ******************* Initialize ***********************
   */
  
  template <class TElastix>
  void
  TransformBendingEnergyPenalty<TElastix>::Initialize(void)
  {
    itk::TimeProbe timer;
    timer.Start();
    this->Superclass1::Initialize();
    timer.Stop();
    elxout << "Initialization of TransformBendingEnergy metric took: " << static_cast<long>(timer.GetMean() * 1000)
           << " ms." << std::endl;
  
  } // end Initialize()
  
  
  /**
   * ***************** BeforeEachResolution ***********************
   */
  
  template <class TElastix>
  void
  TransformBendingEnergyPenalty<TElastix>::BeforeEachResolution(void)
  {
    /** Get the current resolution level. */
    unsigned int level = (this->m_Registration->GetAsITKBaseType())->GetCurrentLevel();
  
    /** 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);
  
  } // end BeforeEachResolution()
  
  
  } // end namespace elastix
  
  #endif // end #ifndef elxTransformBendingEnergyPenaltyTerm_hxx
  