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__