TYRectangularMaillageGraphic.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/color.h"
#include "Tympan/models/common/3d.h"
#include "Tympan/models/business/TYPreferenceManager.h"
#include "Tympan/models/business/TYRectangularMaillage.h"
#include "Tympan/models/business/TYProjet.h"
#include "Tympan/models/business/TYCalcul.h"
#include "Tympan/gui/gl/TYPickingTable.h"
#include "Tympan/gui/gl/TYPaletteGraphic.h"
#include "Tympan/gui/gl/TYPanelGraphic.h"
#include "TYRectangularMaillageGraphic.h"
#include "Tympan/gui/gl/TYRenderContext.h"
#include <math.h>
 
TYRectangularMaillageGraphic::TYRectangularMaillageGraphic(TYRectangularMaillage* pElement)
    : TYElementGraphic(pElement)
{
    _Rendering = RENDER_SURF;
 
    // XXX What are those constants ?
    _scalarX = 0.8;
    _scalarY = 0.25;
 
    _scalarW = 150.0;
    _scalarH = 15.0;
 
    _panelW = 125.0;
    _panelH = 40.0;
 
    TYElementGraphic* pTYElementGraphic = getElement()->getPalette()->getGraphicObject();
    ((TYPaletteGraphic*)pTYElementGraphic)->setWidth(_scalarW);
    ((TYPaletteGraphic*)pTYElementGraphic)->setHeight(_scalarH);
 
    pTYElementGraphic = getElement()->getPanel()->getGraphicObject();
    ((TYPanelGraphic*)pTYElementGraphic)->setWidth(_panelW);
    ((TYPanelGraphic*)pTYElementGraphic)->setHeight(_panelH);
    pTYElementGraphic->setVisible(false);
}
 
TYRectangularMaillageGraphic::~TYRectangularMaillageGraphic() {}
 
void TYRectangularMaillageGraphic::update(bool force /*=false*/)
{
    TYRectangularMaillage* pRectangularMaillage = getElement();
    // TODO We do not want regular scale any longer.
 
    // And now... the mesh
    pRectangularMaillage->computeMesh(_mesh);
 
    // And ... the isocurve
    pRectangularMaillage->computeIsoCurve(_mesh, _isoCurve);
 
    pRectangularMaillage->getPalette()->getGraphicObject()->update(force);
 
    TYElementGraphic::update(force);
}
 
//===================================================================
 
OBox TYRectangularMaillageGraphic::computeBoundingBox() const
{
    OBox boundingBox;
 
    int index1 = 0, index2 = 0;
    unsigned long largeur = 0, hauteur = 0;
 
    TYRectangularMaillage* pMaillage = getElement();
    OVector3D vecX, vecY;
    pMaillage->getDimensionsAndSteps(largeur, hauteur, vecX, vecY);
 
    const unsigned int larg = largeur; // +1;
    const unsigned int haut = hauteur; // +1;
 
    TYTabLPPointCalcul* pPtsCalcul = &pMaillage->getPtsCalcul();
 
    for (int i = 0; i < haut - 1; i++)
    {
        for (int j = 0; j < larg; j++)
        {
 
            index1 = i * larg + j;
            TYPointCalcul* pPtCalcul = pPtsCalcul->at(index1);
            boundingBox.Enlarge((float)(pPtCalcul->_x), (float)(pPtCalcul->_y), (float)(pPtCalcul->_z));
            index2 = (i + 1) * larg + j; // nextPoint(pPtsCalcul, index1, (i+1)*larg+j);
            pPtCalcul = pPtsCalcul->at(index2);
            boundingBox.Enlarge((float)(pPtCalcul->_x), (float)(pPtCalcul->_y), (float)(pPtCalcul->_z));
        }
    }
    return boundingBox;
}
 
