integral_operator.hh
4ecfb2213b97master

integral_operator.hh
View Options

			/**
			* @file
			*
			* @author Lucas Frérot <lucas.frerot@epfl.ch>
			*
			* @section LICENSE
			*
			* Copyright (©) 2017 EPFL (Ecole Polytechnique Fédérale de
			* Lausanne) Laboratory (LSMS - Laboratoire de Simulation en Mécanique des
			* Solides)
			*
			* Tamaas is free software: you can redistribute it and/or modify it under the
			* terms of the GNU Lesser General Public License as published by the Free
			* Software Foundation, either version 3 of the License, or (at your option) any
			* later version.
			*
			* Tamaas 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 Lesser General Public License for more
			* details.
			*
			* You should have received a copy of the GNU Lesser General Public License
			* along with Tamaas. If not, see <http://www.gnu.org/licenses/>.
			*
			*/
			/* -------------------------------------------------------------------------- */
			#ifndef INTEGRAL_OPERATOR_HH
			#define INTEGRAL_OPERATOR_HH
			/* -------------------------------------------------------------------------- */
			#include "grid_base.hh"
			#include "model_type.hh"
			#include <boost/preprocessor/cat.hpp>
			#include <boost/preprocessor/seq.hpp>
			#include <boost/preprocessor/stringize.hpp>
			/* -------------------------------------------------------------------------- */
			namespace tamaas {
			/* -------------------------------------------------------------------------- */

			class Model;

			/**
			* @brief Generic class for integral operators
			*/
			class IntegralOperator {
			public:
			/// Kind of operator
			enum kind { neumann, dirichlet };

			public:
			/// Constructor
			IntegralOperator(Model* model) : model(model) {}
			/// Virtual destructor for subclasses
			virtual ~IntegralOperator() = default;

			/// Apply operator on input
			virtual void apply(GridBase<Real>& input, GridBase<Real>& output) const = 0;

			/// Apply operator on filtered input
			virtual void applyIf(GridBase<Real>& input, GridBase<Real>& output,
			std::function<bool(UInt)>) const {
			apply(input, output);
			}

			/// Update any data dependent on model parameters
			virtual void updateFromModel() = 0;

			/// Get model
			const Model& getModel() const { return *model; }

			/// Kind
			virtual kind getKind() const = 0;

			/// Type
			virtual model_type getType() const = 0;

			/// Norm
			virtual Real getOperatorNorm() {
			TAMAAS_EXCEPTION("operator does not implement norm");
			}

			protected:
			Model* model = nullptr;
			};

			/* -------------------------------------------------------------------------- */
			/* Printing IntegralOperator::kind to string */
			/* -------------------------------------------------------------------------- */
			#define INTEGRAL_KINDS (neumann)(dirichlet)

			inline std::ostream& operator<<(std::ostream& o,
			const IntegralOperator::kind& val) {
			switch (val) {
			#define PRINT_INTEGRAL_KIND(r, data, kind) \
			case data::kind: \
			o << BOOST_PP_STRINGIZE(kind); \
			break;
			BOOST_PP_SEQ_FOR_EACH(PRINT_INTEGRAL_KIND, IntegralOperator,
			INTEGRAL_KINDS);
			#undef PRINT_INTEGRAL_KIND
			}
			return o;
			}

			#undef INTEGRAL_KINDS

			/* -------------------------------------------------------------------------- */
			} // namespace tamaas

			#endif // INTEGRAL_OPERATOR_HH

tamaas-public/src/model/integral_operator.hh4ecfb2213b97master

integral_operator.hhView Options

tamaas-public/src/model/integral_operator.hh
4ecfb2213b97master

integral_operator.hh
View Options