SimplexJacQF.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexJacQF_h__)
#define __geom_SimplexJacQF_h__

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

#include "SimplexJac.h"

BEGIN_NAMESPACE_GEOM


template<int N, typename T = double>
class SimplexJacQF {
private:

  typedef SimplexJac<N,T> Jacobian;

public:

  //
  // Public types.
  //

  typedef T Number;

  typedef typename Jacobian::Vertex Vertex;

  typedef typename Jacobian::Simplex Simplex;

  typedef typename Jacobian::Matrix Matrix;

private:

  //
  // Member data.
  //

  Jacobian _jacobian;

public:

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

  SimplexJacQF() :
    _jacobian()
  {}

  SimplexJacQF(const SimplexJacQF& other) :
    _jacobian(other._jacobian)
  {}

  SimplexJacQF(const Simplex& s) :
    _jacobian(s)
  {}

  SimplexJacQF& 
  operator=(const SimplexJacQF& other) {
    if (&other != this) {
      _jacobian = other._jacobian;
    }
    return *this;
  }
  
  ~SimplexJacQF()
  {}

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

  const Matrix&
  getMatrix() const {
    return _jacobian.getMatrix();
  }

  const ads::FixedArray<N,Matrix>&
  getGradientMatrix() const {
    return _jacobian.getGradientMatrix();
  }

  Number
  getDeterminant() const {
    return _jacobian.getDeterminant();
  }

  const Vertex&
  getGradientDeterminant() const {
    return _jacobian.getGradientDeterminant();
  }

  Number
  computeContent() const {
    return _jacobian.computeContent();
  }

  void
  computeGradientContent(Vertex* grad) const {
    _jacobian.computeGradientContent(grad);
  }

  Vertex
  computeGradientContent() const {
    return _jacobian.computeGradientContent();
  }

  static
  int
  getDimension() {
    return N;
  }

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

  void
  setFunction(const Simplex& s) {
    _jacobian.setFunction(s);
  }

  void
  set(const Simplex& s) {
    _jacobian.set(s);
  }


  void
  setFunction(const geom::Simplex<N,ads::FixedArray<N+1,Number>,Number>& s) {
    Simplex t;
    projectToLowerDimension(s, &t);
    setFunction(t);
  }


  void
  setFunction(const geom::Simplex<N,ads::FixedArray<N+2,Number>,Number>& s) {
    Simplex t;
    projectToLowerDimension(s, &t);
    setFunction(t);
  }


  void
  set(const geom::Simplex<N,ads::FixedArray<N+1,Number>,Number>& s) {
    Simplex t;
    projectToLowerDimension(s, &t);
    set(t);
  }

};

END_NAMESPACE_GEOM

#endif

---