gmTriCellGeometryInfo.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_TRI_CELL_GEOMETRY_INFO_H_
#define _GEMA_TRI_CELL_GEOMETRY_INFO_H_
 
#include "gmCellGeometryInfo.h"
#include "gmCellGeometryIntegrationRuleSet.h"
 
#include "gmGeometryUtils.h"
 
#include "gmTriIntegrationRule.h"
#include "gmTriShape.h"
 
//---------------------------------------------------------
// Tri integration rules
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSet<GmTriGaussIntegrationRule,  GmCellGeometryNullIntegrationRule,
                                         GmTriNewtonIntegrationRule, GmTriNewtonIntegrationRule> GmTriIntegrationRuleSet;
 
typedef GmCellGeometryIntegrationRuleSet<GmTriGaussEdgeIntegrationRule,     GmCellGeometryNullIntegrationRule,
                                         GmCellGeometryNullIntegrationRule, GmCellGeometryNullIntegrationRule> GmTriEdgeIntegrationRuleSet;
 
 
//---------------------------------------------------------
// Tri3 & Tri3D3 rule defaults
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, 1, 1> GmLinearTriIntegrationRuleSetDefaults;
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, -1, -1, -1> GmLinearTriEdgeIntegrationRuleSetDefaults;
 
//---------------------------------------------------------
// Tri6 & Tri3D6 rule defaults
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 6, -1, 2, 2> GmQuadraticTriIntegrationRuleSetDefaults;
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, -1, -1, -1> GmQuadraticTriEdgeIntegrationRuleSetDefaults;
 
 
//---------------------------------------------------------
// GmTri3CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmTri3CellGeometryInfo
  : public GmCellGeometryInfoSurfaceElement<GmTriIntegrationRuleSet,     GmLinearTriIntegrationRuleSetDefaults,
                                            GmTriEdgeIntegrationRuleSet, GmLinearTriEdgeIntegrationRuleSetDefaults>
{
public:
  static const GmTri3CellGeometryInfo* 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 GmTri3Shape; }
 
  virtual double dimension(const GmMatrix& X) const 
  { 
    return GmGeometryUtils::triangleArea(X.col(0), X.col(1), X.col(2));
  }
 
  virtual void centroidCartesian(const GmMatrix& X, GmVector& coord) const 
  { 
    GmGeometryUtils::triangleCentroid(X.col(0), X.col(1), X.col(2), coord);
  }
 
  virtual bool localAxis(const GmMatrix&X, GmMatrix& R) const {Q_UNUSED(X); R = arma::eye(2, 2); return false; }
 
  // ---------  Capabilities API ---------
 
  virtual int isValid(const GmMatrix& X, double tol) const { 
      Q_UNUSED(tol);
      int ret = GM_GEOMETRY_STATE_VALID;
      if (!GmGeometryUtils::polygon2DIsCCW(X)) ret |= GM_GEOMETRY_STATE_BAD_ORIENTATION;
      if (GmGeometryUtils::polygon2DIsSelfIntersecting(X)) ret |= GM_GEOMETRY_STATE_SELF_INTERSECT;
      return ret;
  }
 
  // See comments on the base class
  virtual bool contains(const GmMatrix& X, const GmVector& coord) const { return GmGeometryUtils::pointInConvexPolygon(X, coord); }
 
protected:
  GmTri3CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSurfaceElement<GmTriIntegrationRuleSet,     GmLinearTriIntegrationRuleSetDefaults,
                                       GmTriEdgeIntegrationRuleSet, GmLinearTriEdgeIntegrationRuleSetDefaults>(std::move(metadata)) {}
};
 
//---------------------------------------------------------
// GmTri3D3CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmTri3D3CellGeometryInfo : public GmTri3CellGeometryInfo
{
public:
  static const GmTri3D3CellGeometryInfo* 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 GmTri3D3Shape; }
 
  // ---------  Capabilities API ---------
 
  virtual int isValid(const GmMatrix& X, double tol) const { 
      Q_UNUSED(tol);  
      int ret = GM_GEOMETRY_STATE_VALID;
      if (GmGeometryUtils::tri3DIsColinear(X)) ret |= GM_GEOMETRY_STATE_SELF_INTERSECT;
      return ret; 
  }
  // contains() capability is currently not supported in 3D, so the implementation
  // from GmTri3CellGeometryInfo must be overridden. Quality is ok.
  virtual bool contains(const GmMatrix& X, const GmVector& coord) const { Q_UNUSED(X); Q_UNUSED(coord); return false; }
 
  
  virtual bool localAxis(const GmMatrix&X, GmMatrix& R) const {Q_UNUSED(X); R = arma::eye(3, 3); return false; }
 
