gauss_function_space_view.hpp

#ifndef NIHU_GAUSS_FUNCTION_SPACE_VIEW_HPP_INCLUDED
#define NIHU_GAUSS_FUNCTION_SPACE_VIEW_HPP_INCLUDED
 
#include "gauss_field_view.hpp"
 
#include "core/function_space.hpp"
#include "tmp/algorithm.hpp"
 
namespace NiHu
{
namespace function_space_impl_internal
{
    template <class Mesh, class Dimension>
    struct get_num_dofs<Mesh, field_option::gauss, Dimension>
    {
        template <class ElemType>
        struct count_elem_type_nodes
        {
            struct type 
            {
                void operator()(Mesh const& mesh, size_t &count)
                {
                    typedef typename field_traits::internal::assign_gauss_shape<typename ElemType::domain_t>::type nset_t;
            
                    count += 
                        (mesh.template end<ElemType>() - mesh.template begin<ElemType>()) *     // Number of elements of the given type
                        shape_set_traits::num_nodes<nset_t>::value;                        // Number of nodes per element
                }
            };
        };
 
        typedef function_space_traits< function_space_view<Mesh, field_option::gauss, Dimension> > traits_t;
        static size_t eval(Mesh const& mesh)
        {
            size_t count = 0;
            tmp::call_each<
                typename Mesh::elem_type_vector_t,
                count_elem_type_nodes<tmp::_1>,
                Mesh const&, 
                size_t&                             // Explicit type info required for passing reference
            >(mesh, count);
            return count * traits_t::quantity_dimension;
        }
    };
} // namespace NiHu::function_space_impl_internal
 
template <class Mesh, class Dimension = field_dimension::_1d>
function_space_view<Mesh, field_option::gauss, Dimension> const &
    gauss_view(mesh_base<Mesh> const &msh, Dimension dim = Dimension())
{
    return create_function_space_view(msh.derived(), field_option::gauss(), dim);
}
 
} // namespace NiHu
#endif /* NIHU_GAUSS_FUNCTION_SPACE_VIEW_HPP_INCLUDED*/