TYAcousticModelDefaultSolver.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 __TYACOUSTICMODEL_DEFAULTSOLVER__
#define __TYACOUSTICMODEL_DEFAULTSOLVER__
 
#include "Tympan/core/interfaces.h"
#include "Tympan/solvers/9613/TYAcousticModel.h"
#include "Tympan/solvers/9613/DefaultSolver/TYTrajetDefaultSolver.h"
#include "Tympan/solvers/9613/DefaultSolver/TYSolverDefaultSolver.h"
#include "Tympan/solvers/TYSolverDefines.h"
#include "Tympan/models/common/plan.h"
#include <gtest/gtest_prod.h>
 
class TYTrajetDefaultSolver;
class TYSolver;
 
class TYAcousticModelDefaultSolver : public TYAcousticModel
{
public:
    TYAcousticModelDefaultSolver(TYSolverDefaultSolver& solver); //<! Constructor
 
    virtual void compute(const std::deque<TYSIntersection>& tabIntersect, TYTrajetDefaultSolver& trajet,
                         TabPoint3D& ptsTop, TabPoint3D& ptsLeft, TabPoint3D& ptsRight);
 
    void computeWaveLength() override;
 
    virtual bool computeCheminsAvecEcran(const OSegment3D& rayon, const tympan::AcousticSource& source,
                                         const TabPoint3D& pts, const bool vertical,
                                         TYTabCheminDefaultSolver& TabChemins, double distance,
                                         bool conditionFav = false) const;
 
    void computeCheminReflexion(const std::deque<TYSIntersection>& tabIntersect, const OSegment3D& rayon,
                                const tympan::AcousticSource& source, TYTabCheminDefaultSolver& TabChemins,
                                double distance) const;
 
    void computeCheminSansEcran(const OSegment3D& rayon, const tympan::AcousticSource& source,
                                TYTabCheminDefaultSolver& TabChemins, double distance,
                                bool conditionFav = false) const;
 
    void computeCheminAPlat(const OSegment3D& rayon, const tympan::AcousticSource& source,
                            TYTabCheminDefaultSolver& TabChemins, double distance) const;
 
    OSpectre calculC(const double& epaisseur) const;
 
    OSpectre calculAttDiffraction(const OSegment3D& rayon, const OSegment3D& penteMoyenne, const bool& miroir,
                                  const double& re, const double& epaisseur, const bool& vertical,
                                  const bool& avantApres, bool& bDiffOk, bool conditionFav = false) const;
 
    OSpectre limAttDiffraction(const OSpectre& sNC, const OSpectre& C) const;
 
    bool addGroundSteps(const OPoint3D& ptDebut, const OPoint3D& ptFin, const OSegment3D& penteMoyenne,
                        const tympan::AcousticSource& source, const bool& fromSource, const bool& toRecepteur,
                        TYTabEtapeDefaultSolver& Etapes, double& longueur) const;
 
    bool solve(TYTrajetDefaultSolver& trajet);
 
private:
    OSpectreComplex getReflexionSpectrumAt(const OSegment3D& incident, double length,
                                           const OSegment3D& segPente,
                                           const tympan::AcousticSource& source) const;
 
    TYSolver& getSolver() const override
    {
        return _solver;
    }
 
private:
    FRIEND_TEST(test_TYAcousticModel, calculAttDiffraction);
    FRIEND_TEST(test_TYAcousticModel, calculC);
    FRIEND_TEST(test_TYAcousticModel, limAttDiffraction);
    FRIEND_TEST(test_TYAcousticModel, computeCheminReflexion);
 
    OSpectre _lambda;
    OSpectreComplex _absoNulle;
 
    TYSolverDefaultSolver& _solver;
};
 
#endif // __TYACOUSTICMODEL_DEFAULTSOLVER__