HenckyOps Namespace Reference

Classes
struct	CommonData
struct	OpCalculateEigenVals
struct	OpCalculateHenckyPlasticStress
struct	OpCalculateHenckyStress
struct	OpCalculateLogC
struct	OpCalculateLogC_dC
struct	OpCalculatePiolaStress
struct	OpHenckyTangent
struct	OpPostProcHencky

Functions
auto	is_eq (const double &a, const double &b)
template<int DIM>
auto	get_uniq_nb (double *ptr)
template<int DIM>
auto	sort_eigen_vals (FTensor::Tensor1< double, DIM > &eig, FTensor::Tensor2< double, DIM, DIM > &eigen_vec)

Variables
constexpr double	eps = std::numeric_limits<float>::epsilon()
auto	f = [](double v) { return 0.5 * std::log(static_cast<long double>(v)); }
auto	d_f = [](double v) { return 0.5 / static_cast<long double>(v); }
auto	dd_f

Function Documentation

get_uniq_nb()

template<int DIM>

auto HenckyOps::get_uniq_nb ( double *  ptr )

inline

Definition at line 15 of file HenckyOps.hpp.

                                                        {
  std::array<double, DIM> tmp;
  std::copy(ptr, &ptr[DIM], tmp.begin());
  std::sort(tmp.begin(), tmp.end());
  return std::distance(tmp.begin(), std::unique(tmp.begin(), tmp.end(), is_eq));
};

is_eq()

auto HenckyOps::is_eq ( const double &  a, const double &  b  )

inline

Definition at line 11 of file HenckyOps.hpp.

                                                    {
  return std::abs(a - b) < eps;
};

sort_eigen_vals()

template<int DIM>

auto HenckyOps::sort_eigen_vals ( FTensor::Tensor1< double, DIM > &  eig, FTensor::Tensor2< double, DIM, DIM > &  eigen_vec  )

inline

Definition at line 23 of file HenckyOps.hpp.

                                                                         {
   std::array<std::pair<double, size_t>, DIM> tmp;
  auto is_eq_pair = [](const auto &a, const auto &b) {
    return a.first < b.first;
  };
 
  for (size_t i = 0; i != DIM; ++i)
    tmp[i] = std::make_pair(eig(i), i);
  std::sort(tmp.begin(), tmp.end(), is_eq_pair);
 
  int i, j, k;
  if (is_eq_pair(tmp[0], tmp[1])) {
    i = tmp[0].second;
    j = tmp[2].second;
    k = tmp[1].second;
  } else {
    i = tmp[0].second;
    j = tmp[1].second;
    k = tmp[2].second;
  }
 
  FTensor::Tensor2<double, 3, 3> eigen_vec_c{
      eigen_vec(i, 0), eigen_vec(i, 1), eigen_vec(i, 2),
 
      eigen_vec(j, 0), eigen_vec(j, 1), eigen_vec(j, 2),
 
      eigen_vec(k, 0), eigen_vec(k, 1), eigen_vec(k, 2)};
 
  FTensor::Tensor1<double, 3> eig_c{eig(i), eig(j), eig(k)};
 
  {
    FTensor::Index<'i', DIM> i;
    FTensor::Index<'j', DIM> j;
    eig(i) = eig_c(i);
    eigen_vec(i, j) = eigen_vec_c(i, j);
  }
};

Variable Documentation

d_f

auto HenckyOps::d_f = [](double v) { return 0.5 / static_cast<long double>(v); }

Definition at line 6 of file HenckyOps.hpp.

dd_f

auto HenckyOps::dd_f

Initial value:

= [](double v) {
  return -0.5 / (static_cast<long double>(v) * static_cast<long double>(v));
}

Definition at line 7 of file HenckyOps.hpp.

eps

constexpr double HenckyOps::eps = std::numeric_limits<float>::epsilon()

constexpr

Definition at line 3 of file HenckyOps.hpp.

f

auto HenckyOps::f = [](double v) { return 0.5 * std::log(static_cast<long double>(v)); }

Definition at line 5 of file HenckyOps.hpp.