operations.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_GEOM_OPS_H
#define BEM_GEOM_OPS_H
 
#include <cassert>
 
#include "types.hpp"
#include "constants.hpp"
 
 
namespace bem
{
 
const Float GEOMETRY_DEFAULT_TOL = std::sqrt(float_eps);
 
// Forward declarations
template <uint8_t dim> class Edge;
template <uint8_t dim> class Triangle;
 
template <uint8_t dim>
class GeometryOps
{
 
    static_assert((dim == 1 || dim == 2 || dim == 3), "Dimension must be one, two, or three.");
 
public:
 
    static EigMatNX<Float, dim> polar_to_cartesian(ConstEigRef<EigMatNX<Float, dim>> points);
 
 
    static EigMatNX<Complex, dim> cartesian_to_polar_field(
        ConstEigRef<EigMatNX<Float, dim>> points,
        ConstEigRef<EigMatNX<Complex, dim>> field
        );
 
 
    static EigMatNX<Float, dim> transform_coordinate_system(
        ConstEigRef<EigMatNX<Float, dim>> v_in,
        ConstEigRef<EigColVecN<Float, dim>> new_origin,
        ConstEigRef<EigMatMN<Float, dim, dim>> new_uvw
        );
 
 
    static EigRowVec<Float> angle_between_vectors(
        ConstEigRef<EigMatNX<Float, dim>> v1,
        ConstEigRef<EigMatNX<Float, dim>> v2
        );
 
 
    static EigRowVec<Float> directed_angle_between_vectors(
        ConstEigRef<EigMatNX<Float, dim>> v1,
        ConstEigRef<EigMatNX<Float, dim>> v2
        );
 
 
    static uint8_t common_vertices(
        const Triangle<dim>& triangle1,
        const Triangle<dim>& triangle2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static uint8_t common_vertices(
        EigRef<EigColVecN<Index, 3>> idx1,
        EigRef<EigColVecN<Index, 3>> idx2,
        const Triangle<dim>& triangle1,
        const Triangle<dim>& triangle2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static int8_t check_parallel_triangles(
        const Triangle<3>& triangle1,
        const Triangle<3>& triangle2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static bool check_perpendicular_triangles(
        const Triangle<3>& triangle1,
        const Triangle<3>& triangle2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static bool check_coplanar_triangles(
        const Triangle<3>& triangle1,
        const Triangle<3>& triangle2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static int8_t check_parallel_edges(
        const Edge<dim>& edge1,
        const Edge<dim>& edge2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static bool check_perpendicular_edges(
        const Edge<dim>& edge1,
        const Edge<dim>& edge2,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
 
    static bool point_in_polygon(
        ConstEigRef<EigColVecN<Float, dim>>& point,
        ConstEigRef<EigMatNX<Float, dim>>& polygon,
        const Float tol = GEOMETRY_DEFAULT_TOL
        );
 
};
 
}
 
#ifndef BEM_LINKED
#include "geometry/operations.cpp"
#endif
 
#endif