An Object Finite Element Library

Tetra4.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   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; Version 2 of the License.

   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

  ==============================================================================

      Definition of class 'Tetra4' for shape functions of 4-node tetrahedron

  ==============================================================================*/


#ifndef __TETRA4_H
#define __TETRA4_H

#include <math.h>
#include <float.h>
#include <stdlib.h>
#include <assert.h>

#include "OFELI_Config.h"
#include "FEShape.h"
#include "LocalMatrix.h"


namespace OFELI {

#ifndef DOXYGEN_SHOULD_SKIP_THIS
struct ErrorInTetra4 {
   void Message(const char *file, size_t line, int code, int p1) const;
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

class Tetra4 : public FEShape
{

public :

    Tetra4();

    Tetra4(const Element *element);

    ~Tetra4()
    {
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Tetra4 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    double Sh(size_t i, Point<double> s) const;

    Point<double> DSh(size_t i) const { return _dsh[i-1]; }

    double getVolume() const { return OFELI_SIXTH*_det; }

    double getMaxEdgeLength() const;

    double getMinEdgeLength() const;

    Point<double> getRefCoord(const Point<double> &x) const;

    bool isIn(const Point<double> &x);

    double getInterpolate(const Point<double> &x, const LocalVect<double,4> &v);

    double check() const;

 private :

   void          CalculateShape();
   const Element *_el;
   ErrorInTetra4 _e;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani