single_layer.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_PROJ_SINGLE_LAYER_H
#define BEM_RWG_PROJ_SINGLE_LAYER_H
 
#include "types.hpp"
#include "geometry/primitives/triangle.hpp"
 
#include "rwg/integrators/src/base.hpp"
#include "rwg/integrators/obs/quadrature.hpp"
 
#include "rwg/projectors/base.hpp"
 
 
namespace bem::rwg
{
 
template <typename SrcIntegratorType = SrcStrategic<>>
class VectorSingleLayerProj: public ProjectorBase<3>
{
 
    static_assert(
        std::is_base_of<SrcIntegratorBase, SrcIntegratorType>::value,
        "VectorSingleLayerProj: `SrcIntegratorType` must derive from `SrcIntegratorBase`"
        );
 
public:
 
    VectorSingleLayerProj(const SrcIntegratorType src_integrator = SrcStrategic<>()):
        src_integrator_(src_integrator) {};
 
 
    EigMatXN<Complex, 3> compute(
        const Complex k,
        ConstEigRef<EigMatNX<Float, 3>> obs_points,
        const Triangle<3>& src_tri
        ) override;
 
 
private:
 
    SrcIntegratorType src_integrator_;
 
};
 
 
template <typename SrcIntegratorType = SrcStrategic<>>
class ScalarSingleLayerProj: public ProjectorBase<1>
{
 
    static_assert(
        std::is_base_of<SrcIntegratorBase, SrcIntegratorType>::value,
        "ScalarSingleLayerProj: `SrcIntegratorType` must derive from `SrcIntegratorBase`"
        );
 
public:
 
    ScalarSingleLayerProj(const SrcIntegratorType src_integrator = SrcStrategic<>()):
        src_integrator_(src_integrator) {};
 
 
    EigMatXN<Complex, 1> compute(
        const Complex k,
        ConstEigRef<EigMatNX<Float, 3>> obs_points,
        const Triangle<3>& src_tri
        ) override;
 
 
private:
 
    SrcIntegratorType src_integrator_;
 
};
 
 
template <typename SrcIntegratorType = SrcStrategic<>>
class GradScalarSingleLayerProj: public ProjectorBase<1>
{
 
    static_assert(
        std::is_base_of<SrcIntegratorBase, SrcIntegratorType>::value,
        "GradScalarSingleLayerProj: `SrcIntegratorType` must derive from `SrcIntegratorBase`"
        );
 
public:
 
    GradScalarSingleLayerProj(const SrcIntegratorType src_integrator = SrcStrategic<>()):
        src_integrator_(src_integrator) {};
 
 
    EigMatXN<Complex, 1> compute(
        const Complex k,
        ConstEigRef<EigMatNX<Float, 3>> obs_points,
        const Triangle<3>& src_tri
        ) override;
 
private:
 
    SrcIntegratorType src_integrator_;
 
};
 
 
template <typename SrcIntegratorType = SrcStrategic<>>
class VectorHypersingularProj: public ProjectorBase<3>
{
 
    static_assert(
        std::is_base_of<SrcIntegratorBase, SrcIntegratorType>::value,
        "VectorHypersingularProj: `SrcIntegratorType` must derive from `SrcIntegratorBase`"
        );
 
public:
 
    VectorHypersingularProj(const SrcIntegratorType src_integrator = SrcStrategic<>()):
        proj_g_(src_integrator), proj_gradg_(src_integrator) {};
 
 
    EigMatXN<Complex, 3> compute(
        const Complex k,
        ConstEigRef<EigMatNX<Float, 3>> obs_points,
        const Triangle<3>& src_tri
        ) override;
 
 
private:
 
    VectorSingleLayerProj<SrcIntegratorType> proj_g_;
    GradScalarSingleLayerProj<SrcIntegratorType> proj_gradg_;
 
};
 
}
 
#include "rwg/projectors/single_layer.tpp"
 
#endif