26#ifndef SET_ADVECTION_DIFFUSION_ZERO_GRADIENT_BOUNDARY_3D_HH
27#define SET_ADVECTION_DIFFUSION_ZERO_GRADIENT_BOUNDARY_3D_HH
34template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
37 setZeroGradientBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice, superGeometry.
getMaterialIndicator(material));
41template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
47 for (
int iCloc = 0; iCloc < sLattice.
getLoadBalancer().size(); iCloc++) {
48 setZeroGradientBoundary<T,DESCRIPTOR,MixinDynamics>(sLattice.
getBlock(iCloc), indicator->getBlockIndicatorF(iCloc));
51 addPoints2CommBC(sLattice,std::forward<
decltype(indicator)>(indicator), _overlap);
56template<
typename T,
typename DESCRIPTOR,
typename MixinDynamics>
59 using namespace boundaryhelper;
62 std::vector<int> discreteNormal(4,0);
63 blockGeometryStructure.forSpatialLocations([&](
auto iX,
auto iY,
auto iZ) {
64 if (blockGeometryStructure.getNeighborhoodRadius({iX, iY, iZ}) >= margin
65 && indicator(iX, iY, iZ)) {
66 discreteNormal = blockGeometryStructure.getStatistics().getType(iX, iY, iZ);
67 if (discreteNormal[1]!=0 || discreteNormal[2]!=0 || discreteNormal[3]!=0) {
68 for(
int iPop = 0; iPop < DESCRIPTOR::q; iPop++){
69 auto c = descriptors::c<DESCRIPTOR>(iPop);
70 T k = c[0]*(-discreteNormal[1]) + c[1]*(-discreteNormal[2]) + c[2]*(-discreteNormal[3]);
71 if(k > T{0.} && blockGeometryStructure.getMaterial(iX+c[0], iY+c[1], iZ+c[2]) == 1 ){
72 _block.
get(iX+c[0], iY+c[1], iZ+c[2]).template getFieldPointer<descriptors::NEIGHBOR>()[0] = T{1.};
74 if(k > T{0.} && blockGeometryStructure.getMaterial(iX+2*c[0], iY+2*c[1], iZ+2*c[2]) == 1 ){
75 _block.
get(iX+2*c[0], iY+2*c[1], iZ+2*c[2]).template getFieldPointer<descriptors::NEIGHBOR>()[0] = T{1.};
82 _block.template defineDynamics<NoCollideDynamics>({iX, iY, iZ});
84 _block.template defineDynamics<BounceBack>({iX, iY, iZ});
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.
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.
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 addPoints2CommBC(SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperIndicatorF2D< T > > &&indicator, int _overlap)
Adds needed Cells to the Communicator _commBC in SuperLattice.
void setZeroGradientBoundary(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 3 > &superGeometry, int material)
Initialising the setZeroGradientBoundary function on the superLattice domain This is an AdvectionDiff...
Communication after collision.