1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
// Gmsh - Copyright (C) 1997-2020 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
#ifndef STENSOR63_H_
#define STENSOR63_H_
class STensor63 {
protected:
double _val[3][3][3][3][3][3];
public:
STensor63(const double val = 0);
STensor63(const STensor63 &src);
inline double &operator()(int i, int j, int k, int l, int m, int n)
{
return _val[i][j][k][l][m][n];
};
inline double operator()(int i, int j, int k, int l, int m, int n) const
{
return _val[i][j][k][l][m][n];
};
STensor63 &operator=(const STensor63 &other)
{
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
for(int k = 0; k < 3; k++)
for(int l = 0; l < 3; l++)
for(int m = 0; m < 3; m++)
for(int n = 0; n < 3; n++)
_val[i][j][k][l][m][n] = other._val[i][j][k][l][m][n];
return *this;
};
STensor63 &operator+=(const STensor63 &other)
{
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
for(int k = 0; k < 3; k++)
for(int l = 0; l < 3; l++)
for(int m = 0; m < 3; m++)
for(int n = 0; n < 3; n++)
_val[i][j][k][l][m][n] += other._val[i][j][k][l][m][n];
return *this;
};
STensor63 &operator-=(const STensor63 &other)
{
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
for(int k = 0; k < 3; k++)
for(int l = 0; l < 3; l++)
for(int m = 0; m < 3; m++)
for(int n = 0; n < 3; n++)
_val[i][j][k][l][m][n] -= other._val[i][j][k][l][m][n];
return *this;
};
STensor63 &operator*=(const double &other)
{
for(int i = 0; i < 3; i++)
for(int j = 0; j < 3; j++)
for(int k = 0; k < 3; k++)
for(int l = 0; l < 3; l++)
for(int m = 0; m < 3; m++)
for(int n = 0; n < 3; n++) _val[i][j][k][l][m][n] *= other;
return *this;
};
};
#endif // STENSOR63_H_
|
