25#ifndef LATTICE_PHYS_BOUNDARY_FORCE_3D_HH
26#define LATTICE_PHYS_BOUNDARY_FORCE_3D_HH
44template<
typename T,
typename DESCRIPTOR>
50 _indicatorF(std::move(indicatorF))
52 this->
getName() =
"physBoundaryForce";
53 for (
int iC = 0; iC < this->
_sLattice.getLoadBalancer().size(); ++iC) {
57 _indicatorF->getBlockIndicatorF(iC),
62template<
typename T,
typename DESCRIPTOR>
68 superGeometry.getMaterialIndicator(material),
72template <
typename T,
typename DESCRIPTOR>
78 _indicatorF(indicatorF),
79 _blockGeometry(indicatorF.getBlockGeometry())
81 this->
getName() =
"physBoundaryForce";
84template<
typename T,
typename DESCRIPTOR>
87 for (
int i = 0; i < this->getTargetDim(); ++i) {
91 if (_indicatorF(input)) {
92 for (
int iPop = 1; iPop < DESCRIPTOR::q; ++iPop) {
97 if (_blockGeometry.get({input[0] + c[0], input[1] + c[1], input[2] + c[2]}) == 1) {
99 T f = this->_blockLattice.get(input[0] + c[0], input[1] + c[1], input[2] + c[2])[iPop];
102 f += this->_blockLattice.get(input)[descriptors::opposite<DESCRIPTOR>(iPop)];
104 for (
int i = 0; i < this->getTargetDim(); ++i) {
105 output[i] -= c[i] * f;
109 for (
int i = 0; i < this->getTargetDim(); ++i) {
110 output[i] = this->_converter.getPhysForce(output[i]);
Base block indicator functor.
functor returns pointwise phys force acting on a boundary with a given material on local lattice
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
BlockLatticePhysBoundaryForce3D(BlockLattice< T, DESCRIPTOR > &blockLattice, BlockIndicatorF3D< T > &indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Platform-abstracted block lattice for external access and inter-block interaction.
Smart pointer for managing the various ways of passing functors around.
std::string & getName()
read and write access to name
std::vector< std::unique_ptr< BlockF3D< T > > > _blockF
Super functors may consist of several BlockF3D<W> derived functors.
Representation of a statistic for a parallel 2D geometry.
Base indicator functor (discrete)
SuperLattice< T, DESCRIPTOR > & _sLattice
functor to get pointwise phys force acting on a boundary with a given material on local lattice
SuperLatticePhysBoundaryForce3D(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF3D< T > > &&indicatorF, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
Wrapper functions that simplify the use of MPI.
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.