TYReflectionPathFinder.h

Go to the documentation of this file.

/*
 * Copyright (C) <2012-2024> <EDF-DTG> <FRANCE>
 * This file is part of Code_TYMPAN (R).
 * Code_TYMPAN (R) 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.
 * Code_TYMPAN (R) is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 * See the GNU General Public License for more details.
 * You should have received a copy of the GNU General Public License along
 * with Code_TYMPAN (R). If not, see <https://www.gnu.org/licenses/>.
 */
 
#ifndef __TY_REFLECTION_PATH_FINDER__
#define __TY_REFLECTION_PATH_FINDER__
 
#include "Tympan/solvers/TYSolverDefines.h"
#include "Tympan/geometric_methods/AcousticRaytracer/Engine/Simulation.h"
 
#include <gtest/gtest_prod.h>
 
struct TYReflectionPath
{
    OPoint3D source;
    OPoint3D receptor;
    std::deque<OPoint3D> reflectionPoints;
    std::vector<const TYSIntersection*> reflectingSegments;
};
 
class TYAbstractReflectionPathFinder
{
public:
    virtual ~TYAbstractReflectionPathFinder() = default;
 
    virtual std::unique_ptr<TYAbstractReflectionPathFinder> instanciate() const = 0;
 
    virtual bool setup(size_t maxOrder, const std::deque<TYSIntersection>& tabIntersect,
                       const OPoint3D& source, const OPoint3D& receptor) = 0;
 
    virtual bool remainingPaths() const = 0;
 
    virtual TYReflectionPath popPath() = 0;
 
protected:
    TYAbstractReflectionPathFinder() = default;
};
 
class TYReflectionPathFinderART : public TYAbstractReflectionPathFinder
{
public:
    TYReflectionPathFinderART(unsigned int nbRays, double rcptRadius);
 
    std::unique_ptr<TYAbstractReflectionPathFinder> instanciate() const override;
 
    bool setup(size_t maxOrder, const std::deque<TYSIntersection>& tabIntersect, const OPoint3D& source,
               const OPoint3D& receptor) override;
 
    bool remainingPaths() const override;
 
    TYReflectionPath popPath() override;
 
protected:
    TYReflectionPathFinderART() = default;
 
    void setupConfiguration(size_t maxOrder);
 
    bool setupScene(const std::deque<TYSIntersection>& tabIntersect, const OPoint3D& sourcePoint,
                    const OPoint3D& receptorPoint);
 
    void setupAcceleratingDataStructure();
 
    void setupSolver();
 
    void computeRaytracing();
 
    void buildReflectionPaths(const OPoint3D& source, const OPoint3D& receptor);
 
    bool processIntersectingSegment(const TYSIntersection& inter);
 
    bool processInfraIntersectingSegment(const TYSIntersection& inter);
 
    bool processTopoIntersectingSegment(const TYSIntersection& inter);
 
    void buildQuadranglePointsFromInter(const TYSIntersection& inter, OPoint3D& p1Top, OPoint3D& p1Bottom,
                                        OPoint3D& p2Top, OPoint3D& p2Bottom) const;
 
    void addOrientedQuadrangle(const OPoint3D& p1Top, const OPoint3D& p1Bottom, const OPoint3D& p2Top,
                               const OPoint3D& p2Bottom, const OVector3D vExt,
                               const TYSIntersection& sourceInter);
 
    Shape* addOrientedTriangle(const OPoint3D& p1, const OPoint3D& p2, const OPoint3D& p3,
                               const OVector3D& vExt);
 
    OVector3D computeELPlaneNormal(const OPoint3D& source, const OPoint3D& receptor);
 
    bool isPathValid(const TYReflectionPath& path);
 
    bool isPathAlreadyFound(const TYReflectionPath& path) const;
 
    bool isLastSegmentIntersectingSceneElement(const TYReflectionPath& path);
 
    Simulation _ARTSimulation;
    std::unique_ptr<Solver> _pSolver;
    std::set<std::vector<const TYSIntersection*>> _foundSequences;
    std::deque<TYReflectionPath> _paths;
    std::unordered_map<const Shape*, const TYSIntersection*> _shapeToTYSIntersection;
    unsigned int _nbRays;
    double _rcptRadius;
    double _extrusionLength;
    OVector3D _ELnormal;
 
    FRIEND_TEST(test_ReflectionPathFinderART, computeELPlaneNormal_horizontalPlane1);
    FRIEND_TEST(test_ReflectionPathFinderART, computeELPlaneNormal_horizontalPlane2);
    FRIEND_TEST(test_ReflectionPathFinderART, computeELPlaneNormal_nonHorizontalPlane);
};
 
#endif // __TY_REFLECTION_PATH_FINDER__