src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyMotion/sixDoFRigidBodyMotionIO.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "sixDoFRigidBodyMotion.H"
  28 #include "IOstreams.H"
  29
  30 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  31
  32 void Foam::sixDoFRigidBodyMotion::write(Ostream& os) const
  33 {
  34     motionState_.write(os);
  35
  36     os.writeKeyword("initialCentreOfMass")
  37         << initialCentreOfMass_ << token::END_STATEMENT << nl;
  38     os.writeKeyword("initialOrientation")
  39         << initialQ_ << token::END_STATEMENT << nl;
  40     os.writeKeyword("momentOfInertia")
  41         << momentOfInertia_ << token::END_STATEMENT << nl;
  42     os.writeKeyword("mass")
  43         << mass_ << token::END_STATEMENT << nl;
  44     os.writeKeyword("report")
  45         << report_ << token::END_STATEMENT << nl;
  46
  47     if (!restraints_.empty())
  48     {
  49         os  << indent << "restraints" << nl
  50             << indent << token::BEGIN_BLOCK << incrIndent << nl;
  51
  52         forAll(restraints_, rI)
  53         {
  54             word restraintType = restraints_[rI].type();
  55
  56             os  << indent << restraintNames_[rI] << nl
  57                 << indent << token::BEGIN_BLOCK << incrIndent << endl;
  58
  59             os.writeKeyword("sixDoFRigidBodyMotionRestraint")
  60                 << restraintType << token::END_STATEMENT << nl;
  61
  62             os.writeKeyword(word(restraintType + "Coeffs")) << nl;
  63
  64             os  << indent << token::BEGIN_BLOCK << nl << incrIndent;
  65
  66             restraints_[rI].write(os);
  67
  68             os << decrIndent << indent << token::END_BLOCK << nl;
  69
  70             os  << decrIndent << indent << token::END_BLOCK << endl;
  71         }
  72
  73         os  << decrIndent << indent << token::END_BLOCK << nl;
  74     }
  75
  76     if (!constraints_.empty())
  77     {
  78         os  << indent << "constraints" << nl
  79             << indent << token::BEGIN_BLOCK << incrIndent << nl;
  80
  81         os.writeKeyword("maxIterations")
  82             << maxConstraintIterations_ << token::END_STATEMENT << nl;
  83
  84         forAll(constraints_, rI)
  85         {
  86             word constraintType = constraints_[rI].type();
  87
  88             os  << indent << constraintNames_[rI] << nl
  89                 << indent << token::BEGIN_BLOCK << incrIndent << endl;
  90
  91             os.writeKeyword("sixDoFRigidBodyMotionConstraint")
  92                 << constraintType << token::END_STATEMENT << nl;
  93
  94             constraints_[rI].sixDoFRigidBodyMotionConstraint::write(os);
  95
  96             os.writeKeyword(word(constraintType + "Coeffs")) << nl;
  97
  98             os  << indent << token::BEGIN_BLOCK << nl << incrIndent;
  99
 100             constraints_[rI].write(os);
 101
 102             os << decrIndent << indent << token::END_BLOCK << nl;
 103
 104             os  << decrIndent << indent << token::END_BLOCK << endl;
 105         }
 106
 107         os  << decrIndent << indent << token::END_BLOCK << nl;
 108     }
 109 }
 110
 111
 112 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
 113
 114 Foam::Istream& Foam::operator>>(Istream& is, sixDoFRigidBodyMotion& sDoFRBM)
 115 {
 116     is  >> sDoFRBM.motionState_
 117         >> sDoFRBM.initialCentreOfMass_
 118         >> sDoFRBM.initialQ_
 119         >> sDoFRBM.momentOfInertia_
 120         >> sDoFRBM.mass_;
 121
 122     // Check state of Istream
 123     is.check
 124     (
 125         "Foam::Istream& Foam::operator>>"
 126         "(Foam::Istream&, Foam::sixDoFRigidBodyMotion&)"
 127     );
 128
 129     return is;
 130 }
 131
 132
 133 Foam::Ostream& Foam::operator<<
 134 (
 135     Ostream& os,
 136     const sixDoFRigidBodyMotion& sDoFRBM
 137 )
 138 {
 139     os  << sDoFRBM.motionState()
 140         << token::SPACE << sDoFRBM.initialCentreOfMass()
 141         << token::SPACE << sDoFRBM.initialQ()
 142         << token::SPACE << sDoFRBM.momentOfInertia()
 143         << token::SPACE << sDoFRBM.mass();
 144
 145     // Check state of Ostream
 146     os.check
 147     (
 148         "Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
 149         "const Foam::sixDoFRigidBodyMotion&)"
 150     );
 151
 152     return os;
 153 }
 154
 155
 156 // ************************************************************************* //