domain.hpp

Go to the documentation of this file.

// This file is a part of NiHu, a C++ BEM template library.
//
// Copyright (C) 2012-2019  Peter Fiala <fiala@hit.bme.hu>
// Copyright (C) 2012-2019  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 DOMAIN_HPP_INCLUDED
#define DOMAIN_HPP_INCLUDED
 
#include "space.hpp"
 
#include <array>
#include <string>
 
namespace NiHu
{
 
namespace domain_traits
{
    template <class Derived>
    struct space_type;
 
    template <class Derived>
    struct num_corners;
 
    template <class Derived>
    struct num_edges;
 
    template <class Derived>
    struct volume;
 
    template <class Derived>
    struct id
    {
        enum {
            value = space_type<Derived>::type::dimension * 10 + num_corners<Derived>::value
        };
    };
 
    template <class Derived>
    struct name
    {
        static std::string const value;
    };
} // end of namespace domain_traits
 
template <class Derived>
class domain_base
{
public:
    typedef Derived type;
    typedef typename domain_traits::space_type<Derived>::type space_t;
    enum {
        num_corners = domain_traits::num_corners<Derived>::value,
        num_edges = domain_traits::num_edges<Derived>::value,
        id = domain_traits::id<Derived>::value,
        dimension = space_t::dimension
    };
    typedef typename space_t::scalar_t scalar_t;
    typedef typename space_t::location_t xi_t;
    typedef std::array<xi_t, num_corners> corners_t;
    typedef unsigned domain_edge_t[2];
    typedef domain_edge_t edges_t[num_edges];
 
    static xi_t const &get_center()
    {
        return Derived::get_center_impl();
    }
 
    static corners_t const &get_corners()
    {
        return Derived::get_corners_impl();
    }
 
    static xi_t const &get_corner(unsigned idx)
    {
        return get_corners()[idx];
    }
 
    static edges_t const &get_edges()
    {
        return Derived::get_edges_impl();
    }
 
    static domain_edge_t const &get_edge(unsigned idx)
    {
        return get_edges()[idx];
    }
 
    static constexpr scalar_t get_volume()
    {
        return domain_traits::volume<Derived>::value;
    }
 
    static std::string const &get_name()
    {
        return domain_traits::name<Derived>::value;
    }
 
    static constexpr unsigned get_id()
    {
        return id;
    }
 
    static xi_t constrain_inside(xi_t const &xi)
    {
        return Derived::constrain_inside_impl(xi);
    }
}; // end of class domain_base
 
} // end of namesapce NiHu
 
#endif