TYProjet.cpp

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/>.
 */
 
#include "Tympan/core/logging.h"
#include "Tympan/models/business/TYPluginManager.h"
#include "Tympan/models/business/TYPreferenceManager.h"
#include "Tympan/models/business/OLocalizator.h"
#include "Tympan/models/business/TYProgressManager.h"
#include "Tympan/models/business/TYRectangularMaillage.h"
 
#include "TYProjet.h"
 
#if TY_USE_IHM
    #include "Tympan/gui/widgets/TYProjetWidget.h"
    #include "Tympan/gui/gl/TYProjetGraphic.h"
#endif
 
TY_EXTENSION_INST(TYProjet);
TY_EXT_GRAPHIC_INST(TYProjet);
 
#define TR(id) OLocalizator::getString("OMessageManager", (id))
 
bool TYProjet::gSaveValues = true;
 
TYProjet::TYProjet()
    : _auteur(""), _dateCreation("2001-10-01"), _dateModif("2001-10-01"), _comment(""),
      _pSite(new TYSiteNode()), _pCurrentCalcul(new TYCalcul()), _delaunayTolerance(0.0001),
      _maxDistBetweenPoints(200.0), _useDefaultGeomVal(true), _bStatusSolver(true)
{
    _name = TYNameManager::get()->generateName(getClassName());
 
    _pSite->setParent(this);
    _pSite->setRoot(true);
    _pSite->setUseEmpriseAsCrbNiv(true);
    _pSite->setProjet(this); // Definit le projet courant
 
    _pCurrentCalcul->setParent(this);
    _listCalcul.push_back(_pCurrentCalcul);
 
    // Ajout du site au calcul courant
    _pCurrentCalcul->addToSelection(_pSite, false);
 
    // first step, force spectre form to Octave Form
    _spectreForm = SPECTRE_FORM_TIERS;
 
    forceUseDefault();
}
 
TYProjet::TYProjet(const TYProjet& other)
{
    *this = other;
}
 
TYProjet::~TYProjet()
{
    remAllCalcul();
}
 
TYProjet& TYProjet::operator=(const TYProjet& other)
{
    if (this != &other)
    {
        TYElement::operator=(other);
        _auteur = other._auteur;
        _dateCreation = other._dateCreation;
        _dateModif = other._dateModif;
        _comment = other._comment;
        _pSite = other._pSite;
        _pointsControl = other._pointsControl;
        _pCurrentCalcul = other._pCurrentCalcul;
        _listCalcul = other._listCalcul;
        _delaunayTolerance = other._delaunayTolerance;
        _maxDistBetweenPoints = other._maxDistBetweenPoints;
        _useDefaultGeomVal = other._useDefaultGeomVal;
        _spectreForm = other._spectreForm;
    }
 
    TYCourbeNiveau::setDefaultDistMax(_maxDistBetweenPoints); // MaJ Courbes de niveau
 
    return *this;
}
 
bool TYProjet::operator==(const TYProjet& other) const
{
    if (this != &other)
    {
        if (TYElement::operator!=(other))
        {
            return false;
        }
        if (_auteur != other._auteur)
        {
            return false;
        }
        if (_dateCreation != other._dateCreation)
        {
            return false;
        }
        if (_dateModif != other._dateModif)
        {
            return false;
        }
        if (_comment != other._comment)
        {
            return false;
        }
        if (_pSite != other._pSite)
        {
            return false;
        }
        if (_pointsControl != other._pointsControl)
        {
            return false;
        }
        if (_pCurrentCalcul != other._pCurrentCalcul)
        {
            return false;
        }
        if (!(_listCalcul == other._listCalcul))
        {
            return false;
        }
        if (_delaunayTolerance != other._delaunayTolerance)
        {
            return false;
        }
        if (_maxDistBetweenPoints != other._maxDistBetweenPoints)
        {
            return false;
        }
        if (_useDefaultGeomVal != other._useDefaultGeomVal)
        {
            return false;
        }
        if (_spectreForm != other._spectreForm)
        {
            return false;
        }
    }
    return true;
}
 
bool TYProjet::operator!=(const TYProjet& other) const
{
    return !operator==(other);
}
 
std::string TYProjet::toString() const
{
    return "TYProjet";
}
 
