ENH: decomposePar.C: add shortcircuit to avoid allocating point mappers
[OpenFOAM-2.0.x.git] / applications / utilities / parallelProcessing / decomposePar / dimFieldDecomposerDecomposeFields.C
blobbbe06248895ad71518edb826c796b41d1ac4d47e
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.
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.
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.
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/>.
24 \*---------------------------------------------------------------------------*/
26 #include "dimFieldDecomposer.H"
28 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
30 template<class Type>
31 Foam::tmp<Foam::DimensionedField<Type, Foam::volMesh> >
32 Foam::dimFieldDecomposer::decomposeField
34     const DimensionedField<Type, volMesh>& field
35 ) const
37     // Create and map the internal field values
38     Field<Type> mappedField(field, cellAddressing_);
40     // Create the field for the processor
41     return tmp<DimensionedField<Type, volMesh> >
42     (
43         new DimensionedField<Type, volMesh>
44         (
45             IOobject
46             (
47                 field.name(),
48                 procMesh_.time().timeName(),
49                 procMesh_,
50                 IOobject::NO_READ,
51                 IOobject::NO_WRITE
52             ),
53             procMesh_,
54             field.dimensions(),
55             mappedField
56         )
57     );
61 template<class GeoField>
62 void Foam::dimFieldDecomposer::decomposeFields
64     const PtrList<GeoField>& fields
65 ) const
67     forAll(fields, fieldI)
68     {
69         decomposeField(fields[fieldI])().write();
70     }
74 // ************************************************************************* //