quadrature.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_OPINT_OBS_QUAD_H
#define BEM_RWG_OPINT_OBS_QUAD_H
 
#include "types.hpp"
#include "geometry/primitives/triangle.hpp"
#include "quadrature/triangle/base.hpp"
#include "quadrature/triangle/gauss.hpp"
 
#include "rwg/integrators/obs/base.hpp"
#include "rwg/integrators/src/base.hpp"
#include "rwg/integrators/src/strategic.hpp"
 
 
namespace bem::rwg
{
 
template <typename TriangleQuadratureType = GaussTriangleQuadrature<3>, typename SrcIntegratorType = SrcStrategic<>>
class ObsQuadrature: public ObsIntegratorBase
{
 
    using base = ObsIntegratorBase;
    static_assert(
        std::is_base_of<TriangleQuadratureBase<3>, TriangleQuadratureType>::value,
        "ObsQuadrature: `TriangleQuadratureType` must derive from `TriangleQuadratureBase<3>`"
        );
    static_assert(
        std::is_base_of<SrcIntegratorBase, SrcIntegratorType>::value,
        "ObsQuadrature: `SrcIntegratorType` must derive from `SrcIntegratorBase`"
        );
 
public:
 
    ObsQuadrature(
        const TriangleQuadratureType tri_quad = GaussTriangleQuadrature<3>(),
        const SrcIntegratorType src_integrator = SrcStrategic<>()
        ): tri_quad_(tri_quad), src_integrator_(src_integrator) {};
 
 
    ObsResult integrate(
        const Complex k,
        const Triangle<3>& obs_tri,
        const Triangle<2>& src_tri
        ) override;
 
 
    const TriangleQuadratureType& quadrature_object() const
    { return tri_quad_; };
 
 
    TriangleQuadratureType& quadrature_object()
    { return tri_quad_; };
 
 
    const SrcIntegratorType& src_integrator() const
    { return src_integrator_; };
 
 
    SrcIntegratorType& src_integrator()
    { return src_integrator_; };
 
 
private:
 
    Complex g_term(const SrcResult& src_result);
 
 
    EigRowVecN<Complex, 12> rs_g_terms(const SrcResult& src_result);
 
 
    EigRowVecN<Complex, 9> grad_g_terms(const SrcResult& src_result);
 
 
    EigRowVecN<Complex, 15> rot_grad_g_terms(const SrcResult& src_result);
 
 
    TriangleQuadratureType tri_quad_;
    SrcIntegratorType src_integrator_;
 
};
 
}
 
#include "rwg/integrators/obs/quadrature.tpp"
 
#endif