ChordGenericStandardPlaneComputer.h

#pragma once

#if defined(ChordGenericStandardPlaneComputer_RECURSES)
#error Recursive header files inclusion detected in ChordGenericStandardPlaneComputer.h
#else // defined(ChordGenericStandardPlaneComputer_RECURSES)

#define ChordGenericStandardPlaneComputer_RECURSES

#if !defined ChordGenericStandardPlaneComputer_h

#define ChordGenericStandardPlaneComputer_h

// Inclusions
#include <iostream>
#include <set>
#include "DGtal/base/Common.h"
#include "DGtal/base/IteratorAdapter.h"
#include "DGtal/kernel/CInteger.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/arithmetic/IntegerComputer.h"
#include "DGtal/geometry/surfaces/ChordNaivePlaneComputer.h"

namespace DGtal
{

  // template class ChordGenericStandardPlaneComputer
  template < typename TSpace,
             typename TInputPoint,
             typename TInternalScalar >
  class ChordGenericStandardPlaneComputer
  {
    BOOST_CONCEPT_ASSERT(( concepts::CSpace< TSpace > ));
    BOOST_CONCEPT_ASSERT(( concepts::CSignedNumber< TInternalScalar > ));
    BOOST_STATIC_ASSERT(( TSpace::dimension == 3 ));
    BOOST_STATIC_ASSERT(( TInputPoint::dimension == 3 ));

    // ----------------------- public types ------------------------------
  public:
    typedef TSpace Space;
    typedef TInputPoint InputPoint;
    typedef TInternalScalar InternalScalar;
    typedef InputPoint InputVector;
    typedef ChordNaivePlaneComputer< Space, InputPoint, InternalScalar > ChordComputer;
    typedef typename ChordComputer::Component Component;
    typedef typename ChordComputer::Coordinate Coordinate;
    typedef typename ChordComputer::InternalVector InternalVector;
    typedef typename ChordComputer::Point Point;
    typedef typename ChordComputer::Primitive Primitive;
    typedef typename ChordComputer::InputPointSet InputPointSet;
    typedef typename ChordComputer::Size Size;
    typedef typename ChordComputer::ConstIterator InputPointSetConstIterator;
    typedef typename ChordComputer::Iterator InputPointSetIterator;

    // ----------------------- internal types ------------------------------
  private:
    typedef std::vector<Dimension>::iterator OrthantIterator;
    typedef std::vector<Dimension>::const_iterator OrthantConstIterator;


    // --------- adapter types for standard to naive transformation --------
  public:

    struct Transform {
      typedef InputPoint Value;
      typedef InputPoint value_type;
      bool _posX; 
      bool _posY; 
      inline Transform() {}
      inline Transform( bool posX, bool posY ) : _posX( posX ), _posY( posY ) {}
      inline InputPoint operator()( const InputPoint & p ) const
      { return InputPoint( p[ 0 ] + ( _posX ? p[ 2 ] : -p[ 2 ] ),
                           p[ 1 ] + ( _posY ? p[ 2 ] : -p[ 2 ] ),
                           p[ 2 ] ); }
      inline bool isPlusZOnX() const { return _posX; }
      inline bool isPlusZOnY() const { return _posY; }
      inline InputPoint inverse( const InputPoint & p ) const
      { return InputPoint( p[ 0 ] + ( _posX ? -p[ 2 ] : p[ 2 ] ),
                           p[ 1 ] + ( _posY ? -p[ 2 ] : p[ 2 ] ),
                           p[ 2 ] ); }
      template <typename Vector3D>
      void transformBack( Vector3D & v )
      {
        v[ 2 ] += ( _posX ? v[ 0 ] : -v[ 0 ] );
        v[ 2 ] += ( _posY ? v[ 1 ] : -v[ 1 ] );
      }
      template <typename Vector3D>
      void transform( Vector3D & v )
      {
        v[ 2 ] += ( _posX ? -v[ 0 ] : v[ 0 ] );
        v[ 2 ] += ( _posY ? -v[ 1 ] : v[ 1 ] );
      }
    };

