convected_helmholtz_3d_singular_integrals.hpp

#ifndef NIHU_CONVECTED_HELMHOLTZ_3D_SINGULAR_INTEGRALS_HPP_INCLUDED
#define NIHU_CONVECTED_HELMHOLTZ_3D_SINGULAR_INTEGRALS_HPP_INCLUDED
 
#include "core/element_match.hpp"
#include "core/singular_integral_shortcut.hpp"
 
#include "convected_helmholtz_3d_kernel.hpp"
#include "field_type_helpers.hpp"
 
#define NIHU_GUIGGIANI_DEBUGGING 1
 
namespace NiHu
{
 
template <class WaveNumber, class TestField, class TrialField>
class singular_integral_shortcut<
    convected_helmholtz_3d_DLPt_kernel<WaveNumber>, TestField, TrialField, match::match_2d_type,
    typename std::enable_if<
        is_collocational<TestField, TrialField>::value 
    >::type
>
{
public:
    template <class result_t>
    static result_t &eval(
        result_t &result,
        kernel_base<convected_helmholtz_3d_DLPt_kernel<WaveNumber> > const &kernel,
        field_base<TestField> const &,
        field_base<TrialField> const &trial_field,
        element_match const &)
    {
#if NIHU_GUIGGIANI_DEBUGGING
        static bool printed = false;
        if (!printed)
        {
            std::cout << "Calling Guiggiani for convected_helmholtz_3d_DLPt_kernel now" << std::endl;
            printed = true;
        }
#endif
        typedef guiggiani<TrialField, convected_helmholtz_3d_DLPt_kernel<WaveNumber>, 9, 19> guiggiani_t;
        auto const &elem = trial_field.get_elem();
        guiggiani_t gui(elem, kernel.derived());
        for (unsigned r = 0; r < TestField::num_dofs; ++r)
        {
            auto const &xi0 = TestField::nset_t::corner_at(r);
            gui.integrate(result.row(r), xi0, elem.get_normal(xi0));
        }
        
        return result;
    }
};
 
 
template <class WaveNumber, class TestField, class TrialField>
class singular_integral_shortcut<
    convected_helmholtz_3d_DLPt_n_kernel<WaveNumber>, TestField, TrialField, match::match_2d_type,
    typename std::enable_if<
        is_collocational<TestField, TrialField>::value 
    >::type
>
{
public:
    template <class result_t>
    static result_t &eval(
        result_t &result,
        kernel_base<convected_helmholtz_3d_DLPt_n_kernel<WaveNumber> > const &kernel,
        field_base<TestField> const &,
        field_base<TrialField> const &trial_field,
        element_match const &)
    {
#if NIHU_GUIGGIANI_DEBUGGING
        static bool printed = false;
        if (!printed)
        {
            std::cout << "Calling Guiggiani for convected_helmholtz_3d_DLPt_n_kernel now" << std::endl;
            printed = true;
        }
#endif
        typedef guiggiani<TrialField, convected_helmholtz_3d_DLPt_n_kernel<WaveNumber>, 9, 19> guiggiani_t;
        auto const &elem = trial_field.get_elem();
        guiggiani_t gui(elem, kernel.derived());
        for (unsigned r = 0; r < TestField::num_dofs; ++r)
        {
            auto const &xi0 = TestField::nset_t::corner_at(r);
            gui.integrate(result.row(r), xi0, elem.get_normal(xi0));
        }
        
        return result;
    }
};
 
 
 
template <class WaveNumber, class TestField, class TrialField>
class singular_integral_shortcut<
    convected_helmholtz_3d_HSP_kernel<WaveNumber>, TestField, TrialField, match::match_2d_type,
    typename std::enable_if<
        is_collocational<TestField, TrialField>::value 
    >::type
>
{
public:
    template <class result_t>
    static result_t &eval(
        result_t &result,
        kernel_base<convected_helmholtz_3d_HSP_kernel<WaveNumber> > const &kernel,
        field_base<TestField> const &,
        field_base<TrialField> const &trial_field,
        element_match const &)
    {
#if NIHU_GUIGGIANI_DEBUGGING
        static bool printed = false;
        if (!printed)
        {
            std::cout << "Calling Guiggiani for convected_helmholtz_3d_HSP_kernel now" << std::endl;
            printed = true;
        }
#endif
        typedef guiggiani<TrialField, convected_helmholtz_3d_HSP_kernel<WaveNumber>, 9, 19> guiggiani_t;
        auto const &elem = trial_field.get_elem();
        guiggiani_t gui(elem, kernel.derived());
        for (unsigned r = 0; r < TestField::num_dofs; ++r)
        {
            auto const &xi0 = TestField::nset_t::corner_at(r);
            gui.integrate(result.row(r), xi0, elem.get_normal(xi0));
        }
        
        return result;
    }
};
 
} // end of namespace NiHu
 
#endif /* NIHU_CONVECTED_HELMHOLTZ_3D_SINGULAR_INTEGRALS_HPP_INCLUDED */