#include "itkShrinkImageFilter.h"
#include "itkIdentityTransform.h"
-#include "vnl/vnl_math.h"
-
namespace clitk
{
factors[idim] = this->GetSchedule()[ilevel][idim] /
this->GetSchedule()[ilevel+1][idim];
}
- variance[idim] = vnl_math_sqr( 0.5 *
+ variance[idim] = std::sqr( 0.5 *
static_cast<float>( factors[idim] ) );
if( factors[idim] != 1 )
{
if( factors[idim] > 1 )
{
oper->SetDirection( idim );
- oper->SetVariance( vnl_math_sqr( 0.5 *
+ oper->SetVariance( std::sqr( 0.5 *
static_cast<float>( factors[idim] ) ) );
oper->CreateDirectional();
radius[idim] = oper->GetRadius()[idim];
if( factors[idim] > 1 )
{
oper->SetDirection( idim );
- oper->SetVariance( vnl_math_sqr( 0.5 *
+ oper->SetVariance( std::sqr( 0.5 *
static_cast<float>( factors[idim] ) ) );
oper->CreateDirectional();
radius[idim] = oper->GetRadius()[idim];
for( idim = 0; idim < TInputImage::ImageDimension; idim++ )
{
oper->SetDirection(idim);
- oper->SetVariance( vnl_math_sqr( 0.5 * static_cast<float>(
+ oper->SetVariance( std::sqr( 0.5 * static_cast<float>(
this->GetSchedule()[refLevel][idim] ) ) );
oper->SetMaximumError( this->GetMaximumError() );
oper->CreateDirectional();