private:
  GmTri3D3CellGeometryInfo(const GmTri3CellGeometryInfo& tri3) : GmTri3CellGeometryInfo(tri3) {}
};
 
//---------------------------------------------------------
// GmTri6CellGeometryInfo
//---------------------------------------------------------
 
//tri6 nodes in circular order
static int tri6order[] = {0,3,1,4,2,5};
 
class GMC_API_EXPORT GmTri6CellGeometryInfo
  : public GmCellGeometryInfoSurfaceElement<GmTriIntegrationRuleSet,     GmQuadraticTriIntegrationRuleSetDefaults,
                                            GmTriEdgeIntegrationRuleSet, GmQuadraticTriEdgeIntegrationRuleSetDefaults>
{
public:
  static const GmTri6CellGeometryInfo* 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 GmTri6Shape; }
 
  virtual double dimension(const GmMatrix& X) const { return dimensionByIntegration(X); }
 
  virtual void centroidCartesian(const GmMatrix& X, GmVector& coord) const { centroidByIntegration(X, coord); }
 
  virtual bool localAxis(const GmMatrix&X, GmMatrix& R) const {Q_UNUSED(X); R = arma::eye(2, 2); return false; }
 
  // ---------  Capabilities API ---------
 
  virtual int isValid(const GmMatrix& X, double tol) const {
      Q_UNUSED(tol);
      GmMatrix Xr = GmGeometryUtils::shuffle(X, tri6order, 6);
      int ret = GM_GEOMETRY_STATE_VALID;
      if (!GmGeometryUtils::polygon2DIsCCW(Xr))              ret |= GM_GEOMETRY_STATE_BAD_ORIENTATION;
      if (GmGeometryUtils::polygon2DIsSelfIntersecting(Xr)) ret |= GM_GEOMETRY_STATE_SELF_INTERSECT;
      return ret;
  }
 
  // See comments on the base class
  virtual bool contains(const GmMatrix& X, const GmVector& coord) const { return GmGeometryUtils::pointInTesselatedPolygon_T6(X, coord); }
 
protected:
  GmTri6CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoSurfaceElement<GmTriIntegrationRuleSet,     GmQuadraticTriIntegrationRuleSetDefaults,
                                       GmTriEdgeIntegrationRuleSet, GmQuadraticTriEdgeIntegrationRuleSetDefaults>(std::move(metadata)) {}
};
 
//---------------------------------------------------------
// GmTri3D6CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmTri3D6CellGeometryInfo : public GmTri6CellGeometryInfo
{
public:
  static const GmTri3D6CellGeometryInfo* 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 GmTri3D6Shape; }
 
  
  virtual bool localAxis(const GmMatrix&X, GmMatrix& R) const {Q_UNUSED(X); R = arma::eye(3, 3); return false; }
 
  // ---------  Capabilities API ---------
  
  virtual int isValid(const GmMatrix& X, double tol) const { 
      int ret = GM_GEOMETRY_STATE_VALID;
      if (!GmGeometryUtils::polygon3DIsPlanar(GmGeometryUtils::shuffle(X, tri6order, 6), tol)) ret |= GM_GEOMETRY_STATE_NON_PLANAR;
      if (GmGeometryUtils::tri3DIsColinear(X.cols(0, 2)))                                      ret |= GM_GEOMETRY_STATE_SELF_INTERSECT;
      return ret;
  }
  
  // contains() capability is currently not supported in 3D, so the implementation
  // from GmTri6CellGeometryInfo must be overridden.
  virtual bool contains(const GmMatrix& X, const GmVector& coord) const { Q_UNUSED(X); Q_UNUSED(coord); return false; }
 
private:
  GmTri3D6CellGeometryInfo(const GmTri6CellGeometryInfo& tri6) : GmTri6CellGeometryInfo(tri6) {}
};
 
 
#endif