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_9613__
#define __TYACOUSTICMODEL_9613__
 
#include "Tympan/core/interfaces.h"
#include "Tympan/solvers/9613Solver/TYChemin.h"
#include "Tympan/solvers/9613Solver/TYSolverDefines.h"
#include "Tympan/models/common/plan.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, const bool left) const;
 
    void computeCheminReflexion(const std::deque<TYSIntersection>& tabIntersect, const OSegment3D& ray,
                                const tympan::AcousticSource& source, TYTabChemin& TabChemins,
                                double distance) const;
 
    void computeCheminSansEcran(const std::deque<TYSIntersection>& tabIntersect, 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;
 
    OSpectreOctave calculC3(const double& epaisseur) const;
 
    OSpectreOctave calculAttDiffraction(const OSegment3D& ray, const double& re, const double& dss,
                                        const double& dsr, const double& width, const bool& vertical) const;
 
    OSpectreOctave limAttDiffraction(const OSpectreOctave& sNC, const OSpectreOctave& C) const;
 
    bool addStep(const OPoint3D& ptStart, const OPoint3D& ptEnd, const tympan::AcousticSource& source,
                 const bool& fromSource, TYTabEtape& Etapes) const;
 
    bool getGroundfactors(const std::deque<TYSIntersection>& tabIntersect, const OSegment3D& ray2D, double hs,
                          double hr, double& Gs, double& Gm, double& Gr) const;
 
    bool getGroundfactors(const std::deque<TYSIntersection>& tabIntersectUpSegment,
                          const std::deque<TYSIntersection>& tabIntersectDownSegment, const OSegment3D& SO2D,
                          const OSegment3D& OR2D, double hs, double hr, double& Gs, double& Gm,
                          double& Gr) const;
 
    /*
     * @brief Compute effective barrier attenuation implementing recommandation paragraph 5.6
     *  of ISO TR 17534-3
     * @param [in] Abar_top Barrier attenuation for the vertical path.
     * @param [in] Abar_left Barrier attenuation for the left lateral path.
     * @param [in] Abar_right Barrier attenuation for the right lateral path.
     *
     */
    OSpectreOctave computeEffectiveBarAttenuation(const OSpectreOctave& Abar_top,
                                                  const OSpectreOctave& Abar_left,
                                                  const OSpectreOctave& Abar_right);
 
    /*
     * @brief Compute mean slope plan - Kept for ISO 9613-2 ground effect alternative method.
     *
     * @param penteMoyenne Mean slope used to compute the plan corresponding to this slope.
     *
     * @return OPlan The computed mean slope plan.
     */
    OPlan buildMeanSlopePlan(const OSegment3D& penteMoyenne) const;
 
    bool solve(TYTrajet& trajet);
 
private:
    OSpectreOctave getReflexionSpectrumAt(const OSegment3D& incident, double length,
                                          const OSegment3D& segPente,
                                          const tympan::AcousticSource& source) const;
 
    void meanSlope(const OSegment3D& director, OSegment3D& slope) const;
 
    bool computeSegmentEdgesHeights(double& hauteurA, double& hauteurB, const OSegment3D& meanSlope,
                                    const OSegment3D& ray) const;
 
    bool computeGZone(const OPoint3D& ptDebut, const OPoint3D& ptFin, double& GZone, double& dpZone,
                      const std::deque<TYSIntersection>& tabIntersect) const;
 
private:
    bool _useSol;
    bool _useReflex;
    int _propaCond;
    bool _useAtmo;
    bool _interference;
    double _paramH;
 
    std::unique_ptr<AtmosphericConditions> pSolverAtmos;
 
    OSpectreOctave _lambda;
    OSpectreOctave _absoNulle;
 
    TYSolver& _solver;
};
 
#endif // __TYACOUSTICMODEL_9613__