defineU2D.hh

Go to the documentation of this file.

/*  This file is part of the OpenLB library
 *
 *  Copyright (C) 2020 Alexander Schulz
 *  E-mail contact: info@openlb.net
 *  The most recent release of OpenLB can be downloaded at
 *  <http://www.openlb.net/>
 *
 *  This program is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU General Public License
 *  as published by the Free Software Foundation; either version 2
 *  of the License, or (at your option) any later version.
 *
 *  This program 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 General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public
 *  License along with this program; if not, write to the Free
 *  Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
 *  Boston, MA  02110-1301, USA.
*/
 
#ifndef DEFINE_U_2D_HH
#define DEFINE_U_2D_HH
 
#include "defineU2D.h"
 
namespace olb {
 
namespace legacy {
 
template<typename T, typename DESCRIPTOR>
void defineUBouzidi(SuperLattice<T, DESCRIPTOR>& sLattice, SuperGeometry<T,2>& superGeometry, int material,
                    AnalyticalF2D<T,T>& u,
                    std::vector<int> bulkMaterials)
{
  defineUBouzidi<T,DESCRIPTOR>(sLattice, superGeometry.getMaterialIndicator(material),
                               superGeometry.getMaterialIndicator(std::move(bulkMaterials)),
                               u);
}
 
template<typename T, typename DESCRIPTOR>
void defineUBouzidi(SuperLattice<T, DESCRIPTOR>& sLattice, FunctorPtr<SuperIndicatorF2D<T>>&& indicator,
                    FunctorPtr<SuperIndicatorF2D<T>>&& bulkIndicator,
                    AnalyticalF2D<T,T>& u)
{
 
  for (int iCloc = 0; iCloc < sLattice.getLoadBalancer().size(); ++iCloc) {
    defineUBouzidi<T,DESCRIPTOR>(sLattice.getBlockIndicator(iCloc), indicator->getBlockIndicatorF(iCloc),
                                 bulkIndicator->getBlockIndicatorF(iCloc),
                                 u);
  }
}
 
 
template<typename T, typename DESCRIPTOR>
void defineUBouzidi(BlockLattice<T,DESCRIPTOR>& block, BlockIndicatorF2D<T>& indicator, BlockIndicatorF2D<T>& bulkIndicator, AnalyticalF2D<T,T>& u)
{
  block.forSpatialLocations([&](auto iX, auto iY) {
    if (indicator(iX,iY)) {
      for (int iPop = 1; iPop < DESCRIPTOR::q ; ++iPop) {
        int iXn = iX + descriptors::c<DESCRIPTOR >(iPop,0);
        int iYn = iY + descriptors::c<DESCRIPTOR >(iPop,1);
        if (block.isInside({iXn,iYn}) && bulkIndicator(iXn, iYn)) {
          T intersection[] = { T(), T() }; // coord. of intersection
          int opp = descriptors::opposite<DESCRIPTOR >(iPop);
          if (getBoundaryIntersection<T,DESCRIPTOR>(block, iX, iY, opp, intersection) ) {
            T vel[]= {T(),T()};
            u(vel,intersection);
            defineUBouzidi<T,DESCRIPTOR>(iX, iY, opp, vel);
          }
        }
      }
    }
  });
}
 
 
template<typename T, typename DESCRIPTOR>
void defineUBouzidi(BlockLattice<T,DESCRIPTOR>& block, int iX, int iY, int iPop, const T u[DESCRIPTOR::d])
{
  bool _output = false;
  OstreamManager clout(std::cout, "defineUBouzidi");
  block.getDynamics(iX, iY)->defineU(iPop, u);
  if (_output) {
    clout << "defineUBouzidi(" << iX << ", " << iY << " )" << std::endl;
  }
}
 
 
 
template<typename T, typename DESCRIPTOR>
bool getBoundaryIntersection(BlockLattice<T,DESCRIPTOR>& block, int iX, int iY, int iPop, T point[DESCRIPTOR::d])
{
  return static_cast<legacy::OffDynamics<T,DESCRIPTOR>*>(block.getDynamics(iX,iY))->getBoundaryIntersection(iPop, point);
}
 
template<typename T, typename DESCRIPTOR>
void setBoundaryIntersection(BlockLattice<T,DESCRIPTOR>& block, int iX, int iY, int iPop, T distance)
{
  bool _output = false;
  OstreamManager clout(std::cout, "setBoundaryIntersection");
  static_cast<legacy::OffDynamics<T,DESCRIPTOR>*>(block.getDynamics(iX,iY))->setBoundaryIntersection(iPop, distance);
  if (_output) {
    clout << "setBoundaryIntersection(" << iX << ", " << iY << " )" << std::endl;
  }
}
 
}
 
}//namespace olb
 
#endif