#include <src/ftensor/src/FTensor/Tensor2_symmetric/diffusion_Tensor2_symmetric.hpp>

Collaboration diagram for FTensor::diffusion_Tensor2_symmetric< T, Dim, i, j >:

Public Member Functions
promote< T, double >::V	operator() (const int N1, const int N2) const

	diffusion_Tensor2_symmetric (const Tensor2_symmetric< T *, Dim > &A, const int Di, const int Dj, const int Dk, const double Dx)

Private Attributes
const Tensor2_symmetric< T *, Dim > &	a

const int	di

const int	dj

const int	dk

const double	dx

Detailed Description

template<class T, int Dim, char i, char j>
class FTensor::diffusion_Tensor2_symmetric< T, Dim, i, j >

Definition at line 8 of file diffusion_Tensor2_symmetric.hpp.

Constructor & Destructor Documentation

◆ diffusion_Tensor2_symmetric()

template<class T , int Dim, char i, char j>

FTensor::diffusion_Tensor2_symmetric< T, Dim, i, j >::diffusion_Tensor2_symmetric	(	const Tensor2_symmetric< T *, Dim > &	A,
		const int	Di,
		const int	Dj,
		const int	Dk,
		const double	Dx
	)

inline

Definition at line 32 of file diffusion_Tensor2_symmetric.hpp.

35 : a(A), di(Di), dj(Dj), dk(Dk), dx(Dx)

36 {}

FTensor::diffusion_Tensor2_symmetric::a

const Tensor2_symmetric< T *, Dim > & a

Definition: diffusion_Tensor2_symmetric.hpp:10

FTensor::diffusion_Tensor2_symmetric::di

const int di

Definition: diffusion_Tensor2_symmetric.hpp:11

FTensor::diffusion_Tensor2_symmetric::dx

const double dx

Definition: diffusion_Tensor2_symmetric.hpp:12

FTensor::diffusion_Tensor2_symmetric::dk

const int dk

Definition: diffusion_Tensor2_symmetric.hpp:11

FTensor::diffusion_Tensor2_symmetric::dj

const int dj

Definition: diffusion_Tensor2_symmetric.hpp:11

A

constexpr AssemblyType A

Definition: operators_tests.cpp:30

Member Function Documentation

◆ operator()()

template<class T , int Dim, char i, char j>

promote< T, double >::V FTensor::diffusion_Tensor2_symmetric< T, Dim, i, j >::operator()	(	const int	N1,
		const int	N2
	)		const

inline

Definition at line 15 of file diffusion_Tensor2_symmetric.hpp.

    {
      return ((*(a.ptr(N1, N2) + di) - 2 * a(N1, N2) + *(a.ptr(N1, N2) - di))
              + (*(a.ptr(N1, N2) + dj) - 2 * a(N1, N2) + *(a.ptr(N1, N2) - dj))
              + (*(a.ptr(N1, N2) + dk) - 2 * a(N1, N2) + *(a.ptr(N1, N2) - dk))
              + ((*(a.ptr(N1, N2) + di + dj) + *(a.ptr(N1, N2) + di - dj)
                  + *(a.ptr(N1, N2) - di + dj) + *(a.ptr(N1, N2) - di - dj)
                  - 4 * a(N1, N2))
                 + (*(a.ptr(N1, N2) + di + dk) + *(a.ptr(N1, N2) + di - dk)
                    + *(a.ptr(N1, N2) - di + dk) + *(a.ptr(N1, N2) - di - dk)
                    - 4 * a(N1, N2))
                 + (*(a.ptr(N1, N2) + dj + dk) + *(a.ptr(N1, N2) + dj - dk)
                    + *(a.ptr(N1, N2) - dj + dk) + *(a.ptr(N1, N2) - dj - dk)
                    - 4 * a(N1, N2)))
                  / (std::sqrt(2.0)))
             * dx * dx;
    }