FreeFOAM The Cross-Platform CFD Toolkit
LUscalarMatrixTemplates.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) 1991-2010 OpenCFD Ltd.
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 
27 
28 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
29 
30 template<class Type>
32 {
33  if (Pstream::parRun())
34  {
35  Field<Type> completeSourceSol(n());
36 
37  if (Pstream::master())
38  {
39  typename Field<Type>::subField
40  (
41  completeSourceSol,
42  sourceSol.size()
43  ).assign(sourceSol);
44 
45  for
46  (
47  int slave=Pstream::firstSlave();
48  slave<=Pstream::lastSlave();
49  slave++
50  )
51  {
53  (
55  slave,
56  reinterpret_cast<char*>
57  (
58  &(completeSourceSol[procOffsets_[slave]])
59  ),
60  (procOffsets_[slave + 1] - procOffsets_[slave])*sizeof(Type)
61  );
62  }
63  }
64  else
65  {
67  (
70  reinterpret_cast<const char*>(sourceSol.begin()),
71  sourceSol.byteSize()
72  );
73  }
74 
75  if (Pstream::master())
76  {
77  LUBacksubstitute(*this, pivotIndices_, completeSourceSol);
78 
79  sourceSol = typename Field<Type>::subField
80  (
81  completeSourceSol,
82  sourceSol.size()
83  );
84 
85  for
86  (
87  int slave=Pstream::firstSlave();
88  slave<=Pstream::lastSlave();
89  slave++
90  )
91  {
93  (
95  slave,
96  reinterpret_cast<const char*>
97  (
98  &(completeSourceSol[procOffsets_[slave]])
99  ),
100  (procOffsets_[slave + 1] - procOffsets_[slave])*sizeof(Type)
101  );
102  }
103  }
104  else
105  {
107  (
110  reinterpret_cast<char*>(sourceSol.begin()),
111  sourceSol.byteSize()
112  );
113  }
114  }
115  else
116  {
117  LUBacksubstitute(*this, pivotIndices_, sourceSol);
118  }
119 }
120 
121 
122 // ************************ vim: set sw=4 sts=4 et: ************************ //