partitionspace.h

Go to the documentation of this file.

#ifndef PARTITIONSPACE_H
#define PARTITIONSPACE_H

#include <cassert>
using namespace std;

#include <typedefs.h>

template< typename PT >
class partitionspace
{
public:

  virtual boolc isInside(PT const & x) const = 0; 

  virtual boolc intersection
  (
    PT & x, 
    PT const & poutside, 
    PT const & pinside
  ) const
    { assert(false); return false; }

  virtual boolc isOnBoundary(PT const & w) const
    { assert(false); return false; }

  virtual ~partitionspace() {}

  enum ptclassification { undefined, below, on, above };
  ptclassification classify(PT const & w) const
  { 
    if (this->isOnBoundary(w))
      return on;

    if (this->isInside(w))
      return above;

    return below;
  }
  static stringc classifystring(ptclassification c)
  {
    string s;
    switch (c)
    {
      case undefined: s="undefined"; break;
      case below: s="below"; break;
      case on: s="on"; break;
      case above: s="above"; break;
    }

    return s;
  } 

  void classify
  (
    ptclassification* vc, 
    PT const* beg,
    PT const* end
  ) const;

};

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

template< typename PT >
void partitionspace<PT>::classify
(
  ptclassification* vc, 
  PT const* beg,
  PT const* end
) const
{
  PT const* i(beg);
  for (; i!=end; ++i)
  {
    *vc = classify(*i);
    ++vc;
  }
}


#endif