DOM_Element TYProjet::toXML(DOM_Element& domElement)
{
    DOM_Element domNewElem = TYElement::toXML(domElement);
 
    TYXMLTools::addElementStringValue(domNewElem, "auteur", _auteur);
    TYXMLTools::addElementStringValue(domNewElem, "dateCreation", _dateCreation);
    TYXMLTools::addElementStringValue(domNewElem, "dateModif", _dateModif);
    TYXMLTools::addElementStringValue(domNewElem, "comment", _comment);
    TYXMLTools::addElementDoubleValue(domNewElem, "delaunayTolerence", _delaunayTolerance);
    TYXMLTools::addElementDoubleValue(domNewElem, "distMaxBetweenPoints", _maxDistBetweenPoints);
    TYXMLTools::addElementBoolValue(domNewElem, "useDefGeomVal", _useDefaultGeomVal);
    // Sauvegarde de la forme d'affichage du spectre
    TYXMLTools::addElementIntValue(domNewElem, "spectreForm", static_cast<int>(_spectreForm));
 
    if (_pCurrentCalcul)
    {
        TYXMLTools::addElementStringValue(domNewElem, "currentCalculID", _pCurrentCalcul->getID().toString());
    }
 
    // Ajout du site node sur lequel s'effectue le calcul
    DOM_Document domDoc = domElement.ownerDocument();
 
    // Site
    _pSite->toXML(domNewElem);
 
    // Points de controle
    DOM_Element pointsControlNode = domDoc.createElement("PointsControl");
    domNewElem.appendChild(pointsControlNode);
    for (unsigned int i = 0; i < _pointsControl.size(); i++)
    {
        _pointsControl[i]->toXML(pointsControlNode);
    }
 
    // Noise Maps
    DOM_Element noiseMapsNode = domDoc.createElement("NoiseMaps");
    domNewElem.appendChild(noiseMapsNode);
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        _maillages[i]->toXML(noiseMapsNode);
    }
 
    // Calculs
    for (unsigned int i = 0; i < _listCalcul.size(); i++)
    {
        _listCalcul.at(i)->toXML(domNewElem);
    }
 
    return domNewElem;
}
 
