SimplexAdjJac.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexAdjJac_h__)
#define __geom_SimplexAdjJac_h__

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

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

#include "../../../ads/tensor/SquareMatrix.h"

BEGIN_NAMESPACE_GEOM


template<int N, typename T = double>
class SimplexAdjJac
{
public:

  //
  // Public types.
  //

  typedef T Number;

  typedef ads::SquareMatrix<N,Number> Matrix;

private:

  //
  // Member data.
  //

  // The adjoint (scaled inverse) of the Jacobian matrix.
  Matrix _matrix;

  // The gradient of the scaled inverse of the Jacobian matrix.
  ads::FixedArray<N,Matrix> _gradientMatrix;

public:

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

  SimplexAdjJac() :
    _matrix(),
    _gradientMatrix()
  {}

  SimplexAdjJac(const SimplexAdjJac& other) :
    _matrix(other._matrix),
    _gradientMatrix(other._gradientMatrix)
  {}

  SimplexAdjJac(const Matrix& jacobian) {
    set(jacobian);
  }

  SimplexAdjJac& 
  operator=(const SimplexAdjJac& other) {
    if (&other != this) {
      _matrix = other._matrix;
      _gradientMatrix = other._gradientMatrix;
    }
    return *this;
  }
  
  ~SimplexAdjJac()
  {}

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

  const Matrix&
  getMatrix() const {
    return _matrix;
  }

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

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

  void
  setFunction(const Matrix& jacobian);

  void
  set(const Matrix& jacobian);


};

END_NAMESPACE_GEOM

#define __geom_SimplexAdjJac_ipp__
#include "SimplexAdjJac.ipp"
#undef __geom_SimplexAdjJac_ipp__

#endif

---