gmWedgeCellGeometryInfo.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_WEDGE_CELL_GEOMETRY_INFO_H_
#define _GEMA_WEDGE_CELL_GEOMETRY_INFO_H_

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

#include "gmGeometryUtils.h"

#include "gmWedgeIntegrationRule.h"
#include "gmWedgeShape.h"

//---------------------------------------------------------
// Wedge integration rules
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSet<GmWedgeGaussIntegrationRule,       GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmWedgeIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmWedgeGaussFaceIntegrationRule,   GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmWedgeFaceIntegrationRuleSet;

typedef GmCellGeometryIntegrationRuleSet<GmWedgeGaussEdgeIntegrationRule,   GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmWedgeEdgeIntegrationRuleSet;

//---------------------------------------------------------
// Wedge6 rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDoubleParDefaultRules<GM_GAUSS_RULE_TYPE, 1, 2, -1, -1, -1, -1, -1, -1> GmLinearWedgeIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmLinearWedgeTriFaceIntegrationRuleSetDefaults;
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmLinearWedgeQuadFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetMultiFaceDefaultRules<GmLinearWedgeTriFaceIntegrationRuleSetDefaults, 
                                                              GmLinearWedgeQuadFaceIntegrationRuleSetDefaults> GmLinearWedgeFaceIntegrationRuleSetDefaults;

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

//---------------------------------------------------------
// Wedge15 rule defaults
//---------------------------------------------------------

typedef GmCellGeometryIntegrationRuleSetDoubleParDefaultRules<GM_GAUSS_RULE_TYPE, 3, 2, -1, -1, -1, -1, -1, -1> GmQuadraticWedgeIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmQuadraticWedgeTriFaceIntegrationRuleSetDefaults;
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmQuadraticWedgeQuadFaceIntegrationRuleSetDefaults;

typedef GmCellGeometryIntegrationRuleSetMultiFaceDefaultRules<GmQuadraticWedgeTriFaceIntegrationRuleSetDefaults,
                                                              GmQuadraticWedgeQuadFaceIntegrationRuleSetDefaults> GmQuadraticWedgeFaceIntegrationRuleSetDefaults;

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


//---------------------------------------------------------
// GmWedge6CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmWedge6CellGeometryInfo
  : public GmCellGeometryInfoSolidElement<GmWedgeIntegrationRuleSet,     GmLinearWedgeIntegrationRuleSetDefaults,
                                          GmWedgeEdgeIntegrationRuleSet, GmLinearWedgeEdgeIntegrationRuleSetDefaults,
                                          GmWedgeFaceIntegrationRuleSet, GmLinearWedgeFaceIntegrationRuleSetDefaults>
{
public:
  static const GmWedge6CellGeometryInfo* 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 GmWedge6Shape; }

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

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

  // ---------  Capabilities API ---------
  virtual bool isValid(const GmMatrix& X, double tol) const {
      return GmGeometryUtils::wedge6IsWellDefined(X, tol);
  }
private:
  GmWedge6CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmWedgeIntegrationRuleSet,     GmLinearWedgeIntegrationRuleSetDefaults,
                                     GmWedgeEdgeIntegrationRuleSet, GmLinearWedgeEdgeIntegrationRuleSetDefaults,
                                     GmWedgeFaceIntegrationRuleSet, GmLinearWedgeFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};

//---------------------------------------------------------
// GmWedge15CellGeometryInfo
//---------------------------------------------------------

class GMC_API_EXPORT GmWedge15CellGeometryInfo
  : public GmCellGeometryInfoSolidElement<GmWedgeIntegrationRuleSet,     GmQuadraticWedgeIntegrationRuleSetDefaults,
                                          GmWedgeEdgeIntegrationRuleSet, GmQuadraticWedgeEdgeIntegrationRuleSetDefaults,
                                          GmWedgeFaceIntegrationRuleSet, GmQuadraticWedgeFaceIntegrationRuleSetDefaults>
{
public:
  static const GmWedge15CellGeometryInfo* 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 GmWedge15Shape; }

  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::wedge15ValidEdges(X) && GmGeometryUtils::wedge6IsWellDefined(X, tol);
  }
private:
  GmWedge15CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSolidElement<GmWedgeIntegrationRuleSet,     GmQuadraticWedgeIntegrationRuleSetDefaults,
                                     GmWedgeEdgeIntegrationRuleSet, GmQuadraticWedgeEdgeIntegrationRuleSetDefaults,
                                     GmWedgeFaceIntegrationRuleSet, GmQuadraticWedgeFaceIntegrationRuleSetDefaults>(std::move(metadata)) {}
};


#endif