fe_engine_inline_impl.hh
No OneTemporary
Actions

Subscribers

None

File Metadata

Created: Sun, Nov 3, 17:47

fe_engine_inline_impl.hh
View Options

	/**
	* @file fe_engine_inline_impl.hh
	*
	* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
	* @author Nicolas Richart <nicolas.richart@epfl.ch>
	*
	* @date creation: Tue Jul 20 2010
	* @date last modification: Fri Dec 11 2020
	*
	* @brief Implementation of the inline functions of the FEEngine Class
	*
	*
	* @section LICENSE
	*
	* Copyright (©) 2010-2021 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_class.hh"
	#include "fe_engine.hh"
	#include "mesh.hh"
	/* -------------------------------------------------------------------------- */
	#include "element_type_conversion.hh"
	/* -------------------------------------------------------------------------- */

	#ifndef AKANTU_FE_ENGINE_INLINE_IMPL_HH_
	#define AKANTU_FE_ENGINE_INLINE_IMPL_HH_

	namespace akantu {

	/* -------------------------------------------------------------------------- */
	inline Real FEEngine::getElementInradius(const Matrix<Real> & coord,
	ElementType type) {
	Real inradius = 0;

	#define GET_INRADIUS(type) inradius = ElementClass<type>::getInradius(coord);

	AKANTU_BOOST_ALL_ELEMENT_SWITCH(GET_INRADIUS);
	#undef GET_INRADIUS

	return inradius;
	}

	/* -------------------------------------------------------------------------- */
	inline Real FEEngine::getElementInradius(const Element & element) const {
	auto spatial_dimension = mesh.getSpatialDimension();
	auto positions = make_view(mesh.getNodes(), spatial_dimension).begin();
	auto connectivity = mesh.getConnectivities().get(element);

	Matrix<Real> coords(spatial_dimension, connectivity.size());

	for(auto && data : zip(connectivity, coords)) {
	std::get<1>(data) = positions[std::get<0>(data)];
	}

	return getElementInradius(coords, element.type);
	}


	/* -------------------------------------------------------------------------- */
	inline InterpolationType FEEngine::getInterpolationType(ElementType type) {
	return convertType<ElementType, InterpolationType>(type);
	}

	/* -------------------------------------------------------------------------- */
	/// @todo rewrite this function in order to get the cohesive element
	/// type directly from the facet
	#if defined(AKANTU_COHESIVE_ELEMENT)
	inline ElementType FEEngine::getCohesiveElementType(ElementType type) {
	ElementType ctype;
	#define GET_COHESIVE_TYPE(type) \
	ctype = CohesiveFacetProperty<type>::cohesive_type;

	AKANTU_BOOST_ALL_ELEMENT_SWITCH(GET_COHESIVE_TYPE);
	#undef GET_COHESIVE_TYPE

	return ctype;
	}
	#else
	inline ElementType FEEngine::getCohesiveElementType(__attribute__((unused))
	ElementType type_facet) {
	return _not_defined;
	}
	#endif

	/* -------------------------------------------------------------------------- */
	#if defined(AKANTU_IGFEM)
	} // akantu
	#include "igfem_helper.hh"
	namespace akantu {

	inline Vector<ElementType> FEEngine::getIGFEMElementTypes(ElementType type) {

	#define GET_IGFEM_ELEMENT_TYPES(type) \
	return IGFEMHelper::getIGFEMElementTypes<type>();

	AKANTU_BOOST_REGULAR_ELEMENT_SWITCH(GET_IGFEM_ELEMENT_TYPES);

	#undef GET_IGFEM_ELEMENT_TYPES
	}
	#endif

	/* -------------------------------------------------------------------------- */
	template <typename T>
	void FEEngine::extractNodalToElementField(const Mesh & mesh,
	const Array<T> & nodal_f,
	Array<T> & elemental_f,
	ElementType type,
	GhostType ghost_type,
	const Array<UInt> & filter_elements) {
	AKANTU_DEBUG_IN();

	UInt nb_nodes_per_element = Mesh::getNbNodesPerElement(type);
	UInt nb_degree_of_freedom = nodal_f.getNbComponent();
	UInt nb_element = mesh.getNbElement(type, ghost_type);
	UInt * conn_val = mesh.getConnectivity(type, ghost_type).storage();

	if (filter_elements != empty_filter) {
	nb_element = filter_elements.size();
	}

	elemental_f.resize(nb_element);

	T * nodal_f_val = nodal_f.storage();
	T * f_val = elemental_f.storage();

	UInt * el_conn;
	for (UInt el = 0; el < nb_element; ++el) {
	if (filter_elements != empty_filter) {
	el_conn = conn_val + filter_elements(el) * nb_nodes_per_element;
	} else {
	el_conn = conn_val + el * nb_nodes_per_element;
	}

	for (UInt n = 0; n < nb_nodes_per_element; ++n) {
	UInt node = *(el_conn + n);
	std::copy(nodal_f_val + node * nb_degree_of_freedom,
	nodal_f_val + (node + 1) * nb_degree_of_freedom, f_val);
	f_val += nb_degree_of_freedom;
	}
	}

	AKANTU_DEBUG_OUT();
	}

	/* -------------------------------------------------------------------------- */
	template <typename T>
	void FEEngine::filterElementalData(const Mesh & mesh, const Array<T> & elem_f,
	Array<T> & filtered_f, ElementType type,
	GhostType ghost_type,
	const Array<UInt> & filter_elements) {
	AKANTU_DEBUG_IN();

	UInt nb_element = mesh.getNbElement(type, ghost_type);
	if (nb_element == 0) {
	filtered_f.resize(0);
	return;
	}

	UInt nb_degree_of_freedom = elem_f.getNbComponent();
	UInt nb_data_per_element = elem_f.size() / nb_element;

	if (filter_elements != empty_filter) {
	nb_element = filter_elements.size();
	}

	filtered_f.resize(nb_element * nb_data_per_element);

	T * elem_f_val = elem_f.storage();
	T * f_val = filtered_f.storage();

	UInt el_offset;
	for (UInt el = 0; el < nb_element; ++el) {
	if (filter_elements != empty_filter) {
	el_offset = filter_elements(el);
	} else {
	el_offset = el;
	}

	std::copy(elem_f_val +
	el_offset * nb_data_per_element * nb_degree_of_freedom,
	elem_f_val +
	(el_offset + 1) * nb_data_per_element * nb_degree_of_freedom,
	f_val);
	f_val += nb_degree_of_freedom * nb_data_per_element;
	}

	AKANTU_DEBUG_OUT();
	}

	} // namespace akantu

	#endif /* AKANTU_FE_ENGINE_INLINE_IMPL_HH_ */

fe_engine_inline_impl.hhNo OneTemporaryActions

File Metadata

fe_engine_inline_impl.hhView Options

Event Timeline

fe_engine_inline_impl.hh
No OneTemporary
Actions

fe_engine_inline_impl.hh
View Options