/home/dorival/mechsys/lib/fem/geomelem.h

Go to the documentation of this file.

/************************************************************************
 * MechSys - Open Library for Mechanical Systems                        *
 * Copyright (C) 2005 Dorival M. Pedroso, Raul Durand                   *
 * Copyright (C) 2009 Sergio Galindo                                    *
 *                                                                      *
 * 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 3 of the License, or    *
 * 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, see <http://www.gnu.org/licenses/>  *
 ************************************************************************/

#ifndef MECHSYS_FEM_GEOMELEM_H
#define MECHSYS_FEM_GEOMELEM_H

// MechSys
#include <mechsys/geomtype.h>
#include <mechsys/fem/quadrature.h>
#include <mechsys/fem/node.h>
#include <mechsys/util/string.h>
#include <mechsys/linalg/matvec.h>

namespace FEM
{

class GeomElem
{
public:
    // Constructor & Destructor
    GeomElem (int          NDim,                         
              size_t       NN,                           
              size_t       NFN,                          
              char const * Name="__unnamed_geomelem__"); 
    virtual ~GeomElem () {}

    // Methods to be overloaded
    virtual void   SetIPs     (int TotNIP)                            =0; 
    virtual size_t FNode      (size_t IdxFace, size_t IdxFNode) const =0; 
    virtual void   Shape      (double r, double s, double t)    const =0; 
    virtual void   Derivs     (double r, double s, double t)    const =0; 
    virtual void   FaceShape  (double r, double s)              const =0; 
    virtual void   FaceDerivs (double r, double s)              const =0; 
    virtual void   NatCoords  (Mat_t & C)                       const =0; 

    // Constants
    size_t             NN;   
    size_t             NFN;  
    size_t             NIP;  
    size_t             NFIP; 
    IntegPoint const * IPs;  
    IntegPoint const * FIPs; 
    String             Name; 

    // Mutable data
    mutable Vec_t N;     
    mutable Mat_t dNdR;  
    mutable Vec_t FN;    
    mutable Mat_t FdNdR; 
};




inline GeomElem::GeomElem (int NDim, size_t TheNN, size_t TheNFN, char const * TheName)
    : NN(TheNN), NFN(TheNFN), Name(TheName)
{
    N    .change_dim (NN);
    dNdR .change_dim (NDim,NN);
    FN   .change_dim (NFN);
    FdNdR.change_dim (NDim-1,NFN);
}




SDPair GEOM;

typedef GeomElem * (*GeomElemMakerPtr)(int NDim);

typedef std::map<String, GeomElemMakerPtr> GeomElemFactory_t;

GeomElemFactory_t GeomElemFactory;

GeomElem * AllocGeomElem(String const & Name, int NDim)
{
    GeomElemFactory_t::iterator it = GeomElemFactory.find(Name);
    if (it==GeomElemFactory.end()) throw new Fatal("AllocGeomElem: '%s' is not available", Name.CStr());

    GeomElem * ptr = (*it->second)(NDim);

    return ptr;
}

}; // namespace FEM

#ifdef USE_BOOST_PYTHON
double PyGEOM (BPy::str const & Key) { return FEM::GEOM(BPy::extract<char const *>(Key)()); }
#endif

#endif // MECHSYS_FEM_GEOMELEM