An Object Finite Element Library

FEShape.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 parent class 'FEShape' from which derive all finite 
                               element shape classes

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


#ifndef __FE_SHAPE_H
#define __FE_SHAPE_H

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

#include "OFELI_Config.h"
#include "Element.h"
#include "Side.h"
#include "LocalMatrix.h"
#include "LocalVect.h"
#include "Point.h"
#include "util.h"
#include "Gauss.h"


namespace OFELI {

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

class FEShape
{

 public :

    FEShape() { }

    FEShape(const Element *el);

    FEShape(const Side *sd);

    virtual ~FEShape() { }

// Initialize local point coordinates in element.
// \a s is a point in the reference element
// This function computes jacobian, shape functions and their partial
// derivatives at \a s. Other member functions only return these values.
//    void setLocal(const Point<double> &s);

    double Sh(size_t i) const { return _sh[i-1]; }

    double Sh(size_t i, Point<double> s) const { s = 0; return _sh[i-1]; }

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

    double getDet() const { return _det; }

    Point<double> getCenter() const { return _c; }

    Point<double> getLocalPoint() const
    {
       Point<double> s = 0.;
       for (size_t i=0; i<_x.size(); i++)
          s += _sh[i] * _x[i];
       return s;
    }

    Point<double> getLocalPoint(const Point<double> &s) const
    {
       Point<double> t = 0.;
       for (size_t i=0; i<_x.size(); i++)
          t += Sh(i+1,s) * _x[i];
       return t;
    }

 protected:

#ifndef DOXYGEN_SHOULD_SKIP_THIS
   bool                     _localized;
   double                   _det;
   valarray<double>         _sh;
   valarray<size_t>         _node;
   valarray<Point<double> > _x, _dsh, _dshl;
   Point<double>            _c;
   size_t                   _label;
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani