boundaryDataSurfaceReader.C

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
    Copyright (C) 2022-2023 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/>.

\*---------------------------------------------------------------------------*/

#include "boundaryDataSurfaceReader.H"
#include "rawIOField.H"
#include "Time.H"
#include "addToRunTimeSelectionTable.H"

// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //

namespace Foam
{
    defineTypeNameAndDebug(boundaryDataSurfaceReader, 0);
    addToRunTimeSelectionTable
    (
        surfaceReader,
        boundaryDataSurfaceReader,
        fileName
    );
}


// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //

Foam::pointField Foam::boundaryDataSurfaceReader::readPoints
(
    const Time& runTime,
    const fileName& baseDir,
    const word& pointsName
)
{
    fileName pointsFile
    (
        baseDir / (pointsName.empty() ? "points" : pointsName)
    );
    pointsFile.toAbsolute();

    IOobject io
    (
        pointsFile,             // absolute path
        runTime,
        IOobjectOption::MUST_READ,
        IOobjectOption::NO_WRITE,
        IOobjectOption::NO_REGISTER,
        true                    // global object (currently not used)
    );

    DebugInfo
        << "File: " << io.objectPath() << endl;

    // Read data (no average value!)
    rawIOField<point> rawData(io);

    pointField points(std::move(rawData.field()));

    DebugInfo
        << "File: " << io.objectPath()
        << " " << points.size() << " points" << endl;

    return points;
}


Foam::pointField Foam::boundaryDataSurfaceReader::readPoints
(
    const fileName& dirName,
    const word& pointsName
)
{
    refPtr<Time> timePtr(Time::NewGlobalTime());

    return readPoints(*timePtr, dirName, pointsName);
}


// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //

void Foam::boundaryDataSurfaceReader::readCase()
{
    DebugInFunction << endl;

    timeValues_ = TimePaths::findTimes(baseDir_);
}


void Foam::boundaryDataSurfaceReader::readGeometry
(
    meshedSurface& surf,
    const label timeIndex
)
{
    surf.clear();

    pointField points(this->readPoints(baseDir_, pointsName_));

    surf = meshedSurface(std::move(points), faceList());
}


template<class Type>
Foam::tmp<Foam::Field<Type>>
Foam::boundaryDataSurfaceReader::readFieldTemplate
(
    const label timeIndex,
    const label fieldIndex
) const
{
    Type dummyAvg;

    return readField<Type>
    (
        baseDir_,
        timeValues_[timeIndex],
        fieldNames_[fieldIndex],
        dummyAvg
    );
}


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

Foam::boundaryDataSurfaceReader::boundaryDataSurfaceReader
(
    const fileName& fName,
    const word& pointsName
)
:
    boundaryDataSurfaceReader(fName, dictionary(), pointsName)
{}


Foam::boundaryDataSurfaceReader::boundaryDataSurfaceReader
(
    const fileName& fName,
    const dictionary& options,
    const word& pointsName
)
:
    surfaceReader(fName, options),
    baseDir_(fName.path()),
    pointsName_(pointsName),
    timeValues_(),
    fieldNames_(),
    surfPtr_(nullptr)
{
    options.readIfPresent("points", pointsName_);

    baseDir_.toAbsolute();
    debug = 1;
    DebugInFunction << endl;
    Info<< "create with " << baseDir_ << endl;

    readCase();
}


// * * * * * * * * * * * * * Public Member Functions   * * * * * * * * * * * //

const Foam::meshedSurface& Foam::boundaryDataSurfaceReader::geometry
(
    const label timeIndex
)
{
    DebugInFunction << endl;

    if (!surfPtr_)
    {
        surfPtr_.reset(new meshedSurface);
        readGeometry(*surfPtr_, timeIndex);
    }

    return *surfPtr_;
}


Foam::instantList Foam::boundaryDataSurfaceReader::times() const
{
    return timeValues_;
}


Foam::wordList Foam::boundaryDataSurfaceReader::fieldNames
(
    const label timeIndex
) const
{
    if (timeValues_.empty() || timeIndex >= timeValues_.size())
    {
        fieldNames_.clear();
        return wordList();
    }

    fileNameList items =
        fileHandler().readDir(baseDir_/timeValues_[timeIndex].name());

    fieldNames_.resize_nocopy(items.size());

    std::transform
    (
        items.begin(),
        items.end(),
        fieldNames_.begin(),
        [](const fileName& f) { return word(f); }
    );

    Foam::sort(fieldNames_);

    return fieldNames_;
}


Foam::tmp<Foam::Field<Foam::scalar>> Foam::boundaryDataSurfaceReader::field
(
    const label timeIndex,
    const label fieldIndex,
    const scalar& refValue
) const
{
    return readFieldTemplate<scalar>(timeIndex, fieldIndex);
}


Foam::tmp<Foam::Field<Foam::vector>> Foam::boundaryDataSurfaceReader::field
(
    const label timeIndex,
    const label fieldIndex,
    const vector& refValue
) const
{
    return readFieldTemplate<vector>(timeIndex, fieldIndex);
}


Foam::tmp<Foam::Field<Foam::sphericalTensor>>
Foam::boundaryDataSurfaceReader::field
(
    const label timeIndex,
    const label fieldIndex,
    const sphericalTensor& refValue
) const
{
    return readFieldTemplate<sphericalTensor>(timeIndex, fieldIndex);
}


Foam::tmp<Foam::Field<Foam::symmTensor>> Foam::boundaryDataSurfaceReader::field
(
    const label timeIndex,
    const label fieldIndex,
    const symmTensor& refValue
) const
{
    return readFieldTemplate<symmTensor>(timeIndex, fieldIndex);
}


Foam::tmp<Foam::Field<Foam::tensor>> Foam::boundaryDataSurfaceReader::field
(
    const label timeIndex,
    const label fieldIndex,
    const tensor& refValue
) const
{
    return readFieldTemplate<tensor>(timeIndex, fieldIndex);
}


// ************************************************************************* //