1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkSinRegularizedHeavisideStepFunctionTest1.cxx
Language: C++
Date: $Date$
Version: $Revision$
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.
=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include "itkSinRegularizedHeavisideStepFunction.h"
#include "itkTestingMacros.h"
int itkSinRegularizedHeavisideStepFunctionTest1( int, char* [] )
{
typedef double InputType;
typedef double OutputType;
typedef itk::SinRegularizedHeavisideStepFunction< InputType, OutputType > HeavisideFunctionBaseType;
HeavisideFunctionBaseType::Pointer functionBase0 = HeavisideFunctionBaseType::New();
std::cout << "GetNameOfClass() = " << functionBase0->GetNameOfClass() << std::endl;
functionBase0->Print( std::cout );
const double epsilon0 = 1.0;
const double epsilon1 = 1e-4;
TEST_SET_GET_VALUE( epsilon0, functionBase0->GetEpsilon() );
functionBase0->SetEpsilon( epsilon1 );
TEST_SET_GET_VALUE( epsilon1, functionBase0->GetEpsilon() );
const double epsilon2 = 0.5;
functionBase0->SetEpsilon( epsilon2 );
const signed int minValue = -20;
const signed int maxValue = 20;
const InputType incValue = 0.1;
for( signed int x = minValue; x < maxValue; x++ )
{
const InputType ix = x * incValue;
OutputType f = functionBase0->Evaluate( ix );
OutputType df = functionBase0->EvaluateDerivative( ix );
std::cout << ix << " " << f << " " << df << std::endl;
}
return EXIT_SUCCESS;
}
|