int TYProjet::fromXML(DOM_Element domElement)
{
    TYElement::fromXML(domElement);
 
    _pCurrentCalcul = nullptr;
 
    remAllCalcul();
 
    bool auteurOk = false;
    bool dateCreationOk = false;
    bool dateModifOk = false;
    bool commentOk = false;
    bool currentCalculIDOk = false;
    bool delaunayOk = false;
    bool dMaxOk = false;
    bool useDefaultOk = false;
    bool keepAltiUpTodate = false;
    bool spectreFormOk = false;
    int spectreFormInt = static_cast<int>(_spectreForm); // default current
 
    QString currentCalculID;
    LPTYCalcul pCalcul = new TYCalcul();
    int readOk = 0; // Indicateur de bonne relecture du calcul
    int i = 0;
 
    DOM_Element elemCur;
 
    QDomNodeList childs = domElement.childNodes();
    int childcount = childs.length();
    for (i = 0; i < childcount; i++)
    {
        elemCur = childs.item(i).toElement();
        OMessageManager::get()->info("Charge element du projet %d/%d.", i + 1, childcount);
 
        TYXMLTools::getElementStringValue(elemCur, "auteur", _auteur, auteurOk);
        TYXMLTools::getElementStringValue(elemCur, "dateCreation", _dateCreation, dateCreationOk);
        TYXMLTools::getElementStringValue(elemCur, "dateModif", _dateModif, dateModifOk);
        TYXMLTools::getElementStringValue(elemCur, "comment", _comment, commentOk);
        TYXMLTools::getElementStringValue(elemCur, "currentCalculID", currentCalculID, currentCalculIDOk);
 
        TYXMLTools::getElementDoubleValue(elemCur, "delaunayTolerence", _delaunayTolerance, delaunayOk);
        TYXMLTools::getElementDoubleValue(elemCur, "distMaxBetweenPoints", _maxDistBetweenPoints, dMaxOk);
        TYXMLTools::getElementBoolValue(elemCur, "useDefGeomVal", _useDefaultGeomVal, useDefaultOk);
        TYXMLTools::getElementIntValue(elemCur, "spectreForm", spectreFormInt, spectreFormOk);
 
        if (_pSite->callFromXMLIfEqual(elemCur))
        {
            _pSite->setChildsNotInCurrentCalcul();
        }
        // On sauvegarde l'information de mise à jour de l'altimétrie
        keepAltiUpTodate = _pSite->getAltimetry()->isUpToDate();
 
        // Points de controle
        LPTYPointControl pPointControl = new TYPointControl();
        if (elemCur.nodeName() == "PointsControl")
        {
            DOM_Element elemCur2;
            QDomNodeList childs2 = elemCur.childNodes();
 
            for (unsigned int j = 0; j < childs2.length(); j++)
            {
                elemCur2 = childs2.item(j).toElement();
                if (pPointControl->callFromXMLIfEqual(elemCur2))
                {
                    addPointControl(pPointControl);
                    pPointControl = new TYPointControl();
                }
            }
        }
 
        // Noise maps
        LPTYGeometryNode pGeoNode = new TYGeometryNode();
        if (elemCur.nodeName() == "NoiseMaps")
        {
            DOM_Element elemCur2;
            QDomNodeList childs2 = elemCur.childNodes();
 
            for (unsigned int j = 0; j < childs2.length(); j++)
            {
                elemCur2 = childs2.item(j).toElement();
                if (pGeoNode->callFromXMLIfEqual(elemCur2))
                {
                    addMaillage(pGeoNode);
                    pGeoNode = new TYGeometryNode();
                }
            }
        }
 
        // Calculs
        pCalcul->setParent(this);
        if (pCalcul->callFromXMLIfEqual(elemCur, &readOk))
        {
            // Insure compatibility when loading Code_TYMPAN T310 projects
            // Enforce pCalcul->_mapPointCtrlSpectre consistency
            if (pCalcul->getMapPointCtrlSpectre().empty() && !pCalcul->hasResuCtrlPoints())
            {
                OMessageManager::get()->info("T310 project migration : enforce receptors configuration "
                                             "consistency in calculation settings");
                for (unsigned int i = 0; i < _pointsControl.size(); i++)
                {
                    pCalcul->addPtCtrlToResult(_pointsControl[i]);
                }
            }
            addCalcul(pCalcul);
            pCalcul = new TYCalcul();
        }
    }
 
    // Mise a jour des parametre pour la triangulation du site
    verifGeometricParam();
 
    // Update TYCalcul control point spectrum if needed
    // --> compatibility with previous TYPointControl management paradigm (see CS 0a6aab520874@EDF_RD)
    for (unsigned int i = 0; i < _pointsControl.size(); i++)
    {
        if (_pointsControl[i].getRealPointer()->getAllUses() != nullptr)
        {
            std::map<TYUUID, LPTYSpectre>* compatibilityVector =
                static_cast<std::map<TYUUID, LPTYSpectre>*>(_pointsControl[i].getRealPointer()->getAllUses());
            std::map<TYUUID, LPTYSpectre>::iterator it;
            for (it = compatibilityVector->begin(); it != compatibilityVector->end(); it++)
            {
                TYCalcul* pCalc = dynamic_cast<TYCalcul*>(getInstance((*it).first));
                if (pCalc == nullptr)
                {
                    continue;
                }
 
                TYSpectre* pSpectre = (*it).second;
                pCalc->addPtCtrlToResult(_pointsControl[i]);
                pCalc->setSpectre(_pointsControl[i].getRealPointer(), pSpectre);
            }
 
            // Cleaning map after use
            compatibilityVector->clear();
            _pointsControl[i].getRealPointer()->cleanAllUses();
            compatibilityVector = nullptr;
        }
    }
 
    // update TYCalcul noise map spectrums if needed
    // --> compatibility with previous TYMaillage management paradigm
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        TYMaillage* pMaillage = dynamic_cast<TYMaillage*>(_maillages[i]->getElement());
        if (pMaillage->getAllUses() != nullptr)
        {
            TYTabLPSpectre* compatibilityVector = static_cast<TYTabLPSpectre*>(pMaillage->getAllUses());
 
            // Balayage de tous les calculs
            for (unsigned j = 0; j < _listCalcul.size(); j++)
            {
                if (pMaillage->etat(_listCalcul[j]) == true)
                {
                    _listCalcul[j]->setNoiseMapSpectrums(pMaillage, *compatibilityVector);
                }
            }
 
            // Cleaning map after use
            compatibilityVector->clear();
            pMaillage->cleanAllUses();
            compatibilityVector = nullptr;
        }
    }
 
    // Identification du calcul courant
    bool bCurrentCalculFound = false;
    for (i = 0; i < (int)_listCalcul.size(); i++)
    {
        if (_listCalcul[i]->getID().toString() == currentCalculID)
        {
            setCurrentCalcul(_listCalcul[i]);
            bCurrentCalculFound = true;
            break;
        }
    }
 
    if (!bCurrentCalculFound && (_listCalcul.size() > 0))
    {
        setCurrentCalcul(_listCalcul[0]);
    }
 
    if (spectreFormOk)
    {
        switch (static_cast<TYSpectreForm>(spectreFormInt))
        {
            case SPECTRE_FORM_TIERS:
            case SPECTRE_FORM_OCT:
            case SPECTRE_FORM_CST_DF:
            case SPECTRE_FORM_UNSHAPED:
                _spectreForm = static_cast<TYSpectreForm>(spectreFormInt);
                break;
            default:
                // Rétro-compat / valeur invalide -> garder le défaut (_spectreForm inchangé)
                break;
        }
    }
 
    if (readOk == -1)
    {
        return readOk;
    }
 
    // A REVOIR _pElement ne devrait pas renvoyer le projet dans showsources/TYModelerFrame...
    if (_pSite)
    {
        _pSite->update(_pSite);
        _pSite->getAltimetry()->setIsUpToDate(keepAltiUpTodate);
    }
    return 1;
}
 
