src/meshLibrary/utilities/octrees/meshOctree/refinementControls/objectRefinement/sphereRefinement.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | cfMesh: A library for mesh generation
   4    \\    /   O peration     |
   5     \\  /    A nd           | Author: Franjo Juretic (franjo.juretic@c-fields.com)
   6      \\/     M anipulation  | Copyright (C) Creative Fields, Ltd.
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of cfMesh.
  10
  11     cfMesh is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     cfMesh 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 cfMesh.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "sphereRefinement.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "boundBox.H"
  29
  30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34
  35 defineTypeNameAndDebug(sphereRefinement, 0);
  36 addToRunTimeSelectionTable(objectRefinement, sphereRefinement, dictionary);
  37
  38 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  39
  40 sphereRefinement::sphereRefinement()
  41 :
  42     objectRefinement(),
  43     centre_(),
  44     radius_(-1.0)
  45 {}
  46
  47 sphereRefinement::sphereRefinement
  48 (
  49     const word& name,
  50     const scalar cellSize,
  51     const point& centre,
  52     const scalar radius
  53 )
  54 :
  55     objectRefinement(),
  56     centre_(centre),
  57     radius_(radius)
  58 {
  59     setName(name);
  60     setCellSize(cellSize);
  61 }
  62
  63 sphereRefinement::sphereRefinement
  64 (
  65     const word& name,
  66     const dictionary& dict
  67 )
  68 :
  69     objectRefinement(name, dict)
  70 {
  71     this->operator=(dict);
  72 }
  73
  74 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  75
  76 bool sphereRefinement::intersectsObject(const boundBox& bb) const
  77 {
  78     const point& c = (bb.max() + bb.min()) / 2.0;
  79
  80     if( magSqr(c - centre_) < sqr(radius_) )
  81         return true;
  82
  83     return false;
  84 }
  85
  86 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  87
  88 dictionary sphereRefinement::dict(bool ignoreType) const
  89 {
  90     dictionary dict;
  91
  92     dict.add("cellSize", cellSize());
  93     dict.add("type", type());
  94
  95     dict.add("centre", centre_);
  96     dict.add("radius", radius_);
  97
  98     return dict;
  99 }
 100
 101 void sphereRefinement::write(Ostream& os) const
 102 {
 103     os  << " type:   " << type()
 104         << " centre: " << centre_
 105         << " radius: " << radius_;
 106 }
 107
 108 void sphereRefinement::writeDict(Ostream& os, bool subDict) const
 109 {
 110     if( subDict )
 111     {
 112         os << indent << token::BEGIN_BLOCK << incrIndent << nl;
 113     }
 114
 115     os.writeKeyword("cellSize") << cellSize() << token::END_STATEMENT << nl;
 116
 117     // only write type for derived types
 118     if( type() != typeName_() )
 119     {
 120         os.writeKeyword("type") << type() << token::END_STATEMENT << nl;
 121     }
 122
 123     os.writeKeyword("centre") << centre_ << token::END_STATEMENT << nl;
 124     os.writeKeyword("radius") << radius_ << token::END_STATEMENT << nl;
 125
 126     if( subDict )
 127     {
 128         os << decrIndent << indent << token::END_BLOCK << endl;
 129     }
 130 }
 131
 132 void sphereRefinement::operator=(const dictionary& d)
 133 {
 134     // allow as embedded sub-dictionary "coordinateSystem"
 135     const dictionary& dict =
 136     (
 137         d.found(typeName_())
 138       ? d.subDict(typeName_())
 139       : d
 140     );
 141
 142     // unspecified centre is (0 0 0)
 143     if( dict.found("centre") )
 144     {
 145         dict.lookup("centre") >> centre_;
 146     }
 147     else
 148     {
 149         FatalErrorIn
 150         (
 151             "void sphereRefinement::operator=(const dictionary& d)"
 152         ) << "Entry centre is not specified!" << exit(FatalError);
 153         centre_ = vector::zero;
 154     }
 155
 156     // specify radius
 157     if( dict.found("radius") )
 158     {
 159         radius_ = readScalar(dict.lookup("radius"));
 160     }
 161     else
 162     {
 163         FatalErrorIn
 164         (
 165             "void sphereRefinement::operator=(const dictionary& d)"
 166         ) << "Entry radius is not specified!" << exit(FatalError);
 167         radius_ = -1.0;
 168     }
 169 }
 170
 171 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 172
 173 Ostream& sphereRefinement::operator<<(Ostream& os) const
 174 {
 175     os << "name " << name() << nl;
 176     os << "cell size " << cellSize() << nl;
 177     write(os);
 178     return os;
 179 }
 180
 181 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 182
 183 } // End namespace Foam
 184
 185 // ************************************************************************* //