src/surfMesh/MeshedSurface/MeshedSurfaceIO.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-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 "MeshedSurface.H"
  28 #include "boundBox.H"
  29 #include "Ostream.H"
  30
  31 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  32
  33 template<class Face>
  34 void Foam::MeshedSurface<Face>::writeStats(Ostream& os) const
  35 {
  36     os  << "points      : " << this->points().size() << nl;
  37     if (MeshedSurface<Face>::isTri())
  38     {
  39         os << "triangles   : " << this->size() << nl;
  40     }
  41     else
  42     {
  43         label nTri = 0;
  44         label nQuad = 0;
  45         forAll(*this, i)
  46         {
  47             const label n = this->operator[](i).size();
  48
  49             if (n == 3)
  50             {
  51                 nTri++;
  52             }
  53             else if (n == 4)
  54             {
  55                 nQuad++;
  56             }
  57         }
  58
  59         os  << "faces       : " << this->size()
  60             << "  (tri:" << nTri << " quad:" << nQuad
  61             << " poly:" << (this->size() - nTri - nQuad ) << ")" << nl;
  62     }
  63
  64     os  << "boundingBox : " << boundBox(this->points()) << endl;
  65 }
  66
  67
  68 // ************************************************************************* //