COBAGenericNaivePlane.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

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

//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
~COBAGenericNaivePlane()
{ // Nothing to do.
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
COBAGenericNaivePlane()
{ // Object is invalid
  _axesToErase.reserve( 3 );
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
COBAGenericNaivePlane( const COBAGenericNaivePlane & other )
  : myAxes( other.myAxes )
{
  for ( AxisConstIterator it = myAxes.begin(), itE = myAxes.end();
        it != itE; ++it )
    myComputers[ *it ] = other.myComputers[ *it ];
  _axesToErase.reserve( 3 );
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger> &
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
operator=( const COBAGenericNaivePlane & other )
{
  if ( this != &other )
    {
      myAxes = other.myAxes;
      for ( AxisConstIterator it = myAxes.begin(), itE = myAxes.end();
            it != itE; ++it )
        myComputers[ *it ] = other.myComputers[ *it ];
    }
  return *this;
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
DGtal::Dimension
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
active() const
{
  ASSERT( myAxes.size() > 0 );
  return myAxes[ 0 ];
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::MyIntegerComputer &
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
ic() const
{
  return myComputers[ active() ].ic();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
clear()
{
  myAxes.clear();
  for ( unsigned int i = 0; i < 3; ++i )
    {
      myAxes.push_back( i );
      myComputers[ i ].clear();
    }
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
init( InternalInteger diameter, 
      InternalInteger widthNumerator,
      InternalInteger widthDenominator )
{
  clear();
  for ( unsigned int i = 0; i < 3; ++i )
    myComputers[ i ].init( i, diameter, widthNumerator, widthDenominator );
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Size
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
size() const
{
  return myComputers[ active() ].size();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
empty() const
{
  return myComputers[ active() ].empty();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Size
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
max_size() const
{
  return myComputers[ active() ].max_size();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Size
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
maxSize() const
{
  return max_size();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Size
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
complexity() const
{
  return myComputers[ active() ].complexity();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
operator()( const Point & p ) const
{
  return myComputers[ active() ].operator()( p );
}

//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
extendAsIs( const Point & p )
{ 
  ASSERT( isValid() );
  unsigned int nbok = 0;
  for ( AxisConstIterator it = myAxes.begin(), itE = myAxes.end(); 
        it != itE; ++it )
    {
      nbok += myComputers[ *it ].operator()( p ) ? 1 : 0;
    }
  if ( nbok != 0 ) // at least one is ok.
    {
      for ( AxisIterator it = myAxes.begin(); it != myAxes.end(); )
        // cannot put end() in variable, since end() moves when
        // modifiying a vector.
        {
          bool ok = myComputers[ *it ].extendAsIs( p );
          if ( ! ok ) 
            it = myAxes.erase( *it );
          else
            ++it;
        }
      ASSERT( ! myAxes.empty() );
      return true;
    }
  return false;
}

//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
extend( const Point & p )
{
  ASSERT( isValid() );
  unsigned int nbok = 0;
  _axesToErase.clear();
  for ( AxisConstIterator axIt = myAxes.begin(), axItE = myAxes.end(); 
        axIt != axItE; ++axIt )
    {
      bool ok = myComputers[ *axIt ].extend( p );
      if ( ! ok ) _axesToErase.push_back( *axIt );
      else ++nbok;
    }
  if ( nbok != 0 ) // at least one is ok.
    { // if one is ok, we must remove ko ones from the list of active
      // axes.
      AxisIterator axIt = myAxes.begin();
      for ( unsigned int i = 0; i < _axesToErase.size(); ++i )
        {
          while ( *axIt != _axesToErase[ i ] ) ++axIt;
          axIt = myAxes.erase( axIt );
        }
      return true;
    }
  return false;
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
isExtendable( const Point & p ) const
{
  ASSERT( isValid() );
  unsigned int nbok = 0;
  for ( AxisConstIterator it = myAxes.begin(), itE = myAxes.end(); 
        it != itE; ++it )
    {
      nbok += myComputers[ *it ].isExtendable( p ) ? 1 : 0;
    }
  return nbok != 0;
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
template <typename TInputIterator>
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
extend( TInputIterator it, TInputIterator itE )
{
  BOOST_CONCEPT_ASSERT(( boost::InputIterator<TInputIterator> ));

  ASSERT( isValid() );
  unsigned int nbok = 0;
  _axesToErase.clear();
  for ( AxisConstIterator axIt = myAxes.begin(), axItE = myAxes.end(); 
        axIt != axItE; ++axIt )
    {
      bool ok = myComputers[ *axIt ].extend( it, itE );
      if ( ! ok ) _axesToErase.push_back( *axIt );
      else ++nbok;
    }
  if ( nbok != 0 ) // at least one is ok.
    { // if one is ok, we must remove ko ones from the list of active
      // axes.
      AxisIterator axIt = myAxes.begin();
      for ( unsigned int i = 0; i < _axesToErase.size(); ++i )
        {
          while ( *axIt != _axesToErase[ i ] ) ++axIt;
          axIt = myAxes.erase( axIt );
        }
      return true;
    }
  return false;
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
template <typename TInputIterator>
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
isExtendable( TInputIterator it, TInputIterator itE ) const
{
  BOOST_CONCEPT_ASSERT(( boost::InputIterator<TInputIterator> ));

  ASSERT( isValid() );
  unsigned int nbok = 0;
  for ( AxisConstIterator axIt = myAxes.begin(), axItE = myAxes.end(); 
        axIt != axItE; ++axIt )
    {
      nbok += myComputers[ *axIt ].isExtendable( it, itE ) ? 1 : 0;
    }
  return nbok != 0;
}


//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
template <typename Vector3D>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
getNormal( Vector3D & normal ) const
{
  myComputers[ active() ].getNormal( normal );
}
//-----------------------------------------------------------------------------
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
template <typename Vector3D>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
getUnitNormal( Vector3D & normal ) const
{
  myComputers[ active() ].getUnitNormal( normal );
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
getBounds( double & min, double & max ) const
{
  myComputers[ active() ].getBounds( min, max );
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
const typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Point &
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
minimalPoint() const
{
  return myComputers[ active() ].minimalPoint();
}
//-----------------------------------------------------------------------------
template <typename TSpace, typename TInternalInteger>
inline
const typename DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::Point &
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::
maximalPoint() const
{
  return myComputers[ active() ].maximalPoint();
}



// Interface - public :

template <typename TSpace, typename TInternalInteger>
inline
void
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::selfDisplay ( std::ostream & out ) const
{
  double min, max;
  double N[ 3 ];
  out << "[COBAGenericNaivePlane";
  for ( AxisConstIterator axIt = myAxes.begin(), axItE = myAxes.end(); 
        axIt != axItE; ++axIt )
    out << " " << myComputers[ *axIt ];
  out << " ]"; 
}

template <typename TSpace, typename TInternalInteger>
inline
bool
DGtal::COBAGenericNaivePlane<TSpace, TInternalInteger>::isValid() const
{
  return myComputers[ active() ].isValid();
}


// Internals

// Implementation of inline functions                                        //

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

//                                                                           //