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_OPINT_SRC_BASE_H
#define BEM_RWG_OPINT_SRC_BASE_H
 
#include "types.hpp"
#include "geometry/primitives/triangle.hpp"
 
 
namespace bem::rwg
{
 
struct SrcResult
{
    EigRowVec<Complex> g;
    EigMatNX<Complex, 2> rs_g;
    EigMatNX<Complex, 3> grad_g;
    // EigMatNX<Complex, 3> xs_grad_g, ys_grad_g;
    void resize(Index n) {
        g.resize(1, n); rs_g.resize(2, n); grad_g.resize(3, n);
        g.setZero(); rs_g.setZero(); grad_g.setZero();
        return;
    };
};
 
 
class SrcIntegratorBase
{
public:
 
    virtual SrcResult integrate(
        const Complex k,
        const Triangle<2>& src_tri,
        ConstEigRef<EigMatNX<Float, 3>> r_obs
        ) = 0;
 
 
    virtual void set_compute_terms(
        bool compute_g_terms,
        bool compute_grad_g_terms
        )
    {
        compute_g_terms_ = compute_g_terms;
        compute_grad_g_terms_ = compute_grad_g_terms;
        return;
    };
 
 
    virtual SrcResult zeros(
        const Complex k,
        const Triangle<2>& src_tri,
        ConstEigRef<EigMatNX<Float, 3>> r_obs
        ) { SrcResult result; result.resize(r_obs.cols()); return result; };
 
 
protected:
 
    bool compute_g_terms_ = true;
    bool compute_grad_g_terms_ = true;
 
};
 
}
 
#endif