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

rAKA/extra_packages/igfem/src/interpolation_element_igfem.cc
b56858367c28bugfixes/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
*
*
* 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"
namespace 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]};
} // namespace akantu
c4science · Help