SimplexCenter.h

Go to the documentation of this file.

/* Copyright Benoit Hudson 2008 */
#ifndef SimplexCenter_HEADER
#define SimplexCenter_HEADER

#ifdef HAVE_CONFIG_H
 #include <config.h>
#endif

#include <geometry/Circles.h>

template <class Simplex, unsigned ambient, unsigned topological>
struct SimplexCenter
#ifdef HAVE_LAPACK
{
  typename Geometry::CenterRadius<ambient> operator() (const Simplex& s) const {
    // this is excrutiatingly slow
    boost::array<Point, topological + 1> points;
    for(unsigned i = 0; i <= topological; i++) {
      points[i] = simplex[i]->toPoint();
    }
    return Geometry::circumcenter<ambient, topological>(points);
  }
}
#endif
;

template <class Simplex, unsigned ambient>
struct SimplexCenter<Simplex, ambient, 1>
{
  typename Geometry::CenterRadius<ambient> operator() (const Simplex& s) const {
    return Geometry::circumcenter_segment<ambient>(s[0]->toPoint(),
        s[1]->toPoint());
  }
};

template <class Simplex, unsigned ambient>
struct SimplexCenter<Simplex, ambient, 2>
{
  typename Geometry::CenterRadius<ambient> operator() (const Simplex& s) const {
    return Geometry::circumcenter_triangle<ambient>(s[0]->toPoint(),
        s[1]->toPoint(), s[2]->toPoint());
  }
};

template <class Simplex, unsigned ambient>
struct SimplexCenter<Simplex, ambient, 3>
{
  typename Geometry::CenterRadius<ambient> operator() (const Simplex& s) const {
    return Geometry::circumcenter_tetrahedron<ambient>(s[0]->toPoint(),
        s[1]->toPoint(), s[2]->toPoint(), s[3]->toPoint());
  }
};

#endif

---