/usr/include/palabos/offLattice/offLatticeModel3D.h is in libplb-dev 1.5~r1+repack1-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 | /* This file is part of the Palabos library.
*
* Copyright (C) 2011-2015 FlowKit Sarl
* Route d'Oron 2
* 1010 Lausanne, Switzerland
* E-mail contact: contact@flowkit.com
*
* The most recent release of Palabos can be downloaded at
* <http://www.palabos.org/>
*
* The library Palabos is free software: you can redistribute it and/or
* modify it under the terms of the GNU Affero General Public License as
* published by the Free Software Foundation, either version 3 of the
* License, or (at your option) any later version.
*
* The library 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 Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef OFF_LATTICE_MODEL_3D_H
#define OFF_LATTICE_MODEL_3D_H
#include "core/globalDefs.h"
namespace plb {
template< typename T, class SurfaceData >
class OffLatticeModel3D {
public:
OffLatticeModel3D(BoundaryShape3D<T,SurfaceData>* shape_, int flowType_);
OffLatticeModel3D(OffLatticeModel3D<T,SurfaceData> const& rhs);
OffLatticeModel3D<T,SurfaceData>& operator= (
OffLatticeModel3D<T,SurfaceData> const& rhs );
virtual ~OffLatticeModel3D();
void provideShapeArguments(std::vector<AtomicBlock3D*> args);
plint getTag(plint id) const;
bool pointOnSurface (
Dot3D const& fromPoint, Dot3D const& direction,
Array<T,3>& locatedPoint, T& distance,
Array<T,3>& wallNormal, SurfaceData& surfaceData,
OffBoundary::Type& bdType, plint& id ) const;
Array<T,3> computeContinuousNormal (
Array<T,3> const& p, plint id, bool isAreaWeighted = false ) const;
bool intersectsSurface (
Dot3D const& p1, Dot3D const& p2, plint& id ) const;
bool isFluid(Dot3D const& location) const;
bool velIsJ() const { return velIsJflag; }
void setVelIsJ(bool velIsJflag_) { velIsJflag = velIsJflag_; }
bool getPartialReplace() const { return partialReplaceFlag; }
void setPartialReplace(bool prFlag) { partialReplaceFlag = prFlag; }
virtual OffLatticeModel3D<T,SurfaceData>* clone() const =0;
virtual plint getNumNeighbors() const =0;
virtual void prepareCell (
Dot3D const& cellLocation, AtomicContainerBlock3D& container ) =0;
virtual void boundaryCompletion (
AtomicBlock3D& lattice,
AtomicContainerBlock3D& container,
std::vector<AtomicBlock3D const*> const& args ) =0;
virtual ContainerBlockData* generateOffLatticeInfo() const =0;
virtual Array<T,3> getLocalForce(AtomicContainerBlock3D& container) const =0;
private:
BoundaryShape3D<T,SurfaceData>* shape;
int flowType;
bool velIsJflag;
bool partialReplaceFlag;
};
/// Precompute a list of nodes which are close to the off-lattice boundary
/// (both wet and dry ones).
template<typename T, class SurfaceData>
class OffLatticePatternFunctional3D : public BoxProcessingFunctional3D
{
public:
/// numNeighbors is the number of neighbors a boundary node must be able
/// to access along a given direction.
OffLatticePatternFunctional3D (
OffLatticeModel3D<T,SurfaceData>* offLatticeModel_ );
virtual ~OffLatticePatternFunctional3D();
OffLatticePatternFunctional3D(OffLatticePatternFunctional3D const& rhs);
OffLatticePatternFunctional3D& operator= (
OffLatticePatternFunctional3D const& rhs );
void swap(OffLatticePatternFunctional3D& rhs);
virtual OffLatticePatternFunctional3D<T,SurfaceData>* clone() const;
/// First AtomicBlock: OffLatticeInfo.
/// If there are more atomic-blocks then they are forwarded to the
/// shape function, to provide additional read-only parameters.
virtual void processGenericBlocks(Box3D domain, std::vector<AtomicBlock3D*> fields);
virtual void getTypeOfModification(std::vector<modif::ModifT>& modified) const;
virtual BlockDomain::DomainT appliesTo() const;
private:
OffLatticeModel3D<T,SurfaceData>* offLatticeModel;
};
template<typename T, template<typename U> class Descriptor, class SurfaceData>
class OffLatticeCompletionFunctional3D : public BoxProcessingFunctional3D
{
public:
OffLatticeCompletionFunctional3D (
OffLatticeModel3D<T,SurfaceData>* offLatticeModel_,
plint numShapeArgs_, plint numCompletionArgs_ );
virtual ~OffLatticeCompletionFunctional3D();
OffLatticeCompletionFunctional3D(OffLatticeCompletionFunctional3D<T,Descriptor,SurfaceData> const& rhs);
OffLatticeCompletionFunctional3D& operator= (
OffLatticeCompletionFunctional3D<T,Descriptor,SurfaceData> const& rhs );
void swap(OffLatticeCompletionFunctional3D<T,Descriptor,SurfaceData>& rhs);
/// First AtomicBlock: Lattice; Second AtomicBlock: Container for off-
/// lattice info.
/// If there are more atomic-blocks then they are forwarded to the
/// shape function, to provide additional read-only parameters.
virtual void processGenericBlocks(Box3D domain, std::vector<AtomicBlock3D*> fields);
virtual OffLatticeCompletionFunctional3D<T,Descriptor,SurfaceData>* clone() const;
virtual void getTypeOfModification(std::vector<modif::ModifT>& modified) const;
virtual BlockDomain::DomainT appliesTo() const;
private:
OffLatticeModel3D<T,SurfaceData>* offLatticeModel;
plint numShapeArgs, numCompletionArgs;
};
template< typename T, class SurfaceData >
class GetForceOnObjectFunctional3D : public PlainReductiveBoxProcessingFunctional3D
{
public:
GetForceOnObjectFunctional3D (
OffLatticeModel3D<T,SurfaceData>* offLatticeModel_ );
virtual ~GetForceOnObjectFunctional3D();
GetForceOnObjectFunctional3D(GetForceOnObjectFunctional3D<T,SurfaceData> const& rhs);
GetForceOnObjectFunctional3D<T,SurfaceData>& operator= (
GetForceOnObjectFunctional3D<T,SurfaceData> const& rhs );
/// First AtomicBlock: Container for off-lattice info.
virtual void processGenericBlocks(Box3D domain, std::vector<AtomicBlock3D*> fields);
virtual GetForceOnObjectFunctional3D<T,SurfaceData>* clone() const;
virtual void getTypeOfModification(std::vector<modif::ModifT>& modified) const;
virtual BlockDomain::DomainT appliesTo() const;
Array<T,3> getForce() const;
private:
OffLatticeModel3D<T,Array<T,3> >* offLatticeModel;
Array<plint,3> forceId;
};
} // namespace plb
#endif // OFF_LATTICE_MODEL_3D_H
|