SimplexAdjJacQF.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexAdjJacQF_h__)
#define __geom_SimplexAdjJacQF_h__

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

#include "SimplexJacQF.h"
#include "SimplexAdjJac.h"

BEGIN_NAMESPACE_GEOM


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

  typedef SimplexJacQF<N,T> Base;
  typedef SimplexAdjJac<N,T> Adjoint;

public:

  //
  // Public types.
  //

  typedef T Number;

  typedef typename Base::Vertex Vertex;

  typedef typename Base::Simplex Simplex;

  typedef typename Base::Matrix Matrix;

protected:

  //
  // Member data.
  //

  Adjoint _adjoint;

public:

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

  SimplexAdjJacQF() :
    Base(),
    _adjoint()
  {}

  SimplexAdjJacQF(const SimplexAdjJacQF& other) :
    Base(other),
    _adjoint(other._adjoint)
  {}

  SimplexAdjJacQF(const Simplex& s) :
    Base(s),
    _adjoint(getMatrix())
  {}

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

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


  //
  // Inherited from SimplexJacQF.
  //

  using Base::getMatrix;

  using Base::getGradientMatrix;

  using Base::getDeterminant;

  using Base::getGradientDeterminant;

  using Base::computeContent;

  using Base::computeGradientContent;

  using Base::getDimension;

  //
  // From SimplexAdjJac.
  //

  const Matrix&
  getAdjointMatrix() const {
    return _adjoint.getMatrix();
  }

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

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

  void
  setFunction(const Simplex& s) {
    Base::setFunction(s);
    _adjoint.setFunction(getMatrix());
  }

  void
  set(const Simplex& s) {
    Base::set(s);
    _adjoint.set(getMatrix());
  }


  void
  setFunction(const geom::Simplex<N,ads::FixedArray<N+1,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

---