bool TYProjet::addPointControl(LPTYPointControl pPointControl)
{
    assert(pPointControl);
 
    pPointControl->setParent(this);
    _pointsControl.push_back(pPointControl);
 
    if (getCurrentCalcul()) // Uniquement pour la creation depuis un modeleur
    {
        getCurrentCalcul()->addPtCtrlToResult(
            pPointControl); // Le point de controle est actif dans le calcul courant
    }
 
    setIsGeometryModified(true);
 
    return true;
}
 
bool TYProjet::remPointControl(LPTYPointControl pPointControl)
{
    bool ret = false;
    assert(pPointControl);
    TYTabLPPointControl::iterator ite;
 
    for (ite = _pointsControl.begin(); ite != _pointsControl.end(); ite++)
    {
        if ((*ite) == pPointControl)
        {
            // Nettoyage des calculs concernes
            for (unsigned int i = 0; i < this->_listCalcul.size(); i++)
            {
                if ((*ite)->etat(_listCalcul[i]) == true)
                {
                    _listCalcul[i]->remPtCtrlFromResult((*ite));
                }
            }
            _pointsControl.erase(ite);
            ret = true;
            setIsGeometryModified(true);
            break;
        }
    }
 
    return ret;
}
 
LPTYPointControl TYProjet::duplicatePointControl(const LPTYPointControl& pPoint)
{
    assert(pPoint);
 
    // Duplication du point
    LPTYPointControl pPointCopy = new TYPointControl(*pPoint);
 
    // Ajout a la liste des points du projet
    addPointControl(pPointCopy);
 
    return pPointCopy;
}
 
bool TYProjet::remAllPointControl()
{
    _pointsControl.clear();
    setIsGeometryModified(true);
 
    return true;
}
 
bool TYProjet::addCalcul(LPTYCalcul pCalcul)
{
    assert(pCalcul);
 
    pCalcul->setParent(this);
    _listCalcul.push_back(pCalcul);
 
    return true;
}
 
bool TYProjet::remCalcul(const LPTYCalcul pCalcul)
{
    assert(pCalcul);
    bool ret = false;
 
    // On ne peut detruire le calcul courant
    if (pCalcul == getCurrentCalcul())
    {
        return false;
    }
 
    TYTabLPCalcul::iterator ite;
 
    for (ite = _listCalcul.begin(); ite != _listCalcul.end(); ite++)
    {
        if ((*ite) == pCalcul)
        {
            // remove calcul from noise map and control point status list
            cleanReceptorsStatus((*ite));
 
            // Remove calcul from list
            _listCalcul.erase(ite);
            ret = true;
            break;
        }
    }
 
    return ret;
}
 
bool TYProjet::remCalcul(QString idCalcul)
{
    bool ret = false;
    TYTabLPCalcul::iterator ite;
 
    for (ite = _listCalcul.begin(); ite != _listCalcul.end(); ite++)
    {
        if ((*ite)->getID().toString() == idCalcul)
        {
            return remCalcul((*ite));
        }
    }
 
    return ret;
}
 
void TYProjet::remAllCalcul()
{
    _listCalcul.clear();
}
 
void TYProjet::unsetSite()
{
    _pSite = NULL;
}
 
void TYProjet::setCurrentCalcul(LPTYCalcul pCurCalcul)
{
    assert(pCurCalcul);
 
    // On test si le calcul est present dans la liste
 
    // Statut solveur ok par defaut
    _bStatusSolver = true;
 
    // Calcul courant
    _pCurrentCalcul = pCurCalcul;
 
    // Test si le solveur du calcul est bien disponible
    if (!TYPluginManager::get()->exist(_pCurrentCalcul->getSolverId()))
    {
        _bStatusSolver = false;
    }
 
    if (_pSite)
    {
        // Le site root associe au projet est obligatoirement actif dans le calcul courant
        _pCurrentCalcul->setSite(_pSite);
        _pSite->setInCurrentCalcul(true, false);
        TYListID aTYListID = _pCurrentCalcul->getElementSelection();
        _pSite->updateCurrentCalcul(aTYListID, true);
    }
 
    // Mise  jour des maillages
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        dynamic_cast<TYMaillage*>(_maillages.at(i)->getElement())->updateFromCalcul(_pCurrentCalcul);
    }
 
