material.hh

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MATERIAL_HH_
#define MATERIAL_HH_

#include <dune/common/fmatrix.hh>
#include <dune/common/dynvector.hh>

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