mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
ENH: PointEdgePoint: transport nearest with passive data, using PointEdgeWave
This commit is contained in:
mattijs
55
src/meshTools/algorithms/PointEdgeWave/PointData.C
Normal file
55
src/meshTools/algorithms/PointEdgeWave/PointData.C
Normal file
@ -0,0 +1,55 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
|
||||||
|
\\/ M anipulation | Copyright (C) 2015 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 "PointData.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
Foam::Ostream& Foam::operator<<(Ostream& os, const PointData<DataType>& pd)
|
||||||
|
{
|
||||||
|
if (os.format() == IOstream::ASCII)
|
||||||
|
{
|
||||||
|
return os
|
||||||
|
<< static_cast<const pointEdgePoint&>(pd)
|
||||||
|
<< token::SPACE << pd.data();
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
return os
|
||||||
|
<< static_cast<const pointEdgePoint&>(pd)
|
||||||
|
<< pd.data();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
Foam::Istream& Foam::operator>>(Istream& is, PointData<DataType>& pd)
|
||||||
|
{
|
||||||
|
return is >> static_cast<pointEdgePoint&>(pd) >> pd.data_;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
211
src/meshTools/algorithms/PointEdgeWave/PointData.H
Normal file
211
src/meshTools/algorithms/PointEdgeWave/PointData.H
Normal file
@ -0,0 +1,211 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
|
||||||
|
\\/ M anipulation | Copyright (C) 2015 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/>.
|
||||||
|
|
||||||
|
Class
|
||||||
|
Foam::PointData
|
||||||
|
|
||||||
|
Description
|
||||||
|
Variant of pointEdgePoint with some transported additional data. Templated
|
||||||
|
on the transported data type.
|
||||||
|
|
||||||
|
SourceFiles
|
||||||
|
PointDataI.H
|
||||||
|
PointData.C
|
||||||
|
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
#ifndef PointData_H
|
||||||
|
#define PointData_H
|
||||||
|
|
||||||
|
#include "pointEdgePoint.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
namespace Foam
|
||||||
|
{
|
||||||
|
|
||||||
|
// Forward declaration of classes
|
||||||
|
class Istream;
|
||||||
|
class Ostream;
|
||||||
|
template<class DataType>
|
||||||
|
class PointData;
|
||||||
|
|
||||||
|
// Forward declaration of friend functions and operators
|
||||||
|
template<class DataType>
|
||||||
|
Ostream& operator<<(Ostream&, const PointData<DataType>&);
|
||||||
|
template<class DataType>
|
||||||
|
Istream& operator>>(Istream&, PointData<DataType>&);
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
Class PointData Declaration
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
class PointData
|
||||||
|
:
|
||||||
|
public pointEdgePoint
|
||||||
|
{
|
||||||
|
private:
|
||||||
|
|
||||||
|
// Private data
|
||||||
|
|
||||||
|
//- Additional transported data
|
||||||
|
DataType data_;
|
||||||
|
|
||||||
|
|
||||||
|
public:
|
||||||
|
|
||||||
|
// Constructors
|
||||||
|
|
||||||
|
//- Construct null
|
||||||
|
inline PointData();
|
||||||
|
|
||||||
|
//- Construct from origin, distance and data
|
||||||
|
inline PointData
|
||||||
|
(
|
||||||
|
const point& origin,
|
||||||
|
const scalar distSqr,
|
||||||
|
const DataType& data
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Construct as copy
|
||||||
|
inline PointData(const PointData&);
|
||||||
|
|
||||||
|
|
||||||
|
// Member Functions
|
||||||
|
|
||||||
|
// Access
|
||||||
|
|
||||||
|
//- Const access the data
|
||||||
|
inline const DataType& data() const;
|
||||||
|
|
||||||
|
|
||||||
|
// Needed by meshWave
|
||||||
|
|
||||||
|
//- Apply rotation matrix to the data
|
||||||
|
template<class TrackingData>
|
||||||
|
inline void transform
|
||||||
|
(
|
||||||
|
const tensor& rotTensor,
|
||||||
|
TrackingData& td
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Influence of edge on point
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool updatePoint
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label pointI,
|
||||||
|
const label edgeI,
|
||||||
|
const PointData<DataType>& edgeInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Influence of different value on same point.
|
||||||
|
// Merge new and old info.
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool updatePoint
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label pointI,
|
||||||
|
const PointData<DataType>& newPointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Influence of different value on same point.
|
||||||
|
// No information about current position whatsoever.
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool updatePoint
|
||||||
|
(
|
||||||
|
const PointData<DataType>& newPointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
);
|
||||||
|
|
||||||
|
//- Influence of point on edge.
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool updateEdge
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label edgeI,
|
||||||
|
const label pointI,
|
||||||
|
const PointData<DataType>& pointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
);
|
||||||
|
|
||||||
|
// Member Operators
|
||||||
|
|
||||||
|
inline bool operator==(const PointData<DataType>&) const;
|
||||||
|
inline bool operator!=(const PointData<DataType>&) const;
|
||||||
|
|
||||||
|
|
||||||
|
// IOstream Operators
|
||||||
|
|
||||||
|
friend Ostream& operator<< <DataType>
|
||||||
|
(
|
||||||
|
Ostream&,
|
||||||
|
const PointData<DataType>&
|
||||||
|
);
|
||||||
|
friend Istream& operator>> <DataType>
|
||||||
|
(
|
||||||
|
Istream&,
|
||||||
|
PointData<DataType>&
|
||||||
|
);
|
||||||
|
};
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
//- Data associated with PointData types is contiguous
|
||||||
|
|
||||||
|
template<>
|
||||||
|
inline bool contiguous<PointData<scalar> >()
|
||||||
|
{
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<>
|
||||||
|
inline bool contiguous<PointData<vector> >()
|
||||||
|
{
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
} // End namespace Foam
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#include "PointDataI.H"
|
||||||
|
|
||||||
|
#ifdef NoRepository
|
||||||
|
# include "PointData.C"
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
#endif
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
234
src/meshTools/algorithms/PointEdgeWave/PointDataI.H
Normal file
234
src/meshTools/algorithms/PointEdgeWave/PointDataI.H
Normal file
@ -0,0 +1,234 @@
|
|||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
========= |
|
||||||
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||||
|
\\ / O peration |
|
||||||
|
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
|
||||||
|
\\/ M anipulation | Copyright (C) 2015 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 "polyMesh.H"
|
||||||
|
#include "transform.H"
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline Foam::PointData<DataType>::PointData()
|
||||||
|
:
|
||||||
|
pointEdgePoint()
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline Foam::PointData<DataType>::PointData
|
||||||
|
(
|
||||||
|
const point& origin,
|
||||||
|
const scalar distSqr,
|
||||||
|
const DataType& data
|
||||||
|
)
|
||||||
|
:
|
||||||
|
pointEdgePoint(origin, distSqr),
|
||||||
|
data_(data)
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline Foam::PointData<DataType>::PointData(const PointData<DataType>& wpt)
|
||||||
|
:
|
||||||
|
pointEdgePoint(wpt),
|
||||||
|
data_(wpt.data())
|
||||||
|
{}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline const DataType& Foam::PointData<DataType>::data() const
|
||||||
|
{
|
||||||
|
return data_;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
template<class TrackingData>
|
||||||
|
inline void Foam::PointData<DataType>::transform
|
||||||
|
(
|
||||||
|
const tensor& rotTensor,
|
||||||
|
TrackingData& td
|
||||||
|
)
|
||||||
|
{
|
||||||
|
pointEdgePoint::transform(rotTensor, td);
|
||||||
|
data_ = Foam::transform(rotTensor, data_);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool Foam::PointData<DataType>::updatePoint
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label pointI,
|
||||||
|
const label edgeI,
|
||||||
|
const PointData<DataType>& edgeInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
)
|
||||||
|
{
|
||||||
|
if
|
||||||
|
(
|
||||||
|
pointEdgePoint::updatePoint
|
||||||
|
(
|
||||||
|
mesh,
|
||||||
|
pointI,
|
||||||
|
edgeI,
|
||||||
|
edgeInfo,
|
||||||
|
tol,
|
||||||
|
td
|
||||||
|
)
|
||||||
|
)
|
||||||
|
{
|
||||||
|
data_ = edgeInfo.data_;
|
||||||
|
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool Foam::PointData<DataType>::updatePoint
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label pointI,
|
||||||
|
const PointData<DataType>& newPointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
)
|
||||||
|
{
|
||||||
|
if
|
||||||
|
(
|
||||||
|
pointEdgePoint::updatePoint
|
||||||
|
(
|
||||||
|
mesh,
|
||||||
|
pointI,
|
||||||
|
newPointInfo,
|
||||||
|
tol,
|
||||||
|
td
|
||||||
|
)
|
||||||
|
)
|
||||||
|
{
|
||||||
|
data_ = newPointInfo.data_;
|
||||||
|
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool Foam::PointData<DataType>::updatePoint
|
||||||
|
(
|
||||||
|
const PointData<DataType>& newPointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
)
|
||||||
|
{
|
||||||
|
if (pointEdgePoint::updatePoint(newPointInfo, tol, td))
|
||||||
|
{
|
||||||
|
data_ = newPointInfo.data_;
|
||||||
|
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
template<class TrackingData>
|
||||||
|
inline bool Foam::PointData<DataType>::updateEdge
|
||||||
|
(
|
||||||
|
const polyMesh& mesh,
|
||||||
|
const label edgeI,
|
||||||
|
const label pointI,
|
||||||
|
const PointData<DataType>& pointInfo,
|
||||||
|
const scalar tol,
|
||||||
|
TrackingData& td
|
||||||
|
|
||||||
|
)
|
||||||
|
{
|
||||||
|
if
|
||||||
|
(
|
||||||
|
pointEdgePoint::updateEdge
|
||||||
|
(
|
||||||
|
mesh,
|
||||||
|
edgeI,
|
||||||
|
pointI,
|
||||||
|
pointInfo,
|
||||||
|
tol,
|
||||||
|
td
|
||||||
|
)
|
||||||
|
)
|
||||||
|
{
|
||||||
|
data_ = pointInfo.data_;
|
||||||
|
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
else
|
||||||
|
{
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline bool Foam::PointData<DataType>::operator==
|
||||||
|
(
|
||||||
|
const Foam::PointData<DataType>& rhs
|
||||||
|
)
|
||||||
|
const
|
||||||
|
{
|
||||||
|
return pointEdgePoint::operator==(rhs) && (data() == rhs.data());
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class DataType>
|
||||||
|
inline bool Foam::PointData<DataType>::operator!=
|
||||||
|
(
|
||||||
|
const Foam::PointData<DataType>& rhs
|
||||||
|
)
|
||||||
|
const
|
||||||
|
{
|
||||||
|
return !(*this == rhs);
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
Reference in New Issue
Block a user