void TYRectangularMaillageGraphic::display(TYRenderContext& renderContext)
{
    // Pas affiche si inactif
    if (getElement()->etat() == false)
    {
        return;
    }
 
    if (_mesh.empty())
    {
        update();
    }
 
    // Les points de ce maillage
    TYRectangularMaillage* pMaillage = getElement();
 
    if (_visible)
    {
 
        if (renderContext.type == TYRenderType::Picking)
        {
            TYPickingTable::addElement(getElement());
            glPushName((GLuint)(TYPickingTable::getIndex()));
        }
 
        switch (_Rendering)
        {
            case RENDER_SURF:
                displaySurface(renderContext);
                break;
            case RENDER_LINES:
                displayLines(renderContext);
                break;
            case RENDER_MIXED:
                displayMixed(renderContext);
                break;
            default:
                displaySurface(renderContext);
                break;
        }
 
        if (renderContext.type == TYRenderType::Picking)
        {
            glPopName();
        }
 
        if (_highlight)
        {
            GLint viewport[4];
            glGetIntegerv(GL_VIEWPORT, viewport);
            unsigned int width = viewport[2];
            unsigned int height = viewport[3];
 
            GLfloat pointSize = NAN;
            glGetFloatv(GL_POINT_SIZE, &pointSize);
            float tmpColor[4];
            glGetFloatv(GL_CURRENT_COLOR, tmpColor);
 
            if (renderContext.pass == TYRenderPass::Overlay)
            {
                drawName(renderContext);
            }
            glColor3f(tmpColor[0], tmpColor[1], tmpColor[2]);
 
            glPointSize(pointSize);
 
            // Affichage de la palette graphique
            TYElementGraphic* pTYElementGraphic = pMaillage->getPalette()->getGraphicObject();
            TYPaletteGraphic* pTYPaletteGraphic = (TYPaletteGraphic*)pTYElementGraphic;
            pTYPaletteGraphic->setPosX(width * _scalarX);
            pTYPaletteGraphic->setPosY(height * _scalarY);
            pTYPaletteGraphic->display(renderContext);
        }
 
        // Affichage du panel
        TYElementGraphic* pTYElementGraphic = pMaillage->getPanel()->getGraphicObject();
        TYPanelGraphic* pTYPanelGraphic = (TYPanelGraphic*)pTYElementGraphic;
        pTYPanelGraphic->display(renderContext);
    }
}
 
// Mode de Rendu Surface
 
void TYRectangularMaillageGraphic::displaySurface(TYRenderContext& renderContext)
{
    // Opacite par defaut
    float opacity = 1.0;
 
    // Anticrénelage par défaut
    bool antialiasing = true;
 
#if TY_USE_IHM
    //  static const char prefName[] = "MaillageOpacity";
 
    if (TYPreferenceManager::exists(TYDIRPREFERENCEMANAGER, "MaillageOpacity"))
    {
        opacity = TYPreferenceManager::getFloat(TYDIRPREFERENCEMANAGER, "MaillageOpacity");
    }
    else
    {
        TYPreferenceManager::setFloat(TYDIRPREFERENCEMANAGER, "MaillageOpacity", opacity);
    }
 
    if (TYPreferenceManager::exists(TYDIRPREFERENCEMANAGER, "Antialiasing"))
    {
        antialiasing = TYPreferenceManager::getBool(TYDIRPREFERENCEMANAGER, "Antialiasing");
    }
    else
    {
        TYPreferenceManager::setBool(TYDIRPREFERENCEMANAGER, "Antialiasing", antialiasing);
    }
 
#endif // TY_USE_IHM
 
    // TODO (issue xxx) : instead of a single loop it does
    // several loops (by batch of 1000 triangles)... for no reason at all?
    // Maybe an (old) issue with the number of primitives between the glBegin()
    // and glEnd()... but i've never heard of limit in this case...
    size_t bmin = 0;
    size_t bmax = 0;
    size_t inc = 0;
 
    size_t nbTriangles = _mesh.size();
 
    if (nbTriangles > 1000)
    {
        inc = 1000;
    }
    else
    {
        inc = nbTriangles;
    }
 
    bmin = 0;
    bmax += inc;
    bmax = (bmax <= nbTriangles ? bmax : nbTriangles);
 
    // Palette
    TYRectangularMaillage* pMaillage = getElement();
    LPTYPalette pPalette = pMaillage->getPalette();
 
    // Couleurs des sommets du triangle
    OColor color0 = OColor();
    OColor color1 = OColor();
    OColor color2 = OColor();
 
    if (antialiasing)
    {
        do
        {
            glBegin(GL_TRIANGLES);
            for (size_t i = bmin; i < bmax; i++)
            {
                const MTriangle& vtx = _mesh[i];
                const MPoint& v1 = vtx.pts[0];
                const MPoint& v2 = vtx.pts[1];
                const MPoint& v3 = vtx.pts[2];
 
                color0 = pPalette->getColorFromValue(v1.scalar);
                color1 = pPalette->getColorFromValue(v2.scalar);
                color2 = pPalette->getColorFromValue(v3.scalar);
 
                glColor4f(color0.r, color0.g, color0.b, opacity);
                glVertex3f(v1.pt._x, v1.pt._y, v1.pt._z);
 
                glColor4f(color1.r, color1.g, color1.b, opacity);
                glVertex3f(v2.pt._x, v2.pt._y, v2.pt._z);
 
                glColor4f(color2.r, color2.g, color2.b, opacity);
                glVertex3f(v3.pt._x, v3.pt._y, v3.pt._z);
            }
            glEnd();
 
            bmin = bmax;
            bmax += inc;
            bmax = (bmax <= nbTriangles ? bmax : nbTriangles);
        } while (bmin < nbTriangles);
    }
    else
    {
        glBegin(GL_QUADS);
        for (size_t i = 0; i < nbTriangles; i = i + 2)
        {
            const MTriangle& vtx1 = _mesh[i];
            const MTriangle& vtx2 = _mesh[i + 1];
            const MPoint& v1 = vtx1.pts[0];
            const MPoint& v2 = vtx1.pts[1];
            const MPoint& v3 = vtx2.pts[2];
            const MPoint& v4 = vtx1.pts[2];
 
            color0 = pPalette->getColorFromValue((v1.scalar + v2.scalar + v3.scalar + v4.scalar) * 0.25f);
 
            glColor4f(color0.r, color0.g, color0.b, opacity);
            glVertex3f(v1.pt._x, v1.pt._y, v1.pt._z);
 
            glColor4f(color0.r, color0.g, color0.b, opacity);
            glVertex3f(v2.pt._x, v2.pt._y, v2.pt._z);
 
            glColor4f(color0.r, color0.g, color0.b, opacity);
            glVertex3f(v3.pt._x, v3.pt._y, v3.pt._z);
 
            glColor4f(color0.r, color0.g, color0.b, opacity);
            glVertex3f(v4.pt._x, v4.pt._y, v4.pt._z);
        }
        glEnd();
    }
}
 
