app-cwinkelm/elasticityTensor.hh

Go to the documentation of this file.

// Elasticity Tensor

#ifndef elasticityTensor_hh
#define elasticityTensor_hh

#include "basics/vectorsMatrices.hh"

using concepts::Real;

template <int dim>
class ElasticityTensor {
 public:
  virtual Real entry(uint i, uint j, uint r, uint s) const = 0;
  virtual concepts::Mapping<Real, dim> operator() (uint i, uint j) const;
  virtual concepts::Mapping<Real, dim> operator()
    (concepts::Mapping<Real, dim>& strain) const;
 private:
  concepts::Mapping<Real, dim> rp2map_(Real* temp) const;
};

template <int dim>
class ElasticityTensorIsotropic : public ElasticityTensor<dim> {
 public:
  ElasticityTensorIsotropic(Real mu, Real lambda);
  virtual Real entry(uint i, uint j, uint r, uint s) const;
 private:
  Real mu_;
  Real lambda_;
};

#endif // elasticityTensor_hh