globalIndexI.H 3.5 KB
Newer Older
1 2 3 4
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
Mark Olesen's avatar
Mark Olesen committed
5
    \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

    OpenFOAM is free software; you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by the
    Free Software Foundation; either version 2 of the License, or (at your
    option) any later version.

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
    along with OpenFOAM; if not, write to the Free Software Foundation,
    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA

\*---------------------------------------------------------------------------*/

#include "ListOps.H"

// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

31 32 33 34 35 36 37
//inline const Foam::labelList& Foam::globalIndex::offsets() const
//{
//    return offsets_;
//}


inline Foam::label Foam::globalIndex::offset(const label procI) const
38
{
39
    return (procI == 0 ? 0 : offsets_[procI-1]);
40 41 42
}


mattijs's avatar
mattijs committed
43
inline Foam::label Foam::globalIndex::localSize(const label procI) const
mattijs's avatar
mattijs committed
44 45 46
{
    return
    (
mattijs's avatar
mattijs committed
47 48 49
        procI == 0
      ? offsets_[procI]
      : offsets_[procI] - offsets_[procI-1]
mattijs's avatar
mattijs committed
50 51 52 53
    );
}


mattijs's avatar
mattijs committed
54 55 56 57 58 59
inline Foam::label Foam::globalIndex::localSize() const
{
    return localSize(Pstream::myProcNo());
}


60 61 62 63 64 65
inline Foam::label Foam::globalIndex::size() const
{
    return offsets_[Pstream::nProcs()-1];
}


66 67 68 69 70 71 72 73 74 75
inline Foam::label Foam::globalIndex::toGlobal
(
    const label procI,
    const label i
) const
{
    return(procI == 0 ? i : i + offsets_[procI-1]);
}


76 77
inline Foam::label Foam::globalIndex::toGlobal(const label i) const
{
78
    return toGlobal(Pstream::myProcNo(), i);
79 80
}

81

82
//- Is on local processor
83
inline bool Foam::globalIndex::isLocal(const label procI, const label i) const
84
{
Mattijs Janssens's avatar
Mattijs Janssens committed
85
    return
86 87 88 89 90 91 92 93
        (i < offsets_[procI])
     && (i >= (procI == 0 ? 0 : offsets_[procI-1]));
}


inline bool Foam::globalIndex::isLocal(const label i) const
{
    return isLocal(Pstream::myProcNo(), i);
94 95 96 97
}


inline Foam::label Foam::globalIndex::toLocal(const label procI, const label i)
98
const
99 100 101
{
    label localI = (procI == 0 ? i : i - offsets_[procI-1]);

Mattijs Janssens's avatar
Mattijs Janssens committed
102
    if (localI < 0 || i >= offsets_[procI])
103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135
    {
        FatalErrorIn("globalIndex::toLocal(const label, const label)")
            << "Global " << i << " does not belong on processor "
            << procI << endl << "Offsets:" << offsets_
            << abort(FatalError);
    }
    return localI;
}


inline Foam::label Foam::globalIndex::toLocal(const label i) const
{
    return toLocal(Pstream::myProcNo(), i);
}


inline Foam::label Foam::globalIndex::whichProcID(const label i) const
{
    label index = findLower(offsets_, i+1);

    if (index == Pstream::nProcs()-1)
    {
        FatalErrorIn("globalIndex::whichProcID(const label)")
            << "Global " << i << " does not belong on any processor."
            << " Offsets:" << offsets_
            << abort(FatalError);
    }

    return index+1;
}


// ************************************************************************* //