Initialize at openlb-1-3

1 files changed, 102 insertions, 0 deletions

diff --git a/src/functors/lattice/indicator/superIndicatorBaseF3D.hh b/src/functors/lattice/indicator/superIndicatorBaseF3D.hh
new file mode 100644
index 0000000..9a0bb63
--- /dev/null
+++ b/src/functors/lattice/indicator/superIndicatorBaseF3D.hh
@@ -0,0 +1,102 @@
+/*  This file is part of the OpenLB library
+ *
+ *  Copyright (C) 2016 Benjamin Förster
+ *  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 SUPER_INDICATOR_BASE_F_3D_HH
+#define SUPER_INDICATOR_BASE_F_3D_HH
+
+#include "superIndicatorBaseF3D.h"
+#include "geometry/superGeometry3D.h"
+
+namespace olb {
+
+
+template <typename T>
+SuperIndicatorF3D<T>::SuperIndicatorF3D(SuperGeometry3D<T>& geometry)
+  : SuperF3D<T, bool>(geometry, 1),
+    _superGeometry(geometry)
+{ }
+
+template <typename T>
+BlockIndicatorF3D<T>& SuperIndicatorF3D<T>::getBlockIndicatorF(int iCloc)
+{
+  OLB_ASSERT(iCloc < int(this->_blockF.size()) && iCloc >= 0,
+             "block functor index within bounds");
+  // Note: The type system doesn't guarantee this operation to be valid
+  //       as blockF may contain any implementation of the BlockF3D interface.
+  return *static_cast<BlockIndicatorF3D<T>*>(this->_blockF[iCloc].get());
+}
+
+template <typename T>
+BlockIndicatorF3D<T>& SuperIndicatorF3D<T>::getExtendedBlockIndicatorF(int iCloc)
+{
+  OLB_ASSERT(iCloc < int(this->_extendedBlockF.size()) && iCloc >= 0,
+             "block functor index within bounds");
+  return *(this->_extendedBlockF[iCloc]);
+}
+
+template <typename T>
+SuperGeometry3D<T>& SuperIndicatorF3D<T>::getSuperGeometry()
+{
+  return this->_superGeometry;
+}
+
+template <typename T>
+bool SuperIndicatorF3D<T>::operator() (const int input[])
+{
+  bool output{};
+  if (_cachedData) {
+    output = _cachedData->get(input[0], input[1], input[2], input[3]);
+  }
+  else {
+    this->operator()(&output, input);
+  }
+  return output;
+}
+
+template <typename T>
+bool SuperIndicatorF3D<T>::operator() (int iC, int iX, int iY, int iZ)
+{
+  bool output{};
+  if (_cachedData) {
+    output = _cachedData->get(iC, iX, iY, iZ);
+  }
+  else {
+    this->operator()(&output, iC, iX, iY, iZ);
+  }
+  return output;
+}
+
+template <typename T>
+void SuperIndicatorF3D<T>::cache()
+{
+  _cachedData = std::unique_ptr<SuperData3D<T,bool>>(
+                  new SuperData3D<T,bool>(*this));
+  for (unsigned iC = 0; iC < this->_blockF.size(); ++iC) {
+    getExtendedBlockIndicatorF(iC).setCache(_cachedData->get(iC));
+  }
+}
+
+
+} // namespace olb
+
+#endif