    typedef boost::transform_iterator<Transform,InputPointSetConstIterator,
                                      InputPoint,InputPoint>
    ConstIterator;

    // ----------------------- std public types ------------------------------
  public:
    typedef ConstIterator const_iterator;
    typedef const InputPoint* const_pointer;
    typedef const InputPoint& const_reference;
    typedef InputPoint value_type;
    typedef typename InputPointSet::difference_type difference_type;
    typedef typename InputPointSet::size_type size_type;

    // ----------------------- Standard services ------------------------------
  public:

    ~ChordGenericStandardPlaneComputer();

    ChordGenericStandardPlaneComputer();

    ChordGenericStandardPlaneComputer ( const ChordGenericStandardPlaneComputer & other );

    ChordGenericStandardPlaneComputer & operator= ( const ChordGenericStandardPlaneComputer & other );

    Dimension active() const;

    void clear();

    void init( InternalScalar widthNumerator = NumberTraits< InternalScalar >::ONE,
               InternalScalar widthDenominator = NumberTraits< InternalScalar >::ONE );

    //-------------------- model of ForwardContainer -----------------------------
  public:

    Size size() const;

    bool empty() const;

    ConstIterator begin() const;

    ConstIterator end() const;

    Size max_size() const;

    Size maxSize() const;


    //-------------------- model of concepts::CPointPredicate -----------------------------
  public:

    bool operator()( const Point & p ) const;

    //-------------------- model of CIncrementalPrimitiveComputer -----------------------------
  public:

    bool extendAsIs( const InputPoint & p );

    bool extend( const InputPoint & p );

    bool isExtendable( const InputPoint & p ) const;

    //-------------------- model of CAdditivePrimitiveComputer -----------------------------
  public:

    template <typename TInputIterator>
    bool extend( TInputIterator it, TInputIterator itE );

    template <typename TInputIterator>
    bool isExtendable( TInputIterator it, TInputIterator itE ) const;

    //-------------------- Primitive services -----------------------------
  public:

    Primitive primitive() const;

    //-------------------- Parameters services -----------------------------
  public:

    template <typename Vector3D>
    void getNormal( Vector3D & normal ) const;

    template <typename Vector3D>
    void getUnitNormal( Vector3D & normal ) const;

    void getBounds( double & min, double & max ) const;

    InputPoint minimalPoint() const;

    InputPoint maximalPoint() const;

    void getCharacteristics( InternalVector & n,
                             InternalScalar & imin, InternalScalar & imax,
                             InputPoint & p_min, InputPoint & p_max ) const;

    // ----------------------- Interface --------------------------------------
  public:

    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;

    // ------------------------- Private Datas --------------------------------
  private:
    std::vector<Dimension> myOrthants; 
    ChordComputer myComputers[ 4 ];    
    static Transform myTransforms[ 4 ];
    mutable std::vector<Dimension> _orthantsToErase; 
    // ------------------------- Hidden services ------------------------------
  protected:
    static Transform t( Dimension orthant );
    static Transform invT( Dimension orthant );

    // ------------------------- Internals ------------------------------------
  private:
  }; // end of class ChordGenericStandardPlaneComputer


  template <typename TSpace, typename TInputPoint, typename TInternalScalar>
  std::ostream&
  operator<< ( std::ostream & out,
               const ChordGenericStandardPlaneComputer<TSpace, TInputPoint, TInternalScalar> & object );

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/surfaces/ChordGenericStandardPlaneComputer.ih"

//                                                                           //

#endif // !defined ChordGenericStandardPlaneComputer_h

#undef ChordGenericStandardPlaneComputer_RECURSES
#endif // else defined(ChordGenericStandardPlaneComputer_RECURSES)