SimplexModDet.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexModDet_h__)
#define __geom_SimplexModDet_h__

#if defined(DEBUG_geom) && !defined(DEBUG_SimplexModDet)
#define DEBUG_SimplexModDet
#endif

#include "../../defs.h"

#include <limits>

#include <cmath>
#include <cassert>

BEGIN_NAMESPACE_GEOM


template<typename T = double>
class SimplexModDet {
public:

  //
  // Public types.
  //

  typedef T Number;

private:

  //
  // Static member data.
  //

  // A little bigger than the machine epsilon.
  static Number _epsilon;

  //
  // Not implemented.
  //

  // Default constructor.
  SimplexModDet();

  // Copy constructor.
  SimplexModDet(const SimplexModDet&);

  // Assignment operator.
  SimplexModDet& 
  operator=(const SimplexModDet&);
  
  // Destructor.
  ~SimplexModDet();

public:

  //--------------------------------------------------------------------------

  static
  Number
  getEpsilon() {
    return _epsilon;
  }

  static
  Number
  getDelta(const Number minDeterminant) {
    if (minDeterminant < _epsilon) {
      const Number argument = _epsilon * (_epsilon - minDeterminant);
#ifdef DEBUG_geom
      assert(argument >= 0);
#endif
      return std::sqrt(argument);
    }
    return 0.0;
  }


  static
  Number
  getH(Number determinant, Number minDeterminant);

};

END_NAMESPACE_GEOM

#define __geom_SimplexModDet_ipp__
#include "SimplexModDet.ipp"
#undef __geom_SimplexModDet_ipp__

#endif

---