// Mode de Rendu Courbes
 
void TYRectangularMaillageGraphic::displayLines(TYRenderContext& renderContext, bool invertColors)
{
    // Opacite par defaut
    float opacity = 1.0;
 
#if TY_USE_IHM
    //  static const char prefName[] = "MaillageOpacity";
 
    if (TYPreferenceManager::exists(TYDIRPREFERENCEMANAGER, "MaillageOpacity"))
    {
        opacity = TYPreferenceManager::getFloat(TYDIRPREFERENCEMANAGER, "MaillageOpacity");
    }
    else
    {
        TYPreferenceManager::setFloat(TYDIRPREFERENCEMANAGER, "MaillageOpacity", opacity);
    }
#endif // TY_USE_IHM
 
    float lw = NAN;
    glGetFloatv(GL_LINE_WIDTH, &lw);
    glLineWidth(3.0f);
 
    std::vector<MPoint>::iterator it;
 
    // Palette
    TYRectangularMaillage* pMaillage = getElement();
    LPTYPalette pPalette = pMaillage->getPalette();
 
    glBegin(GL_LINES);
    for (it = _isoCurve.begin(); it != _isoCurve.end(); it++)
    {
        const MPoint& mp = *it;
        const OColor& color = pPalette->getColorFromValue(mp.scalar);
        if (invertColors)
        {
            glColor4f(1.0f - color.r, 1.0f - color.g, 1.0f - color.b, opacity);
        }
        else
        {
            glColor4f(color.r, color.g, color.b, opacity);
        }
        glVertex3f(mp.pt._x, mp.pt._y,
                   mp.pt._z + 0.01); // 0.01 offset so that we can display the lines from top view
    }
    glEnd();
 
    glLineWidth(lw);
}
 
void TYRectangularMaillageGraphic::displayMixed(TYRenderContext& renderContext)
{
    displaySurface(renderContext);
    displayLines(renderContext, true);
}
 
int TYRectangularMaillageGraphic::nextPoint(const TYTabLPPointCalcul* pPtsCalcul, const int& index1,
                                            const int& index2)
{
 
    int indice = index2;
 
    while ((pPtsCalcul->at(indice)->etat() == false) && (indice > 0))
    {
        indice--;
    }
 
    return indice;
}