SimplexCondNum.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexCondNum_h__)
#define __geom_SimplexCondNum_h__

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

#include "SimplexAdjJacQF.h"

#include <limits>

BEGIN_NAMESPACE_GEOM


template<int N, typename T = double>
class SimplexCondNum :
  public SimplexAdjJacQF<N,T> {
private:

  typedef SimplexAdjJacQF<N,T> Base;

public:

  //
  // Public types.
  //

  typedef T Number;

  typedef typename Base::Vertex Vertex;

  typedef typename Base::Simplex Simplex;

  typedef typename Base::Matrix Matrix;

public:

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

  SimplexCondNum() :
    Base()
  {}

  SimplexCondNum(const SimplexCondNum& other) :
    Base(other)
  {}

  SimplexCondNum(const Simplex& s) :
    Base(s)
  {}

  SimplexCondNum& 
  operator=(const SimplexCondNum& other) {
    if (&other != this) {
      Base::operator=(other);
    }
    return *this;
  }
  
  ~SimplexCondNum()
  {}

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


  using Base::getMatrix;

  using Base::getGradientMatrix;

  using Base::getDeterminant;

  using Base::getGradientDeterminant;

  using Base::computeContent;

  using Base::computeGradientContent;

  using Base::getDimension;

  using Base::getAdjointMatrix;

  using Base::getAdjointGradientMatrix;

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

  using Base::setFunction;

  using Base::set;

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


  Number
  operator()() const;


  Number
  operator()(const Simplex& simplex) const {
    setFunction(simplex);
    return operator()();
  }


  void
  computeGradient(Vertex* gradient) const;


protected:

  Number
  computeFunction(Number snj, Number sna) const;

};

END_NAMESPACE_GEOM

#define __geom_SimplexCondNum_ipp__
#include "SimplexCondNum.ipp"
#undef __geom_SimplexCondNum_ipp__

#endif

---