delaunaysimpleD2.h

Go to the documentation of this file.

#ifndef DELAUNAYSIMPLED2_H
#define DELAUNAYSIMPLED2_H

#include <vector>
using namespace std;

#include <point.h>
#include <simplexD2indexed.h>
#include <triangle.h>
#include <typedefs.h>

template< typename PT, typename PD, typename INDX >
class delaunaysimpleD2
{
  void evaluateTriangle
  (
    INDX const i, 
    INDX const j, 
    INDX const k
  );
public:

  delaunaysimpleD2();

  vector< PT > pts;
  vector< simplexD2indexed<INDX> > vi;

  void get( PT & a, PT & b, PT & c, INDX const i) const;

  void eval();

};


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

template< typename PT, typename PD, typename INDX >
delaunaysimpleD2<PT,PD,INDX>::delaunaysimpleD2()
{
  pts.push_back(PT());
  vi.push_back( simplexD2indexed<INDX>() );
}

template< typename PT, typename PD, typename INDX >
void delaunaysimpleD2<PT,PD,INDX>::evaluateTriangle
(
  INDX const i, 
  INDX const j, 
  INDX const k
)
{
  triangle<PT,PD> t;
  t.constructUnordered( pts[i],pts[j],pts[k]);

  PT p0;
  PD radius;
  t.outercircle(radius,p0);

//cout << SHOW(i) << " " << SHOW(j) << " " << SHOW(k) << endl;

  INDX const n=pts.size();
  for (uint w=1; w<n; ++w)
  {
    if (w==i)
      continue;
    if (w==j)
      continue;
    if (w==k)
      continue;

    if ((p0-pts[w]).dot() < radius*radius )
      return;
  }

  vi.push_back( simplexD2indexed<INDX>(i,j,k) );
}

template< typename PT, typename PD, typename INDX >
void delaunaysimpleD2<PT,PD,INDX>::eval()
{
  INDX const n=pts.size();
  INDX i,j,k;
  for (i=1; i<n-2; ++i)
  {
    for (j=i+1; j<n-1; ++j)
    {
      for (k=j+1; k<n; ++k)
      {
//cout << SHOW(i) << " " << SHOW(j) << " " << SHOW(k) << endl;
        evaluateTriangle(i,j,k);
      }
    }
  }

  //cout << print(vi,"\n") << endl;
}

template< typename PT, typename PD, typename INDX >
void delaunaysimpleD2<PT,PD,INDX>::get
( 
  PT & a, 
  PT & b, 
  PT & c, 
  INDX const i
) const
{
  simplexD2indexed<INDX> const & t(vi[i]);
  
  a = pts[t.pi[0]];
  b = pts[t.pi[1]];
  c = pts[t.pi[2]];
}



#endif