BlockElmMats.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _BLOCK_ELM_MATS_H
#define _BLOCK_ELM_MATS_H
 
#include "ElmMats.h"
 
 
class BlockElmMats : public ElmMats
{
public:
  BlockElmMats(size_t nBlk, size_t nBas) :
    neldof(0), blockInfo(nBlk), basisInfo(nBas) {}
  virtual ~BlockElmMats() {}
 
  using ElmMats::redim;
  bool redim(size_t blkIndex, size_t nen, size_t ncmp = 1, char basis = 1);
  bool redimOffDiag(size_t blkIndex, char symmetric = 1);
  bool finalize();
 
  virtual const Matrix& getNewtonMatrix() const;
  virtual const Vector& getRHSVector() const;
 
private:
  struct Block
  {
    char   basis; 
    size_t ncmp;  
    size_t idof;  
 
    Block(size_t c = 0, char b = 0) : basis(b), ncmp(c), idof(0) {}
  };
 
  struct Basis
  {
    size_t nen;  
    size_t ncmp; 
 
    Basis() : nen(0), ncmp(0) {}
  };
 
  size_t neldof; 
  std::vector<Block> blockInfo; 
  std::vector<Basis> basisInfo; 
  std::vector<char>  symmFlag;  
};
 
#endif