FreeFOAM The Cross-Platform CFD Toolkit
writeAddSubtractField.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) 2008-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 
26 template<class Type>
28 (
29  const IOobject& baseHeader,
30  const IOobject& addHeader,
31  const fvMesh& mesh,
32  bool& processed
33 )
34 {
36 
37  if
38  (
39  baseHeader.headerClassName() == fieldType::typeName
40  && baseHeader.headerClassName() == addHeader.headerClassName()
41  )
42  {
43  if (resultName_ == "")
44  {
45  if (calcMode_ == ADD)
46  {
47  resultName_ = baseHeader.name() + "_add_" + addHeader.name();
48  }
49  else
50  {
51  resultName_ = baseHeader.name() + "_subtract_"
52  + addHeader.name();
53  }
54  }
55 
56  Info<< " Reading " << baseHeader.name() << endl;
57  fieldType baseField(baseHeader, mesh);
58 
59  Info<< " Reading " << addHeader.name() << endl;
60  fieldType addField(addHeader, mesh);
61 
62  if (baseField.dimensions() == addField.dimensions())
63  {
64  Info<< " Calculating " << resultName_ << endl;
65 
66  fieldType newField
67  (
68  IOobject
69  (
70  resultName_,
71  mesh.time().timeName(),
72  mesh,
73  IOobject::NO_READ
74  ),
75  calcMode_ == ADD ? baseField + addField : baseField - addField
76  );
77  newField.write();
78  }
79  else
80  {
81  Info<< " Cannot calculate " << resultName_ << nl
82  << " - inconsistent dimensions: "
83  << baseField.dimensions() << " - " << addField.dimensions()
84  << endl;
85  }
86 
87  processed = true;
88  }
89 }
90 
91 
92 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //