gmShapeH.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_SHAPEH_H_
#define _GEMA_SHAPEH_H_

#include "gmShape.h"

class GMC_API_EXPORT GmHShape : public GmShape
{
public:

  GmHShape(int P, int Q) : _POrder(P), _QOrder(Q) { assert(P > 0 && Q > 0); }

  int POrder() const { return _POrder; }

  int QOrder() const { return _QOrder; }

  // See comments on the base class
  virtual int numFunctions()  const { return numPFunctions() + numQFunctions(); }
  virtual int numPFunctions() const { return _nPFunctions; }
  virtual int numQFunctions() const { return _nQFunctions; }

  virtual int numPVertexFunctions() const { return _nPVertexFunctions; }
  virtual int numPEdgeFunctions()   const { return _nPEdgeFunctions;   }
  virtual int numPFaceFunctions()   const { return _nPFaceFunctions;   }
  virtual int numPBubbleFunctions() const { return _nPBubbleFunctions; }

  virtual int numQEdgeFunctions()   const { return _nQEdgeFunctions;   }
  virtual int numQFaceFunctions()   const { return _nQFaceFunctions;   }
  virtual int numQBubbleFunctions() const { return _nQBubbleFunctions; }

  virtual int numPEdgeFunctionsEdge() const { return _nPEdgeFunctionsEdge;};         //<! Number of Pressure Edge Functions associated for edge
  virtual int numPFaceFunctionsFace() const { return _nPFaceFunctionsFace;};         //<! Number of Pressure Face Functions associated for face.
  virtual int numQEdgeFunctionsEdge() const { return _nQEdgeFunctionsEdge;};             //<! Number of Flux Edge Functions associated for edge
  virtual int numQFaceFunctionsFace() const { return _nQFaceFunctionsFace;};             //<! Number of Flux Face Functions associated for face.

  virtual void shapeValuesVec(const GmVector& ncoord, GmMatrix &N, bool transposed = false) const = 0;
  virtual void shapeValuesVecBorder(const int border, const GmVector& ncoord, GmMatrix& N, bool transposed = false) const = 0;

  virtual void shapePartialsVec(const GmVector& ncoord, GmMatrix& N) const = 0;

  // See comments on the base class
  virtual double scaledJacobianDet(const GmMatrix& J) const { return arma::det(J); }

  virtual void jacobianAndPartials(const GmVector& ncoord, const GmMatrix& X, GmMatrix& J, GmMatrix& P, bool transposed = false) const;

  void fillNaturalCoordinates(GmMatrix& coord, bool transposed = false) const;
  virtual void naturalToCartesian(const GmVector& ncoord, const GmMatrix& X, GmVector& coord) const;

protected:
  // See comments on the base class. At the moment we don't have a working version for H elements
  // What do we need to change from standard quad / hex methods?
  virtual bool gradientBasedCartesianToNatural(const GmVector& coord, const GmMatrix& X, double tol,
                                               int maxIter, double natTol, GmVector& ncoord, bool* inside) const
  {
    Q_UNUSED(coord); Q_UNUSED(X); Q_UNUSED(tol); Q_UNUSED(maxIter); Q_UNUSED(natTol); Q_UNUSED(ncoord); Q_UNUSED(inside);
    return false;
  }

  int _POrder; 
  int _QOrder; 

  int _nPVertexFunctions;
  int _nPEdgeFunctions;
  int _nPFaceFunctions;
  int _nPBubbleFunctions;
  int _nPEdgeFunctionsEdge;
  int _nPFaceFunctionsFace;
  int _nPFunctions;

  int _nQEdgeFunctions;
  int _nQFaceFunctions;
  int _nQBubbleFunctions;
  int _nQEdgeFunctionsEdge;
  int _nQFaceFunctionsFace;
  int _nQFunctions;
};


#endif