function_space.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_FUNCSP_H
#define BEM_RWG_FUNCSP_H
 
#include <functional>
 
#include "types.hpp"
#include "geometry/primitives/triangle.hpp"
#include "quadrature/triangle/base.hpp"
 
 
namespace bem
{
// Forward declarations
template <typename T> class MatrixBase;
template <uint8_t dim> class TriangleMesh;
}
 
 
namespace bem::rwg
{
 
class Rwg
{
public:
 
    static EigRowVecN<Float, 3> normalization(const Triangle<3>& tri)
    { return tri.edge_polarities() / tri.area() / two; };
 
 
    static EigMatNX<Float, 3> value(
        const Triangle<3>& tri,
        uint8_t edge,
        ConstEigRef<EigMatNX<Float, 3>> points,
        const bool rotated = false
        );
 
 
    static EigColVecN<Complex, 3> test_field(
        const Triangle<3>& tri,
        std::function<EigMatNX<Complex, 3> (ConstEigRef<EigMatNX<Float, 3>>)> field_eval,
        TriangleQuadratureBase<3>& tri_quad,
        const bool rotated = false
        );
 
 
    static EigMatNX<Complex, 3> reconstruct_field(
        const TriangleMesh<3>& mesh,
        const MatrixBase<Complex>& coeffs,
        ConstEigRef<EigMatNX<Float, 3>> points,
        const bool rotated = false
        );
 
};
 
 
class Pulse
{
public:
 
    static Float normalization(const Triangle<3>& tri)
    { return one / tri.area(); };
 
 
    static EigRowVec<Float> value(const Triangle<3>& tri, ConstEigRef<EigMatNX<Float, 3>> points)
    {
        return EigRowVec<Float>::Constant(1, points.cols(), normalization(tri));
    }
 
 
    static Complex test_field(
        const Triangle<3>& tri,
        std::function<EigRowVec<Complex> (ConstEigRef<EigMatNX<Float, 3>>)> field_eval,
        TriangleQuadratureBase<3>& tri_quad
        );
 
 
    static EigRowVec<Complex> reconstruct_field(
        const TriangleMesh<3>& mesh,
        const MatrixBase<Complex>& coeffs,
        ConstEigRef<EigMatNX<Float, 3>> points
        );
 
};
 
}
 
#ifndef BEM_LINKED
#include "rwg/function_space.cpp"
#endif
 
#endif