v0.14.0
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 
6 namespace 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 }
FTensor
JSON compatible output.
Definition: Christof_constructor.hpp:6
FTensor::Riemann_plus_Riemann::Riemann_plus_Riemann
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)
Definition: Riemann_plus_Riemann.hpp:22
A
constexpr AssemblyType A
Definition: operators_tests.cpp:30
a
constexpr double a
Definition: approx_sphere.cpp:30
FTensor::Riemann_Expr
Definition: Riemann_Expr.hpp:15
FTensor::promote::V
T1 V
Definition: promote.hpp:17
FTensor::Riemann_plus_Riemann::iterB
Riemann_Expr< B, U, Dim, i, j, k, l > iterB
Definition: Riemann_plus_Riemann.hpp:13
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
FTensor::Riemann_plus_Riemann
Definition: Riemann_plus_Riemann.hpp:10
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
FTensor::operator+
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)
Definition: Ddg_plus_Ddg.hpp:66
FTensor::Riemann_plus_Riemann::iterA
Riemann_Expr< A, T, Dim, i, j, k, l > iterA
Definition: Riemann_plus_Riemann.hpp:12
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
FTensor::Riemann_plus_Riemann::operator()
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Definition: Riemann_plus_Riemann.hpp:17
EshelbianPlasticity::U
@ U
Definition: EshelbianContact.cpp:197
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21