tria_1_gauss_shape_set.hpp

#ifndef TRIA_1_GAUSS_SHAPE_SET_HPP_INCLUDED
#define TRIA_1_GAUSS_SHAPE_SET_HPP_INCLUDED
 
#include "../core/shapeset.hpp"
#include "lib_domain.hpp"
 
// define the new shape set's traits
namespace NiHu
{
    
// predefine the new shape set
class tria_1_gauss_shape_set;
 
namespace shape_set_traits
{
    // the shape set is defined over the tria domain
    template <>
    struct domain<tria_1_gauss_shape_set> : tria_domain {};
 
    // it has 4 shapeset nodes
    template <>
    struct num_nodes<tria_1_gauss_shape_set> { enum { value = 3 }; };
 
    // if used as geometrical shape set, Jacobian is 1-st order in xi/eta
    template <>
    struct jacobian_order<tria_1_gauss_shape_set> { enum { value = 0 }; };
 
    // N(xi,eta) is first order
    template <>
    struct polynomial_order<tria_1_gauss_shape_set> { enum { value = 1 }; };
 
    // N(xi, eta) is computed on the fly
    template <>
    struct shape_complexity<tria_1_gauss_shape_set, 0>
    {
        typedef matrix_function_complexity::general type;
    };
 
    // N'(xi, eta) is constant
    template <>
    struct shape_complexity<tria_1_gauss_shape_set, 1>
    {
        typedef matrix_function_complexity::constant type;
    };
 
    // N''(xi, eta) is zero
    template <>
    struct shape_complexity<tria_1_gauss_shape_set, 2>
    {
        typedef matrix_function_complexity::zero type;
    };
 
    // DOF locations are inside the element (2DOF)
    template <>
    struct position_dof_vector<tria_1_gauss_shape_set>
    {
        typedef tmp::vector<dof2, dof2, dof2> type;
    };
 
    template <>
    struct corner_index_vector<tria_1_gauss_shape_set>
    {
        typedef tmp::vector<not_corner, not_corner, not_corner> type;
    };
 
} // end of namespace shape_set_traits
 
 
// define the new shape set class
class tria_1_gauss_shape_set
    : public shape_set_base<tria_1_gauss_shape_set>
{
    typedef shape_set_base<tria_1_gauss_shape_set> base_t;
public:
 
    using base_t::corner_at;
    
    // return iterator to corners
    static corner_iterator_t corner_begin_impl()
    {
        return m_corners.cbegin();
    }
 
protected:
    static const corners_t m_corners;
};
 
// define expression for computing the zeroth order derivative N(xi, eta)
template <>
class shape_function<tria_1_gauss_shape_set, 0>
{
    typedef shape_set_traits::shape_value_type<tria_1_gauss_shape_set, 0>::type shape_t;
    typedef shape_set_traits::domain<tria_1_gauss_shape_set>::type::xi_t xi_t;
public:
    static shape_t eval(xi_t const &_xi)
    {
        auto xi = _xi[0], eta = _xi[1];
        return ( tria_1_gauss_shape_set::shape_t() <<
            5./3. - 2.*xi - 2.*eta,
            2. * xi - 1./3.,
            2. * eta - 1./3.
        ).finished();
    }
};
 
 
// define expression for computing the first order derivative N(xi, eta)
template <>
class shape_function<tria_1_gauss_shape_set, 1>
{
    typedef shape_set_traits::shape_value_type<tria_1_gauss_shape_set, 1>::type shape_t;
    typedef shape_set_traits::domain<tria_1_gauss_shape_set>::type::xi_t xi_t;
public:
    static shape_t eval(xi_t const &)
    {
        return ( tria_1_gauss_shape_set::dshape_t() <<
            -2., -2.,
            2., 0.,
            0., 2.
        ).finished();
    }
};
 
} // end of namespace NiHu
 
#endif // TRIA_1_GAUSS_SHAPE_SET_HPP_INCLUDED