gmTetCellGeometryInfo.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GEMA_TET_CELL_GEOMETRY_INFO_H_
#define _GEMA_TET_CELL_GEOMETRY_INFO_H_

#include "gmCellGeometryInfo.h"
#include "gmCellGeometryIntegrationRuleSet.h"

#include "gmGeometryUtils.h"

#include "gmTetIntegrationRule.h"
#include "gmTetShape.h"

//---------------------------------------------------------
// Tet integration rules
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSet<GmTetGaussIntegrationRule,  GmCellGeometryNullIntegrationRule,
                                         GmTetNewtonIntegrationRule, GmTetNewtonIntegrationRule> GmTetIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmTetGaussFaceIntegrationRule,     GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmTetFaceIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmTetGaussEdgeIntegrationRule,     GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmTetEdgeIntegrationRuleSet;

//---------------------------------------------------------
// Tet4 rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 4, -1, 1, 1> GmLinearTetIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmLinearTetFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmLinearTetEdgeIntegrationRuleSetDefaults;

//---------------------------------------------------------
// Tet10  rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 4, -1, 2, 2> GmQuadraticTetIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmQuadraticTetFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmQuadraticTetEdgeIntegrationRuleSetDefaults;


//---------------------------------------------------------
// GmTet4CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmTet4CellGeometryInfo 
  : public GmCellGeometryInfoSolidElement<GmTetIntegrationRuleSet,     GmLinearTetIntegrationRuleSetDefaults,
                                          GmTetEdgeIntegrationRuleSet, GmLinearTetEdgeIntegrationRuleSetDefaults,
                                          GmTetFaceIntegrationRuleSet, GmLinearTetFaceIntegrationRuleSetDefaults>
{
public:
  static const GmTet4CellGeometryInfo* instance();

  // Shape factory. See comments on the base class.
  virtual GmShape* shapeInstance(int P, int Q) const { Q_UNUSED(P); Q_UNUSED(Q); return new GmTet4Shape; }

  virtual double dimension(const GmMatrix& X) const 
  { 
    return GmGeometryUtils::tetrahedronVolume(X.col(0), X.col(1), X.col(2), X.col(3));
  }

  virtual void centroidCartesian(const GmMatrix& X, GmVector& coord) const 
  { 
    GmGeometryUtils::tetrahedronCentroid(X.col(0), X.col(1), X.col(2), X.col(3), coord);
  }

  // ---------  Capabilities API ---------

  // See comments on the base class
  virtual bool isValid(const GmMatrix& X, double tol) const
  {
    Q_UNUSED(tol);
    return GmGeometryUtils::tet4IsWellDefined(X.col(0), X.col(1), X.col(2), X.col(3));
  }

  // See comments on the base class
  virtual bool contains(const GmMatrix& X, const GmVector& coord) const { return GmGeometryUtils::pointInTetrahedron(X, coord); }

private:
  GmTet4CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmTetIntegrationRuleSet,     GmLinearTetIntegrationRuleSetDefaults,
                                     GmTetEdgeIntegrationRuleSet, GmLinearTetEdgeIntegrationRuleSetDefaults,
                                     GmTetFaceIntegrationRuleSet, GmLinearTetFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};

//---------------------------------------------------------
// GmTet10CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmTet10CellGeometryInfo
  : public GmCellGeometryInfoSolidElement<GmTetIntegrationRuleSet,     GmQuadraticTetIntegrationRuleSetDefaults,
                                          GmTetEdgeIntegrationRuleSet, GmQuadraticTetEdgeIntegrationRuleSetDefaults,
                                          GmTetFaceIntegrationRuleSet, GmQuadraticTetFaceIntegrationRuleSetDefaults>
{
public:
  static const GmTet10CellGeometryInfo* instance();

  // Shape factory. See comments on the base class.
  virtual GmShape* shapeInstance(int P, int Q) const { Q_UNUSED(P); Q_UNUSED(Q); return new GmTet10Shape; }

  virtual double dimension(const GmMatrix& X) const { return dimensionByIntegration(X); }

  virtual void centroidCartesian(const GmMatrix& X, GmVector& coord) const { centroidByIntegration(X, coord); }

  // ---------  Capabilities API ---------
  //For now same check as Tet4
  virtual bool isValid(const GmMatrix& X, double tol) const
  {
      Q_UNUSED(tol);
      return GmGeometryUtils::tet10ValidEdges(X) &&
             GmGeometryUtils::tet4IsWellDefined(X.col(0), X.col(1), X.col(2), X.col(3));
  }

private:
  GmTet10CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmTetIntegrationRuleSet,     GmQuadraticTetIntegrationRuleSetDefaults,
                                     GmTetEdgeIntegrationRuleSet, GmQuadraticTetEdgeIntegrationRuleSetDefaults,
                                     GmTetFaceIntegrationRuleSet, GmQuadraticTetFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};

#endif