src/OpenFOAM/meshes/polyMesh/polyMeshTetDecomposition/tetIndices.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  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
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "tetIndices.H"
  27
  28 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  29
  30 Foam::tetIndices::tetIndices()
  31 :
  32     cellI_(-1),
  33     faceI_(-1),
  34     faceBasePtI_(-1),
  35     facePtAI_(-1),
  36     facePtBI_(-1),
  37     tetPtI_(-1)
  38 {}
  39
  40
  41 Foam::tetIndices::tetIndices
  42 (
  43     label cellI,
  44     label faceI,
  45     label faceBasePtI,
  46     label facePtAI,
  47     label facePtBI,
  48     label tetPtI
  49 )
  50 :
  51     cellI_(cellI),
  52     faceI_(faceI),
  53     faceBasePtI_(faceBasePtI),
  54     facePtAI_(facePtAI),
  55     facePtBI_(facePtBI),
  56     tetPtI_(tetPtI)
  57 {}
  58
  59
  60 Foam::tetIndices::tetIndices
  61 (
  62     label cellI,
  63     label faceI,
  64     label tetPtI,
  65     const polyMesh& mesh
  66 )
  67 :
  68     cellI_(cellI),
  69     faceI_(faceI),
  70     faceBasePtI_(-1),
  71     facePtAI_(-1),
  72     facePtBI_(-1),
  73     tetPtI_(tetPtI)
  74 {
  75     const faceList& pFaces = mesh.faces();
  76     const labelList& pOwner = mesh.faceOwner();
  77
  78     const Foam::face& f = pFaces[faceI_];
  79
  80     bool own = (pOwner[faceI_] == cellI_);
  81
  82     faceBasePtI_ = mesh.tetBasePtIs()[faceI_];
  83
  84     label facePtI = (tetPtI_ + faceBasePtI_) % f.size();
  85     label otherFacePtI = f.fcIndex(facePtI);
  86
  87     if (own)
  88     {
  89         facePtAI_ = facePtI;
  90         facePtBI_ = otherFacePtI;
  91     }
  92     else
  93     {
  94         facePtAI_ = otherFacePtI;
  95         facePtBI_ = facePtI;
  96     }
  97 }
  98
  99
 100 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 101
 102 Foam::tetIndices::~tetIndices()
 103 {}
 104
 105
 106 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
 107
 108 Foam::Istream& Foam::operator>>(Istream& is, tetIndices& tI)
 109 {
 110     is  >> tI.cell()
 111         >> tI.face()
 112         >> tI.faceBasePt()
 113         >> tI.facePtA()
 114         >> tI.facePtB()
 115         >> tI.tetPt();
 116
 117     // Check state of Istream
 118     is.check
 119     (
 120         "Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::tetIndices&)"
 121     );
 122
 123     return is;
 124 }
 125
 126
 127 Foam::Ostream& Foam::operator<<(Ostream& os, const tetIndices& tI)
 128 {
 129     os  << tI.cell() << token::SPACE
 130         << tI.face() << token::SPACE
 131         << tI.faceBasePt() << token::SPACE
 132         << tI.facePtA() << token::SPACE
 133         << tI.facePtB() << token::SPACE
 134         << tI.tetPt() << token::SPACE
 135         << endl;
 136
 137     // Check state of Ostream
 138     os.check
 139     (
 140         "Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
 141         "const Foam::tetIndices&)"
 142     );
 143
 144     return os;
 145 }
 146
 147
 148 // ************************************************************************* //