Homec4science
Diffusion akantu (bugfixes/39-the-python-example-for-phasefield-are-broken-due-to-a-change-to-setintegrationscheme) Browse

rAKA/src/interpolation_element_igfem.cc
266687a791a4bugfixes/39-the-python…

Last Change
Raw File

interpolation_element_igfem.cc
View Options

/**
* @file element_class_igfem.hh
*
* @author Aurelia Isabel Cuba Ramos <aurelia.cubaramos@epfl.ch>
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
*
* @brief Specialization for interface-enriched finite elements
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
*/
/* -------------------------------------------------------------------------- */
#include "element_class.hh"
__BEGIN_AKANTU__
/* -------------------------------------------------------------------------- */
template<> UInt InterpolationElement<_itp_igfem_segment_3>::nb_sub_elements = 2;
template<> UInt InterpolationElement<_itp_igfem_triangle_4>::nb_sub_elements = 2;
template<> UInt InterpolationElement<_itp_igfem_triangle_5>::nb_sub_elements = 2;
/* !!! stored as a matrix nb_subelements X nb_nodes_per_subelement in COL MAJOR */
/* -------------------------------------------------------------------------- */
template<> UInt InterpolationElement<_itp_igfem_segment_3>::sub_element_connectivity_vect[] = { 0, // first type
2,
2, // second type
1};
template<> UInt InterpolationElement<_itp_igfem_triangle_4>::sub_element_connectivity_vect[] = { 0, // first type
1,
3,
0, // second type
3,
2};
template<> UInt InterpolationElement<_itp_igfem_triangle_5>::sub_element_connectivity_vect[] = { 0, // first type
3,
4,
3, // second type
1,
2,
4};
template<> UInt * InterpolationElement<_itp_igfem_segment_3>::sub_element_connectivity[] = { &sub_element_connectivity_vect[0],
&sub_element_connectivity_vect[2] };
template<> UInt * InterpolationElement<_itp_igfem_triangle_4>::sub_element_connectivity[] = { &sub_element_connectivity_vect[0],
&sub_element_connectivity_vect[3] };
template<> UInt * InterpolationElement<_itp_igfem_triangle_5>::sub_element_connectivity[] = { &sub_element_connectivity_vect[0],
&sub_element_connectivity_vect[3] };
/* -------------------------------------------------------------------------- */
template<> UInt InterpolationElement<_itp_igfem_segment_3>::nb_enrichments = 1;
template<> UInt InterpolationElement<_itp_igfem_triangle_4>::nb_enrichments = 1;
template<> UInt InterpolationElement<_itp_igfem_triangle_5>::nb_enrichments = 2;
template<> UInt InterpolationElement<_itp_igfem_segment_3>::enrichment_vect[] = { 1, // on first subelement
0}; // on second subelement
template<> UInt InterpolationElement<_itp_igfem_triangle_4>::enrichment_vect[] = { 2, // on first subelement
1}; // on second subelement
template<> UInt InterpolationElement<_itp_igfem_triangle_5>::enrichment_vect[] = { 1, // on first subelement
2,
0,
3}; // on second subelement
template<> UInt * InterpolationElement<_itp_igfem_segment_3>::enrichments[] = { &enrichment_vect[0],
&enrichment_vect[1]};
template<> UInt * InterpolationElement<_itp_igfem_triangle_4>::enrichments[] = { &enrichment_vect[0],
&enrichment_vect[1]};
template<> UInt * InterpolationElement<_itp_igfem_triangle_5>::enrichments[] = { &enrichment_vect[0],
&enrichment_vect[2]};
__END_AKANTU__
c4science · Help