complexity_estimator.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 COMPLEXITY_ESTIMATOR_HPP_INCLUDED
#define COMPLEXITY_ESTIMATOR_HPP_INCLUDED
 
#include "../tmp/integer.hpp"
#include "field.hpp"
 
namespace NiHu
{
 
template <class TestField, class TrialField, class KernelEstimator>
class complexity_estimator
{
public:
    typedef TestField test_field_t;
    typedef TrialField trial_field_t;
 
    enum {
        test_field_complexity =
            shape_set_traits::polynomial_order<typename test_field_t::nset_t>::value +
            shape_set_traits::jacobian_order<typename test_field_t::elem_t::lset_t>::value,
        trial_field_complexity =
            shape_set_traits::polynomial_order<typename trial_field_t::nset_t>::value +
            shape_set_traits::jacobian_order<typename trial_field_t::elem_t::lset_t>::value
    };
 
    static unsigned const total_field_complexity = tmp::max_<
        tmp::integer<unsigned, test_field_complexity>,
        tmp::integer<unsigned, trial_field_complexity>
    >::type::value;
 
    static unsigned eval(
        field_base<test_field_t> const &test_field,
        field_base<trial_field_t> const &trial_field
    )
    {
        return total_field_complexity + KernelEstimator::eval(test_field, trial_field);
    }
};
 
 
template <class Estim1, class...Estims>
struct merge_kernel_complexity_estimators :
    merge_kernel_complexity_estimators<
        Estim1,
        typename merge_kernel_complexity_estimators<Estims...>::type
    >
{
};
 
 
template <class Estim1, class Estim2>
struct merge_kernel_complexity_estimators<Estim1, Estim2>
{
    struct type
    {
        template <class test_field_t, class trial_field_t>
        static unsigned eval(
            field_base<test_field_t> const &test_field,
            field_base<trial_field_t> const &trial_field
        )
        {
            return std::max(
                Estim1::eval(test_field, trial_field),
                Estim2::eval(test_field, trial_field)
            );
        }
    };
};
 
}
 
#endif