ElmMats.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ELM_MATS_H
#define _ELM_MATS_H
 
#include "LocalIntegral.h"
#include "MatVec.h"
 
 
class ElmMats : public LocalIntegral
{
public:
  explicit ElmMats(bool lhs = true) : rhsOnly(false), withLHS(lhs) {}
  virtual ~ElmMats() {}
 
  void resize(size_t nA, size_t nB, size_t nC = 0);
 
  void redim(size_t ndim);
 
  virtual bool empty() const { return A.empty() && b.empty(); }
 
  virtual void setStepSize(double, int) {}
 
  virtual const Matrix& getNewtonMatrix() const;
  virtual const Vector& getRHSVector() const;
 
  void printMat(std::ostream& os, size_t idx = 0,
                const char* prefix = nullptr) const;
  void printVec(std::ostream& os, size_t idx = 0,
                const char* prefix = nullptr) const;
  void printScl(std::ostream& os, size_t idx = 0,
                const char* prefix = nullptr) const;
 
  std::vector<Matrix> A; 
  std::vector<Vector> b; 
  std::vector<double> c; 
 
  std::vector<const char*> Aname; 
  std::vector<const char*> Bname; 
  std::vector<const char*> Cname; 
 
  bool rhsOnly; 
  bool withLHS; 
};
 
#endif