EigenMatrix::Fdd4MImpl< E, C > Struct Template Reference

#include "src/ftensor/src/MatrixFunctionTemplate.hpp"

Collaboration diagram for EigenMatrix::Fdd4MImpl< E, C >:

Public Types
using	Val = typename E::Val
using	Vec = typename E::Vec
using	Fun = typename E::Fun
using	NumberNb = typename E::NumberNb
using	NumberDim = typename E::NumberDim

Public Member Functions
	Fdd4MImpl (E &e)
template<int A, int a, int b, int I, int J, int K, int L>
auto	fd2M () const
template<int a, int b, int i, int j, int k, int l, int m, int n>
auto	term_fd2S (const Number< a > &, const Number< b > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int NB, int a, int b, int i, int j, int k, int l, int m, int n>
C	term_SM (const Number< a > &, const Number< b > &, const Number< NB > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int nb, int a, int i, int j, int k, int l, int m, int n>
C	eval_fdS2 (const Number< nb > &, const Number< a > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int a, int i, int j, int k, int l, int m, int n>
C	eval_fdS2 (const Number< 1 > &, const Number< a > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int nb, int a, int i, int j, int k, int l, int m, int n>
C	eval_SM (const Number< nb > &, const Number< a > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int a, int i, int j, int k, int l, int m, int n>
C	eval_SM (const Number< 1 > &, const Number< a > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const
template<int nb, int a, int i, int j, int k, int l, int m, int n>
C	eval (const Number< nb > &, const Number< a > &, const Number< i > &, const Number< j > &, const Number< k > &, const Number< l > &, const Number< m > &, const Number< n > &) const

Public Attributes
E &	e

Detailed Description

template<typename E, typename C>
struct EigenMatrix::Fdd4MImpl< E, C >

Definition at line 163 of file MatrixFunctionTemplate.hpp.

Member Typedef Documentation

Fun

template<typename E , typename C >

using EigenMatrix::Fdd4MImpl< E, C >::Fun = typename E::Fun

Definition at line 166 of file MatrixFunctionTemplate.hpp.

NumberDim

template<typename E , typename C >

using EigenMatrix::Fdd4MImpl< E, C >::NumberDim = typename E::NumberDim

Definition at line 172 of file MatrixFunctionTemplate.hpp.

NumberNb

template<typename E , typename C >

using EigenMatrix::Fdd4MImpl< E, C >::NumberNb = typename E::NumberNb

Definition at line 171 of file MatrixFunctionTemplate.hpp.

Val

template<typename E , typename C >

using EigenMatrix::Fdd4MImpl< E, C >::Val = typename E::Val

Definition at line 164 of file MatrixFunctionTemplate.hpp.

Vec

template<typename E , typename C >

using EigenMatrix::Fdd4MImpl< E, C >::Vec = typename E::Vec

Definition at line 165 of file MatrixFunctionTemplate.hpp.

Constructor & Destructor Documentation

Fdd4MImpl()

template<typename E , typename C >

EigenMatrix::Fdd4MImpl< E, C >::Fdd4MImpl ( E &  e )

inline

Definition at line 168 of file MatrixFunctionTemplate.hpp.

: e(e) {}

Member Function Documentation

eval()

template<typename E , typename C >

template<int nb, int a, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::eval ( const Number< nb > &  , const Number< a > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 337 of file MatrixFunctionTemplate.hpp.

                                                            {
    return
 
        (2 * e.fVal(a)) * eval_SM(NumberDim(), Number<a>(), Number<i>(),
                                  Number<j>(), Number<k>(), Number<l>(),
                                  Number<m>(), Number<n>())
 
        +
 
        eval_fdS2(NumberDim(), Number<a>(), Number<i>(), Number<j>(),
                  Number<k>(), Number<l>(), Number<m>(), Number<n>());
  }

eval_fdS2() [1/2]

template<typename E , typename C >

template<int a, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::eval_fdS2 ( const Number< 1 > &  , const Number< a > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 298 of file MatrixFunctionTemplate.hpp.

                                                                 {
    if constexpr (a != 0)
      return term_fd2S(Number<a>(), Number<0>(), Number<i>(), Number<j>(),
                       Number<k>(), Number<l>(), Number<m>(), Number<n>());
    else
      return 0;
  }

eval_fdS2() [2/2]

template<typename E , typename C >

template<int nb, int a, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::eval_fdS2 ( const Number< nb > &  , const Number< a > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 284 of file MatrixFunctionTemplate.hpp.

                                                                 {
    if constexpr (a != nb - 1)
      return term_fd2S(Number<a>(), Number<nb - 1>(), Number<i>(), Number<j>(),
                       Number<k>(), Number<l>(), Number<m>(), Number<n>()) +
             eval_fdS2(Number<nb - 1>(), Number<a>(), Number<i>(), Number<j>(),
                       Number<k>(), Number<l>(), Number<m>(), Number<n>());
    else
      return eval_fdS2(Number<nb - 1>(), Number<a>(), Number<i>(), Number<j>(),
                       Number<k>(), Number<l>(), Number<m>(), Number<n>());
  }

eval_SM() [1/2]

template<typename E , typename C >

template<int a, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::eval_SM ( const Number< 1 > &  , const Number< a > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 325 of file MatrixFunctionTemplate.hpp.

                                                               {
    if constexpr (a != 0)
      return term_SM(Number<a>(), Number<0>(), NumberNb(), Number<i>(),
                     Number<j>(), Number<k>(), Number<l>(), Number<m>(),
                     Number<n>());
    else
      return 0;
  }

eval_SM() [2/2]

template<typename E , typename C >

template<int nb, int a, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::eval_SM ( const Number< nb > &  , const Number< a > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 309 of file MatrixFunctionTemplate.hpp.

                                                               {
    if constexpr (a != nb - 1)
      return term_SM(Number<a>(), Number<nb - 1>(), NumberNb(), Number<i>(),
                     Number<j>(), Number<k>(), Number<l>(), Number<m>(),
                     Number<n>()) +
             eval_SM(Number<nb - 1>(), Number<a>(), Number<i>(), Number<j>(),
                     Number<k>(), Number<l>(), Number<m>(), Number<n>());
 
    else
      return eval_SM(Number<nb - 1>(), Number<a>(), Number<i>(), Number<j>(),
                     Number<k>(), Number<l>(), Number<m>(), Number<n>());
  }

fd2M()

template<typename E , typename C >

template<int A, int a, int b, int I, int J, int K, int L>

auto EigenMatrix::Fdd4MImpl< E, C >::fd2M ( ) const

inline

Definition at line 175 of file MatrixFunctionTemplate.hpp.

                           {
    return e.d2MType1[b][get_sym_index(Number<A>(), Number<a>(), NumberDim())](
        Number<I>(), Number<J>(), Number<K>(), Number<L>());
  }

term_fd2S()

template<typename E , typename C >

template<int a, int b, int i, int j, int k, int l, int m, int n>

auto EigenMatrix::Fdd4MImpl< E, C >::term_fd2S ( const Number< a > &  , const Number< b > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 181 of file MatrixFunctionTemplate.hpp.

                                                                    {
    if constexpr (i == j && k == l) {
      return 4 *
             (
 
                 fd2M<a, a, b, i, k, m, n>() * e.aM[b](Number<j>(), Number<l>())
 
                 +
 
                 fd2M<b, a, b, i, k, m, n>() * e.aM[a](Number<j>(), Number<l>())
 
             );
 
    } else if constexpr (i == j)
      return 2 *
             (
 
                 fd2M<a, a, b, i, k, m, n>() *
                     e.aM[b](Number<j>(), Number<l>()) +
                 fd2M<b, a, b, j, l, m, n>() * e.aM[a](Number<i>(), Number<k>())
 
                 +
 
                 fd2M<b, a, b, i, k, m, n>() *
                     e.aM[a](Number<j>(), Number<l>()) +
                 fd2M<a, a, b, j, l, m, n>() * e.aM[b](Number<i>(), Number<k>())
 
             );
    else if constexpr (k == l)
      return 2 *
             (
 
                 fd2M<a, a, b, i, k, m, n>() *
                     e.aM[b](Number<j>(), Number<l>()) +
                 fd2M<b, a, b, j, l, m, n>() *
                     e.aM[a](Number<i>(), Number<k>()) +
 
                 +
 
                     fd2M<b, a, b, i, k, m, n>() *
                     e.aM[a](Number<j>(), Number<l>()) +
                 fd2M<a, a, b, j, l, m, n>() * e.aM[b](Number<i>(), Number<k>())
 
             );
    else
      return fd2M<a, a, b, i, k, m, n>() * e.aM[b](Number<j>(), Number<l>()) +
             fd2M<b, a, b, j, l, m, n>() * e.aM[a](Number<i>(), Number<k>()) +
             fd2M<a, a, b, i, l, m, n>() * e.aM[b](Number<j>(), Number<k>()) +
             fd2M<b, a, b, j, k, m, n>() * e.aM[a](Number<i>(), Number<l>())
 
             +
 
             fd2M<b, a, b, i, k, m, n>() * e.aM[a](Number<j>(), Number<l>()) +
             fd2M<a, a, b, j, l, m, n>() * e.aM[b](Number<i>(), Number<k>()) +
             fd2M<b, a, b, i, l, m, n>() * e.aM[a](Number<j>(), Number<k>()) +
             fd2M<a, a, b, j, k, m, n>() * e.aM[b](Number<i>(), Number<l>());
  }

term_SM()

template<typename E , typename C >

template<int NB, int a, int b, int i, int j, int k, int l, int m, int n>

C EigenMatrix::Fdd4MImpl< E, C >::term_SM ( const Number< a > &  , const Number< b > &  , const Number< NB > &  , const Number< i > &  , const Number< j > &  , const Number< k > &  , const Number< l > &  , const Number< m > &  , const Number< n > &    ) const

inline

Definition at line 242 of file MatrixFunctionTemplate.hpp.

                                            {
 
    if constexpr (NB == 1) {
      return 0;
 
    } else if constexpr (NB == 2) {
 
      if constexpr (a == 1 || b == 1) {
        return
 
            e.aF2(Number<a>(), Number<b>()) *
            (e.aSM[get_nodiag_index(Number<b>(), Number<a>(), NumberDim())]
                  [get_sym_index(Number<m>(), Number<n>(), NumberDim())](
                      Number<i>(), Number<j>(), Number<k>(), Number<l>()) -
             e.aSM[get_nodiag_index(Number<a>(), Number<b>(), NumberDim())]
                  [get_sym_index(Number<m>(), Number<n>(), NumberDim())](
                      Number<i>(), Number<j>(), Number<k>(), Number<l>()));
 
      } else {
        return 0;
      }
 
    } else {
 
      return
 
          e.aF2(Number<a>(), Number<b>()) *
          (e.aSM[get_nodiag_index(Number<b>(), Number<a>(), NumberDim())]
                [get_sym_index(Number<m>(), Number<n>(), NumberDim())](
                    Number<i>(), Number<j>(), Number<k>(), Number<l>()) -
           e.aSM[get_nodiag_index(Number<a>(), Number<b>(), NumberDim())]
                [get_sym_index(Number<m>(), Number<n>(), NumberDim())](
                    Number<i>(), Number<j>(), Number<k>(), Number<l>()));
    }
 
    return 0;
  }

Member Data Documentation

e

template<typename E , typename C >

E& EigenMatrix::Fdd4MImpl< E, C >::e

Definition at line 169 of file MatrixFunctionTemplate.hpp.

---

The documentation for this struct was generated from the following file:

• src/ftensor/src/MatrixFunctionTemplate.hpp