#include <src/ftensor/src/FTensor/Tensor0/dd_boundary_Tensor0.hpp>

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

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

	dd_boundary_Tensor0 (const Tensor0< T * > &A, const Tensor1< int, Dim > &D_ijk, const Tensor1< double, Dim > &D_xyz, const Tensor2< bool, Dim, 2 > &Boundary)

Private Attributes
const Tensor0< T * > &	a

const Tensor1< int, Dim > &	d_ijk

const Tensor1< double, Dim > &	d_xyz

const Tensor2< bool, Dim, 2 > &	boundary

Detailed Description

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

Definition at line 9 of file dd_boundary_Tensor0.hpp.

Constructor & Destructor Documentation

◆ dd_boundary_Tensor0()

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

FTensor::dd_boundary_Tensor0< T, Dim, i, j >::dd_boundary_Tensor0	(	const Tensor0< T * > &	A,
		const Tensor1< int, Dim > &	D_ijk,
		const Tensor1< double, Dim > &	D_xyz,
		const Tensor2< bool, Dim, 2 > &	Boundary
	)

inline

Definition at line 71 of file dd_boundary_Tensor0.hpp.

74 : a(A), d_ijk(D_ijk), d_xyz(D_xyz), boundary(Boundary)

75 {}

FTensor::dd_boundary_Tensor0::d_ijk

const Tensor1< int, Dim > & d_ijk

Definition: dd_boundary_Tensor0.hpp:12

FTensor::dd_boundary_Tensor0::a

const Tensor0< T * > & a

Definition: dd_boundary_Tensor0.hpp:11

FTensor::dd_boundary_Tensor0::d_xyz

const Tensor1< double, Dim > & d_xyz

Definition: dd_boundary_Tensor0.hpp:13

FTensor::dd_boundary_Tensor0::boundary

const Tensor2< bool, Dim, 2 > & boundary

Definition: dd_boundary_Tensor0.hpp:14

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::dd_boundary_Tensor0< T, Dim, i, j >::operator()	(	const int	N1,
		const int	N2
	)		const

inline

Definition at line 17 of file dd_boundary_Tensor0.hpp.

    {
      return N1 == N2
               ? (boundary(N1, 0)
                    ? (*(&a + 2 * d_ijk(N1)) - 2 * *(&a + d_ijk(N1)) + a)
                    : (boundary(N1, 1)
                         ? (a - 2 * *(&a - d_ijk(N1)) + *(&a - 2 * d_ijk(N1)))
                         : (*(&a + d_ijk(N1)) - 2 * a + *(&a - d_ijk(N1)))))
                   * d_xyz(N1) * d_xyz(N1)
               : (boundary(N1, 0)
                    ? (boundary(N2, 0)
                         ? ((*(&a + d_ijk(N1) + d_ijk(N2)) - *(&a + d_ijk(N2))
                             - *(&a + d_ijk(N1)) + a)
                            * d_xyz(N1) * d_xyz(N2))
                         : (boundary(N2, 1) ? ((*(&a + d_ijk(N1)) - a
                                                - *(&a + d_ijk(N1) - d_ijk(N2))
                                                + *(&a - d_ijk(N2)))
                                               * d_xyz(N1) * d_xyz(N2))
                                            : (*(&a + d_ijk(N1) + d_ijk(N2))
                                               - *(&a + d_ijk(N2))
                                               - *(&a + d_ijk(N1) - d_ijk(N2))
                                               + *(&a - d_ijk(N2)))
                                                * d_xyz(N1) * d_xyz(N2) * 0.5))
                    : (boundary(N1, 1)
                         ? (boundary(N2, 0)
                              ? ((*(&a + d_ijk(N2))
                                  - *(&a - d_ijk(N1) + d_ijk(N2)) - a
                                  + *(&a - d_ijk(N1)))
                                 * d_xyz(N1) * d_xyz(N2))
                              : (boundary(N2, 1)
                                   ? ((a - *(&a - d_ijk(N1)) - *(&a - d_ijk(N2))
                                       + *(&a - d_ijk(N1) - d_ijk(N2)))
                                      * d_xyz(N1) * d_xyz(N2))
                                   : (*(&a + d_ijk(N2))
                                      - *(&a - d_ijk(N1) + d_ijk(N2))
                                      - *(&a - d_ijk(N2))
                                      + *(&a - d_ijk(N1) - d_ijk(N2)))
                                       * d_xyz(N1) * d_xyz(N2) * 0.5))
                         : (boundary(N2, 0)
                              ? ((*(&a + d_ijk(N1) + d_ijk(N2))
                                  - *(&a - d_ijk(N1) + d_ijk(N2))
                                  - *(&a + d_ijk(N1)) + *(&a - d_ijk(N1)))
                                 * d_xyz(N1) * d_xyz(N2) * 0.5)
                              : (boundary(N2, 1)
                                   ? ((*(&a + d_ijk(N1)) - *(&a - d_ijk(N1))
                                       - *(&a + d_ijk(N1) - d_ijk(N2))
                                       + *(&a - d_ijk(N1) - d_ijk(N2)))
                                      * d_xyz(N1) * d_xyz(N2) * 0.5)
                                   : ((*(&a + d_ijk(N1) + d_ijk(N2))
                                       - *(&a - d_ijk(N1) + d_ijk(N2))
                                       - *(&a + d_ijk(N1) - d_ijk(N2))
                                       + *(&a - d_ijk(N1) - d_ijk(N2)))
                                      * d_xyz(N1) * d_xyz(N2) * 0.25)))));
    }