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;
}


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