lib_shape.cpp

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/>.
 
#include "lib_shape.hpp"
 
namespace NiHu
{
 
namespace shape_set_traits
{
    template <>
    const std::string name<line_0_shape_set>::value = "Line 0 shape set";
 
    template <>
    const std::string name<line_1_shape_set>::value = "Line 1 shape set";
 
    template <>
    const std::string name<line_2_shape_set>::value = "Line 2 shape set";
 
    template <>
    const std::string name<tria_0_shape_set>::value = "Tria 0 shape set";
 
    template <>
    const std::string name<tria_1_shape_set>::value = "Tria 1 shape set";
 
    template <>
    const std::string name<tria_2_shape_set>::value = "Tria 2 shape set";
 
    template <>
    const std::string name<quad_0_shape_set>::value = "Quad 0 shape set";
 
    template <>
    const std::string name<quad_1_shape_set>::value = "Quad 1 shape set";
 
    template <>
    const std::string name<quad_2_shape_set>::value = "Quad 2 shape set";
 
    template <>
    const std::string name<quad_28_shape_set>::value = "Quad 28 shape set";
 
    template <>
    const std::string name<brick_0_shape_set>::value = "Brick 0 shape set";
 
    template <>
    const std::string name<brick_1_shape_set>::value = "Brick 1 shape set";
 
    template <>
    const std::string name<line_1_gauss_shape_set>::value = "Line 1 Gauss shape set";
 
    template <>
    const std::string name<quad_1_gauss_shape_set>::value = "Quad 1 Gauss shape set";
 
    template <>
    const std::string name<tria_1_gauss_shape_set>::value = "Tria 1 Gauss shape set";
}
 
 
line_2_shape_set::corners_t
const line_2_shape_set::m_corners = {
    line_2_shape_set::xi_t::Constant(-1.0),
    line_2_shape_set::xi_t::Constant(0.0),
    line_2_shape_set::xi_t::Constant(1.0),
};
 
 
tria_2_shape_set::corners_t
const tria_2_shape_set::m_corners = {
    tria_2_shape_set::xi_t(0.0, 0.0),
    tria_2_shape_set::xi_t(0.5, 0.0),
    tria_2_shape_set::xi_t(1.0, 0.0),
    tria_2_shape_set::xi_t(0.5, 0.5),
    tria_2_shape_set::xi_t(0.0, 1.0),
    tria_2_shape_set::xi_t(0.0, 0.5)
};
 
 
quad_2_shape_set::corners_t
const quad_2_shape_set::m_corners = {
    quad_2_shape_set::xi_t(-1.0, -1.0),
    quad_2_shape_set::xi_t(0.0, -1.0),
    quad_2_shape_set::xi_t(+1.0, -1.0),
    quad_2_shape_set::xi_t(+1.0, 0.0),
    quad_2_shape_set::xi_t(+1.0, +1.0),
    quad_2_shape_set::xi_t(0.0, +1.0),
    quad_2_shape_set::xi_t(-1.0, +1.0),
    quad_2_shape_set::xi_t(-1.0, 0.0),
    quad_2_shape_set::xi_t(0.0, 0.0)
};
 
 
quad_28_shape_set::corners_t
const quad_28_shape_set::m_corners = {
    quad_28_shape_set::xi_t(-1.0, -1.0),
    quad_28_shape_set::xi_t(0.0, -1.0),
    quad_28_shape_set::xi_t(+1.0, -1.0),
    quad_28_shape_set::xi_t(+1.0, 0.0),
    quad_28_shape_set::xi_t(+1.0, +1.0),
    quad_28_shape_set::xi_t(0.0, +1.0),
    quad_28_shape_set::xi_t(-1.0, +1.0),
    quad_28_shape_set::xi_t(-1.0, 0.0)
};
 
 
line_1_gauss_shape_set::corners_t
    const line_1_gauss_shape_set::m_corners = {
        line_1_gauss_shape_set::xi_t::Constant(-std::sqrt(3.0)/3.0),
        line_1_gauss_shape_set::xi_t::Constant(+std::sqrt(3.0)/3.0)
};
 
quad_1_gauss_shape_set::corners_t
    const quad_1_gauss_shape_set::m_corners = {
        quad_1_gauss_shape_set::xi_t(-std::sqrt(3.0)/3.0, -std::sqrt(3.0)/3.0),
        quad_1_gauss_shape_set::xi_t(+std::sqrt(3.0)/3.0, -std::sqrt(3.0)/3.0),
        quad_1_gauss_shape_set::xi_t(-std::sqrt(3.0)/3.0, +std::sqrt(3.0)/3.0),
        quad_1_gauss_shape_set::xi_t(+std::sqrt(3.0)/3.0, +std::sqrt(3.0)/3.0)
};
 
tria_1_gauss_shape_set::corners_t
    const tria_1_gauss_shape_set::m_corners = {
        tria_1_gauss_shape_set::xi_t(1.0/6.0, 1.0/6.0),
        tria_1_gauss_shape_set::xi_t(4.0/6.0, 1.0/6.0),
        tria_1_gauss_shape_set::xi_t(1.0/6.0, 4.0/6.0)
};
 
 
} // end of namespace NiHu