for( unsigned int d = 0; d < I::ImageDimension; d++ )
dist += std::abs( long( a[ d ] ) - long( b[ d ] ) );
if( this->m_NeighborhoodOrder == 1 )
- return( dist == 0 || dist == 1 );
+ return( dist <= 1 );
else
- return( dist == 0 || dist == 1 || dist == 2 );
+ return( dist <= I::ImageDimension );
}
// -------------------------------------------------------------------------