TYAcousticModel.h

Go to the documentation of this file.

/*
 * Copyright (C) <2012-2014> <EDF-R&D> <FRANCE>
 * This program 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 2 of the License, or
 * (at your option) any later version.
 * This program 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 this program; if not, write to the Free Software Foundation, Inc.,
 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 */
 
#ifndef __TYACOUSTICMODEL__
#define __TYACOUSTICMODEL__
 
#include "Tympan/core/interfaces.h"
#include "Tympan/solvers/DefaultSolver/TYChemin.h"
#include "Tympan/solvers/DefaultSolver/TYSolverDefines.h"
#include <gtest/gtest_prod.h>
 
class TYTrajet;
class TYSolver;
 
class TYAcousticModel : public AcousticModelInterface
{
public:
    TYAcousticModel(TYSolver& solver); //<! Constructor
    virtual ~TYAcousticModel();        //<! Destructor
 
    virtual void compute(const std::deque<TYSIntersection>& tabIntersect, TYTrajet& trajet,
                         TabPoint3D& ptsTop, TabPoint3D& ptsLeft, TabPoint3D& ptsRight);
 
    void init();
 
    virtual bool computeCheminsAvecEcran(const OSegment3D& rayon, const tympan::AcousticSource& source,
                                         const TabPoint3D& pts, const bool vertical, TYTabChemin& TabChemins,
                                         double distance, bool conditionFav = false) const;
 
    void computeCheminReflexion(const std::deque<TYSIntersection>& tabIntersect, const OSegment3D& rayon,
                                const tympan::AcousticSource& source, TYTabChemin& TabChemins,
                                double distance) const;
 
    // * \fn getSpectreAttVegetation
    // * @brief Calcule les spectres d'attenuation caracteristiques de la vegetation
    // * pour chaque terrain survole par rayon.
    // *
    // * @param rayon Un segment representant le rayon acoustique.
    // * @param tabSegSp Le tableau dans lequel sont stockes les spectres d'attenuation et segments associes.
    // * @param penteMoyenne La pente moyenne sur le rayon.
    // */
    // void getSpectreAttVegetation(const OSegment3D& rayon, const OSegment3D& penteMoyenne,
    // TYTabStructSegLPSp& tabSegSp) const;
 
    // * \fn computeCheminAvecVeg
    // * @brief Calcule le chemin direct en tenant compte de la vegetation.
    // *
    // * @param rayon Un segment representant le rayon acoustique.
    // * @param Atmo L'atmosphere courante.
    // * @param pSrcGeoNode La source ponctuelle.
    // * @param penteMoyenne La pente moyenne sur le rayon.
    // *
    // * @return Le chemin determine.
    // */
    // void computeCheminAvecVeg(const OSegment3D& rayon, const tympan::AcousticSource& source, const
    // OSegment3D& penteMoyenne, TYTabChemin& TabChemin, double distance) const;
 
    void computeCheminSansEcran(const OSegment3D& rayon, const tympan::AcousticSource& source,
                                TYTabChemin& TabChemins, double distance, bool conditionFav = false) const;
 
    void computeCheminAPlat(const OSegment3D& rayon, const tympan::AcousticSource& source,
                            TYTabChemin& 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 addEtapesSol(const OPoint3D& ptDebut, const OPoint3D& ptFin, const OSegment3D& penteMoyenne,
                      const tympan::AcousticSource& source, const bool& fromSource, const bool& toRecepteur,
                      TYTabEtape& Etapes, double& longueur) const;
 
    void addEtapeDirectToChemins(const TYEtape& EtapeDirecte, TYTabChemin& TabChemin) const;
 
    bool solve(TYTrajet& trajet);
 
private:
    OSpectreComplex getReflexionSpectrumAt(const OSegment3D& incident, double length,
                                           const OSegment3D& segPente,
                                           const tympan::AcousticSource& source) const;
 
    void meanSlope(const OSegment3D& director, OSegment3D& slope) const;
 
private:
    FRIEND_TEST(test_TYAcousticModel, calculAttDiffraction);
    FRIEND_TEST(test_TYAcousticModel, calculC);
    FRIEND_TEST(test_TYAcousticModel, limAttDiffraction);
    FRIEND_TEST(test_TYAcousticModel, computeCheminReflexion);
    bool _useSol;
    bool _useReflex;
    int _propaCond;
    bool _useAtmo;
    bool _interference;
    double _paramH;
 
    std::unique_ptr<AtmosphericConditions> pSolverAtmos;
 
    OSpectre _lambda;
    OSpectreComplex _absoNulle;
 
    TYSolver& _solver;
};
 
#endif // __TYACOUSTICMODEL__