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/lattice/integral/superLpNorm3D.h | 106 ++++++++++++++++++++++++++
 1 file changed, 106 insertions(+)
 create mode 100644 src/functors/lattice/integral/superLpNorm3D.h

(limited to 'src/functors/lattice/integral/superLpNorm3D.h')

diff --git a/src/functors/lattice/integral/superLpNorm3D.h b/src/functors/lattice/integral/superLpNorm3D.h
new file mode 100644
index 0000000..c526ec1
--- /dev/null
+++ b/src/functors/lattice/integral/superLpNorm3D.h
@@ -0,0 +1,106 @@
+/*  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_LP_NORM_3D_H
+#define SUPER_LP_NORM_3D_H
+
+#include <memory>
+#include <vector>
+
+#include "utilities/functorPtr.h"
+
+namespace olb {
+
+template<typename T, typename W> class SuperF3D;
+template<typename T> class SuperIndicatorF3D;
+template<typename T> class SuperGeometry3D;
+
+/// Functor that returns the Lp norm over omega of the the euklid norm of the input functor
+/**
+ * Maintains block level BlockLpNorm3D functors as required.
+ *
+ * P == 0: inf norm
+ * P >= 1: p norm
+ */
+template <typename T, typename W, int P>
+class SuperLpNorm3D : public SuperF3D<T,W> {
+private:
+  FunctorPtr<SuperF3D<T,W>>        _f;
+  FunctorPtr<SuperIndicatorF3D<T>> _indicatorF;
+public:
+  /**
+   * \param f          (non-)owning pointer or reference to SuperF3D<T,W>
+   * \param indicatorF (non-)owning pointer or reference to SuperIndicatorF3D<T>.
+   *                   Describes the subset to be integrated.
+   **/
+  SuperLpNorm3D(FunctorPtr<SuperF3D<T,W>>&&        f,
+                FunctorPtr<SuperIndicatorF3D<T>>&& indicatorF);
+
+  /**
+   * Legacy constructor accepting super geometry reference.
+   *
+   * \param f          (non-)owning pointer or reference to SuperF3D<T,W>
+   * \param indicatorF (non-)owning pointer or reference to SuperIndicatorF3D<T>.
+   *                   Describes the subset to be integrated.
+   **/
+  SuperLpNorm3D(FunctorPtr<SuperF3D<T,W>>&& f,
+                SuperGeometry3D<T>&,
+                FunctorPtr<SuperIndicatorF3D<T>>&& indicatorF);
+
+  /**
+   * \param f         (non-)owning pointer or reference to SuperF3D<T,W>
+   * \param geometry  super geometry required to construct SuperIndicatorMaterial3D using materials
+   * \param materials vector of material numbers to be included in the Lp norm
+   **/
+  SuperLpNorm3D(FunctorPtr<SuperF3D<T,W>>&& f,
+                SuperGeometry3D<T>& geometry,
+                std::vector<int>    materials);
+
+  /**
+   * \param f        (non-)owning pointer or reference to SuperF3D<T,W>
+   * \param geometry super geometry required to construct SuperIndicatorMaterial3D using material
+   * \param material single material number to be included in the Lp norm
+   **/
+  SuperLpNorm3D(FunctorPtr<SuperF3D<T,W>>&& f,
+                SuperGeometry3D<T>& geometry,
+                int                 material);
+
+  bool operator() (W output[], const int input[]) override;
+};
+
+
+/// Functor that returns the L1 norm over omega of the the euklid norm of the input functor
+template <typename T, typename W = T>
+using SuperL1Norm3D = SuperLpNorm3D<T,W,1>;
+
+/// Functor that returns the L2 norm over omega of the the euklid norm of the input functor
+template <typename T, typename W = T>
+using SuperL2Norm3D = SuperLpNorm3D<T,W,2>;
+
+/// Functor that returns the Linf norm over omega of the the euklid norm of the input functor
+template <typename T, typename W = T>
+using SuperLinfNorm3D = SuperLpNorm3D<T,W,0>;
+
+}
+
+#endif
-- 
cgit v1.2.3