v0.15.0
Loading...
Searching...
No Matches
Ddg_mod_Ddg.hpp
Go to the documentation of this file.
1/* Subtracts Ddg-Ddg -> Riemann */
2
3#pragma once
4
5namespace FTensor
6{
7 /* A(i,j,k,l) - B(i,l,k,j) */
8
9 template <class A, class B, class T, class U, int Dim, char i, char j,
10 char k, char l>
12 {
15
16 public:
17 typename promote<T, U>::V
18 operator()(const int N1, const int N2, const int N3, const int N4) const
19 {
20 return iterA(N1, N2, N3, N4) - iterB(N1, N4, N3, N2);
21 }
22
27 };
28
29 template <class A, class B, class T, class U, int Dim, char i, char j,
30 char k, char l>
31 Riemann_Expr<Ddg_mod_Ddg<A, B, T, U, Dim, i, j, k, l>,
32 typename promote<T, U>::V, Dim, i, j, k, l>
40}
constexpr double a
promote< T, U >::V operator()(const int N1, const int N2, const int N3, const int N4) const
Ddg_Expr< B, U, Dim, Dim, i, l, k, j > iterB
Ddg_mod_Ddg(const Ddg_Expr< A, T, Dim, Dim, i, j, k, l > &a, const Ddg_Expr< B, U, Dim, Dim, i, l, k, j > &b)
Ddg_Expr< A, T, Dim, Dim, i, j, k, l > iterA
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'l', 3 > l
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
Tensors class implemented by Walter Landry.
Definition FTensor.hpp:51
Riemann_Expr< Ddg_mod_Ddg< A, B, T, U, Dim, i, j, k, l >, typename promote< T, U >::V, Dim, i, j, k, l > operator%(const Ddg_Expr< A, T, Dim, Dim, i, j, k, l > &a, const Ddg_Expr< B, U, Dim, Dim, i, l, k, j > &b)
constexpr AssemblyType A