26#ifndef SET_EXT_FIELD_BOUNDARY_3D_HH
27#define SET_EXT_FIELD_BOUNDARY_3D_HH
34template <
typename T,
typename DESCRIPTOR,
typename FIELD_A,
typename FIELD_B>
37 setExtFieldBoundary<T,DESCRIPTOR,FIELD_A,FIELD_B>(sLattice,
42template <
typename T,
typename DESCRIPTOR,
typename FIELD_A,
typename FIELD_B>
47 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); ++iCloc) {
48 setExtFieldBoundary<T,DESCRIPTOR,FIELD_A,FIELD_B>(
50 indicator->getBlockIndicatorF(iCloc));
53 addPoints2CommBC(sLattice,std::forward<
decltype(indicator)>(indicator), _overlap);
58template <
typename T,
typename DESCRIPTOR,
typename FIELD_A,
typename FIELD_B>
64 std::vector<int> discreteNormal(4, 0);
65 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
66 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
67 && indicator(iX, iY, iZ)) {
68 discreteNormal = blockGeometryStructure.getStatistics().getType(iX, iY, iZ);
69 if (discreteNormal[1]!=0 || discreteNormal[2]!=0 || discreteNormal[3]!=0) {
70 auto postProcessor = std::unique_ptr<PostProcessorGenerator3D<T, DESCRIPTOR>>{
71 new ExtFieldBoundaryProcessorGenerator3D<T,DESCRIPTOR,FIELD_A,FIELD_B>(iX, iX, iY, iY, iZ, iZ, -discreteNormal[1], -discreteNormal[2], -discreteNormal[3]) };
77 clout <<
"Warning: Could not setExternalFieldBoundary (" << iX <<
", " << iY <<
", " << iZ <<
"), discreteNormal=(" << discreteNormal[0] <<
","<< discreteNormal[1] <<
","<< discreteNormal[2] <<
"," << discreteNormal[3] <<
")" << std::endl;
Base block indicator functor.
BlockGeometry< T, 3 > & 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.
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.
Base indicator functor (discrete)
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 setExtFieldBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 3 > &superGeometry, int material)
Initialising the ExternalFieldBoundary on the superLattice domain This is an AdvectionDiffusionBounda...
void addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.