v0.14.0
Searching...
No Matches
Riemann_plus_Riemann.hpp
Go to the documentation of this file.
1/* Adds a Riemann to a Riemann, yielding a
2 Riemann. */
3
4#pragma once
5
6namespace FTensor
7{
8 template <class A, class B, class T, class U, int Dim, char i, char j,
9 char k, char l>
11 {
14
15 public:
16 typename promote<T, U>::V
17 operator()(const int N1, const int N2, const int N3, const int N4) const
18 {
19 return iterA(N1, N2, N3, N4) + iterB(N1, N2, N3, N4);
20 }
21
24 : iterA(a), iterB(b)
25 {}
26 };
27
28 template <class A, class B, class T, class U, int Dim, char i, char j,
29 char k, char l>
30 Riemann_Expr<Riemann_plus_Riemann<A, B, T, U, Dim, i, j, k, l>,
31 typename promote<T, U>::V, Dim, i, j, k, l>
34 {
37 TensorExpr(a, b));
38 }
39}
static Number< 2 > N2
static Number< 1 > N1
constexpr double a
Riemann_Expr< B, U, Dim, i, j, k, l > iterB
Riemann_Expr< A, T, Dim, i, j, k, l > iterA
Riemann_plus_Riemann(const Riemann_Expr< A, T, Dim, i, j, k, l > &a, const Riemann_Expr< B, U, Dim, i, j, k, l > &b)
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
const double T
Tensors class implemented by Walter Landry.
Definition: FTensor.hpp:51
Ddg_Expr< Ddg_plus_Ddg< A, B, T, U, Dim01_0, Dim23_0, Dim01_1, Dim23_1, i0, j0, k0, l0, i1, j1, k1, l1 >, typename promote< T, U >::V, Dim01_0, Dim23_0, i0, j0, k0, l0 > operator+(const Ddg_Expr< A, T, Dim01_0, Dim23_0, i0, j0, k0, l0 > &a, const Ddg_Expr< B, U, Dim01_1, Dim23_1, i1, j1, k1, l1 > &b)
constexpr AssemblyType A