RegularPointEmbedder.ih

#include <cstdlib>
#include <cmath>
#include "DGtal/kernel/NumberTraits.h"

// IMPLEMENTATION of inline methods.

// ----------------------- Standard services ------------------------------

//-----------------------------------------------------------------------------
template <typename TSpace>
inline
DGtal::RegularPointEmbedder<TSpace>::~RegularPointEmbedder()
{
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
DGtal::RegularPointEmbedder<TSpace>::RegularPointEmbedder()
{}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
DGtal::RegularPointEmbedder<TSpace> &
DGtal::RegularPointEmbedder<TSpace>::
operator=( const RegularPointEmbedder & other )
{
  if ( this != &other )
    {
      myGridSteps = other.myGridSteps;
    }
  return *this;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
void 
DGtal::RegularPointEmbedder<TSpace>
::init( typename RealVector::Component gridStep )
{
  myGridSteps = RealVector::diagonal( gridStep );
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
void 
DGtal::RegularPointEmbedder<TSpace>
::init( const RealVector & aGridSteps )
{
  myGridSteps = aGridSteps;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::Point
DGtal::RegularPointEmbedder<TSpace>
::floor( const RealPoint & p ) const
{
  Point aPoint;
  for ( Dimension i = 0; i < Space::dimension; ++i )
    aPoint[ i ] = (Integer) ::floor( p[ i ] / myGridSteps[ i ] );
  return aPoint;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::Point
DGtal::RegularPointEmbedder<TSpace>
::ceil( const RealPoint & p ) const
{
  Point aPoint;
  for ( Dimension i = 0; i < Space::dimension; ++i )
    aPoint[ i ] = (Integer) ::ceil( p[ i ] / myGridSteps[ i ] );
  return aPoint;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::Point
DGtal::RegularPointEmbedder<TSpace>
::round( const RealPoint & p ) const
{
  Point aPoint;
  for ( Dimension i = 0; i < Space::dimension; ++i )
    aPoint[ i ] = (Integer) ::round( p[ i ] / myGridSteps[ i ] );
  return aPoint;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::RealPoint
DGtal::RegularPointEmbedder<TSpace>
::embed( const Point & p ) const
{
  return this->operator()( p );
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::RealPoint
DGtal::RegularPointEmbedder<TSpace>
::operator()( const Point & p ) const
{
  RealPoint aRealPoint;
  for ( Dimension i = 0; i < Space::dimension; ++i )
    aRealPoint[ i ] = NumberTraits<Integer>::castToDouble( p[ i ] ) 
      * myGridSteps[ i ];
  return aRealPoint;
}
//-----------------------------------------------------------------------------
template <typename TSpace>
inline
typename DGtal::RegularPointEmbedder<TSpace>::RealVector
DGtal::RegularPointEmbedder<TSpace>
::gridSteps() const
{
  return myGridSteps;
}


// Interface - public :

template <typename TSpace>
inline
void
DGtal::RegularPointEmbedder<TSpace>::selfDisplay ( std::ostream & out ) const
{
  out << "[RegularPointEmbedder]";
}

template <typename TSpace>
inline
bool
DGtal::RegularPointEmbedder<TSpace>::isValid() const
{
    return true;
}



// Implementation of inline functions                                        //

template <typename TSpace>
inline
std::ostream&
DGtal::operator<< ( std::ostream & out, 
      const RegularPointEmbedder<TSpace> & object )
{
  object.selfDisplay( out );
  return out;
}

//                                                                           //