location_normal.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 LOCATION_NORMAL_HPP_INCLUDED
#define LOCATION_NORMAL_HPP_INCLUDED
 
#include "../core/element.hpp"
 
namespace NiHu
{
 
template <class Space>
class location_input
{
public:
    typedef Space space_t;
    typedef typename space_t::location_t x_t;
    typedef typename space_t::scalar_t scalar_t;
 
    template <class elem_t>
    location_input(elem_t const &elem, typename elem_t::xi_t const &xi) :
        m_x(elem.get_x(xi))
    {
    }
 
    location_input(x_t const &x)
        : m_x(x)
    {
    }
 
    x_t const &get_x(void) const
    {
        return m_x;
    }
 
private:
    x_t m_x;
};
 
 
template <class Space>
class location_normal_jacobian_input
    : public location_input<Space>
{
public:
    typedef Space space_t;
    typedef typename space_t::location_t x_t;
    typedef typename space_t::scalar_t scalar_t;
 
    template <class elem_t>
    location_normal_jacobian_input(elem_t const &elem, typename elem_t::xi_t const &xi)
        : location_input<Space>(elem, xi)
        , m_jac_vector(elem.get_normal(xi))
        , m_jac(m_jac_vector.norm())
        , m_unit_normal(m_jac_vector / m_jac)
    {
    }
 
    location_normal_jacobian_input(x_t const &x, x_t const &jac_vector)
        : location_input<Space>(x)
        , m_jac_vector(jac_vector)
        , m_jac(m_jac_vector.norm())
        , m_unit_normal(m_jac_vector / m_jac)
    {
    }
 
    x_t const &get_unit_normal(void) const
    {
        return m_unit_normal;
    }
 
    scalar_t const &get_jacobian(void) const
    {
        return m_jac;
    }
    
    x_t const &get_jacobian_vector() const
    {
        return m_jac_vector;
    }
 
private:
    x_t m_jac_vector;
    scalar_t m_jac;
    x_t m_unit_normal;
};
 
 
template <class Space>
class location_volume_jacobian_input
    : public location_input<Space>
{
public:
    typedef Space space_t;
    typedef typename space_t::scalar_t scalar_t;
    
    typedef typename space_t::location_t location_t;
 
    template <class elem_t>
    location_volume_jacobian_input(elem_t const &elem, typename elem_t::xi_t const &xi) :
        location_input<Space>(elem, xi),
        m_jac(elem.get_dx(xi).determinant())
    {
    }
 
    location_volume_jacobian_input(location_t const &x, scalar_t const &jac) :
        location_input<Space>(x),
        m_jac(jac)
    {
    }
 
    scalar_t const &get_jacobian(void) const
    {
        return m_jac;
    }
 
private:
    scalar_t m_jac;
};
 
 
typedef location_input<space_1d<> > location_input_1d;
typedef location_input<space_2d<> > location_input_2d;
typedef location_input<space_3d<> > location_input_3d;
 
typedef location_normal_jacobian_input<space_2d<> > location_normal_input_2d;
typedef location_normal_jacobian_input<space_3d<> > location_normal_input_3d;
 
 
} // end of namespace NiHu
 
#endif // LOCATION_NORMAL_HPP_INCLUDED