#if TY_USE_IHM
    updateGraphic();
#endif
    setIsGeometryModified(true);
    setIsAcousticModified(true);
}
 
bool TYProjet::updateAltiPointControle(TYPointControl* pPtControl, const TYAltimetrie* pAlti)
{
    assert(pPtControl);
    assert(pAlti);
 
    bool modified = true;
 
    // On va modifier le point de contrele (l'altitude seulement)
    // (Utile pour l'update du GraphicObject du maillage)
    pPtControl->setIsGeometryModified(true);
 
    // Init
    OPoint3D pt;
 
    pt._x = pPtControl->_x;
    pt._y = pPtControl->_y;
    pt._z = 0.0;
 
    // Recherche de l'altitude
    modified = pAlti->updateAltitude(pt);
 
    // Ajout de l'offset en Z
    pt._z += pPtControl->getHauteur();
 
    pPtControl->_z = pt._z;
 
    // Le point est maintenant up to date
    pPtControl->setIsGeometryModified(false);
 
    return modified;
}
 
bool TYProjet::updateAltiRecepteurs()
{
    TYAltimetrie* pAlti = getSite()->getAltimetry()._pObj;
    return updateAltiRecepteurs(pAlti);
}
 
bool TYProjet::updateAltiRecepteurs(const TYAltimetrie* pAlti)
{
    assert(pAlti);
    bool bNoPbAlti = true;
    bool modified = false;
 
    // Mise a jour de l'altitude pour les points de controle
    unsigned int i = 0;
    for (i = 0; i < _pointsControl.size(); i++)
    {
        // XXX See ticket https://extranet.logilab.fr/ticket/1484180
        // The coordinates of the problematic point need to be properly reported
        // or idealy an exception should be thrown.
        modified |= updateAltiPointControle(getPointControl(i), pAlti);
        bNoPbAlti &= modified;
    }
 
    // XXX See ticket https://extranet.logilab.fr/ticket/1484180:
    // Why handling of PointControle and Recepteur are so different ?
 
    // Mise a jour de l'altitude pour les recepteurs du calcul
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        updateAltiMaillage(_maillages.at(i), pAlti);
    }
 
    // Done
    if (!bNoPbAlti)
    {
        OMessageManager::get()->info(TR("msg_pbalti"));
    }
 
    setIsGeometryModified(modified);
    return bNoPbAlti;
}
 
LPTYCalcul TYProjet::duplicateCalcul(const LPTYCalcul pCalcul)
{
    assert(pCalcul);
 
    LPTYCalcul pCalculCopy = new TYCalcul();
 
    // Duplication
    // On dit qu'on copie aussi l'ID pour que les ID des elements
    // composant le calcul soit copiees
    pCalculCopy->deepCopy(pCalcul, false);
 
    // On incremente aussi le numero du calcul
    pCalculCopy->setNumero(pCalculCopy->getNumero() + 1);
 
    // Duplication de l'etat des points de controle
    duplicatePtCalcState(pCalcul, pCalculCopy);
 
    // Duplication de l'tat des cartographies
    duplicateNoiseMapState(pCalcul, pCalculCopy);
 
    // Ajout
    addCalcul(pCalculCopy);
 
    return pCalculCopy;
}
 
void TYProjet::duplicatePtCalcState(const TYCalcul* pCalculRef, TYCalcul* pCalculNew)
{
    QString idCalculRef = pCalculRef->getID().toString();
    QString idCalculNew = pCalculNew->getID().toString();
 
    for (unsigned int i = 0; i < _pointsControl.size(); i++)
    {
        _pointsControl[i]->duplicateEtat(idCalculRef, idCalculNew);
    }
}
 
void TYProjet::duplicateNoiseMapState(const TYCalcul* pCalculRef, TYCalcul* pCalculNew)
{
    QString idCalculRef = pCalculRef->getID().toString();
    QString idCalculNew = pCalculNew->getID().toString();
 
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        dynamic_cast<TYMaillage*>(_maillages[i]->getElement())->duplicateEtat(idCalculRef, idCalculNew);
    }
}
 
void TYProjet::remElmtFromCalculs(TYElement* pElement)
{
    for (unsigned int i = 0; i < _listCalcul.size(); i++)
    {
        _listCalcul[i]->remToSelection(pElement);
    }
}
 
void TYProjet::remTabElmtFromCalculs(std::vector<LPTYGeometryNode> tabGeoNode)
{
    for (unsigned int i = 0; i < tabGeoNode.size(); i++)
    {
        remElmtFromCalculs(tabGeoNode[i]->getElement());
    }
}
 
