gmPyraCellGeometryInfo.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_PYRA_CELL_GEOMETRY_INFO_H_
#define _GEMA_PYRA_CELL_GEOMETRY_INFO_H_

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

#include "gmGeometryUtils.h"

#include "gmPyraIntegrationRule.h"
#include "gmPyraShape.h"

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

typedef GmCellGeometryIntegrationRuleSet<GmPyraGaussIntegrationRule,        GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmPyraIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmPyraGaussFaceIntegrationRule,    GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmPyraFaceIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmPyraGaussEdgeIntegrationRule,    GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmPyraEdgeIntegrationRuleSet;

//---------------------------------------------------------
// Pyra5 rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 5, -1, -1, -1> GmLinearPyraIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmLinearPyraTriFaceIntegrationRuleSetDefaults;
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmLinearPyraQuadFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetMultiFaceDefaultRules<GmLinearPyraTriFaceIntegrationRuleSetDefaults,
                                                              GmLinearPyraQuadFaceIntegrationRuleSetDefaults> GmLinearPyraFaceIntegrationRuleSetDefaults;

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

//---------------------------------------------------------
// Pyra13  rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 5, -1, -1, -1> GmQuadraticPyraIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmQuadraticPyraTriFaceIntegrationRuleSetDefaults;
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmQuadraticPyraQuadFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetMultiFaceDefaultRules<GmQuadraticPyraTriFaceIntegrationRuleSetDefaults,
                                                              GmQuadraticPyraQuadFaceIntegrationRuleSetDefaults> GmQuadraticPyraFaceIntegrationRuleSetDefaults;

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


//---------------------------------------------------------
// GmPyra5CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmPyra5CellGeometryInfo 
  : public GmCellGeometryInfoSolidElement<GmPyraIntegrationRuleSet,     GmLinearPyraIntegrationRuleSetDefaults,
                                          GmPyraEdgeIntegrationRuleSet, GmLinearPyraEdgeIntegrationRuleSetDefaults,
                                          GmPyraFaceIntegrationRuleSet, GmLinearPyraFaceIntegrationRuleSetDefaults>
{
public:
  static const GmPyra5CellGeometryInfo* 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 GmPyra5Shape; }

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

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

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

  virtual bool isValid(const GmMatrix& X, double tol) const
  {
      return GmGeometryUtils::pyra5IsWellDefined(X, tol);
  }

private:
  GmPyra5CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmPyraIntegrationRuleSet,     GmLinearPyraIntegrationRuleSetDefaults,
                                     GmPyraEdgeIntegrationRuleSet, GmLinearPyraEdgeIntegrationRuleSetDefaults,
                                     GmPyraFaceIntegrationRuleSet, GmLinearPyraFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};

//---------------------------------------------------------
// GmPyra13CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmPyra13CellGeometryInfo
  : public GmCellGeometryInfoSolidElement<GmPyraIntegrationRuleSet,     GmQuadraticPyraIntegrationRuleSetDefaults,
                                          GmPyraEdgeIntegrationRuleSet, GmQuadraticPyraEdgeIntegrationRuleSetDefaults,
                                          GmPyraFaceIntegrationRuleSet, GmQuadraticPyraFaceIntegrationRuleSetDefaults>
{
public:
  static const GmPyra13CellGeometryInfo* 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 GmPyra13Shape; }

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

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

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

  virtual bool isValid(const GmMatrix& X, double tol) const
  {
      return GmGeometryUtils::pyra13ValidEdges(X) &&
          GmGeometryUtils::pyra5IsWellDefined(X, tol);
  }

private:
  GmPyra13CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmPyraIntegrationRuleSet,     GmQuadraticPyraIntegrationRuleSetDefaults,
                                     GmPyraEdgeIntegrationRuleSet, GmQuadraticPyraEdgeIntegrationRuleSetDefaults,
                                     GmPyraFaceIntegrationRuleSet, GmQuadraticPyraFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};

#endif