procLduMatrix.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011 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 "procLduMatrix.H"
27 #include "procLduInterface.H"
28 #include "lduMatrix.H"
29 
30 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
31 
32 Foam::procLduMatrix::procLduMatrix
33 (
34  const lduMatrix& ldum,
35  const FieldField<Field, scalar>& interfaceCoeffs,
36  const lduInterfaceFieldPtrsList& interfaces
37 )
38 :
39  upperAddr_(ldum.lduAddr().upperAddr()),
40  lowerAddr_(ldum.lduAddr().lowerAddr()),
41  diag_(ldum.diag()),
42  upper_(ldum.upper()),
43  lower_(ldum.lower())
44 {
45  label nInterfaces = 0;
46 
47  forAll(interfaces, i)
48  {
49  if (interfaces.set(i))
50  {
51  nInterfaces++;
52  }
53  }
54 
55  interfaces_.setSize(nInterfaces);
56 
57  nInterfaces = 0;
58 
59  forAll(interfaces, i)
60  {
61  if (interfaces.set(i))
62  {
63  interfaces_.set
64  (
65  nInterfaces++,
67  (
68  interfaces[i],
69  interfaceCoeffs[i]
70  )
71  );
72  }
73  }
74 
75 }
76 
77 
78 Foam::procLduMatrix::procLduMatrix(Istream& is)
79 :
80  upperAddr_(is),
81  lowerAddr_(is),
82  diag_(is),
83  upper_(is),
84  lower_(is),
85  interfaces_(is)
86 {}
87 
88 
89 // * * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * //
90 
92 {
93  os << cldum.upperAddr_
94  << cldum.lowerAddr_
95  << cldum.diag_
96  << cldum.upper_
97  << cldum.lower_
98  << cldum.interfaces_;
99 
100  return os;
101 }
102 
103 
104 // ************************************************************************* //
virtual const labelUList & upperAddr() const =0
Return upper addressing.
Generic field type.
Definition: FieldField.H:51
scalarField & lower()
Definition: lduMatrix.C:165
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
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
I/O for lduMatrix and interface values.
Definition: procLduMatrix.H:61
Ostream & operator<<(Ostream &, const edgeMesh &)
Definition: edgeMeshIO.C:133
const lduAddressing & lduAddr() const
Return the LDU addressing.
Definition: lduMatrix.H:550
IO interface for processorLduInterface.
lduMatrix is a general matrix class in which the coefficients are stored as three arrays...
Definition: lduMatrix.H:77
#define forAll(list, i)
Definition: UList.H:421
virtual const labelUList & lowerAddr() const =0
Return lower addressing.
scalarField & diag()
Definition: lduMatrix.C:183
bool set(const label) const
Is element set.
Definition: UPtrListI.H:80
An Ostream is an abstract base class for all output systems (streams, files, token lists...
Definition: Ostream.H:53
scalarField & upper()
Definition: lduMatrix.C:194