unit_kernel.hpp

Go to the documentation of this file.

// This file is a part of NiHu, a C++ BEM template library.
//
// Copyright (C) 2012-2014  Peter Fiala <fiala@hit.bme.hu>
// Copyright (C) 2012-2014  Peter Rucz <rucz@hit.bme.hu>
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
#ifndef UNIT_KERNEL_HPP_INCLUDED
#define UNIT_KERNEL_HPP_INCLUDED
 
#include "../core/kernel.hpp"
#include "../core/gaussian_quadrature.hpp"
 
namespace NiHu
{
 
template <class Scalar>
class unit_kernel;
 
template<class Scalar>
struct kernel_traits<unit_kernel<Scalar> >
{
    typedef empty_wall test_input_t;
    typedef empty_wall trial_input_t;
    typedef Scalar result_t;
    typedef gauss_family_tag quadrature_family_t;
    static bool const is_symmetric = true;
    static unsigned const singularity_order = 0;
    static unsigned const singular_quadrature_order = 0;
};
 
template <class Scalar>
class unit_kernel :
    public kernel_base<unit_kernel<Scalar> >
{
public:
    typedef kernel_base<unit_kernel<Scalar> > base_t;
    typedef typename base_t::scalar_t scalar_t;
    typedef typename base_t::result_t result_t;
    typedef typename base_t::test_input_t test_input_t;
    typedef typename base_t::trial_input_t trial_input_t;
 
    constexpr result_t operator()(
        test_input_t const &, trial_input_t const &) const
    {
        return result_t(1.0);
    }
 
    constexpr unsigned estimate_complexity(
        test_input_t const &, trial_input_t const &, scalar_t const &) const
    {
        return 0;
    }
};
 
}
 
#endif // UNIT_KERNEL_HPP_INCLUDED