Files
openfoam/src/OpenFOAM/fields/pointPatchFields/basic/value/valuePointPatchField.C
Mark Olesen fc86e52451 ENH: extrapolateInternal() method for patch fields
- provides a more succinct way of writing

      {fa,fv}PatchField<Type>::patchInternalField(*this)

  as well as a consistent naming that can be used for patches derived
  from valuePointPatchField

ENH: readGradientEntry helper method for fixedGradient conditions

- simplifies coding and logic.
- support different read construct modes for fixedGradient
2023-06-19 17:47:03 +02:00

305 lines
6.6 KiB
C

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2011-2016 OpenFOAM Foundation
Copyright (C) 2019-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 "valuePointPatchField.H"
#include "pointPatchFieldMapper.H"
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
template<class Type>
bool Foam::valuePointPatchField<Type>::readValueEntry
(
const dictionary& dict,
IOobjectOption::readOption readOpt
)
{
if (!IOobjectOption::isAnyRead(readOpt)) return false;
const auto& p = pointPatchFieldBase::patch();
const auto* eptr = dict.findEntry("value", keyType::LITERAL);
if (eptr)
{
Field<Type>::assign(*eptr, p.size());
return true;
}
if (IOobjectOption::isReadRequired(readOpt))
{
FatalIOErrorInFunction(dict)
<< "Required entry 'value' : missing for patch " << p.name()
<< " in dictionary " << dict.relativeName() << nl
<< exit(FatalIOError);
}
return false;
}
template<class Type>
void Foam::valuePointPatchField<Type>::extrapolateInternal()
{
const labelUList& meshPoints = pointPatchFieldBase::patch().meshPoints();
const Field<Type>& iF = this->primitiveField();
Field<Type>& pfld = *this;
pfld.resize_nocopy(meshPoints.size()); // In general this is a no-op
forAll(meshPoints, i)
{
pfld[i] = iF[meshPoints[i]];
}
}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * * //
template<class Type>
Foam::valuePointPatchField<Type>::valuePointPatchField
(
const pointPatch& p,
const DimensionedField<Type, pointMesh>& iF
)
:
pointPatchField<Type>(p, iF),
Field<Type>(p.size())
{}
template<class Type>
Foam::valuePointPatchField<Type>::valuePointPatchField
(
const pointPatch& p,
const DimensionedField<Type, pointMesh>& iF,
const Type& value
)
:
pointPatchField<Type>(p, iF),
Field<Type>(p.size(), value)
{}
template<class Type>
Foam::valuePointPatchField<Type>::valuePointPatchField
(
const pointPatch& p,
const DimensionedField<Type, pointMesh>& iF,
const dictionary& dict,
IOobjectOption::readOption requireValue
)
:
pointPatchField<Type>(p, iF, dict),
Field<Type>(p.size())
{
if (!readValueEntry(dict, requireValue))
{
// Not read (eg, optional and missing): define zero
Field<Type>::operator=(Zero);
}
}
template<class Type>
Foam::valuePointPatchField<Type>::valuePointPatchField
(
const valuePointPatchField<Type>& ptf,
const pointPatch& p,
const DimensionedField<Type, pointMesh>& iF,
const pointPatchFieldMapper& mapper
)
:
pointPatchField<Type>(ptf, p, iF, mapper),
Field<Type>(ptf, mapper)
{}
template<class Type>
Foam::valuePointPatchField<Type>::valuePointPatchField
(
const valuePointPatchField<Type>& ptf,
const DimensionedField<Type, pointMesh>& iF
)
:
pointPatchField<Type>(ptf, iF),
Field<Type>(ptf)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
template<class Type>
void Foam::valuePointPatchField<Type>::autoMap
(
const pointPatchFieldMapper& m
)
{
Field<Type>::autoMap(m);
}
template<class Type>
void Foam::valuePointPatchField<Type>::rmap
(
const pointPatchField<Type>& ptf,
const labelList& addr
)
{
Field<Type>::rmap
(
refCast<const valuePointPatchField<Type>>
(
ptf
),
addr
);
}
template<class Type>
void Foam::valuePointPatchField<Type>::updateCoeffs()
{
if (this->updated())
{
return;
}
// Get internal field to insert values into
Field<Type>& iF = const_cast<Field<Type>&>(this->primitiveField());
this->setInInternalField(iF, *this);
pointPatchField<Type>::updateCoeffs();
}
template<class Type>
void Foam::valuePointPatchField<Type>::evaluate(const Pstream::commsTypes)
{
// Get internal field to insert values into
Field<Type>& iF = const_cast<Field<Type>&>(this->primitiveField());
this->setInInternalField(iF, *this);
pointPatchField<Type>::evaluate();
}
template<class Type>
void Foam::valuePointPatchField<Type>::write(Ostream& os) const
{
pointPatchField<Type>::write(os);
this->writeValueEntry(os);
}
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
template<class Type>
void Foam::valuePointPatchField<Type>::operator=
(
const valuePointPatchField<Type>& ptf
)
{
Field<Type>::operator=(ptf);
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator=
(
const pointPatchField<Type>& ptf
)
{
Field<Type>::operator=(this->patchInternalField());
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator=
(
const Field<Type>& tf
)
{
Field<Type>::operator=(tf);
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator=
(
const Type& t
)
{
Field<Type>::operator=(t);
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator==
(
const valuePointPatchField<Type>& ptf
)
{
Field<Type>::operator=(ptf);
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator==
(
const pointPatchField<Type>& ptf
)
{
Field<Type>::operator=(this->patchInternalField());
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator==
(
const Field<Type>& tf
)
{
Field<Type>::operator=(tf);
}
template<class Type>
void Foam::valuePointPatchField<Type>::operator==
(
const Type& t
)
{
Field<Type>::operator=(t);
}
// ************************************************************************* //