iterative_trapz.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_TRAPZ_LINE_QUAD_H
#define ITER_TRAPZ_LINE_QUAD_H
 
#include <functional>
#include <stdexcept>
 
#include "types.hpp"
#include "quadrature/line/base.hpp"
#include "quadrature/line/trapz.hpp"
 
 
namespace bem
{
 
const Float ITER_TRAPZ_LINE_DEFAULT_TOL = 1e-3;
 
template <uint8_t dim>
class IterativeTrapzLineQuadrature: public LineQuadratureBase<dim>
{
 
    using base = LineQuadratureBase<dim>;
 
public:
 
    void set_starting_num_segments(const uint16_t starting_num_segments)
    {
        if (starting_num_segments > TRAPZ_LINE_MAX_NUM_SEGMENTS)
            throw std::domain_error(
                std::string("IterativeTrapzLineQuadrature::set_starting_num_segments(): number of ") +
                std::string("segments must be less than ") +
                std::to_string(TRAPZ_LINE_MAX_NUM_SEGMENTS) + ".");
 
        starting_num_segments_ = starting_num_segments;
        base::points_weights_computed_ = false;
        converged_ = false;
        return;
    };
 
 
    uint16_t starting_num_segments() const { return starting_num_segments_; };
 
 
    void set_tol(const Float tol)
    {
        tol_ = tol;
        base::points_weights_computed_ = false;
        converged_ = false;
        return;
    };
 
 
    Float tol() const { return tol_; };
 
 
    void set_max_iters(const uint16_t max_iters)
    {
        max_iters_ = max_iters;
        base::points_weights_computed_ = false;
        converged_ = false;
        return;
    };
 
 
    uint16_t max_iters() const { return max_iters_; };
 
 
    void compute_points_weights(
        ConstEigRef<EigColVecN<Float, dim>> p1,
        ConstEigRef<EigColVecN<Float, dim>> p2,
        std::function<EigRowVec<Complex> (ConstEigRef<EigMatNX<Float, dim>>)> eval = {}
        ) override;
 
 
    bool converged() const
    {
        if (!base::points_weights_computed_)
            throw std::runtime_error(
                "IterativeTrapzLineQuadrature::converged(): must call `compute_points_weights()` first.");
        return converged_;
    };
 
 
    uint16_t converged_num_segments() const { return converged_num_segments_; };
 
 
private:
 
    bool converged_ = false;
    uint16_t starting_num_segments_ = 1;
    uint16_t converged_num_segments_ = 0;
    uint16_t max_iters_ = LINE_MAX_ORDER;
    Float tol_ = ITER_TRAPZ_LINE_DEFAULT_TOL;
 
};
 
}
 
#ifndef BEM_LINKED
#include "quadrature/line/iterative_trapz.cpp"
#endif
 
#endif