material.hh

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MATERIAL_HH_
#define MATERIAL_HH_
 
#include <opm/common/utility/platform_dependent/disable_warnings.h>
 
#include <dune/common/fmatrix.hh>
#include <dune/common/dynvector.hh>
 
#include <opm/common/utility/platform_dependent/reenable_warnings.h>
 
namespace Opm {
namespace Elasticity {
 
 
class Material
{
protected:
  Material(int ID = 0, double density = 0.0)
    : id(ID), rho(density)
  {
  }
 
  virtual std::ostream& write(std::ostream& os) const
  { 
    return os;
  }
public:
  virtual ~Material() {}
 
  int num() const
  {
    return id;
  }
 
  virtual int numPar() const = 0;
  virtual double getPar(int /* ipar */ = 1) const
  {
    return double(0);
  }
 
  virtual bool getConstitutiveMatrix(Dune::FieldMatrix<double,6,6>& C,
                     bool invers = false) const = 0;
 
  virtual bool getConstitutiveMatrix(Dune::FieldMatrix<double,3,3>& C,
                     bool invers = false) const = 0;
 
  double getMassDensity() const
  {
    return rho;
  }
 
  friend std::ostream& operator<<(std::ostream& os, const Material& m)
  {
    return m.write(os);
  }
 
  static Material* create(int ID, const Dune::DynamicVector<double>& params);
 
  static Material* create(int ID, const std::string& file);
private:
  int  id;  
  double rho; 
};
 
}
}
 
#endif