LocalIntegral.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _LOCAL_INTEGRAL_H
#define _LOCAL_INTEGRAL_H
 
#include "MatVec.h"
 
 
class LocalIntegral
{
protected:
  LocalIntegral() {}
 
public:
  virtual ~LocalIntegral() {}
  virtual void destruct() { delete this; }
  virtual const LocalIntegral* ref() const { return this; }
 
  void getSolution(size_t nsd, size_t nen, Matrix* u = nullptr,
                   Matrix* v = nullptr, Matrix* a = nullptr,
                   bool forceCurrent = false) const
  {
    int nvec = vec.size();
    if (nvec > 5 && forceCurrent) nvec = 3;
    int idis = nvec > 5 ? 3 : 0; // index to element displacement vector (u)
    if (nvec == 2 && !forceCurrent) idis = 1; // probably a quasi-static problem
    int ivel = idis + 1; // index to element velocity vector (v)
    int iacc = ivel + 1; // index to element acceleration vector (a)
    if (u && idis < nvec) u->fill(vec[idis],nsd,nen);
    if (v && ivel < nvec) v->fill(vec[ivel],nsd,nen);
    if (a && iacc < nvec) a->fill(vec[iacc],nsd,nen);
  }
 
  Vectors vec; 
};
 
#endif