From 94d3e79a8617f88dc0219cfdeedfa3147833719d Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Mon, 24 Jun 2019 14:43:36 +0200
Subject: Initialize at openlb-1-3

---
 .../lattice/indicator/superIndicatorF2D.hh         | 164 +++++++++++++++++++++
 1 file changed, 164 insertions(+)
 create mode 100644 src/functors/lattice/indicator/superIndicatorF2D.hh

(limited to 'src/functors/lattice/indicator/superIndicatorF2D.hh')

diff --git a/src/functors/lattice/indicator/superIndicatorF2D.hh b/src/functors/lattice/indicator/superIndicatorF2D.hh
new file mode 100644
index 0000000..d83c745
--- /dev/null
+++ b/src/functors/lattice/indicator/superIndicatorF2D.hh
@@ -0,0 +1,164 @@
+/*  This file is part of the OpenLB library
+ *
+ *  Copyright (C) 2017 Adrian Kummerlaender
+ *  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_F_2D_HH
+#define SUPER_INDICATOR_F_2D_HH
+
+#include <numeric>
+#include <algorithm>
+
+#include "core/util.h"
+#include "superIndicatorF2D.h"
+#include "blockIndicatorF2D.h"
+
+namespace olb {
+
+
+template <typename T>
+SuperIndicatorFfromIndicatorF2D<T>::SuperIndicatorFfromIndicatorF2D(
+  IndicatorF2D<T>& indicatorF, SuperGeometry2D<T>& geometry)
+  : SuperIndicatorF2D<T>(geometry),
+    _indicatorF(indicatorF)
+{
+  this->getName() = "SuperIndicator_from_" + _indicatorF.getName();
+
+  for (int iC = 0; iC < this->getSuperStructure().getLoadBalancer().size(); ++iC) {
+    this->_blockF.emplace_back(
+      new BlockIndicatorFfromIndicatorF2D<T>(
+        _indicatorF, geometry.getBlockGeometry(iC))
+    );
+    this->_extendedBlockF.emplace_back(
+      new BlockIndicatorFfromIndicatorF2D<T>(
+        _indicatorF, geometry.getExtendedBlockGeometry(iC))
+    );
+  }
+}
+
+template <typename T>
+bool SuperIndicatorFfromIndicatorF2D<T>::operator() (bool output[], const int input[])
+{
+  T physR[2];
+  this->_superStructure.getCuboidGeometry().getPhysR(physR, input);
+  return _indicatorF(output, physR);
+}
+
+
+template <typename T, bool HLBM>
+SuperIndicatorFfromSmoothIndicatorF2D<T,HLBM>::SuperIndicatorFfromSmoothIndicatorF2D(
+  FunctorPtr<SmoothIndicatorF2D<T,T,HLBM>>&& indicatorF,
+  SuperGeometry2D<T>&                     geometry)
+  : SuperIndicatorF2D<T>(geometry),
+    _indicatorF(std::move(indicatorF))
+{
+  this->getName() = "SuperIndicator_from_" + _indicatorF->getName();
+
+  LoadBalancer<T>& load = this->getSuperStructure().getLoadBalancer();
+
+  for (int iC = 0; iC < load.size(); ++iC) {
+    this->_blockF.emplace_back(
+      new BlockIndicatorFfromSmoothIndicatorF2D<T,HLBM>(
+        *_indicatorF, geometry.getBlockGeometry(iC))
+    );
+    this->_extendedBlockF.emplace_back(
+      new BlockIndicatorFfromSmoothIndicatorF2D<T,HLBM>(
+        *_indicatorF, geometry.getExtendedBlockGeometry(iC))
+    );
+  }
+}
+
+template <typename T, bool HLBM>
+bool SuperIndicatorFfromSmoothIndicatorF2D<T,HLBM>::operator() (bool output[], const int input[])
+{
+  T physR[2];
+  T inside[1];
+  this->_superStructure.getCuboidGeometry().getPhysR(physR, input);
+  _indicatorF(inside, physR);
+  return !util::nearZero(inside[0]);
+}
+
+
+template <typename T>
+SuperIndicatorMaterial2D<T>::SuperIndicatorMaterial2D(
+  SuperGeometry2D<T>& geometry, std::vector<int> materials)
+  : SuperIndicatorF2D<T>(geometry)
+{
+  const std::string matString = std::accumulate(
+                                  materials.begin()+1,
+                                  materials.end(),
+                                  std::to_string(materials[0]),
+  [](const std::string& a, int b) {
+    return a + '_' + std::to_string(b);
+  });
+  this->getName() = "SuperIndicator_on_Material_" + matString;
+
+  for (int iC = 0; iC < this->_superGeometry.getLoadBalancer().size(); ++iC) {
+    this->_blockF.emplace_back(
+      new BlockIndicatorMaterial2D<T>(this->_superGeometry.getBlockGeometry(iC),
+                                      materials)
+    );
+    this->_extendedBlockF.emplace_back(
+      new BlockIndicatorMaterial2D<T>(this->_superGeometry.getExtendedBlockGeometry(iC),
+                                      materials)
+    );
+  }
+}
+
+template <typename T>
+bool SuperIndicatorMaterial2D<T>::operator() (bool output[], const int input[])
+{
+  output[0] = false;
+
+  LoadBalancer<T>& load = this->_superGeometry.getLoadBalancer();
+
+  if (!this->_blockF.empty() && load.isLocal(input[0])) {
+    // query material number of appropriate block indicator
+    return this->getBlockF(load.loc(input[0]))(output, &input[1]);
+  } else {
+    return false;
+  }
+}
+
+template <typename T>
+SuperIndicatorIdentity2D<T>::SuperIndicatorIdentity2D(FunctorPtr<SuperIndicatorF2D<T>>&& indicatorF)
+  : SuperIndicatorF2D<T>(indicatorF->getSuperGeometry()),
+    _indicatorF(std::move(indicatorF))
+{
+  this->getName() = _indicatorF->getName();
+
+  for (int iC = 0; iC < _indicatorF->getBlockFSize(); ++iC) {
+    this->_blockF.emplace_back(
+      new BlockIndicatorIdentity2D<T>(_indicatorF->getBlockIndicatorF(iC)));
+    this->_extendedBlockF.emplace_back(
+      new BlockIndicatorIdentity2D<T>(_indicatorF->getExtendedBlockIndicatorF(iC)));
+  }
+}
+
+template <typename T>
+bool SuperIndicatorIdentity2D<T>::operator()(bool output[], const int input[])
+{
+  return _indicatorF(output, input);
+}
+
+} // namespace olb
+
+#endif
-- 
cgit v1.2.3