Diffusion akantu (bugfixes/contact_mechanics_model) Browse

contact_mechan…

dumper_generic_elemental_field.hh
View Options

			/**
			* @file dumper_generic_elemental_field.hh
			*
			* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
			* @author Nicolas Richart <nicolas.richart@epfl.ch>
			*
			* @date creation: Tue Sep 02 2014
			* @date last modification: Tue Jan 19 2016
			*
			* @brief Generic interface for elemental fields
			*
			* @section LICENSE
			*
			* Copyright (©) 2014, 2015 EPFL (Ecole Polytechnique Fédérale de Lausanne)
			* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
			*
			* Akantu 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.
			*
			* Akantu 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 Akantu. If not, see <http://www.gnu.org/licenses/>.
			*
			*/

			#ifndef __AKANTU_DUMPER_GENERIC_ELEMENTAL_FIELD_HH__
			#define __AKANTU_DUMPER_GENERIC_ELEMENTAL_FIELD_HH__
			/* -------------------------------------------------------------------------- */
			#include "dumper_field.hh"
			#include "element_type_map_filter.hh"
			#include "dumper_element_iterator.hh"
			#include "dumper_homogenizing_field.hh"
			/* -------------------------------------------------------------------------- */
			namespace akantu {
			__BEGIN_AKANTU_DUMPER__
			/* -------------------------------------------------------------------------- */

			template <class _types, template <class> class iterator_type>
			class GenericElementalField : public Field {
			/* ------------------------------------------------------------------------ */
			/* Typedefs */
			/* ------------------------------------------------------------------------ */
			public:
			// check dumper_type_traits.hh for additional information over these types
			typedef _types types;
			typedef typename types::data_type data_type;
			typedef typename types::it_type it_type;
			typedef typename types::field_type field_type;
			typedef typename types::array_type array_type;
			typedef typename types::array_iterator array_iterator;
			typedef typename field_type::type_iterator field_type_iterator;
			typedef iterator_type<types> iterator;

			/* ------------------------------------------------------------------------ */
			/* Constructors/Destructors */
			/* ------------------------------------------------------------------------ */
			public:
			GenericElementalField(const field_type & field,
			UInt spatial_dimension = _all_dimensions,
			GhostType ghost_type = _not_ghost,
			ElementKind element_kind = _ek_not_defined)
			: field(field), spatial_dimension(spatial_dimension),
			ghost_type(ghost_type), element_kind(element_kind) {
			this->checkHomogeneity();
			}

			/* ------------------------------------------------------------------------ */
			/* Methods */
			/* ------------------------------------------------------------------------ */
			public:
			/// get the number of components of the hosted field
			ElementTypeMap<UInt>
			getNbComponents(UInt dim = _all_dimensions, GhostType ghost_type = _not_ghost,
			ElementKind kind = _ek_not_defined) override {
			return this->field.getNbComponents(dim, ghost_type, kind);
			};

			/// return the size of the contained data: i.e. the number of elements ?
			virtual UInt size() {
			checkHomogeneity();
			return this->nb_total_element;
			}

			/// return the iohelper datatype to be dumped
			iohelper::DataType getDataType() {
			return iohelper::getDataType<data_type>();
			}

			protected:
			/// return the number of entries per element
			UInt getNbDataPerElem(const ElementType & type,
			const GhostType & ghost_type = _not_ghost) const {
			if (!nb_data_per_elem.exists(type, ghost_type))
			return field(type, ghost_type).getNbComponent();

			return nb_data_per_elem(type, this->ghost_type);
			}

			/// check if the same quantity of data for all element types
			void checkHomogeneity() override;

			public:
			void registerToDumper(const std::string & id,
			iohelper::Dumper & dumper) override {
			dumper.addElemDataField(id, *this);
			};

			/// for connection to a FieldCompute
			inline Field * connect(FieldComputeProxy & proxy) override {
			return proxy.connectToField(this);
			}

			/// for connection to a Homogenizer
			inline ComputeFunctorInterface * connect(HomogenizerProxy & proxy) override {
			return proxy.connectToField(this);
			};

			virtual iterator begin() {
			field_type_iterator tit;
			field_type_iterator end;

			/// type iterators on the elemental field
			tit = this->field.firstType(this->spatial_dimension, this->ghost_type,
			this->element_kind);
			end = this->field.lastType(this->spatial_dimension, this->ghost_type,
			this->element_kind);

			/// skip all types without data
			ElementType type = *tit;
			for (; tit != end && this->field(*tit, this->ghost_type).size() == 0;
			++tit) {
			}

			type = *tit;

			if (tit == end)
			return this->end();

			/// getting information for the field of the given type
			const array_type & vect = this->field(type, this->ghost_type);
			UInt nb_data_per_elem = this->getNbDataPerElem(type);
			UInt nb_component = vect.getNbComponent();
			UInt size = (vect.size() * nb_component) / nb_data_per_elem;

			/// define element-wise iterator
			array_iterator it = vect.begin_reinterpret(nb_data_per_elem, size);
			array_iterator it_end = vect.end_reinterpret(nb_data_per_elem, size);
			/// define data iterator
			iterator rit =
			iterator(this->field, tit, end, it, it_end, this->ghost_type);
			rit.setNbDataPerElem(this->nb_data_per_elem);
			return rit;
			}

			virtual iterator end() {
			field_type_iterator tit;
			field_type_iterator end;

			tit = this->field.firstType(this->spatial_dimension, this->ghost_type,
			this->element_kind);
			end = this->field.lastType(this->spatial_dimension, this->ghost_type,
			this->element_kind);

			ElementType type = *tit;
			for (; tit != end; ++tit)
			type = *tit;

			const array_type & vect = this->field(type, this->ghost_type);
			UInt nb_data = this->getNbDataPerElem(type);
			UInt nb_component = vect.getNbComponent();
			UInt size = (vect.size() * nb_component) / nb_data;
			array_iterator it = vect.end_reinterpret(nb_data, size);

			iterator rit = iterator(this->field, end, end, it, it, this->ghost_type);
			rit.setNbDataPerElem(this->nb_data_per_elem);
			return rit;
			}

			virtual UInt getDim() {
			if (this->homogeneous) {
			field_type_iterator tit = this->field.firstType(
			this->spatial_dimension, this->ghost_type, this->element_kind);
			return this->getNbDataPerElem(*tit);
			}

			throw;
			return 0;
			}

			void setNbDataPerElem(const ElementTypeMap<UInt> & nb_data) override {
			nb_data_per_elem = nb_data;
			}

			/* ------------------------------------------------------------------------ */
			/* Class Members */
			/* ------------------------------------------------------------------------ */
			protected:
			/// the ElementTypeMapArray embedded in the field
			const field_type & field;
			/// total number of elements
			UInt nb_total_element;
			/// the spatial dimension of the problem
			UInt spatial_dimension;
			/// whether this is a ghost field or not (for type selection)
			GhostType ghost_type;
			/// The element kind to operate on
			ElementKind element_kind;
			/// The number of data per element type
			ElementTypeMap<UInt> nb_data_per_elem;
			};

			/* -------------------------------------------------------------------------- */
			#include "dumper_generic_elemental_field_tmpl.hh"
			/* -------------------------------------------------------------------------- */

			__END_AKANTU_DUMPER__ } // akantu

			#endif /* __AKANTU_DUMPER_GENERIC_ELEMENTAL_FIELD_HH__ */

rAKA/src/io/dumper/dumper_generic_elemental_field.hh3cb848d814e7bugfixes/contact_mechan…

dumper_generic_elemental_field.hhView Options

rAKA/src/io/dumper/dumper_generic_elemental_field.hh
3cb848d814e7bugfixes/contact_mechan…

dumper_generic_elemental_field.hh
View Options