void TYProjet::verifGeometricParam()
{
 
    _delaunayTolerance = _delaunayTolerance <= 0.0 ? 0.0001 : _delaunayTolerance;
    _delaunayTolerance = _delaunayTolerance > 0.05 ? 0.05 : _delaunayTolerance;
    _maxDistBetweenPoints = _maxDistBetweenPoints < 1.0 ? 1.0 : _maxDistBetweenPoints;
 
    TYCourbeNiveau::setDefaultDistMax(_maxDistBetweenPoints); // MaJ Courbes de niveau
}
 
void TYProjet::forceUseDefault()
{
    _delaunayTolerance = getDefaultDelaunayTolerence();
    _maxDistBetweenPoints = getDefaultDMax();
    verifGeometricParam();
}
 
const double TYProjet::getDefaultDelaunayTolerence() const
{
    double delaunay = 0.0001;
 
#if TY_USE_IHM
    if (TYPreferenceManager::exists(TYDIRPREFERENCEMANAGER, "DelaunayTolerance"))
    {
        delaunay = TYPreferenceManager::getDouble(TYDIRPREFERENCEMANAGER, "DelaunayTolerance");
    }
    else
    {
        TYPreferenceManager::setDouble(TYDIRPREFERENCEMANAGER, "DelaunayTolerance", delaunay);
    }
#endif
 
    return delaunay;
}
 
const double TYProjet::getDefaultDMax() const
{
    double DMax = 100.0;
 
#if TY_USE_IHM
    if (TYPreferenceManager::exists(TYDIRPREFERENCEMANAGER, "DistMinPtCrbNiv"))
    {
        DMax = TYPreferenceManager::getDouble(TYDIRPREFERENCEMANAGER, "DistMinPtCrbNiv");
    }
    else
    {
        TYPreferenceManager::setDouble(TYDIRPREFERENCEMANAGER, "DistMinPtCrbNiv", DMax);
    }
#endif
 
    return DMax;
}
 
void TYProjet::setDMax(const double& val)
{
    _maxDistBetweenPoints = val;
    TYCourbeNiveau::setDefaultDistMax(_maxDistBetweenPoints); // MaJ Courbes de niveau
}
 
void TYProjet::setSite(const LPTYSiteNode pSite)
{
    assert(pSite);
    _pSite = pSite;
    _pSite->setParent(this);
    _pSite->setProjet(this);
    _pSite->setRoot(true); // Ce site est forcemment le site racine
    _pSite->addToCalcul();
    setIsGeometryModified(true);
}
 
void TYProjet::updateCalculsWithMaillage(TYMaillage* pMaillage)
{
    for (unsigned int i = 0; i < _listCalcul.size(); i++)
    {
        _listCalcul[i]->updateMaillage(pMaillage);
    }
}
 
bool TYProjet::addMaillage(LPTYMaillageGeoNode pMaillageGeoNode)
{
    assert(pMaillageGeoNode);
    assert(pMaillageGeoNode->getElement());
 
    pMaillageGeoNode->setParent(this);
    pMaillageGeoNode->getElement()->setParent(this);
 
    _maillages.push_back(pMaillageGeoNode);
 
    // Set in current calcul
    if (getCurrentCalcul() != nullptr)
    {
        getCurrentCalcul()->addMaillage(dynamic_cast<TYMaillage*>(pMaillageGeoNode->getElement()));
    }
 
    setIsGeometryModified(true);
 
    return true;
}
 
bool TYProjet::addMaillage(LPTYMaillage pMaillage)
{
    return addMaillage(new TYMaillageGeoNode((LPTYElement)pMaillage));
}
 
bool TYProjet::remMaillage(const LPTYMaillage pMaillage)
{
    assert(pMaillage);
 
    TYGeometryNode* pNode = TYGeometryNode::GetGeoNode(dynamic_cast<TYElement*>(pMaillage._pObj));
 
    return remMaillage(pNode);
}
 
bool TYProjet::remMaillage(const LPTYMaillageGeoNode pMaillageGeoNode)
{
    assert(pMaillageGeoNode);
    bool ret = false;
 
    // Suppression dans tous les calculs
    for (unsigned int i = 0; i < _listCalcul.size(); i++)
    {
        _listCalcul.at(i)->remMaillage(dynamic_cast<TYMaillage*>(pMaillageGeoNode->getElement()));
    }
 
    // Suppression de la liste des maillages du projet
    TYTabMaillageGeoNode::iterator ite;
 
    for (ite = _maillages.begin(); ite != _maillages.end(); ite++)
    {
        if ((*ite) == pMaillageGeoNode)
        {
            _maillages.erase(ite);
            ret = true;
            setIsGeometryModified(true);
            break;
        }
    }
 
    return ret;
}
 
