26#ifndef SET_FD_BOUNDARY_2D_DEV03_HH
27#define SET_FD_BOUNDARY_2D_DEV03_HH
34template<
typename T,
typename DESCRIPTOR,
typename MODEL,
typename SCHEME_BOUND,
typename PARAMETERS,
typename FIELD,
typename SOURCE>
37 setFdBoundary2D<T,DESCRIPTOR,MODEL,SCHEME_BOUND,PARAMETERS,FIELD,SOURCE>(sLattice, superGeometry.
getMaterialIndicator(material));
41template<
typename T,
typename DESCRIPTOR,
typename MODEL,
typename SCHEME_BOUND,
typename PARAMETERS,
typename FIELD,
typename SOURCE>
44 bool includeOuterCells =
false;
45 int _overlap = indicator->getSuperGeometry().getOverlap();
47 if (indicator->getSuperGeometry().getOverlap() == 1) {
48 includeOuterCells =
true;
49 clout <<
"WARNING: overlap == 1, boundary conditions set on overlap despite unknown neighbor materials" << std::endl;
51 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
52 setFdBoundary2D<T,DESCRIPTOR,MODEL,SCHEME_BOUND,PARAMETERS,FIELD,SOURCE>(sLattice.
getBlock(iCloc),
53 indicator->getBlockIndicatorF(iCloc), includeOuterCells);
57 addPoints2CommBC<T,DESCRIPTOR>(sLattice, std::forward<
decltype(indicator)>(indicator), _overlap);
64template<
typename T,
typename DESCRIPTOR,
typename MODEL,
typename SCHEME_BOUND,
typename PARAMETERS,
typename FIELD,
typename SOURCE>
69 const int margin = includeOuterCells ? 0 : 1;
70 std::vector<int> discreteNormal(3, 0);
71 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY) {
72 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY}) >= margin
73 && indicator(iX, iY)) {
74 discreteNormal = indicator.
getBlockGeometry().getStatistics().getType(iX, iY);
75 block.
get(iX, iY).
template setField<descriptors::NORMAL_X>(-discreteNormal[1]);
76 block.
get(iX, iY).
template setField<descriptors::NORMAL_Y>(-discreteNormal[2]);
Base block indicator functor (discrete)
BlockGeometry< T, 2 > & getBlockGeometry()
Get underlying block geometry structure.
Platform-abstracted block lattice for external access and inter-block interaction.
virtual void addPostProcessor(std::type_index stage, LatticeR< DESCRIPTOR::d > latticeR, PostProcessorPromise< T, DESCRIPTOR > &&promise)=0
Schedule post processor for application to latticeR in stage.
Cell< T, DESCRIPTOR > get(CellID iCell)
Get Cell interface for index iCell.
Smart pointer for managing the various ways of passing functors around.
class for marking output with some text
Representation of a statistic for a parallel 2D geometry.
std::unique_ptr< SuperIndicatorF< T, D > > getMaterialIndicator(std::vector< int > &&materials)
Returns a material indicator using the given vector of materials.
Super class maintaining block lattices for a cuboid decomposition.
BlockLattice< T, DESCRIPTOR > & getBlock(int locC)
Return BlockLattice with local index locC.
LoadBalancer< T > & getLoadBalancer()
Read and write access to the load balancer.
Top level namespace for all of OpenLB.
void setFdBoundary2D(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 2 > &superGeometry, int material)
Initialising the setFdBoundary2D function on the superLattice domain Interpolated Boundaries use the ...
Communication after propagation.