TensorFunction.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _TENSOR_FUNCTION_H
#define _TENSOR_FUNCTION_H
 
#include "Function.h"
#include "matrixnd.h"
#include "Tensor.h"
 
 
class TensorFunc : public utl::SpatialFunction<Tensor>, public FunctionBase
{
protected:
  explicit TensorFunc(size_t n = 0) : utl::SpatialFunction<Tensor>(Tensor(n))
  {
    ncmp = zero.size();
  }
 
public:
  virtual ~TensorFunc() {}
 
  virtual unsigned char getType() const { return 3; }
 
  virtual std::vector<Real> getValue(const Vec3& X) const
  {
    return this->evaluate(X);
  }
 
  virtual Real getScalarValue(const Vec3& X) const
  {
    return this->evaluate(X).trace();
  }
 
  utl::matrix3d<Real> gradient(const Vec3& X) const;
 
  utl::matrix4d<Real> hessian(const Vec3& X) const;
 
  Tensor timeDerivative(const Vec3& X) const;
};
 
 
class STensorFunc : public utl::SpatialFunction<SymmTensor>,
                    public FunctionBase
{
  size_t index(size_t nsd, size_t i, size_t j) const;
 
protected:
  STensorFunc(size_t n = 0, bool with33 = false)
    : utl::SpatialFunction<SymmTensor>(SymmTensor(n,with33))
  {
    ncmp = zero.size();
  }
 
public:
  virtual ~STensorFunc() {}
 
  virtual unsigned char getType() const { return 3; }
 
  virtual std::vector<Real> getValue(const Vec3& X) const
  {
    return this->evaluate(X);
  }
 
  virtual Real getScalarValue(const Vec3& X) const
  {
    return this->evaluate(X).trace();
  }
 
  utl::matrix3d<Real> gradient(const Vec3& X) const;
 
  utl::matrix4d<Real> hessian(const Vec3& X) const;
 
  SymmTensor timeDerivative(const Vec3& X) const;
};
 
#endif