v0.14.0
Loading...
Searching...
No Matches
Functions
test_T3as_14.cpp File Reference
#include "../../../src/FTensor.hpp"
#include "../test_for_zero.hpp"
#include <iostream>

Go to the source code of this file.

Functions

void test_T3as_14 (const Tensor1< double, 3 > &t1_2, const Tensor2< double, 3, 3 > &t2_2, const Tensor3_antisymmetric< double, 3, 3 > &t3as_2)
 

Function Documentation

◆ test_T3as_14()

void test_T3as_14 ( const Tensor1< double, 3 > & t1_2,
const Tensor2< double, 3, 3 > & t2_2,
const Tensor3_antisymmetric< double, 3, 3 > & t3as_2 )

Definition at line 7 of file test_T3as_14.cpp.

10{
11 Index<'i', 3> i;
12 Index<'j', 3> j;
13 Index<'k', 3> k;
14
18
19 test_for_zero(t3as_2(i, j, k) * (t2_2(j, k) * t1_2(i))
20 - (t3as_2(0, 1, 0) * t2_2(1, 0) * t1_2(0)
21 + t3as_2(1, 1, 0) * t2_2(1, 0) * t1_2(1)
22 + t3as_2(2, 1, 0) * t2_2(1, 0) * t1_2(2)
23 + t3as_2(0, 2, 0) * t2_2(2, 0) * t1_2(0)
24 + t3as_2(1, 2, 0) * t2_2(2, 0) * t1_2(1)
25 + t3as_2(2, 2, 0) * t2_2(2, 0) * t1_2(2)
26 + t3as_2(0, 0, 1) * t2_2(0, 1) * t1_2(0)
27 + t3as_2(1, 0, 1) * t2_2(0, 1) * t1_2(1)
28 + t3as_2(2, 0, 1) * t2_2(0, 1) * t1_2(2)
29 + t3as_2(0, 2, 1) * t2_2(2, 1) * t1_2(0)
30 + t3as_2(1, 2, 1) * t2_2(2, 1) * t1_2(1)
31 + t3as_2(2, 2, 1) * t2_2(2, 1) * t1_2(2)
32 + t3as_2(0, 0, 2) * t2_2(0, 2) * t1_2(0)
33 + t3as_2(1, 0, 2) * t2_2(0, 2) * t1_2(1)
34 + t3as_2(2, 0, 2) * t2_2(0, 2) * t1_2(2)
35 + t3as_2(0, 1, 2) * t2_2(1, 2) * t1_2(0)
36 + t3as_2(1, 1, 2) * t2_2(1, 2) * t1_2(1)
37 + t3as_2(2, 1, 2) * t2_2(1, 2) * t1_2(2)),
38 "T3as(i,j,k)*T3dg(j,k,i)");
39 test_for_zero((t2_2(j, k) * t1_2(i) * t3as_2(i, j, k))
40 - (t3as_2(0, 1, 0) * t2_2(1, 0) * t1_2(0)
41 + t3as_2(1, 1, 0) * t2_2(1, 0) * t1_2(1)
42 + t3as_2(2, 1, 0) * t2_2(1, 0) * t1_2(2)
43 + t3as_2(0, 2, 0) * t2_2(2, 0) * t1_2(0)
44 + t3as_2(1, 2, 0) * t2_2(2, 0) * t1_2(1)
45 + t3as_2(2, 2, 0) * t2_2(2, 0) * t1_2(2)
46 + t3as_2(0, 0, 1) * t2_2(0, 1) * t1_2(0)
47 + t3as_2(1, 0, 1) * t2_2(0, 1) * t1_2(1)
48 + t3as_2(2, 0, 1) * t2_2(0, 1) * t1_2(2)
49 + t3as_2(0, 2, 1) * t2_2(2, 1) * t1_2(0)
50 + t3as_2(1, 2, 1) * t2_2(2, 1) * t1_2(1)
51 + t3as_2(2, 2, 1) * t2_2(2, 1) * t1_2(2)
52 + t3as_2(0, 0, 2) * t2_2(0, 2) * t1_2(0)
53 + t3as_2(1, 0, 2) * t2_2(0, 2) * t1_2(1)
54 + t3as_2(2, 0, 2) * t2_2(0, 2) * t1_2(2)
55 + t3as_2(0, 1, 2) * t2_2(1, 2) * t1_2(0)
56 + t3as_2(1, 1, 2) * t2_2(1, 2) * t1_2(1)
57 + t3as_2(2, 1, 2) * t2_2(1, 2) * t1_2(2)),
58 "T3dg(j,k,i)*T3as(i,j,k)");
59}
static Number< 2 > N2
static Number< 1 > N1
static Number< 0 > N0
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
void test_for_zero(const T &t, const std::string &s)