halfspaceContainer.h

Go to the documentation of this file.

#ifndef HALFSPACECONTAINER_H
#define HALFSPACECONTAINER_H

#include <list>
#include <iostream>
using namespace std;

#include <typedefs.h>

template<typename HS, typename PT >
class halfspaceContainer 
{
public:

  HS halfspace;

  list<uint> index;

  vector< PT > const & pts;

  halfspaceContainer
  (
    HS const & halfspace_,
    vector< PT > const & pts_
  ) : halfspace(halfspace_), pts(pts_) {} 

  void isInsideOrOnBoundary(list<uint> const & target);

  void subtractfrom( list<uint> & target );

};

//---------------------------------------------------------
//  Implementation

template< typename HS, typename PT >
void halfspaceContainer<HS,PT>::subtractfrom
(
  list<uint> & target
)
{
  list<uint>::iterator i=target.begin();
  for ( ; i!=target.end(); ++i)
  {
    if (halfspace.isInsideOrOnBoundary(pts[*i]))
    {
      index.push_back(*i);
      i=target.erase(i);
    }
  }
}

template< typename HS, typename PT >
void halfspaceContainer<HS,PT>::isInsideOrOnBoundary
(
  list<uint> const & target
)
{
  list< uint >::const_iterator i=target.begin();
  list< uint >::const_iterator iend=target.end();
  for ( ; i!=iend; ++i)
  {
    if (halfspace.isInsideOrOnBoundary(pts[*i]))
      index.push_back(*i);
  }
}

#endif