SimplexAdjJac.h

Go to the documentation of this file.

// -*- C++ -*-

#if !defined(__geom_SimplexAdjJac_h__)
#define __geom_SimplexAdjJac_h__

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

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

namespace geom {


template < std::size_t 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.
   std::tr1::array<Matrix, N> _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 std::tr1::array<Matrix, N>&
   getGradientMatrix() const {
      return _gradientMatrix;
   }

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

   void
   setFunction(const Matrix& jacobian);

   void
   set(const Matrix& jacobian);


};

} // namespace geom

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

#endif

---