FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

GAMGAgglomerationTemplates.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 \*---------------------------------------------------------------------------*/
00025 
00026 #include <OpenFOAM/GAMGAgglomeration.H>
00027 
00028 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00029 
00030 template<class Type>
00031 void Foam::GAMGAgglomeration::restrictField
00032 (
00033     Field<Type>& cf,
00034     const Field<Type>& ff,
00035     const label fineLevelIndex
00036 ) const
00037 {
00038     const labelList& fineToCoarse = restrictAddressing_[fineLevelIndex];
00039 
00040     if (ff.size() != fineToCoarse.size())
00041     {
00042         FatalErrorIn
00043         (
00044             "void GAMGAgglomeration::restrictField"
00045             "(Field<Type>& cf, const Field<Type>& ff, "
00046             "const label fineLevelIndex) const"
00047         )   << "field does not correspond to level " << fineLevelIndex
00048             << " sizes: field = " << ff.size()
00049             << " level = " << fineToCoarse.size()
00050             << abort(FatalError);
00051     }
00052 
00053     cf = pTraits<Type>::zero;
00054 
00055     forAll(ff, i)
00056     {
00057         cf[fineToCoarse[i]] += ff[i];
00058     }
00059 }
00060 
00061 
00062 template<class Type>
00063 void Foam::GAMGAgglomeration::restrictFaceField
00064 (
00065     Field<Type>& cf,
00066     const Field<Type>& ff,
00067     const label fineLevelIndex
00068 ) const
00069 {
00070     const labelList& fineToCoarse = faceRestrictAddressing_[fineLevelIndex];
00071 
00072     cf = pTraits<Type>::zero;
00073 
00074     forAll(fineToCoarse, ffacei)
00075     {
00076         label cFace = fineToCoarse[ffacei];
00077 
00078         if (cFace >= 0)
00079         {
00080             cf[cFace] += ff[ffacei];
00081         }
00082     }
00083 }
00084 
00085 
00086 template<class Type>
00087 void Foam::GAMGAgglomeration::prolongField
00088 (
00089     Field<Type>& ff,
00090     const Field<Type>& cf,
00091     const label coarseLevelIndex
00092 ) const
00093 {
00094     const labelList& fineToCoarse = restrictAddressing_[coarseLevelIndex];
00095 
00096     forAll(fineToCoarse, i)
00097     {
00098         ff[i] = cf[fineToCoarse[i]];
00099     }
00100 }
00101 
00102 
00103 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines