iterative_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 ITER_GAUSS_TRI_QUAD_H
#define ITER_GAUSS_TRI_QUAD_H
 
#include <functional>
 
#include "types.hpp"
#include "geometry/primitives/triangle.hpp"
#include "quadrature/triangle/base.hpp"
#include "quadrature/triangle/gauss.hpp"
 
 
namespace bem
{
 
const Float ITER_GAUSS_TRI_DEFAULT_TOL = 1e-3;
 
template <uint8_t dim>
class IterativeGaussTriangleQuadrature: public TriangleQuadratureBase<dim>
{
 
    using base = TriangleQuadratureBase<dim>;
 
public:
 
    void set_starting_order(const uint8_t starting_order)
    {
        starting_order_ = starting_order;
        converged_ = false;
        return;
    };
 
 
    uint8_t starting_order() const { return starting_order_; };
 
 
    void set_tol(const Float tol)
    {
        tol_ = tol;
        converged_ = false;
        return;
    };
 
 
    Float tol() const { return tol_; };
 
 
    void set_max_iters(const uint16_t max_iters)
    {
        max_iters_ = max_iters;
        converged_ = false;
        return;
    };
 
 
    uint16_t max_iters() const { return max_iters_; };
 
 
    void compute_points_weights(
        const Triangle<dim>& tri,
        std::function<EigRowVec<Complex> (ConstEigRef<EigMatNX<Float, dim>>)> eval = {}
        ) override;
 
 
    bool converged() const
    {
        if (!base::points_weights_computed_)
            throw std::runtime_error(
                "IterativeGaussTriangleQuadrature::converged(): must call `compute_points_weights()` first.");
        return converged_;
    };
 
 
    uint8_t converged_order() const
    {
        if (!base::points_weights_computed_)
            throw std::runtime_error(
                "IterativeGaussTriangleQuadrature::converged_order(): must call `compute_points_weights()` first.");
        return converged_order_;
    };
 
 
private:
 
    bool converged_ = false;
    uint8_t starting_order_ = 1;
    uint8_t converged_order_ = 0;
    uint16_t max_iters_ = TRI_MAX_ORDER;
    Float tol_ = ITER_GAUSS_TRI_DEFAULT_TOL;
    GaussTriangleQuadrature<dim> gauss_quad_ = GaussTriangleQuadrature<dim> (starting_order_);
 
};
 
}
 
#ifndef BEM_LINKED
#include "quadrature/triangle/iterative_gauss.cpp"
#endif
 
#endif