#include <src/ftensor/src/FTensor/Tensor1/diffusion_Tensor1.hpp>

Collaboration diagram for FTensor::diffusion_Tensor1< T, Dim, i >:

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

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

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

const int	di

const int	dj

const int	dk

const double	dx

Detailed Description

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

Definition at line 8 of file diffusion_Tensor1.hpp.

Constructor & Destructor Documentation

◆ diffusion_Tensor1()

template<class T , int Dim, char i>

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

inline

Definition at line 32 of file diffusion_Tensor1.hpp.

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

35 {}

FTensor::diffusion_Tensor1::dx

const double dx

Definition: diffusion_Tensor1.hpp:12

FTensor::diffusion_Tensor1::di

const int di

Definition: diffusion_Tensor1.hpp:11

FTensor::diffusion_Tensor1::a

const Tensor1< T *, Dim > & a

Definition: diffusion_Tensor1.hpp:10

FTensor::diffusion_Tensor1::dk

const int dk

Definition: diffusion_Tensor1.hpp:11

FTensor::diffusion_Tensor1::dj

const int dj

Definition: diffusion_Tensor1.hpp:11

A

constexpr AssemblyType A

Definition: operators_tests.cpp:30

Member Function Documentation

◆ operator()()

template<class T , int Dim, char i>

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

inline

Definition at line 15 of file diffusion_Tensor1.hpp.

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