inf_gap.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 BEM_RWG_EXC_INF_GAP_H
#define BEM_RWG_EXC_INF_GAP_H
 
#include <stdexcept>
#include <vector>
 
#include "types.hpp"
#include "rwg/excitations/base.hpp"
 
 
namespace bem
{
// Forward declarations
template <uint8_t dim> class Triangle;
}
 
 
namespace bem::rwg
{
 
const Float INF_GAP_DEFAULT_TOL = 1.0e-6;
 
class InfinitesimalGap: public ExcitationBase<3>
{
public:
 
    InfinitesimalGap(
        const std::vector<EigMatNX<Float, 3>>& segments,
        ConstEigRef<EigRowVec<Complex>> amp,
        const Float tol = INF_GAP_DEFAULT_TOL
        ):
            segments_(segments),
            amp_(amp),
            tol_(tol)
    {
        if (segments_.size() != amp_.cols())
        {
            throw std::invalid_argument(
                "InfinitesimalGap: `segments` and `amp` must have the same number of columns."
                );
        }
        return;
    };
 
 
    Index num_excitations() const override { return amp_.cols(); };
 
 
    EigMatNX<Complex, 3> compute(const Complex k, const Triangle<3>& obs_tri) override;
 
 
private:
 
    const std::vector<EigMatNX<Float, 3>> segments_;
    const EigRowVec<Complex> amp_;
    const Float tol_ = INF_GAP_DEFAULT_TOL;
 
};
 
}
 
#ifndef BEM_LINKED
#include "rwg/excitations/inf_gap.cpp"
#endif
 
#endif