gmGaussCoordinates.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_GAUSSCOORDINATES_H_
#define _GEMA_GAUSSCOORDINATES_H_

#include "gmMatrix.h"
#include "gmIntegrationRule.h"
#include "gmElement.h"

class GmElementMesh;
class GmValueAccessor;

class GMC_API_EXPORT GmGaussCoordinates
{
public:
  int numPoints() const { assert(_ir);  return _ir->numPoints(); }

  void coord(int ip, GmVector& c) const
  {
    assert(_ir && _shape);
    GmVector ncrd;
    double   dummyWeight;
    _ir->integrationPoint(ip, ncrd, &dummyWeight);
    _shape->naturalToCartesian(ncrd, _X, c);
  }

  const GmShape* shape() const { return _shape; }

  const GmIntegrationRule* integrationRule() const { return _ir; }

  GmMatrix& coordMatrix() { return _X; }

protected:
  GmGaussCoordinates(const GmShape* shape, const GmIntegrationRule* ir)
    : _shape(shape), _ir(ir) {}

  const GmShape*           _shape; 
  const GmIntegrationRule* _ir;    
  GmMatrix                 _X;     
};

class GMC_API_EXPORT GmSingleGaussCoordinates : public GmGaussCoordinates
{
public:
  GmSingleGaussCoordinates(const GmIntegrationRule* ir, const GmElement* e, const GmValueAccessor* coordAc);
  GmSingleGaussCoordinates(const GmElementMesh* mesh, int ruleSet, const GmElement* e, const GmValueAccessor* coordAc);
};

class GMC_API_EXPORT GmMultiGaussCoordinates : public GmGaussCoordinates
{
public:
  GmMultiGaussCoordinates(const GmElementMesh* mesh, int ruleSet, const GmValueAccessor* coordAc);

  void setElement(const GmElement* e);

private:
  const GmElementMesh*   _mesh;     
  int                    _ruleSet;  
  const GmValueAccessor* _coordAc;  
  GmCellType             _lastType; 
};


#endif