src/surfMesh/MeshedSurface/MeshedSurfaceIO.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "MeshedSurface.H"
  27 #include "boundBox.H"
  28 #include "Ostream.H"
  29
  30 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  31
  32 template<class Face>
  33 void Foam::MeshedSurface<Face>::writeStats(Ostream& os) const
  34 {
  35     os  << "points      : " << this->points().size() << nl;
  36     if (MeshedSurface<Face>::isTri())
  37     {
  38         os << "triangles   : " << this->size() << nl;
  39     }
  40     else
  41     {
  42         label nTri = 0;
  43         label nQuad = 0;
  44         forAll(*this, i)
  45         {
  46             const label n = this->operator[](i).size();
  47
  48             if (n == 3)
  49             {
  50                 nTri++;
  51             }
  52             else if (n == 4)
  53             {
  54                 nQuad++;
  55             }
  56         }
  57
  58         os  << "faces       : " << this->size()
  59             << "  (tri:" << nTri << " quad:" << nQuad
  60             << " poly:" << (this->size() - nTri - nQuad ) << ")" << nl;
  61     }
  62
  63     os  << "boundingBox : " << boundBox(this->points()) << endl;
  64 }
  65
  66
  67 // ************************************************************************* //