ElementBase.h

//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis Claudio Pérez Tato
//
//  This program derives from OpenSees <http://opensees.berkeley.edu>
//  developed by the  «Pacific earthquake engineering research center».
//
//  Except for the restrictions that may arise from the copyright
//  of the original program (see copyright_opensees.txt)
//  XC 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 
//  (at your option) any later version.
//
//  This software 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/>.
//----------------------------------------------------------------------------
//ElementBase.h

#ifndef ElementBase_h
#define ElementBase_h

#include "Element.h"
#include "domain/mesh/element/utils/NodePtrsWithIDs.h"
#include "material/Material.h"
#include "domain/domain/Domain.h"
#include "xc_utils/src/geom/pos_vec/Pos3d.h"

namespace XC {
class Node;

//
template <int NNODOS>
class ElementBase: public Element
  {
  protected:
    NodePtrsWithIDs theNodes;  

    template <class TIPOMAT>
    TIPOMAT *cast_material(const Material *ptr_mat);

    int sendData(CommParameters &cp);
    int recvData(const CommParameters &cp);

  public:
    ElementBase(int tag, int classTag);
    ElementBase(const ElementBase &otro);
    ElementBase &operator=(const ElementBase &otro);

    static inline const int numNodos(void)
      { return NNODOS; }
    // public methods to obtain inforrmation about dof & connectivity    
    int getNumExternalNodes(void) const;
    NodePtrsWithIDs &getNodePtrs(void);
    const NodePtrsWithIDs &getNodePtrs(void) const;
    Pos3d getPosCdg(bool initialGeometry= true) const;
  };



template <int NNODOS>
XC::ElementBase<NNODOS>::ElementBase(int tag, int classTag)
  : Element(tag,classTag), theNodes(this,NNODOS) {}

template <int NNODOS>
XC::ElementBase<NNODOS>::ElementBase(const ElementBase<NNODOS> &otro)
  : Element(otro), theNodes(otro.theNodes) 
  { theNodes.set_owner(this); }

template <int NNODOS>
XC::ElementBase<NNODOS> &XC::ElementBase<NNODOS>::operator=(const ElementBase &otro)
  {
    Element::operator=(otro);
    theNodes= otro.theNodes;
    theNodes.set_owner(this);
    return *this;
  }

template <int NNODOS>
int XC::ElementBase<NNODOS>::getNumExternalNodes(void) const
  { return theNodes.size(); }

template <int NNODOS>
const XC::NodePtrsWithIDs &XC::ElementBase<NNODOS>::getNodePtrs(void) const
  { return theNodes; }

template <int NNODOS>
XC::NodePtrsWithIDs &XC::ElementBase<NNODOS>::getNodePtrs(void)
  { return theNodes; }

template <int NNODOS> template <class TIPOMAT>
TIPOMAT *ElementBase<NNODOS>::cast_material(const Material *ptr_mat)
  {
    TIPOMAT *retval= nullptr;
    const TIPOMAT *tmp = dynamic_cast<const TIPOMAT *>(ptr_mat);
    if(tmp)
      retval= tmp->getCopy();
    else
      {
        std::cerr << nombre_clase() << "::" << __FUNCTION__
              << "; on element: " << getTag() 
                  << " the material " << ptr_mat->getTag()
                  << " has not a suitable type." << std::endl;
        abort();
      }
    if(!retval)
      {
    std::cerr << nombre_clase() << "::" << __FUNCTION__
              << "; on element: " << getTag()
          << "can't get a copy of the material with tag: "
          << ptr_mat->getTag() << std::endl;
        abort();
      }
    return retval;
  }

template <int NNODOS>
int XC::ElementBase<NNODOS>::sendData(CommParameters &cp)
  {
    int res= Element::sendData(cp);
    res+= cp.sendMovable(theNodes,getDbTagData(),CommMetaData(6));
    return res;
  }

template <int NNODOS>
int XC::ElementBase<NNODOS>::recvData(const CommParameters &cp)
  {
    int res= Element::recvData(cp);
    res+= cp.receiveMovable(theNodes,getDbTagData(),CommMetaData(6));
    return res;
  }

template <int NNODOS>
Pos3d XC::ElementBase<NNODOS>::getPosCdg(bool initialGeometry) const
  { return theNodes.getPosCdg(initialGeometry); }

} //end of XC namespace
#endif