gauss.hpp

Go to the documentation of this file.

// OpenBEM - Copyright (C) 2026 Shashwat Sharma
 
// This file is part of OpenBEM.
 
// OpenBEM 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.
 
// You should have received a copy of the GNU General Public License along with OpenBEM.
// If not, see <https://www.gnu.org/licenses/>.
 
 
#ifndef GAUSS_TRI_QUAD_H
#define GAUSS_TRI_QUAD_H
 
#include <functional>
 
#include "types.hpp"
#include "constants.hpp"
#include "geometry/primitives/triangle.hpp"
#include "quadrature/utility.hpp"
#include "quadrature/triangle/base.hpp"
 
 
namespace bem
{
 
template <uint8_t dim>
class GaussTriangleQuadrature: public TriangleQuadratureBase<dim>
{
 
    using base = TriangleQuadratureBase<dim>;
 
public:
 
    GaussTriangleQuadrature(const uint8_t order = TRI_DEFAULT_ORDER)
    {
        EigColVec<Index> orders = EigColVec<Index>::LinSpaced(TRI_MAX_ORDER, 1, TRI_MAX_ORDER);
        rules_ = load_rules<2> (rule_file_, orders);
        set_order(order);
        return;
    };
 
 
    void set_order(const uint8_t order) override;
 
 
    void compute_points_weights(
        const Triangle<dim>& tri,
        std::function<EigRowVec<Complex> (ConstEigRef<EigMatNX<Float, dim>>)> eval = {}
        ) override;
 
 
    const EigMatNX<Float, 2>& ref_points() const { return rules_[base::order_ - 1].nodes;};
 
 
    const EigRowVec<Float>& ref_weights() const { return rules_[base::order_ - 1].weights; };
 
 
private:
 
    std::string rule_file_ = "tri_quad_xiao_gimbutas.json";
    std::vector<QuadratureRule<2>> rules_;
 
};
 
}
 
#ifndef BEM_LINKED
#include "quadrature/triangle/gauss.cpp"
#endif
 
#endif