SymmetricSquareMatrix.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) 2011-2025 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 Class
25  Foam::SymmetricSquareMatrix
26 
27 Description
28  A templated 2D square symmetric matrix of objects of <T>, where the
29  n x n matrix dimension is known and used for subscript bounds checking, etc.
30 
31 SourceFiles
32  SymmetricSquareMatrixI.H
33  SymmetricSquareMatrix.C
34 
35 \*---------------------------------------------------------------------------*/
36 
37 #ifndef SymmetricSquareMatrix_H
38 #define SymmetricSquareMatrix_H
39 
40 #include "SquareMatrix.H"
41 #include "Identity.H"
42 
43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
44 
45 namespace Foam
46 {
47 
48 /*---------------------------------------------------------------------------*\
49  Class SymmetricSquareMatrix Declaration
50 \*---------------------------------------------------------------------------*/
51 
52 template<class Type>
54 :
55  public Matrix<SymmetricSquareMatrix<Type>, Type>
56 {
57 
58 public:
59 
60  // Constructors
61 
62  //- Null constructor.
63  inline SymmetricSquareMatrix();
64 
65  //- Construct given number of rows/columns.
66  inline SymmetricSquareMatrix(const label n);
67 
68  //- Construct given number of rows/columns, initialising to zero
69  inline SymmetricSquareMatrix(const label n, const zero);
70 
71  //- Construct given number of rows/columns,
72  inline SymmetricSquareMatrix(const label n, const Identity<Type>);
73 
74  //- Construct with given number of rows/columns
75  // initialising all elements to the given value
76  inline SymmetricSquareMatrix(const label n, const Type&);
77 
78  //- Construct with given number of rows/columns
79  // and start and end iterators
80  template<class InputIterator>
82  (
83  const label n,
84  InputIterator first,
85  InputIterator last
86  );
87 
88  //- Construct with given number of rows/columns and initialiser list
90  (
91  const label n,
92  std::initializer_list<Type>
93  );
94 
95  //- Construct from an initialiser list
96  inline SymmetricSquareMatrix(std::initializer_list<Type>);
97 
98  //- Construct from initialiser list list
100  (
101  std::initializer_list<std::initializer_list<Type>>
102  );
103 
104  //- Construct from Istream.
106 
107  //- Clone
109 
110 
111  // Member Functions
112 
113  //- Return the instantiated type name
114  static word typeName()
115  {
116  return
117  word("SymmetricSquareMatrix<") + pTraits<Type>::typeName + '>';
118  }
119 
120 
121  // Member Operators
122 
124 };
125 
126 
127 // Global functions
128 
129 //- Return the LU decomposed SymmetricSquareMatrix inverse
130 template<class Type>
132 
133 //- Return the SymmetricSquareMatrix inverse
134 template<class Type>
136 
137 //- Return the LU decomposed SymmetricSquareMatrix det
138 template<class Type>
140 
141 //- Return the SymmetricSquareMatrix det
142 template<class Type>
143 Type det(const SymmetricSquareMatrix<Type>&);
144 
145 
146 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
147 
148 } // End namespace Foam
149 
150 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
151 
152 #include "SymmetricSquareMatrixI.H"
153 
154 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
155 
156 #ifdef NoRepository
157  #include "SymmetricSquareMatrix.C"
158 #endif
159 
160 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
161 
162 #endif
163 
164 // ************************************************************************* //
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
A templated (m x n) matrix of objects of <T>.
Definition: Matrix.H:83
label n() const
Return the number of columns.
Definition: MatrixI.H:64
friend Ostream & operator(Ostream &, const mType &)
Write Matrix to Ostream.
A templated 2D square symmetric matrix of objects of <T>, where the n x n matrix dimension is known a...
static word typeName()
Return the instantiated type name.
autoPtr< SymmetricSquareMatrix< Type > > clone() const
Clone.
SymmetricSquareMatrix()
Null constructor.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: autoPtr.H:51
Traits class for primitives.
Definition: pTraits.H:53
A class for handling words, derived from string.
Definition: word.H:62
A class representing the concept of 0 used to avoid unnecessary manipulations for objects that are kn...
Definition: zero.H:50
Namespace for OpenFOAM.
SymmetricSquareMatrix< Type > invDecomposed(const SymmetricSquareMatrix< Type > &)
Return the LU decomposed SymmetricSquareMatrix inverse.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
void det(LagrangianPatchField< scalar > &f, const LagrangianPatchField< tensor > &f1)
labelList first(const UList< labelPair > &p)
Definition: patchToPatch.C:39
void inv(LagrangianPatchField< tensor > &f, const LagrangianPatchField< tensor > &f1)
scalar detDecomposed(const SquareMatrix< Type > &, const label sign)
Return the LU decomposed SquareMatrix det.