SpatialTensorI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2016-2018 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "Identity.H"
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class Cmpt>
32 {}
33 
34 
35 template<class Cmpt>
37 :
39 {}
40 
41 
42 template<class Cmpt>
44 (
45  const typename SpatialTensor::msType& ms
46 )
47 :
49 {}
50 
51 
52 template<class Cmpt>
54 (
55  const Tensor<Cmpt>& t00, const Tensor<Cmpt>& t01,
56  const Tensor<Cmpt>& t10, const Tensor<Cmpt>& t11
57 )
58 {
59  // Block (0, 0)
60  this->v_[0] = t00.xx(); this->v_[1] = t00.xy(); this->v_[2] = t00.xz();
61  this->v_[6] = t00.yx(); this->v_[7] = t00.yy(); this->v_[8] = t00.yz();
62  this->v_[12] = t00.zx(); this->v_[13] = t00.zy(); this->v_[14] = t00.zz();
63 
64  // Block (0, 1)
65  this->v_[3] = t01.xx(); this->v_[4] = t01.xy(); this->v_[5] = t01.xz();
66  this->v_[9] = t01.yx(); this->v_[10] = t01.yy(); this->v_[11] = t01.yz();
67  this->v_[15] = t01.zx(); this->v_[16] = t01.zy(); this->v_[17] = t01.zz();
68 
69  // Block (1, 0)
70  this->v_[18] = t10.xx(); this->v_[19] = t10.xy(); this->v_[20] = t10.xz();
71  this->v_[24] = t10.yx(); this->v_[25] = t10.yy(); this->v_[26] = t10.yz();
72  this->v_[30] = t10.zx(); this->v_[31] = t10.zy(); this->v_[32] = t10.zz();
73 
74  // Block (1, 1)
75  this->v_[21] = t11.xx(); this->v_[22] = t11.xy(); this->v_[23] = t11.xz();
76  this->v_[27] = t11.yx(); this->v_[28] = t11.yy(); this->v_[29] = t11.yz();
77  this->v_[33] = t11.zx(); this->v_[34] = t11.zy(); this->v_[35] = t11.zz();
78 }
79 
80 
81 template<class Cmpt>
83 (
84  const Cmpt& t00, const Cmpt& t01, const Cmpt& t02,
85  const Cmpt& t03, const Cmpt& t04, const Cmpt& t05,
86 
87  const Cmpt& t10, const Cmpt& t11, const Cmpt& t12,
88  const Cmpt& t13, const Cmpt& t14, const Cmpt& t15,
89 
90  const Cmpt& t20, const Cmpt& t21, const Cmpt& t22,
91  const Cmpt& t23, const Cmpt& t24, const Cmpt& t25,
92 
93  const Cmpt& t30, const Cmpt& t31, const Cmpt& t32,
94  const Cmpt& t33, const Cmpt& t34, const Cmpt& t35,
95 
96  const Cmpt& t40, const Cmpt& t41, const Cmpt& t42,
97  const Cmpt& t43, const Cmpt& t44, const Cmpt& t45,
98 
99  const Cmpt& t50, const Cmpt& t51, const Cmpt& t52,
100  const Cmpt& t53, const Cmpt& t54, const Cmpt& t55
101 )
102 {
103  // Row 0
104  this->v_[0] = t00; this->v_[1] = t01; this->v_[2] = t02;
105  this->v_[3] = t03; this->v_[4] = t04; this->v_[5] = t05;
106 
107  // Row 1
108  this->v_[6] = t10; this->v_[7] = t11; this->v_[8] = t12;
109  this->v_[9] = t13; this->v_[10] = t14; this->v_[11] = t15;
110 
111  // Row 2
112  this->v_[12] = t20; this->v_[13] = t21; this->v_[14] = t22;
113  this->v_[15] = t23; this->v_[16] = t24; this->v_[17] = t25;
114 
115  // Row 3
116  this->v_[18] = t30; this->v_[19] = t31; this->v_[20] = t32;
117  this->v_[21] = t33; this->v_[22] = t34; this->v_[23] = t35;
118 
119  // Row 4
120  this->v_[24] = t40; this->v_[25] = t41; this->v_[26] = t42;
121  this->v_[27] = t43; this->v_[28] = t44; this->v_[29] = t45;
122 
123  // Row 5
124  this->v_[30] = t50; this->v_[31] = t51; this->v_[32] = t52;
125  this->v_[33] = t53; this->v_[34] = t54; this->v_[35] = t55;
126 }
127 
128 
129 template<class Cmpt>
131 :
132  SpatialTensor::msType(is)
133 {}
134 
135 
136 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
137 
138 namespace Foam
139 {
140 
141 // * * * * * * * * * * * * * * * Global Operators * * * * * * * * * * * * * //
142 
143 //- Return the cross-product tensor
144 template<class Cmpt>
146 (
147  const SpatialVector<Cmpt>& v,
148  const Identity<Cmpt>&
149 )
150 {
151  return SpatialTensor<Cmpt>
152  (
153  0, -v.wz(), v.wy(), 0, 0, 0,
154  v.wz(), 0, -v.wx(), 0, 0, 0,
155  -v.wy(), v.wx(), 0, 0, 0, 0,
156  0, -v.lz(), v.ly(), 0, -v.wz(), v.wy(),
157  v.lz(), 0, -v.lx(), v.wz(), 0, -v.wx(),
158  -v.ly(), v.lx(), 0, -v.wy(), v.wx(), 0
159  );
160 }
161 
162 
163 //- Return the dual cross-product tensor
164 template<class Cmpt>
166 (
167  const SpatialVector<Cmpt>& f,
168  const typename Identity<Cmpt>::dual&
169 )
170 {
171  return SpatialTensor<Cmpt>
172  (
173  0, -f.wz(), f.wy(), 0, -f.lz(), f.ly(),
174  f.wz(), 0, -f.wx(), f.lz(), 0, -f.lx(),
175  -f.wy(), f.wx(), 0, -f.ly(), f.lx(), 0,
176  0, 0, 0, 0, -f.wz(), f.wy(),
177  0, 0, 0, f.wz(), 0, -f.wx(),
178  0, 0, 0, -f.wy(), f.wx(), 0
179  );
180 }
181 
182 
183 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
184 
185 } // End namespace Foam
186 
187 // ************************************************************************* //
The identity type in the dual space.
Definition: Identity.H:64
Templated identity and dual space identity tensors derived from SphericalTensor.
Definition: Identity.H:50
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
Definition: Istream.H:60
Templated matrix space.
Definition: MatrixSpace.H:58
MatrixSpace< SpatialTensor< Cmpt >, Cmpt, Mrows, Ncols > msType
MatrixSpace type.
Definition: MatrixSpace.H:63
Templated 3D spatial tensor derived from MatrixSpace used to represent transformations of spatial vec...
Definition: SpatialTensor.H:67
SpatialTensor()
Construct null.
Templated 3D spatial vector derived from VectorSpace used to represent the anglular and linear compon...
Definition: SpatialVector.H:65
Templated 3D tensor derived from MatrixSpace adding construction from 9 components,...
Definition: Tensor.H:67
const Cmpt & xx() const
Definition: TensorI.H:163
const Cmpt & yx() const
Definition: TensorI.H:184
const Cmpt & yz() const
Definition: TensorI.H:198
const Cmpt & xz() const
Definition: TensorI.H:177
const Cmpt & zz() const
Definition: TensorI.H:219
const Cmpt & xy() const
Definition: TensorI.H:170
const Cmpt & zx() const
Definition: TensorI.H:205
const Cmpt & zy() const
Definition: TensorI.H:212
const Cmpt & yy() const
Definition: TensorI.H:191
A class representing the concept of 0 used to avoid unnecessary manipulations for objects that are kn...
Definition: zero.H:50
Namespace for OpenFOAM.
static const zero Zero
Definition: zero.H:97
labelList f(nPoints)