continuity.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_CONT_H
#define BEM_RWG_INTEQ_CONT_H
 
#include "types.hpp"
#include "constants.hpp"
#include "materials.hpp"
 
#include "matrix/base.hpp"
#include "matrix/eigen_dense.hpp"
 
#include "rwg/integral_equations/base.hpp"
#include "rwg/operators/incidence.hpp"
#include "rwg/assemblers/operator_matrix.hpp"
 
 
namespace bem::rwg
{
 
template <typename MatrixType = EigenDenseMatrix<Complex>>
class Continuity: public IntegralEquationBase<MatrixType>
{
 
    using base = IntegralEquationBase<MatrixType>;
    using base::base;
 
public:
 
    MatrixType j_matrix(
        const Float f,
        const Material& material
        )
    {
        DivRwgOp op_D;
        MatrixType D;
        assm_div_.assemble(D, op_D, 0);
        return D;
    }
 
 
    MatrixType rho_matrix(
        const Float f,
        const Material& material
        )
    {
        MatrixType I (base::obs_mesh_.num_elems(), base::src_mesh_.num_elems());
        I.set_identity();
        I.scale(J * two_pi * f);
        return I;
    }
 
 
protected:
 
    FaceEdgeOperatorAssembler assm_div_ = FaceEdgeOperatorAssembler(
        base::obs_mesh_, base::src_mesh_, base::elem_pairs_
        );
 
};
 
}
 
#endif