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
---
src/functors/analytical/analyticCalcF.hh | 544 +++++++++++++++++++++++++++++++
1 file changed, 544 insertions(+)
create mode 100644 src/functors/analytical/analyticCalcF.hh
(limited to 'src/functors/analytical/analyticCalcF.hh')
diff --git a/src/functors/analytical/analyticCalcF.hh b/src/functors/analytical/analyticCalcF.hh
new file mode 100644
index 0000000..bca8f6f
--- /dev/null
+++ b/src/functors/analytical/analyticCalcF.hh
@@ -0,0 +1,544 @@
+/* This file is part of the OpenLB library
+ *
+ * Copyright (C) 2012-2018 Lukas Baron, Tim Dornieden, Mathias J. Krause,
+ * Albert Mink, Adrian Kummerlaender
+ * E-mail contact: info@openlb.net
+ * The most recent release of OpenLB can be downloaded at
+ *
+ *
+ * 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 ANALYTICAL_CALC_F_HH
+#define ANALYTICAL_CALC_F_HH
+
+
+#include "analyticCalcF.h"
+#include "analyticalF.h"
+#include "core/olbDebug.h"
+
+namespace olb {
+
+
+
+//////////////////////////////// AnalyticCalcF1D ////////////////////////////////
+template class F>
+AnalyticCalcF1D::AnalyticCalcF1D(FunctorPtr>&& f,
+ FunctorPtr>&& g)
+ : AnalyticalF1D(f->getTargetDim()),
+ _f(std::move(f)),
+ _g(std::move(g))
+{
+ OLB_ASSERT(g->getTargetDim() == f->getTargetDim(),
+ "the dimensions of both functors need to be equal");
+ std::swap(f->_ptrCalcC, this->_ptrCalcC);
+ this->getName() = "(" + _f->getName() + F::symbol + _g->getName() + ")";
+}
+
+template class F>
+AnalyticCalcF1D::AnalyticCalcF1D(T scalar, FunctorPtr>&& g)
+ : AnalyticCalcF1D(
+ std::unique_ptr>(
+ new AnalyticalConst1D(std::vector(g->getTargetDim(), scalar))),
+ std::forward(g))
+{ }
+
+template class F>
+AnalyticCalcF1D::AnalyticCalcF1D(FunctorPtr>&& f, T scalar)
+ : AnalyticCalcF1D(
+ std::forward(f),
+ std::unique_ptr>(
+ new AnalyticalConst1D(std::vector(f->getTargetDim(), scalar))))
+{ }
+
+template class F>
+bool AnalyticCalcF1D::operator()(T output[], const S input[])
+{
+ T outputTmp[this->_g->getTargetDim()];
+ this->_g(outputTmp, input);
+ this->_f(output, input);
+ for (int i = 0; i < this->_f->getTargetDim(); ++i) {
+ output[i] = F()(output[i], outputTmp[i]);
+ }
+ return true;
+}
+
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus1D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus1D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator+(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus1D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator-(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus1D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator-(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus1D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator-(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus1D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator*(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication1D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator*(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication1D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator*(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication1D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator/(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision1D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator/(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision1D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator/(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision1D(lhs, std::move(rhs)));
+}
+
+/////////////////////////////////operator()/// ////////////////////////////////
+template
+AnalyticalF1D& AnalyticalF1D::operator+(AnalyticalF1D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcPlus1D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+template
+AnalyticalF1D& AnalyticalF1D::operator-(AnalyticalF1D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcMinus1D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+template
+AnalyticalF1D& AnalyticalF1D::operator*(AnalyticalF1D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcMultiplication1D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+template
+AnalyticalF1D& AnalyticalF1D::operator/(AnalyticalF1D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcDivision1D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+
+//////////////////////////////// AnalyticCalcF2D ////////////////////////////////
+template class F>
+AnalyticCalcF2D::AnalyticCalcF2D(FunctorPtr>&& f,
+ FunctorPtr>&& g)
+ : AnalyticalF2D(f->getTargetDim()),
+ _f(std::move(f)),
+ _g(std::move(g))
+{
+ OLB_ASSERT(g->getTargetDim() == f->getTargetDim(),
+ "the dimensions of both functors need to be equal");
+ // pass through the shared_ptr from the first argument f to the arithmetic class itself.
+ // used by secsessive calls: e.g. (functorA + functor B) followed by (functorA + functorC)
+ // the result of the first operation is overwritten by the second.
+
+ // equivalent operations
+ // std::swap(f._ptrCalcC, this->_ptrCalcC);
+ // this->_ptrCalcC = f._ptrCalcC;
+ this->_ptrCalcC.swap(f->_ptrCalcC);
+
+ this->getName() = "(" + _f->getName() + F::symbol + _g->getName() + ")";
+}
+
+template class F>
+AnalyticCalcF2D::AnalyticCalcF2D(T scalar, FunctorPtr>&& g)
+ : AnalyticCalcF2D(
+ std::unique_ptr>(
+ new AnalyticalConst2D(std::vector(g->getTargetDim(), scalar))),
+ std::forward(g))
+{ }
+
+template class F>
+AnalyticCalcF2D::AnalyticCalcF2D(FunctorPtr>&& f, T scalar)
+ : AnalyticCalcF2D(
+ std::forward(f),
+ std::unique_ptr>(
+ new AnalyticalConst2D(std::vector(f->getTargetDim(), scalar))))
+{ }
+
+template class F>
+bool AnalyticCalcF2D::operator()(T output[], const S input[])
+{
+ T outputTmp[this->_g->getTargetDim()];
+ this->_g(outputTmp, input);
+ this->_f(output, input);
+ for (int i = 0; i < this->_f->getTargetDim(); ++i) {
+ output[i] = F()(output[i], outputTmp[i]);
+ }
+ return true;
+}
+
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus2D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus2D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator+(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus2D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator-(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus2D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator-(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus2D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator-(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMinus2D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator*(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication2D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator*(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication2D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator*(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcMultiplication2D(lhs, std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator/(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision2D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator/(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision2D(std::move(lhs), rhs));
+}
+
+template
+std::shared_ptr> operator/(T lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcDivision2D(lhs, std::move(rhs)));
+}
+
+/////////////////////////////////operator()////////////////////////////////////
+template
+AnalyticalF2D& AnalyticalF2D::operator+(AnalyticalF2D& rhs)
+{
+ // version 1
+ // AnalyticalF2D* tmp = new AnalyticCalcPlus2D(*this,rhs);
+ // std::shared_ptr< GenericF > ptr( tmp );
+ // this->_ptrCalcC = ptr;
+
+ // version 2
+ // std::shared_ptr< AnalyticalF2D > tmp = std::make_shared< AnalyticCalcPlus2D >(*this,rhs);
+
+ // version 2.5
+ // std::shared_ptr< AnalyticCalcPlus2D > tmp( new AnalyticCalcPlus2D(*this,rhs) );
+
+ // version 3
+ auto tmp = std::make_shared< AnalyticCalcPlus2D >(*this,rhs);
+
+ this->_ptrCalcC = tmp;
+
+ // std::cout << "operator+(): " << this->_ptrCalcC.get()->getName() << std::endl;
+ return *tmp;
+}
+
+template
+AnalyticalF2D& AnalyticalF2D::operator-(AnalyticalF2D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcMinus2D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+template
+AnalyticalF2D& AnalyticalF2D::operator*(AnalyticalF2D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcMultiplication2D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+template
+AnalyticalF2D& AnalyticalF2D::operator/(AnalyticalF2D& rhs)
+{
+ auto tmp = std::make_shared< AnalyticCalcDivision2D >(*this,rhs);
+ this->_ptrCalcC = tmp;
+ return *tmp;
+}
+
+
+//////////////////////////////// AnalyticCalcF3D ////////////////////////////////
+template class F>
+AnalyticCalcF3D::AnalyticCalcF3D(FunctorPtr>&& f,
+ FunctorPtr>&& g):
+ AnalyticalF3D(f->getTargetDim()),
+ _f(std::move(f)),
+ _g(std::move(g))
+{
+ OLB_ASSERT(g->getTargetDim() == f->getTargetDim(),
+ "the dimensions of both functors need to be equal");
+ std::swap(f->_ptrCalcC, this->_ptrCalcC);
+ this->getName() = "(" + _f->getName() + F::symbol + _g->getName() + ")";
+}
+
+template class F>
+AnalyticCalcF3D::AnalyticCalcF3D(T scalar, FunctorPtr>&& g)
+ : AnalyticCalcF3D(
+ std::unique_ptr>(
+ new AnalyticalConst3D(std::vector(g->getTargetDim(), scalar))),
+ std::forward(g))
+{ }
+
+template class F>
+AnalyticCalcF3D::AnalyticCalcF3D(FunctorPtr>&& f, T scalar)
+ : AnalyticCalcF3D(
+ std::forward(f),
+ std::unique_ptr>(
+ new AnalyticalConst3D(std::vector(f->getTargetDim(), scalar))))
+{ }
+
+template class F>
+bool AnalyticCalcF3D::operator()(T output[], const S input[])
+{
+ T outputTmp[this->_g->getTargetDim()];
+ this->_g(outputTmp, input);
+ this->_f(output, input);
+ for (int i = 0; i < this->_f->getTargetDim(); ++i) {
+ output[i] = F()(output[i], outputTmp[i]);
+ }
+ return true;
+}
+
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, std::shared_ptr> rhs)
+{
+ return std::shared_ptr>(
+ new AnalyticCalcPlus3D(std::move(lhs), std::move(rhs)));
+}
+
+template
+std::shared_ptr> operator+(std::shared_ptr> lhs, T rhs)
+{
+ return std::shared_ptr