operator_matrix.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_OP_ASSEMBLER_H
#define BEM_RWG_OP_ASSEMBLER_H
 
#include "types.hpp"
#include "rwg/assemblers/base.hpp"
 
 
namespace bem
{
// Forward declarations
template <typename T> class MatrixBase;
}
 
 
namespace bem::rwg
{
 
class EdgeOperatorAssembler: public OperatorAssemblerBase<3, 3>
{
 
    using base = OperatorAssemblerBase<3, 3>;
    using base::base;
 
public:
 
    void prep_matrix(MatrixBase<Complex>& mat) override;
 
 
    void fill_matrix(
        MatrixBase<Complex>& mat,
        ConstEigRef<EigColVecN<Index, 2>> elem_pair,
        ConstEigRef<EigMatMN<Complex, 3, 3>> values
        ) override;
 
};
 
 
class FaceOperatorAssembler: public OperatorAssemblerBase<1, 1>
{
 
    using base = OperatorAssemblerBase<1, 1>;
    using base::base;
 
public:
 
    void prep_matrix(MatrixBase<Complex>& mat) override;
 
 
    void fill_matrix(
        MatrixBase<Complex>& mat,
        ConstEigRef<EigColVecN<Index, 2>> elem_pair,
        ConstEigRef<EigMatMN<Complex, 1, 1>> values
        ) override;
 
};
 
 
class FaceEdgeOperatorAssembler: public OperatorAssemblerBase<1, 3>
{
 
    using base = OperatorAssemblerBase<1, 3>;
    using base::base;
 
public:
 
    void prep_matrix(MatrixBase<Complex>& mat) override;
 
 
    void fill_matrix(
        MatrixBase<Complex>& mat,
        ConstEigRef<EigColVecN<Index, 2>> elem_pair,
        ConstEigRef<EigMatMN<Complex, 1, 3>> values
        ) override;
 
};
 
 
class EdgeFaceOperatorAssembler: public OperatorAssemblerBase<3, 1>
{
 
    using base = OperatorAssemblerBase<3, 1>;
    using base::base;
 
public:
 
    void prep_matrix(MatrixBase<Complex>& mat) override;
 
 
    void fill_matrix(
        MatrixBase<Complex>& mat,
        ConstEigRef<EigColVecN<Index, 2>> elem_pair,
        ConstEigRef<EigMatMN<Complex, 3, 1>> values
        ) override;
 
};
 
 
class VectorOperatorsAssembler: public OperatorAssemblerBase<12, 3>
{
 
    using base = OperatorAssemblerBase<12, 3>;
    using base::base;
 
public:
 
    void prep_matrix(MatrixBase<Complex>& mat) override;
 
 
    void fill_matrix(
        MatrixBase<Complex>& mat,
        ConstEigRef<EigColVecN<Index, 2>> elem_pair,
        ConstEigRef<EigMatMN<Complex, 12, 3>> values
        ) override;
 
};
 
}
 
#ifndef BEM_LINKED
#include "rwg/assemblers/operator_matrix.cpp"
#endif
 
#endif