base.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_INTEQ_BASE_H
#define BEM_RWG_INTEQ_BASE_H
 
#include "types.hpp"
#include "geometry/mesh/triangle_mesh.hpp"
 
#include "matrix/base.hpp"
#include "matrix/eigen_dense.hpp"
 
#include "rwg/assemblers/operator_matrix.hpp"
 
 
namespace bem::rwg
{
 
template <typename MatrixType = EigenDenseMatrix<Complex>>
class IntegralEquationBase
{
 
static_assert(
    std::is_base_of<MatrixBase<Complex>, MatrixType>::value,
    "IntegralEquationBase: MatrixType must be derived from MatrixBase<Complex>."
    );
 
public:
 
    IntegralEquationBase(
        const TriangleMesh<3>& obs_mesh,
        const TriangleMesh<3>& src_mesh,
        ConstEigRef<EigMatNX<Index, 2>> elem_pairs,
        const bool flip_normals = false
        ):
            obs_mesh_(obs_mesh),
            src_mesh_(src_mesh),
            elem_pairs_(elem_pairs),
            flip_normals_(flip_normals) {};
 
 
    IntegralEquationBase(
        const TriangleMesh<3>& obs_mesh,
        const TriangleMesh<3>& src_mesh,
        const bool flip_normals = false
        ):
            obs_mesh_(obs_mesh),
            src_mesh_(src_mesh),
            elem_pairs_(OperatorAssemblerBase<3, 3>::make_pairs(obs_mesh, src_mesh)),
            flip_normals_(flip_normals) {};
 
 
    IntegralEquationBase(
        const TriangleMesh<3>& mesh,
        ConstEigRef<EigMatNX<Index, 2>> elem_pairs,
        const bool flip_normals = false
        ):
            obs_mesh_(mesh),
            src_mesh_(mesh),
            elem_pairs_(elem_pairs),
            flip_normals_(flip_normals) {};
 
 
    IntegralEquationBase(
        const TriangleMesh<3>& mesh,
        const bool flip_normals = false
        ):
            obs_mesh_(mesh),
            src_mesh_(mesh),
            elem_pairs_(OperatorAssemblerBase<3, 3>::make_pairs(mesh, mesh)),
            flip_normals_(flip_normals) {};
 
 
protected:
 
    const TriangleMesh<3>& obs_mesh_;
    const TriangleMesh<3>& src_mesh_;
    const EigMatNX<Index, 2> elem_pairs_;
    const bool flip_normals_ = false;
 
};
 
}
 
#endif