embedded_interface_intersector.cc
No OneTemporary
Actions

Subscribers

None

File Metadata

Created: Sun, Nov 10, 22:04

embedded_interface_intersector.cc
View Options

	/**
	* @file embedded_interface_intersector.cc
	*
	* @author Lucas Frerot <lucas.frerot@epfl.ch>
	*
	* @date creation: Wed Apr 29 2015
	* @date last modification: Wed Apr 29 2015
	*
	* @brief Class that loads the interface from mesh and computes intersections
	*
	* @section LICENSE
	*
	* Copyright (©) 2010-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/>.
	*
	*/

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

	#include "embedded_interface_intersector.hh"
	#include "mesh_segment_intersector.hh"

	/// Helper macro for types in the mesh. Creates an intersector and computes intersection queries
	#define INTERFACE_INTERSECTOR_CASE(dim, type) do { \
	MeshSegmentIntersector<dim, type> intersector(this->mesh, interface_mesh); \
	name_to_primitives_it = name_to_primitives_map.begin(); \
	for (; name_to_primitives_it != name_to_primitives_end ; ++name_to_primitives_it) { \
	intersector.setPhysicalName(name_to_primitives_it->first); \
	intersector.buildResultFromQueryList(name_to_primitives_it->second); \
	} } while(0)

	#define INTERFACE_INTERSECTOR_CASE_2D(type) INTERFACE_INTERSECTOR_CASE(2, type)
	#define INTERFACE_INTERSECTOR_CASE_3D(type) INTERFACE_INTERSECTOR_CASE(3, type)

	__BEGIN_AKANTU__

	EmbeddedInterfaceIntersector::EmbeddedInterfaceIntersector(Mesh & mesh, const Mesh & primitive_mesh) :
	MeshGeomAbstract(mesh),
	interface_mesh(mesh.getSpatialDimension(), "interface_mesh"),
	primitive_mesh(primitive_mesh)
	{
	// Initiating mesh connectivity and data
	interface_mesh.addConnectivityType(_segment_2, _not_ghost);
	interface_mesh.addConnectivityType(_segment_2, _ghost);
	interface_mesh.registerData<Element>("associated_element").alloc(0, 1, _segment_2);
	interface_mesh.registerData<std::string>("physical_names").alloc(0, 1, _segment_2);
	}

	EmbeddedInterfaceIntersector::~EmbeddedInterfaceIntersector()
	{}

	void EmbeddedInterfaceIntersector::constructData() {
	AKANTU_DEBUG_IN();

	const UInt dim = this->mesh.getSpatialDimension();

	if (dim == 1)
	AKANTU_DEBUG_ERROR("No embedded model in 1D. Deactivate intersection initialization");

	Array<std::string> * physical_names = NULL;

	try {
	physical_names = &const_cast<Array<std::string> &>(this->primitive_mesh.getData<std::string>("physical_names", _segment_2));
	} catch (debug::Exception & e) {
	AKANTU_DEBUG_ERROR("You must define physical names to reinforcements in order to use the embedded model");
	throw e;
	}

	const UInt nb_nodes_per_element = Mesh::getNbNodesPerElement(_segment_2);
	Array<UInt>::const_vector_iterator connectivity = primitive_mesh.getConnectivity(_segment_2).begin(nb_nodes_per_element);

	Array<std::string>::scalar_iterator
	names_it = physical_names->begin(),
	names_end = physical_names->end();

	std::map<std::string, std::list<K::Segment_3> > name_to_primitives_map;

	// Loop over the physical names and register segment lists in name_to_primitives_map
	for (; names_it != names_end ; ++names_it) {
	UInt element_id = names_it - physical_names->begin();
	const Vector<UInt> el_connectivity = connectivity[element_id];

	K::Segment_3 segment = this->createSegment(el_connectivity);
	name_to_primitives_map[*names_it].push_back(segment);
	}

	// Loop over the background types of the mesh
	Mesh::type_iterator
	type_it = this->mesh.firstType(dim, _not_ghost),
	type_end = this->mesh.lastType(dim, _not_ghost);

	std::map<std::string, std::list<K::Segment_3> >::iterator
	name_to_primitives_it,
	name_to_primitives_end = name_to_primitives_map.end();

	for (; type_it != type_end ; ++type_it) {
	// Used in AKANTU_BOOST_ELEMENT_SWITCH
	ElementType type = *type_it;

	AKANTU_DEBUG_INFO("Computing intersections with background element type " << type);

	switch(dim) {
	case 1:
	break;

	case 2:
	// Compute intersections for supported 2D elements
	AKANTU_BOOST_ELEMENT_SWITCH(INTERFACE_INTERSECTOR_CASE_2D, (_triangle_3));
	break;

	case 3:
	// Compute intersections for supported 3D elements
	AKANTU_BOOST_ELEMENT_SWITCH(INTERFACE_INTERSECTOR_CASE_3D, (_tetrahedron_4));
	break;
	}
	}

	AKANTU_DEBUG_OUT();
	}

	K::Segment_3 EmbeddedInterfaceIntersector::createSegment(const Vector<UInt> & connectivity) {
	AKANTU_DEBUG_IN();

	K::Point_3 * source = NULL, * target = NULL;
	const Array<Real> & nodes = this->primitive_mesh.getNodes();

	if (this->mesh.getSpatialDimension() == 2) {
	source = new K::Point_3(nodes(connectivity(0), 0), nodes(connectivity(0), 1), 0.);
	target = new K::Point_3(nodes(connectivity(1), 0), nodes(connectivity(1), 1), 0.);
	} else if (this->mesh.getSpatialDimension() == 3) {
	source = new K::Point_3(nodes(connectivity(0), 0), nodes(connectivity(0), 1), nodes(connectivity(0), 2));
	target = new K::Point_3(nodes(connectivity(1), 0), nodes(connectivity(1), 1), nodes(connectivity(1), 2));
	}

	K::Segment_3 segment(source, target);
	delete source;
	delete target;

	AKANTU_DEBUG_OUT();
	return segment;
	}

	__END_AKANTU__

	#undef INTERFACE_INTERSECTOR_CASE
	#undef INTERFACE_INTERSECTOR_CASE_2D
	#undef INTERFACE_INTERSECTOR_CASE_3D

embedded_interface_intersector.ccNo OneTemporaryActions

File Metadata

embedded_interface_intersector.ccView Options

Event Timeline

embedded_interface_intersector.cc
No OneTemporary
Actions

embedded_interface_intersector.cc
View Options