gmBarCellGeometryInfo.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_BAR_CELL_GEOMETRY_INFO_H_
#define _GEMA_BAR_CELL_GEOMETRY_INFO_H_
 
#include "gmCellGeometryInfo.h"
#include "gmCellGeometryIntegrationRuleSet.h"
 
#include "gmGeometryUtils.h"
 
#include "gmBarIntegrationRule.h"
#include "gmBarShape.h"
 
//---------------------------------------------------------
// Bar integration rules
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSet<GmLineGaussIntegrationRule,  GmLineLobattoIntegrationRule,
                                         GmLineNewtonIntegrationRule, GmLineNewtonIntegrationRule> GmBarIntegrationRuleSet;
 
 
//---------------------------------------------------------
// Bar2 & Bar3D2 rule defaults
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 2, 2, 2, 2> GmLinearBarIntegrationRuleSetDefaults;
 
 
//---------------------------------------------------------
// Bar3 & Bar3D3 rule defaults
//---------------------------------------------------------
 
typedef GmCellGeometryIntegrationRuleSetDefaultRules<GM_GAUSS_RULE_TYPE, 3, 3, 3, 3> GmQuadraticBarIntegrationRuleSetDefaults;
 
 
//---------------------------------------------------------
// GmBar2CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmBar2CellGeometryInfo
  : public GmCellGeometryInfoLineElement<GmBarIntegrationRuleSet, GmLinearBarIntegrationRuleSetDefaults>
{
public:
  static const GmBar2CellGeometryInfo* 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 GmBar2Shape; }
 
  virtual double dimension(const GmMatrix& X) const { return GmGeometryUtils::edgeLength(X.col(0), X.col(1)); }
 
  virtual bool localAxis(const GmMatrix&X, GmMatrix& R) const {Q_UNUSED(X); R = arma::eye(1,1); return false;}
 
  virtual void centroidCartesian(const GmMatrix& X, GmVector& coord) const { GmGeometryUtils::edgeCentroid(X.col(0), X.col(1), coord); }
 
  // ---------  Capabilities API ---------
 
protected:
  GmBar2CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoLineElement<GmBarIntegrationRuleSet, GmLinearBarIntegrationRuleSetDefaults>(std::move(metadata)) {}
};
 
//---------------------------------------------------------
// GmBar3D2CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmBar3D2CellGeometryInfo : public GmBar2CellGeometryInfo
{
public:
  static const GmBar3D2CellGeometryInfo* 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 GmBar3D2Shape; }
 
private:
  GmBar3D2CellGeometryInfo(const GmBar2CellGeometryInfo& bar2) : GmBar2CellGeometryInfo(bar2) {}
};
 
 
//---------------------------------------------------------
// GmBar3CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmBar3CellGeometryInfo
  : public GmCellGeometryInfoLineElement<GmBarIntegrationRuleSet, GmQuadraticBarIntegrationRuleSetDefaults>
{
public:
  static const GmBar3CellGeometryInfo* 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 GmBar3Shape; }
 
  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(1,1); return false;}
 
  // ---------  Capabilities API ---------
 
protected:
  GmBar3CellGeometryInfo(GmCellGeometryMetadata&& metadata)
    : GmCellGeometryInfoLineElement<GmBarIntegrationRuleSet, GmQuadraticBarIntegrationRuleSetDefaults>(std::move(metadata)) {}
};
 
//---------------------------------------------------------
// GmBar3D3CellGeometryInfo
//---------------------------------------------------------
 
class GMC_API_EXPORT GmBar3D3CellGeometryInfo : public GmBar3CellGeometryInfo
{
public:
  static const GmBar3D3CellGeometryInfo* 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 GmBar3D3Shape; }
 
private:
  GmBar3D3CellGeometryInfo(const GmBar3CellGeometryInfo& bar3) : GmBar3CellGeometryInfo(bar3) {}
};
 
 
#endif