Diffusion akantu (features/node_to_node) Browse

node_to_node

element_info_per_processor_tmpl.hh
View Options

			/**
			* @file element_info_per_processor_tmpl.hh
			*
			* @author Nicolas Richart <nicolas.richart@epfl.ch>
			*
			* @date Fri Mar 11 15:03:12 2016
			*
			* @brief
			*
			* @section LICENSE
			*
			* Copyright (©) 2010-2011 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/>.
			*
			*/

			/* -------------------------------------------------------------------------- */
			#include "element_info_per_processor.hh"
			#include "element_group.hh"
			#include "mesh.hh"
			/* -------------------------------------------------------------------------- */

			#ifndef __AKANTU_ELEMENT_INFO_PER_PROCESSOR_TMPL_HH__
			#define __AKANTU_ELEMENT_INFO_PER_PROCESSOR_TMPL_HH__

			__BEGIN_AKANTU__

			/* -------------------------------------------------------------------------- */
			template <typename T, typename BufferType>
			void ElementInfoPerProc::fillMeshDataTemplated(BufferType & buffer,
			const std::string & tag_name) {

			AKANTU_DEBUG_ASSERT(this->mesh.getNbElement(this->type) == nb_local_element,
			"Did not got enought informations for the tag "
			<< tag_name << " and the element type " << this->type
			<< ":"
			<< "_not_ghost."
			<< " Got " << nb_local_element << " values, expected "
			<< mesh.getNbElement(this->type));
			MeshData & mesh_data = this->getMeshData();

			UInt nb_component = mesh_data.getNbComponent(tag_name, this->type);
			Array<T> & data = mesh_data.getElementalDataArrayAlloc<T>(
			tag_name, this->type, _not_ghost, nb_component);
			data.resize(nb_local_element);
			/// unpacking the data, element by element
			for (UInt i(0); i < nb_local_element; ++i) {
			for (UInt j(0); j < nb_component; ++j) {
			buffer >> data(i, j);
			}
			}

			AKANTU_DEBUG_ASSERT(mesh.getNbElement(this->type, _ghost) == nb_ghost_element,
			"Did not got enought informations for the tag "
			<< tag_name << " and the element type " << this->type
			<< ":"
			<< "_ghost."
			<< " Got " << nb_ghost_element << " values, expected "
			<< mesh.getNbElement(this->type, _ghost));

			mesh_data.registerElementalData<T>(tag_name);
			Array<T> & data_ghost = mesh_data.getElementalDataArrayAlloc<T>(
			tag_name, this->type, _ghost, nb_component);
			data_ghost.resize(nb_ghost_element);

			/// unpacking the ghost data, element by element
			for (UInt j(0); j < nb_ghost_element; ++j) {
			for (UInt k(0); k < nb_component; ++k) {
			buffer >> data_ghost(j, k);
			}
			}
			}

			/* -------------------------------------------------------------------------- */
			template <typename BufferType>
			void ElementInfoPerProc::fillMeshData(BufferType & buffer,
			const std::string & tag_name) {
			#define AKANTU_DISTRIBUTED_SYNHRONIZER_TAG_DATA(r, extra_param, elem) \
			case BOOST_PP_TUPLE_ELEM(2, 0, elem): { \
			fillMeshDataTemplated<BOOST_PP_TUPLE_ELEM(2, 1, elem)>(buffer, tag_name); \
			break; \
			}

			MeshData & mesh_data = this->getMeshData();
			MeshDataTypeCode type_code = mesh_data.getTypeCode(tag_name);

			switch (type_code) {
			BOOST_PP_SEQ_FOR_EACH(AKANTU_DISTRIBUTED_SYNHRONIZER_TAG_DATA, ,
			AKANTU_MESH_DATA_TYPES)
			default:
			AKANTU_DEBUG_ERROR("Could not determine the type of tag" << tag_name
			<< "!");
			break;
			}
			#undef AKANTU_DISTRIBUTED_SYNHRONIZER_TAG_DATA
			}

			/* -------------------------------------------------------------------------- */
			template <class CommunicationBuffer>
			void ElementInfoPerProc::fillElementGroupsFromBuffer(
			CommunicationBuffer & buffer) {
			AKANTU_DEBUG_IN();

			Element el;
			el.type = type;

			for (ghost_type_t::iterator gt = ghost_type_t::begin();
			gt != ghost_type_t::end(); ++gt) {
			UInt nb_element = mesh.getNbElement(type, *gt);
			el.ghost_type = *gt;

			for (UInt e = 0; e < nb_element; ++e) {
			el.element = e;

			std::vector<std::string> element_to_group;
			buffer >> element_to_group;

			AKANTU_DEBUG_ASSERT(e < mesh.getNbElement(type, *gt),
			"The mesh does not have the element " << e);

			std::vector<std::string>::iterator it = element_to_group.begin();
			std::vector<std::string>::iterator end = element_to_group.end();
			for (; it != end; ++it) {
			mesh.getElementGroup(*it).add(el, false, false);
			}
			}
			}

			AKANTU_DEBUG_OUT();
			}

			/* -------------------------------------------------------------------------- */

			__END_AKANTU__

			#endif /* __AKANTU_ELEMENT_INFO_PER_PROCESSOR_TMPL_HH__ */

rAKA/src/synchronizer/element_info_per_processor_tmpl.hhcafceb59140efeatures/node_to_node

element_info_per_processor_tmpl.hhView Options

rAKA/src/synchronizer/element_info_per_processor_tmpl.hh
cafceb59140efeatures/node_to_node

element_info_per_processor_tmpl.hh
View Options