src/fileFormats/starcd/STARCDCore.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 "STARCDCore.H"
  27 #include "ListOps.H"
  28 #include "clock.H"
  29 #include "PackedBoolList.H"
  30 #include "IStringStream.H"
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 Foam::fileFormats::STARCDCore::STARCDCore()
  35 {}
  36
  37
  38 // * * * * * * * * * * *  Protected Member Functions * * * * * * * * * * * * //
  39
  40 bool Foam::fileFormats::STARCDCore::readHeader
  41 (
  42     IFstream& is,
  43     const word& signature
  44 )
  45 {
  46     if (!is.good())
  47     {
  48         FatalErrorIn
  49         (
  50             "fileFormats::STARCDCore::readHeader(...)"
  51         )
  52             << "cannot read " << signature  << "  " << is.name()
  53             << abort(FatalError);
  54     }
  55
  56     word header;
  57     label majorVersion;
  58
  59     string line;
  60
  61     is.getLine(line);
  62     IStringStream(line)() >> header;
  63
  64     is.getLine(line);
  65     IStringStream(line)() >> majorVersion;
  66
  67     // add other checks ...
  68     if (header != signature)
  69     {
  70         Info<< "header mismatch " << signature << "  " << is.name()
  71             << endl;
  72     }
  73
  74     return true;
  75 }
  76
  77
  78 void Foam::fileFormats::STARCDCore::writeHeader
  79 (
  80     Ostream& os,
  81     const word& filetype
  82 )
  83 {
  84     os  << "PROSTAR_" << filetype << nl
  85         << 4000
  86         << " " << 0
  87         << " " << 0
  88         << " " << 0
  89         << " " << 0
  90         << " " << 0
  91         << " " << 0
  92         << " " << 0
  93         << endl;
  94 }
  95
  96
  97 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  98
  99 bool Foam::fileFormats::STARCDCore::readPoints
 100 (
 101     IFstream& is,
 102     pointField& points,
 103     labelList& ids
 104 )
 105 {
 106     if (!is.good())
 107     {
 108         FatalErrorIn
 109         (
 110             "fileFormats::STARCDedgeFormat::readPoints(...)"
 111         )
 112             << "Cannot read file " << is.name()
 113             << exit(FatalError);
 114     }
 115
 116     readHeader(is, "PROSTAR_VERTEX");
 117
 118
 119     // reuse memory if possible
 120     DynamicList<point> dynPoints(points.xfer());
 121     DynamicList<label> dynPointId(ids.xfer());    // STAR-CD index of points
 122
 123     dynPoints.clear();
 124     dynPointId.clear();
 125
 126     label lineLabel;
 127     while ((is >> lineLabel).good())
 128     {
 129         scalar x, y, z;
 130
 131         is >> x >> y >> z;
 132
 133         dynPoints.append(point(x, y, z));
 134         dynPointId.append(lineLabel);
 135     }
 136
 137     points.transfer(dynPoints);
 138     ids.transfer(dynPointId);
 139
 140     return true;
 141 }
 142
 143
 144 void Foam::fileFormats::STARCDCore::writePoints
 145 (
 146     Ostream& os,
 147     const pointField& pointLst
 148 )
 149 {
 150     writeHeader(os, "VERTEX");
 151
 152     // Set the precision of the points data to 10
 153     os.precision(10);
 154
 155     // force decimal point for Fortran input
 156     os.setf(std::ios::showpoint);
 157
 158     forAll(pointLst, ptI)
 159     {
 160         os
 161             << ptI + 1 << " "
 162             << pointLst[ptI].x() << " "
 163             << pointLst[ptI].y() << " "
 164             << pointLst[ptI].z() << nl;
 165     }
 166     os.flush();
 167 }
 168
 169
 170
 171
 172 // ************************************************************************* //