Initialize at openlb-1-3

1 files changed, 198 insertions, 0 deletions

diff --git a/src/geometry/superGeometry2D.h b/src/geometry/superGeometry2D.h
new file mode 100644
index 0000000..6ddb49c
--- /dev/null
+++ b/src/geometry/superGeometry2D.h
@@ -0,0 +1,198 @@
+/*  This file is part of the OpenLB library
+ *
+ *  Copyright (C) 2013, 2014 Mathias J. Krause, Peter Weisbrod
+ *  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.
+ */
+
+/** \file
+ * Representation of a parallel 2D geometry -- header file.
+ */
+
+/// A super geometry represents a parrallel voxel mesh
+/** A super geometry consits of a number of block geometries,
+ * where the material numbers are stored. It is constructed
+ * from a cuboid geometry. All coboids of the cuboid geometry
+ * are asigned to block geometries which are extended by an
+ * overlap in order to enable efficient parallelisation.
+ *
+ * By the class access is provied to the material numbers of
+ * the mesh. Methods for renaming materials are provided as
+ * well as a statistic class.
+ *
+ * This class is not intended to be derived from.
+ */
+
+
+#ifndef SUPER_GEOMETRY_2D_H
+#define SUPER_GEOMETRY_2D_H
+
+#include <vector>
+#include <iostream>
+#include <string>
+
+#include "geometry/cuboidGeometry2D.h"
+#include "geometry/superGeometryStatistics2D.h"
+#include "geometry/blockGeometry2D.h"
+#include "geometry/blockGeometryView2D.h"
+#include "communication/superStructure2D.h"
+#include "communication/loadBalancer.h"
+#include "functors/analytical/indicator/indicatorF2D.h"
+#include "io/ostreamManager.h"
+
+
+/// All OpenLB code is contained in this namespace.
+namespace olb {
+
+template<typename T> class CuboidGeometry2D;
+template<typename T> class BlockGeometry2D;
+template<typename T> class BlockGeometryView2D;
+template<typename T> class IndicatorF2D;
+template<typename T> class SuperIndicatorF2D;
+template<typename T> class SuperStructure2D;
+template<typename T> class SuperGeometryStatistics2D;
+
+template<typename T>
+class SuperGeometry2D : public SuperStructure2D<T> {
+
+private:
+  /// Vector of block geometries with overlap
+  std::vector<BlockGeometry2D<T> > _extendedBlockGeometries;
+  /// Vector of block geometries without overlap
+  std::vector<BlockGeometryView2D<T> > _blockGeometries;
+  /// Statistic class
+  SuperGeometryStatistics2D<T> _statistics;
+  /// class specific output stream
+  mutable OstreamManager clout;
+
+public:
+  /// Constructor
+  SuperGeometry2D(CuboidGeometry2D<T>& cuboidGeometry,
+                  LoadBalancer<T>& lb, int overlap = 2);
+  /// Copy constructor
+  SuperGeometry2D(SuperGeometry2D const& rhs);
+  /// Copy assignment
+  SuperGeometry2D& operator=(SuperGeometry2D const& rhs);
+
+  /// Interface for the communicator class: Write access to the memory of the data of the super structure
+  bool* operator() (int iCloc, int iX, int iY, int iData) override;
+  /// Interface for the communicator class: Read only access to the dim of the data of the super structure
+  int getDataSize() const override;
+  /// Interface for the communicator class: Read only access to the data type dim of the data of the super structure
+  int getDataTypeSize() const override;
+
+  /// Write access to the material numbers, error handling: stops the program if data is not available
+  int& set(int iCglob, int iXloc, int iYloc); //TODO to be removed set->get, problem: with get calling wrong function
+  /// Read only access to the material numbers, error handling: returns 0 if data is not available
+  int const&  get(int iCglob, int iXloc, int iYloc) const;
+  /// Read only access to the material numbers with global communication to all ranks
+  int getAndCommunicate(int iCglob, int iXloc, int iYloc) const;
+  /// Write access to the material numbers, error handling: stops the program if data is not available
+  int& set(std::vector<int> latticeR); //TODO to be removed set->get, problem: with get calling wrong function
+  /// Read only access to the material numbers, error handling: returns 0 if data is not available
+  int const& get(std::vector<int> latticeR) const;
+  /// Read only access to the material numbers with global communication to all ranks
+  int getAndCommunicate(std::vector<int> latticeR) const;
+
+  /// Transforms a lattice to physical position (SI unites)
+  std::vector<T> getPhysR(int iCglob, int iX, int iY) const;
+  /// Transforms a lattice to physical position (SI unites)
+  std::vector<T> getPhysR(std::vector<int> latticeR) const;
+  /// Transforms a lattice to physical position (SI unites)
+  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;
+
+  /// Read and write access to a single extended block geometry
+  BlockGeometry2D<T>& getExtendedBlockGeometry(int locIC);
+  /// Read only access to a single extended block geometry
+  BlockGeometry2D<T> const& getExtendedBlockGeometry(int locIC) const;
+  /// Read and write access to a single block geometry
+  BlockGeometryView2D<T>& getBlockGeometry(int locIC);
+  /// Read only access to a single block geometry
+  BlockGeometryView2D<T> const& getBlockGeometry(int locIC) const;
+
+  /// Returns the statistics object
+  SuperGeometryStatistics2D<T>& getStatistics();
+  /// Read and write access to the statistic status flag, update needed = true
+  bool& getStatisticsStatus();
+  /// Read only access to the statistic status flag, update needed = true
+  bool const& getStatisticsStatus() const;
+  /// Updates the super geometry at the boundaries if needed and afterwards the statisics if needed
+  void updateStatistics(bool verbose=true);
+
+  /// Executes an outer cleaning
+  int clean(bool verbose=true);
+  /// Removes not needed material fluids from the outer domain
+  int outerClean(bool verbose=true);
+  /// inner cleaning for all boundary types
+  int innerClean(bool verbose=true);
+  /// inner cleaning for specific boundary types
+  int innerClean(int material, bool verbose=true);
+  /// check for errors (searches for all outer voxels (=0) with an inner voxel (=1) as a direct neighbour)
+  bool checkForErrors(bool verbose=true);
+
+  /// reset all cell materials inside of a domain to 0
+  void reset(IndicatorF2D<T>& domain);
+
+  /// replace one material with another
+  void rename(int fromM, int toM);
+  /// replace one material that fulfills an indicator functor condition with another
+  void rename(int fromM, int toM, IndicatorF2D<T>& condition);
+  /// replace one material with another respecting an offset (overlap)
+  void rename(int fromM, int toM, unsigned offsetX, unsigned offsetY);
+  /// renames all voxels of material fromM to toM if the number of voxels given by testDirection is of material testM
+  void rename(int fromM, int toM, int testM, std::vector<int> testDirection);
+  /// renames all boundary voxels of material fromBcMat to toBcMat if two neighbour voxel in the direction of the discrete normal are fluid voxel with material fluidM in the region where the indicator function is fulfilled
+  void rename(int fromBcMat, int toBcMat, int fluidMat, IndicatorF2D<T>& condition);
+
+  /// Prints some information about the super geometry
+  void print();
+
+  template<template<typename> class Indicator, typename... Args>
+  std::unique_ptr<SuperIndicatorF2D<T>> getIndicator(Args&&... args)
+  {
+    static_assert(std::is_base_of<SuperIndicatorF2D<T>, Indicator<T>>::value,
+                  "Indicator to be constructed is SuperIndicatorF2D implementation");
+
+    return std::unique_ptr<SuperIndicatorF2D<T>>(
+             new Indicator<T>(*this, std::forward<Args>(args)...)
+           );
+  }
+
+  /**
+   * Returns a material indicator using the given vector of materials
+   *
+   * \param  materials Materials to be indicated
+   * \returns          Unique ownership of the constructed indicator.
+   *                   May be stored or passed directly to e.g. defineDynamics
+   **/
+  std::unique_ptr<SuperIndicatorF2D<T>> getMaterialIndicator(std::vector<int>&& materials);
+  /**
+   * Returns a material indicator using a single material number
+   *
+   * \param material Material to be indicated
+   * \returns        Unique ownership of the constructed indicator.
+   *                 May be stored or passed directly to e.g. defineDynamics
+   **/
+  std::unique_ptr<SuperIndicatorF2D<T>> getMaterialIndicator(int material);
+};
+
+} // namespace olb
+
+#endif