FaceRemoval.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_mesh_simplicial__FaceRemoval_h__)
#define __geom_mesh_simplicial__FaceRemoval_h__

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

#include "../simplex/SimplexModCondNum.h"

#include "../../../ads/array/Array.h"

#include <cassert>

BEGIN_NAMESPACE_GEOM


template<class _QualityMetric = SimplexModCondNum<3>,
         class _Point = ads::FixedArray<3>,
         typename _Number = typename _Point::value_type>
class FaceRemoval {

public:

  //
  // Public types.
  //

  typedef _QualityMetric QualityMetric;
  typedef _Point Point;
  typedef _Number Number;

private:
      
  //
  // Private types.
  //

  typedef typename QualityMetric::Simplex Simplex;

  //
  // Data
  //

  // The source and target of the edge.
  Point _source, _target;
  // The three vertices in the face shared by the two initial tetrahedra.
  ads::FixedArray<3,Point> _face;
  // The quality function.
  mutable QualityMetric _qualityFunction;
  
  //
  // Not implemented.
  //

  // Copy constructor not implemented.
  FaceRemoval(const FaceRemoval&);

  // Assignment operator.
  FaceRemoval& 
  operator=(const FaceRemoval&);

public:

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

  FaceRemoval() :
    _source(),
    _target(),
    _face(),
    _qualityFunction()
  {}

  ~FaceRemoval()
  {}

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

  void
  setSource(const Point& src) {
    _source = src;
  }

  void
  setTarget(const Point& tgt) {
    _target = tgt;
  }


  void
  setFace(const Point& a, const Point& b, const Point& c) {
    _face[0] = a;
    _face[1] = b;
    _face[2] = c;
  }

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

  bool
  flip23() {
    return computeQuality3() > computeQuality2();
  }


  //
  // Private member functions.
  //
  // CONTINUE: Should these be private?
  
  // Return the worse quality of the 2 tetrahedra: 
  // _face[0], _face[1], _face[2], _target
  // and
  // _source, _face[0], _face[1], _face[2]
  Number
  computeQuality2() const;

  // Return the worst quality of the 3 tetrahedra: 
  // _source, _target, _face[0], _face[1]
  // _source, _target, _face[1], _face[2]
  // _source, _target, _face[2], _face[0]
  Number
  computeQuality3() const;
};


END_NAMESPACE_GEOM


#define __geom_mesh_simplicial_FaceRemoval_ipp__
#include "FaceRemoval.ipp"
#undef __geom_mesh_simplicial_FaceRemoval_ipp__

#endif

---