bool TYProjet::remMaillage(QString idMaillage)
{
    bool ret = false;
    TYTabMaillageGeoNode::iterator ite;
 
    for (ite = _maillages.begin(); ite != _maillages.end(); ite++)
    {
        if ((*ite)->getElement()->getID().toString() == idMaillage)
        {
            _maillages.erase(ite);
            ret = true;
            setIsGeometryModified(true);
            break;
        }
    }
 
    return ret;
}
 
bool TYProjet::remAllMaillage()
{
    _maillages.clear();
    setIsGeometryModified(true);
 
    return true;
}
 
LPTYMaillageGeoNode TYProjet::findMaillage(const LPTYMaillage pMaillage)
{
    assert(pMaillage);
    TYTabMaillageGeoNode::iterator ite;
 
    for (ite = _maillages.begin(); ite != _maillages.end(); ite++)
    {
        if (TYMaillage::safeDownCast((*ite)->getElement()) == pMaillage)
        {
            return (*ite);
        }
    }
 
    return NULL;
}
 
void TYProjet::updateGraphicMaillage()
{
#if TY_USE_IHM
    for (unsigned int i = 0; i < getMaillages().size(); i++)
    {
        TYMaillage* pMaillage = getMaillage(i);
        pMaillage->getGraphicObject()->update();
    }
#endif
}
 
bool TYProjet::updateAltiMaillage(TYMaillageGeoNode* pMaillageGeoNode)
{
    TYAltimetrie* pAlti = getSite()->getAltimetry()._pObj;
    return updateAltiMaillage(pMaillageGeoNode, pAlti);
}
 
bool TYProjet::updateAltiMaillage(TYMaillageGeoNode* pMaillageGeoNode, const TYAltimetrie* pAlti)
{
    assert(pMaillageGeoNode);
    assert(pAlti);
 
    bool modified = true;
 
    if (pMaillageGeoNode == NULL)
    {
        return false;
    }
 
    TYPoint pt;
    TYMaillage* pMaillage = TYMaillage::safeDownCast(pMaillageGeoNode->getElement());
 
    if (pMaillage == NULL)
    {
        return false;
    }
 
    OMatrix matrix = pMaillageGeoNode->getMatrix();
    OMatrix matrixinv = matrix.getInvert();
    TYTabLPPointCalcul& tabpoint = pMaillage->getPtsCalcul();
 
    bool bNoPbAlti = true; // Permet de tester si tous les points sont altimtriss correctement.
 
    if (pMaillage->getComputeAlti()) // Cas des maillages rectangulaires et lineaires horizontaux
    {
 
#if TY_USE_IHM
        TYProgressManager::setMessage("Calcul des altitudes des points de maillage");
        TYProgressManager::set(static_cast<uint32>(tabpoint.size()));
#endif // TY_USE_IHM
 
        for (unsigned int i = 0; i < tabpoint.size(); i++)
        {
 
#if TY_USE_IHM
            bool cancel = false;
            TYProgressManager::step(cancel);
            if (cancel)
            {
                break;
            }
#endif // TY_USE_IHM
 
            pt = *tabpoint[i];
            // Passage au repere du site
            pt = matrix * pt;
            pt._z = 0;
 
            // Recherche de l'altitude
            bNoPbAlti &= pAlti->updateAltitude(pt);
 
            // Retour au repere d'origine
            pt = matrixinv * pt;
 
            // Ajout de l'offset en Z
            pt._z += pMaillage->getHauteur();
 
            // Application du calcul
            tabpoint[i]->_x = pt._x;
            tabpoint[i]->_y = pt._y;
            tabpoint[i]->_z = pt._z;
        }
 
#if TY_USE_IHM
        TYProgressManager::stepToEnd();
#endif // TY_USE_IHM
    }
    else // Cas des maillages verticaux
    {
        // Init
        pt = pMaillageGeoNode->getORepere3D()._origin;
        pt._z = 0.0;
 
        // Recherche de l'altitude
        bNoPbAlti &= pAlti->updateAltitude(pt);
 
        // Ajout de l'offset en Z
        pt._z += pMaillage->getHauteur();
 
        if (pMaillage->isA("TYRectangularMaillage"))
        {
            // Ajout d'un offset de la demi-hauteur (maillages verticaux)
            pt._z += ((TYRectangularMaillage*)pMaillage)->getRectangle()->getSizeY() / 2.0;
        }
 
        ORepere3D repere3D = pMaillageGeoNode->getORepere3D();
        repere3D._origin._z = pt._z;
        pMaillageGeoNode->setRepere(repere3D);
 
        modified = true;
    }
 
    if (!bNoPbAlti) // Certains point pas altimetrises
    {
        OMessageManager::get()->info(TR("msg_pbalti"));
    }
 
    // Done
    if (modified)
    {
#if TY_USE_IHM
        pMaillage->updateGraphicTree();
#endif
        setIsGeometryModified(true);
    }
 
    return modified;
}
 
