openfoam/applications/test/plotFunction1/Test-plotFunction1.C
Mark Olesen c9333a5ac8 ENH: improve flexibility of error, messageStream output
- provide a plain stream() method on messageStream to reduce reliance
  on casting operators and slightly opaque operator()() calls etc

- support alternative stream for messageStream serial output.
  This can be used to support local redirection of output.
  For example,

     refPtr<OFstream> logging;   // or autoPtr, unique_ptr etc

     // Later...
     Info.stream(logging.get())
        << "Detailed output ..." << endl;

  This will use the stdout semantics in the normal case, or allow
  redirection to an output file if a target output stream is defined,
  but still effectively use /dev/null on non-master processes.

  This is mostly the same as this ternary

      (logging ? *logging : Info())

  except that the ternary could be incorrect on sub-processes,
  requires more typing etc.

ENH: use case-relative names of dictionary, IOstream for FatalIOError

- normally yields more easily understandable information
2021-11-03 11:46:13 +01:00

106 lines
3.3 KiB
C

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2020-2021 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Application
Test-plotFunction1
Description
Plot scalar Function1 entries
\*---------------------------------------------------------------------------*/
#include "argList.H"
#include "fieldTypes.H"
#include "Function1.H"
#include "PtrList.H"
#include "Fstream.H"
using namespace Foam;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[])
{
argList::noBanner();
argList::noParallel();
argList::setAdvanced("case"); // Hide -case : has no meaning here
argList::addOption("begin", "scalar", "The start time (default: 0)");
argList::addOption("end", "scalar", "The end time (default: 1)");
argList::addOption("incr", "scalar", "The time increment (default: 0.1)");
argList::addOption("timeBase", "scalar", "The time base (default: 1)");
argList::addNote
(
"Read scalar functions from each file and produces"
" time/value output for each"
);
argList::noMandatoryArgs();
argList::addArgument("file1");
argList::addArgument("...");
argList::addArgument("fileN");
#include "setRootCase.H"
const scalar begTime = args.getOrDefault<scalar>("begin", 0);
const scalar endTime = args.getOrDefault<scalar>("end", 1);
const scalar incrTime = args.getOrDefault<scalar>("incr", 0.1);
const scalar timeBase = args.getOrDefault<scalar>("timeBase", 1);
Info.stream().precision(10);
for (label argi=1; argi < args.size(); ++argi)
{
IFstream is(args.get<fileName>(argi));
dictionary dict(is);
for (const entry& dEntry : dict)
{
autoPtr<Function1<scalar>> funPtr =
Function1<scalar>::New(dEntry.keyword(), dict);
auto& fun = *funPtr;
InfoErr<< nl;
fun.writeData(InfoErr);
InfoErr<< nl;
Info<< nl << "# " << fun.type() << nl;
for (scalar t = begTime; t < endTime; t += incrTime)
{
Info<< t << tab << fun.value(t*timeBase) << nl;
}
Info<< nl;
}
}
return 0;
}
// ************************************************************************* //