cuboidGeometry2D.h

Go to the documentation of this file.

/*  This file is part of the OpenLB library
 *
 *  Copyright (C) 2007, 2014 Mathias J. Krause
 *  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 CUBOID_GEOMETRY_2D_H
#define CUBOID_GEOMETRY_2D_H
 
#include <vector>
#include "geometry/cuboid2D.h"
#include "io/ostreamManager.h"
 
 
// All OpenLB code is contained in this namespace.
namespace olb {
 
template<typename T> class IndicatorF2D;
 
 
template<typename T>
class CuboidGeometry2D {
 
private:
  // TODO quick and dirty solution for 3d split to be removed after 2d clearence
  bool _oldApproach;
 
  std::vector<Cuboid2D<T> > _cuboids;
  Cuboid2D<T> _motherCuboid;
  std::vector<bool> _periodicityOn;
  mutable OstreamManager clout;
 
public:
 
  CuboidGeometry2D();
  CuboidGeometry2D(T originX, T originY, T deltaR, int nX, int nY, int nC=1);
  CuboidGeometry2D(IndicatorF2D<T>& indicatorF, T voxelSize, int nC=1);
 
  void reInit(T globPosX, T globPosY,
              T delta, int nX, int nY, int nC=1);
  Cuboid2D<T>& get(int i);
  Cuboid2D<T> const& get(int i) const;
  void setPeriodicity(bool periodicityX, bool periodicityY);
 
  bool getC(std::vector<T> physR, int& iC) const;
  bool getC(const Vector<T,2>& physR, int& iC) const;
  bool getLatticeR(std::vector<T> physR, std::vector<int>& latticeR) const;
  bool getLatticeR(int latticeR[], const T physR[]) const;
  bool getLatticeR(const Vector<T,2>& physR, Vector<int,3>& latticeR) const;
  bool getFloorLatticeR(std::vector<T> physR, std::vector<int>& latticeR) const;
  bool getFloorLatticeR(const Vector<T,2>& physR, Vector<int,3>& latticeR) const;
  std::vector<T> getPhysR(int iCglob, int iX, int iY) const;
  std::vector<T> getPhysR(std::vector<int> latticeR) const;
  void getPhysR(T output[2], const int latticeR[3]) const;
  void getPhysR(T output[2], const int iCglob, const int iX, const int iY) const;
  void getPhysR(T output[2], const int iCglob, LatticeR<2> latticeR) const;
  void getPhysR(T output[3], LatticeR<3> latticeR) const;
 
  int getNc() const;
  T getMinRatio() const;
  T getMaxRatio() const;
  std::vector<T> getMinPhysR() const;
  std::vector<T> getMaxPhysR() const;
  T getMinPhysVolume() const;
  T getMaxPhysVolume() const;
  size_t getMinLatticeVolume() const;
  size_t getMaxLatticeVolume() const;
  T getMinDeltaR() const;
  T getMaxDeltaR() const;
  Cuboid2D<T> getMotherCuboid() const;
 
  int get_iC(T globX, T globY, int offset = 0) const;
  int get_iC(T globX, T globY, int orientationX, int orientationY) const;
 
  void add(Cuboid2D<T> cuboid);
  void remove(int iC);
  //void remove(olb::ScalarField2D<int>* geometryData);
  void split(int iC, int p);
  void shrink(IndicatorF2D<T>& indicatorF);
 
  void getNeighbourhood(int cuboid, std::vector<int>& neighbours, int offset = 0);
 
  void print() const;
  void printExtended();
};
 
 
template<typename T>
void evaluateCuboidGeometryNeighbourhood(
  CuboidGeometry2D<T>& cuboidGeometry,
  std::map<int,std::vector<int>>& neighbourhood,
  int offset=0 )
{
  for (int iC=0; iC<cuboidGeometry.getNc(); ++iC){
    //Retrieve neighbours of iC
    std::vector<int> neighbours;
    cuboidGeometry.getNeighbourhood(iC, neighbours, offset);
    neighbourhood[iC]=neighbours;
  }
}
 
}  // namespace olb
 
#endif