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

globalIndex.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/globalIndex.H>
00027 
00028 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00029 
00030 Foam::globalIndex::globalIndex(const label localSize)
00031 :
00032     offsets_(Pstream::nProcs())
00033 {
00034     labelList localSizes(Pstream::nProcs());
00035     localSizes[Pstream::myProcNo()] = localSize;
00036     Pstream::gatherList(localSizes);
00037     Pstream::scatterList(localSizes);   // just to balance out comms
00038 
00039     label offset = 0;
00040     forAll(offsets_, procI)
00041     {
00042         label oldOffset = offset;
00043         offset += localSizes[procI];
00044 
00045         if (offset < oldOffset)
00046         {
00047             FatalErrorIn("globalIndex::globalIndex(const label)")
00048                 << "Overflow : sum of sizes " << localSizes
00049                 << " exceeds capability of label (" << labelMax
00050                 << "). Please recompile with larger datatype for label."
00051                 << exit(FatalError);
00052         }
00053         offsets_[procI] = offset;
00054     }
00055 }
00056 
00057 
00058 Foam::globalIndex::globalIndex(Istream& is)
00059 {
00060     is >> offsets_;
00061 }
00062 
00063 
00064 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
00065 
00066 Foam::Istream& Foam::operator>>(Istream& is, globalIndex& gi)
00067 {
00068     return is >> gi.offsets_;
00069 }
00070 
00071 
00072 Foam::Ostream& Foam::operator<<(Ostream& os, const globalIndex& gi)
00073 {
00074     return os << gi.offsets_;
00075 }
00076 
00077 
00078 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines