EshelbianAux.hpp

Go to the documentation of this file.

/**
 * @file EshelbianAux.hpp
 * @brief Auxilary functions for Eshelbian plasticity
 * @date 2024-08-30
 *
 * @copyright Copyright (c) 2024
 *
 */
 
#include <Lie.hpp>
 
namespace EshelbianPlasticity {
 
constexpr static auto size_symm = (3 * (3 + 1)) / 2;
 
inline auto voigt_to_symm() {
  FTensor::Index<'L', size_symm> L;
  FTensor::Index<'M', size_symm> M;
  FTensor::Tensor2<double, size_symm, size_symm> t_L;
  t_L(L, M) = 0;
  t_L(0, 0) = 1;
  t_L(1, 3) = 1;
  t_L(2, 5) = 1;
  t_L(3, 1) = 1;
  t_L(4, 2) = 1;
  t_L(5, 4) = 1;
  return t_L;
}
 
template <class T> struct TensorTypeExtractor {
  typedef typename std::remove_pointer<T>::type Type;
};
template <class T, int I> struct TensorTypeExtractor<FTensor::PackPtr<T, I>> {
  typedef typename std::remove_pointer<T>::type Type;
};
 
}