P2TetrahedronFiniteElement.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stéphane Del Pino

//  This program is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2, or (at your option)
//  any later version.

//  This program is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.

//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  

//  $Id: P2TetrahedronFiniteElement.hpp,v 1.2 2006/05/20 15:22:46 delpinux Exp $

#ifndef P2_TETRAHEDRON_FINITE_ELEMENT_HPP
#define P2_TETRAHEDRON_FINITE_ELEMENT_HPP

#include <TinyVector.hpp>
#include <TinyMatrix.hpp>

#include <QuadratureFormula.hpp>
#include <ConformTransformation.hpp>

#include <LagrangianFiniteElement.hpp>
#include <ThreadStaticBase.hpp>

class P2TetrahedronFiniteElement
  : public ThreadStaticBase<P2TetrahedronFiniteElement>,
    public LagrangianFiniteElement<10, P2TetrahedronFiniteElement,
                                   QuadratureFormulaP2Tetrahedron>
{
private:
  static TinyVector<3, real_t> __massCenter; 
public:
  enum {
    numberOfDegreesOfFreedom = 10,
    numberOfVertexDegreesOfFreedom = 1,
    numberOfEdgeDegreesOfFreedom = 1,
    numberOfFaceDegreesOfFreedom = 0,
    numberOfVolumeDegreesOfFreedom = 0,
    numberOfFaceLivingDegreesOfFreedom = 6 // degrees of freedom carried by a face
  };

  static const size_t facesDOF[Tetrahedron::NumberOfFaces][numberOfFaceLivingDegreesOfFreedom];

  static const TinyVector<3, real_t>& massCenter()
  {
    return __massCenter;
  }

  real_t W  (const size_t& i, const TinyVector<3>& X) const;
  real_t dxW(const size_t& i, const TinyVector<3>& X) const;
  real_t dyW(const size_t& i, const TinyVector<3>& X) const;
  real_t dzW(const size_t& i, const TinyVector<3>& X) const;

  inline const
  TinyVector<QuadratureType::numberOfQuadraturePoints,
             TinyVector<3> >&
  integrationVertices() const
  {
    return QuadratureType::instance().vertices();
  }

  P2TetrahedronFiniteElement()
  {
    ;
  }

  ~P2TetrahedronFiniteElement()
  {
    ;
  }
};

#endif // P2_TETRAHEDRON_FINITE_ELEMENT_HPP

---