void TYProjet::selectActivePoint(const LPTYSiteNode pSite)
{
    int i = 0, j = 0;
 
    // Recuperation des volumes
    TYInfrastructure* pInfra = pSite->getInfrastructure();
    TYTabAcousticVolumeNodeGeoNode tabVolNodeGeoNode;
 
    // On commence par recuperer la liste des machines et des btiments
    // Batiments
    for (j = 0; j < static_cast<int>(pInfra->getListBatiment().size()); j++)
    {
        // Si ce batiment est actif pour le calcul
        if (TYBatiment::safeDownCast(pInfra->getBatiment(j)->getElement())->isInCurrentCalcul())
        {
            tabVolNodeGeoNode.push_back((LPTYAcousticVolumeNodeGeoNode&)pInfra->getListBatiment()[j]);
        }
    }
 
    // Machines
    for (j = 0; j < static_cast<int>(pInfra->getListMachine().size()); j++)
    {
        // Si cette machine est active pour le calcul
        if (TYMachine::safeDownCast(pInfra->getMachine(j)->getElement())->isInCurrentCalcul())
        {
            tabVolNodeGeoNode.push_back((LPTYAcousticVolumeNodeGeoNode&)pInfra->getListMachine()[j]);
        }
    }
 
    // Ensuite, on verifie si des points du maillage ne sont pas a l'interieur
    std::vector<OPoint3D> tabPts;
 
    for (i = 0; i < static_cast<int>(getMaillages().size()); i++)
    {
        TYMaillage* pMaillage = getMaillage(i);
        TYMaillageGeoNode* pMaillageGeoNode = getMaillages()[i];
 
        // Matrice pour la position de ce maillage
        OMatrix matrixMaillage = pMaillageGeoNode->getMatrix();
        // Nb de points de calcul dans ce maillage
        size_t nbPtsCalcul = pMaillage->getPtsCalcul().size();
 
        tabPts.reserve(nbPtsCalcul); // Reservation du nombre de points pour eviter les mouvements memoire
 
        for (j = 0; (int)j < nbPtsCalcul; j++)
        {
            pMaillage->getPtsCalcul()[j]->setEtat(true);
            tabPts.push_back(matrixMaillage * (*(pMaillage->getPtsCalcul()[j])));
        }
 
        // Detection des points a l'interieur de volumes
        OPoint3D pt;
        for (j = 0; j < tabVolNodeGeoNode.size(); j++)
        {
            OMatrix matrixVolNode = tabVolNodeGeoNode[j]->getMatrix();
            matrixVolNode.invert();
 
            LPTYAcousticVolumeNode pVolumeNode =
                TYAcousticVolumeNode::safeDownCast(tabVolNodeGeoNode[j]->getElement());
            for (int k = 0; k < nbPtsCalcul; k++)
            {
                // Passage au repere du volume
                pt = matrixVolNode * tabPts[k];
 
                // Si le point se trouve a l'interieur de ce volume
                if (pVolumeNode->isInside(pt))
                {
                    // Desactivation de ce point de calcul
                    pMaillage->getPtsCalcul()[k]->setEtat(false);
                }
            }
        }
 
        // Le maillage et son GeoNode sont maintenant up to date
        pMaillage->setIsGeometryModified(false);
        pMaillageGeoNode->setIsGeometryModified(false);
 
        // Vidage du tableau de points
        tabPts.clear();
    }
 
    tabVolNodeGeoNode.clear();
}
 
void TYProjet::cleanReceptorsStatus(TYCalcul* pCalcul)
{
    for (unsigned int i = 0; i < _pointsControl.size(); i++)
    {
        _pointsControl[i]->remEtat(pCalcul);
    }
 
    for (unsigned int i = 0; i < _maillages.size(); i++)
    {
        dynamic_cast<TYMaillage*>(_maillages[i]->getElement())->remEtat(pCalcul);
    }
}
 
const bool TYProjet::displaySpectrumInOctaveForm()
{
    return _spectreForm == SPECTRE_FORM_OCT;
}
 
void TYProjet::updateGraphic()
{
    updateGraphicMaillage();
    TYElement::updateGraphic();
}