src/dynamicMesh/meshCut/refineCell/refineCell.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 "refineCell.H"
  27 #include "error.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 Foam::refineCell::refineCell()
  32 :
  33     cellNo_(-1),
  34     direction_(vector::zero)
  35 {}
  36
  37
  38 Foam::refineCell::refineCell(const label cellI, const vector& direction)
  39 :
  40     cellNo_(cellI),
  41     direction_(direction)
  42 {
  43     scalar magDir = mag(direction_);
  44
  45     if (magDir < SMALL)
  46     {
  47         FatalErrorIn("refineCell(const label, const vector&)")
  48             << "(almost)zero vector as direction for cell " << cellNo_
  49             << abort(FatalError);
  50     }
  51     else if (mag(magDir - 1) > SMALL)
  52     {
  53         // Normalize
  54         direction_ /= mag(direction_);
  55     }
  56 }
  57
  58
  59 Foam::refineCell::refineCell(Istream& is)
  60 :
  61     cellNo_(readLabel(is)),
  62     direction_(is)
  63 {
  64     scalar magDir = mag(direction_);
  65
  66     if (magDir < SMALL)
  67     {
  68         FatalErrorIn("refineCell(Istream&)")
  69             << "(almost)zero vector as direction for cell " << cellNo_
  70             << abort(FatalError);
  71     }
  72     else if (mag(magDir - 1) > SMALL)
  73     {
  74         // Normalize
  75         direction_ /= mag(direction_);
  76     }
  77 }
  78
  79
  80 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
  81
  82 Foam::Ostream& Foam::operator<<(Ostream& os, const refineCell& r)
  83 {
  84     if (os.format() == IOstream::ASCII)
  85     {
  86         os << r.cellNo() << token::SPACE << r.direction();
  87     }
  88     else
  89     {
  90         os << r.cellNo() << r.direction();
  91     }
  92     return os;
  93 }
  94
  95
  96 // ************************************************************************* //
  97