1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2008 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
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 the
13 Free Software Foundation; either version 2 of the License, or (at your
14 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
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
27 \*---------------------------------------------------------------------------*/
29 #include "GatherBase.H"
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
41 Type GatherBase::flatten(const List<Type> lst)
47 sum += lst[lstI].size();
56 const Type& sub = lst[lstI];
60 result[index++] = sub[subI];
68 template <class DataType, class IndexType, class AddOp>
69 IndexType GatherBase::offset
71 const List<DataType>& values,
72 const List<IndexType>& indices,
76 if (values.size() != indices.size())
80 "GatherBase::offset(const List<DataType>&, "
81 "const List<IndexType>&, AddOp)"
82 ) << "Input data and indices lists not equal size." << endl
83 << "data size:" << values.size()
84 << " indices:" << indices.size()
93 sum += indices[lstI].size();
96 IndexType result(sum);
102 forAll(indices, lstI)
104 const IndexType& sub = indices[lstI];
108 result[index++] = aop(sub[subI], offset);
110 offset += values[lstI].size();
117 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
119 } // End namespace Foam
121 // ************************************************************************* //