ASMsupel.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ASM_SUPEL_H
#define _ASM_SUPEL_H
 
#include "ASMbase.h"
#include "ASMutils.h"
#include "ElmMats.h"
#include "Vec3.h"
 
 
class ASMsupel : public ASMbase
{
public:
  explicit ASMsupel(unsigned char n_f = 6) : ASMbase(0,3,n_f) {}
  ASMsupel(const ASMsupel& patch, unsigned char n_f) : ASMbase(patch,n_f) {}
  ASMsupel(const ASMsupel& patch) : ASMbase(patch) {}
  virtual ~ASMsupel() {}
 
  virtual bool read(std::istream& is);
  virtual bool write(std::ostream&, int) const;
  virtual bool generateFEMTopology();
  virtual bool empty() const { return myElmMat.empty(); }
 
  virtual int getNodeSetIdx(const std::string& setName) const;
  virtual const IntVec& getNodeSet(int iset) const;
  virtual int parseNodeSet(const std::string& setName, const char* cset);
 
  virtual Vec3 getCoord(size_t inod) const;
  virtual void getNodalCoordinates(Matrix& X, bool = false) const;
  virtual bool getElementCoordinates(Matrix& X, int iel, bool = false) const;
 
  virtual void getBoundaryNodes(int lIndex, IntVec& nodes,
                                int, int, int, bool local) const;
 
  virtual void findBoundaryElms(IntVec&, int, int) const {}
  virtual bool getParameterDomain(Real2DMat&, IntVec*) const { return false; }
  virtual void getElmConnectivities(IntMat&, int) const {}
  virtual bool updateCoords(const Vector&) { return false; }
  virtual Vec3 getElementCenter(int) const { return Vec3(); }
 
  virtual bool integrate(Integrand& integrand, GlobalIntegral& glbInt,
                         const TimeDomain&);
 
  virtual bool integrate(Integrand&, int, GlobalIntegral&,
                         const TimeDomain&) { return false; }
 
  virtual int evalPoint(const double*, double*, Vec3&) const { return 0; }
  virtual int findElementContaining(const double*) const { return 0; }
  virtual double findPoint(Vec3&, double*) const { return -1.0; }
 
  virtual bool tesselate(ElementBlock& grid, const int*) const;
 
  using ASMbase::evalSolution;
  virtual bool evalSolution(Matrix& sField, const Vector& locSol,
                            const int*, int, bool) const;
 
  virtual bool transform(const Matrix& Tlg);
 
private:
  Vec3Vec myNodes;  
  ElmMats myElmMat; 
 
  std::vector<ASM::NodeSet> nodeSets; 
};
 
#endif