incidence.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_OPS_INCIDENCE_H
#define BEM_RWG_OPS_INCIDENCE_H
 
#include "types.hpp"
#include "geometry/operations.hpp"
#include "rwg/operators/base.hpp"
 
 
namespace bem::rwg
{
 
class DivRwgOp: public OperatorBase<1, 3>
{
public:
 
    EigMatMN<Complex, 1, 3> compute(
        const Complex k,
        const Triangle<3>& obs_tri,
        const Triangle<3>& src_tri
        ) override
    {
        EigMatMN<Complex, 1, 3> result = EigMatMN<Complex, 1, 3>::Zero(1, 3);
 
        if (GeometryOps<3>::common_vertices(obs_tri, src_tri) < 3)
            return result;
 
        result = src_tri.edge_polarities();
        return result;
    };
 
};
 
}
 
#endif