ASMu2DLag.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ASM_U2D_LAG_H
#define _ASM_U2D_LAG_H
 
#include "ASMs2DLag.h"
#include "ASMutils.h"
 
 
class ASMu2DLag : public ASMs2DLag
{
  class BasisFunctionCache : public ASMs2DLag::BasisFunctionCache
  {
  public:
    BasisFunctionCache(const ASMu2DLag& pch)
      : ASMs2DLag::BasisFunctionCache(pch) {}
 
    virtual ~BasisFunctionCache() = default;
 
    double getParam(int, size_t, size_t, bool) const override { return 0.0; }
 
  protected:
    void setupParameters() override {}
  };
 
public:
  ASMu2DLag(unsigned char n = 2, unsigned char n_f = 2, char fType = 'm');
  ASMu2DLag(const ASMu2DLag& pch, unsigned char n_f);
  ASMu2DLag(const ASMu2DLag& pch);
  virtual ~ASMu2DLag() {}
 
  // Methods for model generation
  // ============================
 
  virtual bool read(std::istream& is);
  virtual bool generateFEMTopology();
  virtual bool empty() const { return nel == 0; }
 
  virtual int getNodeSetIdx(const std::string& setName) const;
  virtual const IntVec& getNodeSet(int iset) const;
  virtual bool isInNodeSet(int iset, int inod) const;
  virtual int parseNodeSet(const std::string& setName, const char* cset);
  virtual int parseNodeBox(const std::string& setName, const char* bbox);
  int addToNodeSet(const std::string& setName, int inod, bool extId = false);
  bool getNodeSet(int iset, std::string& name) const;
 
  virtual int getElementSetIdx(const std::string& setName) const;
  virtual const IntVec& getElementSet(int iset) const;
  virtual bool getElementSet(int iset, std::string& name) const;
  virtual bool isInElementSet(int iset, int iel) const;
  virtual int parseElemSet(const std::string& setName, const char* cset);
  virtual int parseElemBox(const std::string& setName,
                           const std::string& unionSet, const char* bbox);
  int addToElemSet(const std::string& setName, int iel, bool extId = false);
 
  virtual void getBoundaryNodes(int lIndex, IntVec& nodes,
                                int, int, int, bool local) const;
 
  using ASMs2DLag::integrate;
  virtual bool integrate(Integrand& integrand,
                         GlobalIntegral& glbInt, const TimeDomain& time);
 
  using ASMs2DLag::generateThreadGroups;
  virtual void generateThreadGroups(const Integrand&, bool silence,
                                    bool separateGroup1noded);
 
  // Post-processing methods
  // =======================
 
  virtual bool tesselate(ElementBlock& grid, const int*) const;
 
  bool writeXML(const char* fname) const;
 
protected:
  virtual bool evalSolPt(int iel, double xi, double eta,
                         size_t nCmp, const Vector& pchSol,
                         RealArray& ptSol, RealArray& N) const;
 
  void generateThreadGroupsMultiColored(bool silence, bool separateGroup1Noded);
 
  bool swapNode34; 
 
  std::vector<ASM::NodeSet> nodeSets; 
  std::vector<ASM::NodeSet> elemSets; 
 
private:
  char fileType; 
};
 
#endif