glscope.cpp

// SPDX-License-Identifier: GPL-2.0+

#include <cmath>
#include <iostream>

#include <QColor>
#include <QCoreApplication>
#include <QDebug>
#include <QMatrix4x4>
#include <QMessageBox>
#include <QMouseEvent>
#include <QOpenGLShaderProgram>
#include <QPainter>

#include <QOffscreenSurface>
#include <QOpenGLFunctions>

#include "glscope.h"

#include "post/graphgenerator.h"
#include "post/ppresult.h"
#include "scopesettings.h"
#include "viewconstants.h"
#include "viewsettings.h"

GlScope *GlScope::createNormal( DsoSettingsScope *scope, DsoSettingsView *view, QWidget *parent ) {
    GlScope *s = new GlScope( scope, view, parent );
    s->zoomed = false;
    return s;
}

GlScope *GlScope::createZoomed( DsoSettingsScope *scope, DsoSettingsView *view, QWidget *parent ) {
    GlScope *s = new GlScope( scope, view, parent );
    s->zoomed = true;
    return s;
}

void GlScope::useQSurfaceFormat( QSurfaceFormat::RenderableType t ) {
    QCoreApplication::setAttribute( Qt::AA_ShareOpenGLContexts, true );

    // Prefer full desktop OpenGL without fixed pipeline
    QSurfaceFormat format;
    // TODO: check from time to time if warning has gone, comment next line for now
    format.setSamples( 4 ); // antialiasing gives warning with some HW, Qt & OpenGL versions.
    format.setProfile( QSurfaceFormat::CoreProfile );
    if ( t == QSurfaceFormat::OpenGLES ) {
        format.setVersion( 2, 0 );
        format.setRenderableType( QSurfaceFormat::OpenGLES );
        QCoreApplication::setAttribute( Qt::AA_UseOpenGLES, true );
    } else {
        format.setVersion( 3, 2 );
        format.setRenderableType( QSurfaceFormat::OpenGL );
    }
    QSurfaceFormat::setDefaultFormat( format );
}

GlScope::GlScope( DsoSettingsScope *scope, DsoSettingsView *view, QWidget *parent )
    : QOpenGLWidget( parent ), scope( scope ), view( view ) {
    // get OpenGL version to define appropriate OpenGLSL version
    // reason:
    // some not so new intel graphic driver report a very conservative version
    // e.g. debian buster -> "2.1 Mesa 18.3.4"
    QOffscreenSurface surface;
    surface.create();
    QOpenGLContext context;
    context.create();
    context.makeCurrent( &surface );
    QString glVersion = reinterpret_cast< const char * >( context.functions()->glGetString( GL_VERSION ) );
    GLSLversion = glVersion >= "3.2" ? 150 : 120; // version string "3.2 xxxx" > "3.2" is true
    // qDebug() << "OpenGL version" << glVersion << GLSLversion;
    surface.destroy();

    cursorInfo.clear();
    cursorInfo.push_back( &scope->horizontal.cursor );
    selectedCursor = 0;
    for ( ChannelID channel = 0; channel < scope->voltage.size(); ++channel ) {
        cursorInfo.push_back( &scope->voltage[ channel ].cursor );
    }
    for ( ChannelID channel = 0; channel < scope->spectrum.size(); ++channel ) {
        cursorInfo.push_back( &scope->spectrum[ channel ].cursor );
    }
    vaMarker.resize( cursorInfo.size() );
}

GlScope::~GlScope() { /* virtual destructor necessary */
}

QPointF GlScope::posToPosition( QPoint pos ) {
    QPointF position( double( pos.x() - width() / 2 ) * DIVS_TIME / double( width() ),
                      double( height() / 2 - pos.y() ) * DIVS_VOLTAGE / double( height() ) );
    if ( zoomed ) {
        double m1 = scope->getMarker( 0 );
        double m2 = scope->getMarker( 1 );
        if ( m1 > m2 )
            std::swap( m1, m2 );
        position.setX( m1 + ( 0.5 + ( position.x() / DIVS_TIME ) ) * ( m2 - m1 ) );
    }
    return position;
}

void GlScope::mousePressEvent( QMouseEvent *event ) {
    if ( !( zoomed && selectedCursor == 0 ) && event->button() == Qt::LeftButton ) {
        QPointF position = posToPosition( event->pos() );
        selectedMarker = NO_MARKER;
        DsoSettingsScopeCursor *cursor = cursorInfo[ selectedCursor ];
        // Capture nearest marker located within snap area (+/- 1% of full scale).
        double dX0 = fabs( cursor->pos[ 0 ].x() - position.x() );
        double dX1 = fabs( cursor->pos[ 1 ].x() - position.x() );
        double dY0 = fabs( cursor->pos[ 0 ].y() - position.y() );
        double dY1 = fabs( cursor->pos[ 1 ].y() - position.y() );

        switch ( cursor->shape ) {
        case DsoSettingsScopeCursor::RECTANGULAR:
            if ( std::min( dX0, dX1 ) < 1.0 / DIVS_SUB && std::min( dY0, dY1 ) < 1.0 / DIVS_SUB ) {
                // Do we need to swap Y-coords?
                if ( ( dX0 < dX1 && dY0 > dY1 ) || ( dX0 > dX1 && dY0 < dY1 ) ) {
                    std::swap( cursor->pos[ 0 ].ry(), cursor->pos[ 1 ].ry() );
                }
                selectedMarker = ( dX0 < dX1 ) ? 0 : 1;
            }
            break;
        case DsoSettingsScopeCursor::VERTICAL:
            if ( dX0 < dX1 ) {
                if ( dX0 < 1.0 / DIVS_SUB )
                    selectedMarker = 0;
            } else {
                if ( dX1 < 1.0 / DIVS_SUB )
                    selectedMarker = 1;
            }
            break;
        case DsoSettingsScopeCursor::HORIZONTAL:
            if ( dY0 < dY1 ) {
                if ( dY0 < 1.0 / DIVS_SUB )
                    selectedMarker = 0;
            } else {
                if ( dY1 < 1.0 / DIVS_SUB )
                    selectedMarker = 1;
            }
            break;
        case DsoSettingsScopeCursor::NONE:
            break;
        }
        if ( selectedMarker != NO_MARKER ) {
            cursorInfo[ selectedCursor ]->pos[ selectedMarker ] = position;
            if ( selectedCursor == 0 )
                emit markerMoved( selectedCursor, selectedMarker );
        }
    }
    event->accept();
}

void GlScope::mouseMoveEvent( QMouseEvent *event ) {
    if ( !( zoomed && selectedCursor == 0 ) && ( event->buttons() & Qt::LeftButton ) != 0 ) {
        QPointF position = posToPosition( event->pos() );
        if ( selectedMarker == NO_MARKER ) {
            // qDebug() << "mouseMoveEvent";
            // User started draging outside the snap area of any marker:
            // move all markers to current position and select last marker in the array.
            for ( unsigned marker = 0; marker < 2; ++marker ) {
                cursorInfo[ selectedCursor ]->pos[ marker ] = position;
                emit markerMoved( selectedCursor, marker );
                selectedMarker = marker;
            }
        } else if ( selectedMarker < 2 ) {
            cursorInfo[ selectedCursor ]->pos[ selectedMarker ] = position;
            emit markerMoved( selectedCursor, selectedMarker );
        }
    }
    event->accept();
}

void GlScope::mouseReleaseEvent( QMouseEvent *event ) {
    if ( !( zoomed && selectedCursor == 0 ) && event->button() == Qt::LeftButton ) {
        QPointF position = posToPosition( event->pos() );
        if ( selectedMarker < 2 ) {
            // qDebug() << "mouseReleaseEvent";
            cursorInfo[ selectedCursor ]->pos[ selectedMarker ] = position;
            emit markerMoved( selectedCursor, selectedMarker );
        }
        selectedMarker = NO_MARKER;
    }
    event->accept();
}

void GlScope::mouseDoubleClickEvent( QMouseEvent *event ) {
    if ( !( zoomed && selectedCursor == 0 ) && ( event->buttons() & Qt::LeftButton ) != 0 ) {
        // left double click positions two markers left and right of clicked pos with zoom=100
        QPointF position = posToPosition( event->pos() );
        if ( selectedMarker == NO_MARKER ) {
            // User double clicked outside the snap area of any marker
            QPointF p = QPointF( 0.5, 0 );       // 10x zoom
            if ( event->modifiers() & Qt::CTRL ) // 100x zoom
                p /= 10;
            if ( event->modifiers() & Qt::SHIFT ) // center at trigger position
                position = QPointF( 10 * scope->trigger.offset - 5, 0 );
            // move 1st marker left of current position.
            cursorInfo[ selectedCursor ]->pos[ 0 ] = position - p;
            emit markerMoved( selectedCursor, 0 );
            // move 2nd marker right of current position to make zoom=10 or 100.
            cursorInfo[ selectedCursor ]->pos[ 1 ] = position + p;
            emit markerMoved( selectedCursor, 1 );
            //  select no marker
            selectedMarker = NO_MARKER;
        }
    } else if ( !( zoomed && selectedCursor == 0 ) && ( event->buttons() & Qt::RightButton ) != 0 ) {
        // right double click moves all markers out of the way
        cursorInfo[ selectedCursor ]->pos[ 0 ] = QPointF( MARGIN_LEFT, 0 );
        cursorInfo[ selectedCursor ]->pos[ 1 ] = QPointF( MARGIN_RIGHT, 0 );
        emit markerMoved( selectedCursor, 0 );
        emit markerMoved( selectedCursor, 1 );
    }
    event->accept();
}


void GlScope::wheelEvent( QWheelEvent *event ) {
    static std::vector< int > zoomList = {1, 2, 5, 10, 20, 50, 100, 200, 500};
    if ( !( zoomed && selectedCursor == 0 ) ) {
        if ( selectedMarker == NO_MARKER ) {
            double step = event->angleDelta().y() / 1200.0; // one click = 0.1
            // qDebug() << "wheeelEvent" << selectedCursor << event->globalPos() << step;
            double &m1 = cursorInfo[ selectedCursor ]->pos[ 0 ].rx();
            double &m2 = cursorInfo[ selectedCursor ]->pos[ 1 ].rx();
            if ( m1 > m2 )
                std::swap( m1, m2 );
            double dm = m2 - m1;
            if ( event->modifiers() & Qt::CTRL ) { // zoom in/out
                if ( ( step > 0 && dm <= 1 ) || ( step < 0 && dm < 1 ) )
                    step *= 0.1;
                if ( dm >= 3 * step ) {
                    m1 += step;
                    m2 -= step;
                }
            } else {
                if ( step < 0 ) {                           // shift zoom range left ..
                    step = qMax( step, MARGIN_LEFT - m1 );  // .. until m1 == MARGIN_LEFT
                } else {                                    // shift zoom range right ..
                    step = qMin( step, MARGIN_RIGHT - m2 ); // .. until m2 == MARGIN_RIGHT
                }
                m1 += step;
                m2 += step;
            }
        }
        emit markerMoved( selectedCursor, 0 );
        emit markerMoved( selectedCursor, 1 );
        selectedMarker = NO_MARKER;
    }
    event->accept();
}


void GlScope::paintEvent( QPaintEvent *event ) {
    if ( shaderCompileSuccess ) {
        QOpenGLWidget::paintEvent( event );
    } else if ( !zoomed ) { // draw error message on normal view if OpenGL failed
        QPainter painter( this );
        painter.setRenderHint( QPainter::Antialiasing, true );
        QFont font = painter.font();
        font.setPointSize( 18 );
        painter.setFont( font );
        painter.drawText( rect(), Qt::AlignCenter | Qt::TextWordWrap, errorMessage );
    }
    event->accept(); // consume the event
}

unsigned GlScope::forceGLSLversion = 0;

void GlScope::initializeGL() {
    if ( !QOpenGLShaderProgram::hasOpenGLShaderPrograms( context() ) ) {
        errorMessage = tr( "System does not support OpenGL Shading Language (GLSL)" );
        return;
    }
    if ( m_program ) {
        qWarning() << tr( "OpenGL init called twice!" );
        return;
    }

    auto program = std::unique_ptr< QOpenGLShaderProgram >( new QOpenGLShaderProgram( context() ) );

    const char *vshaderES = R"(
          #version 100
          attribute highp vec3 vertex;
          uniform mat4 matrix;
          void main()
          {
              gl_Position = matrix * vec4(vertex, 1.0);
              gl_PointSize = 1.0;
          }
    )";
    const char *fshaderES = R"(
          #version 100
          uniform highp vec4 colour;
          void main() { gl_FragColor = colour; }
    )";

    const char *vshaderDesktop120 = R"(
          #version 120
          attribute highp vec3 vertex;
          uniform mat4 matrix;
          void main()
          {
              gl_Position = matrix * vec4(vertex, 1.0);
              gl_PointSize = 1.0;
          }
    )";
    const char *fshaderDesktop120 = R"(
          #version 120
          uniform highp vec4 colour;
          void main() { gl_FragColor = colour; }
    )";

    const char *vshaderDesktop150 = R"(
          #version 150
          in highp vec3 vertex;
          uniform mat4 matrix;
          void main()
          {
              gl_Position = matrix * vec4(vertex, 1.0);
              gl_PointSize = 1.0;
          }
    )";
    const char *fshaderDesktop150 = R"(
          #version 150
          uniform highp vec4 colour;
          out vec4 flatColor;
          void main() { flatColor = colour; }
    )";

    if ( GlScope::forceGLSLversion )
        GLSLversion = GlScope::forceGLSLversion;
    // qDebug() << "compile shaders" << GlScope::forceGLSLversion << GLSLversion;

    const char *vshaderDesktop = GLSLversion == 120 ? vshaderDesktop120 : vshaderDesktop150;
    const char *fshaderDesktop = GLSLversion == 120 ? fshaderDesktop120 : fshaderDesktop150;

    // Compile vertex shader
    bool usesOpenGL = QSurfaceFormat::defaultFormat().renderableType() == QSurfaceFormat::OpenGL;
    if ( !program->addShaderFromSourceCode( QOpenGLShader::Vertex, usesOpenGL ? vshaderDesktop : vshaderES ) ||
         !program->addShaderFromSourceCode( QOpenGLShader::Fragment, usesOpenGL ? fshaderDesktop : fshaderES ) ) {
        errorMessage = tr( "Failed to compile OpenGL shader programs.\n" ) + program->log();
        return;
    }

    // Link shader pipeline
    if ( !program->link() || !program->bind() ) {
        errorMessage = tr( "Failed to link/bind OpenGL shader programs.\n" ) + program->log();
        return;
    }

    vertexLocation = program->attributeLocation( "vertex" );
    matrixLocation = program->uniformLocation( "matrix" );
    colorLocation = program->uniformLocation( "colour" );

    if ( vertexLocation == -1 || colorLocation == -1 || matrixLocation == -1 ) {
        qWarning() << tr( "Failed to locate shader variable." );
        return;
    }

    program->bind();

    auto *gl = context()->functions();
    gl->glDisable( GL_DEPTH_TEST );
    gl->glEnable( GL_BLEND );
    // Enable depth buffer
    gl->glEnable( GL_DEPTH_TEST );

    // Enable back face culling
    gl->glEnable( GL_CULL_FACE );
    gl->glBlendFunc( GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA );

    QColor bg = view->colors->background;
    gl->glClearColor( GLfloat( bg.redF() ), GLfloat( bg.greenF() ), GLfloat( bg.blueF() ), GLfloat( bg.alphaF() ) );

    generateGrid( program.get() );

    {
        m_vaoMarker.create();
        QOpenGLVertexArrayObject::Binder b( &m_vaoMarker );
        m_marker.create();
        m_marker.bind();
        m_marker.setUsagePattern( QOpenGLBuffer::StaticDraw );
        m_marker.allocate( int( vaMarker.size() * sizeof( Vertices ) ) );
        program->enableAttributeArray( vertexLocation );
        program->setAttributeBuffer( vertexLocation, GL_FLOAT, 0, 3, 0 );
    }
    updateCursor();

    m_program = std::move( program );
    shaderCompileSuccess = true;
}

void GlScope::showData( std::shared_ptr< PPresult > newData ) {
    if ( !shaderCompileSuccess )
        return;
    makeCurrent();
    // Remove too much entries
    while ( view->digitalPhosphorDraws() < m_GraphHistory.size() )
        m_GraphHistory.pop_back();

    // Add if missing
    if ( view->digitalPhosphorDraws() > m_GraphHistory.size() ) {
        m_GraphHistory.resize( m_GraphHistory.size() + 1 );
    }

    // Move last item to front
    m_GraphHistory.splice( m_GraphHistory.begin(), m_GraphHistory, std::prev( m_GraphHistory.end() ) );

    // Add new entry
    m_GraphHistory.front().writeData( newData.get(), m_program.get(), vertexLocation );
    // doneCurrent();

    update();
}

void GlScope::generateVertices( unsigned marker, const DsoSettingsScopeCursor &cursor ) {
    const float Z_ORDER = 1.0f;
    switch ( cursor.shape ) {
    case DsoSettingsScopeCursor::NONE:
        vaMarker[ marker ] = {QVector3D( -DIVS_TIME, -DIVS_VOLTAGE, Z_ORDER ), QVector3D( -DIVS_TIME, DIVS_VOLTAGE, Z_ORDER ),
                              QVector3D( DIVS_TIME, DIVS_VOLTAGE, Z_ORDER ), QVector3D( DIVS_TIME, -DIVS_VOLTAGE, Z_ORDER )};
        break;
    case DsoSettingsScopeCursor::VERTICAL:
        vaMarker[ marker ] = {QVector3D( GLfloat( cursor.pos[ 0 ].x() ), -GLfloat( DIVS_VOLTAGE ), Z_ORDER ),
                              QVector3D( GLfloat( cursor.pos[ 0 ].x() ), GLfloat( DIVS_VOLTAGE ), Z_ORDER ),
                              QVector3D( GLfloat( cursor.pos[ 1 ].x() ), GLfloat( DIVS_VOLTAGE ), Z_ORDER ),
                              QVector3D( GLfloat( cursor.pos[ 1 ].x() ), -GLfloat( DIVS_VOLTAGE ), Z_ORDER )};
        break;
    case DsoSettingsScopeCursor::HORIZONTAL:
        vaMarker[ marker ] = {QVector3D( -GLfloat( DIVS_TIME ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER ),
                              QVector3D( GLfloat( DIVS_TIME ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER ),
                              QVector3D( GLfloat( DIVS_TIME ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER ),
                              QVector3D( -GLfloat( DIVS_TIME ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER )};
        break;
    case DsoSettingsScopeCursor::RECTANGULAR:
        if ( ( cursor.pos[ 1 ].x() - cursor.pos[ 0 ].x() ) * ( cursor.pos[ 1 ].y() - cursor.pos[ 0 ].y() ) > 0.0 ) {
            vaMarker[ marker ] = {QVector3D( GLfloat( cursor.pos[ 0 ].x() ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 1 ].x() ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 1 ].x() ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 0 ].x() ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER )};
        } else {
            vaMarker[ marker ] = {QVector3D( GLfloat( cursor.pos[ 0 ].x() ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 0 ].x() ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 1 ].x() ), GLfloat( cursor.pos[ 1 ].y() ), Z_ORDER ),
                                  QVector3D( GLfloat( cursor.pos[ 1 ].x() ), GLfloat( cursor.pos[ 0 ].y() ), Z_ORDER )};
        }
        break;
    }
}

void GlScope::updateCursor( unsigned index ) {
    if ( index > 0 ) {
        generateVertices( index, *cursorInfo[ index ] );
    } else
        for ( index = 0; index < cursorInfo.size(); ++index ) {
            generateVertices( index, *cursorInfo[ index ] );
        }
    // Write coordinates to GPU
    makeCurrent();
    m_marker.bind();
    m_marker.write( 0, vaMarker.data(), int( vaMarker.size() * sizeof( Vertices ) ) );
}

void GlScope::paintGL() {
    if ( !shaderCompileSuccess )
        return;

    auto *gl = context()->functions();

    QColor bg = view->colors->background;
    gl->glClearColor( GLfloat( bg.redF() ), GLfloat( bg.greenF() ), GLfloat( bg.blueF() ), GLfloat( bg.alphaF() ) );

    // Clear OpenGL buffer and configure settings
    // TODO Don't clear if view->digitalPhosphorDraws()>1
    gl->glClear( GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT );
    gl->glLineWidth( 1 );

    m_program->bind();

    // Apply zoom settings via matrix transformation
    if ( zoomed ) {
        QMatrix4x4 m;
        m.scale( QVector3D( GLfloat( DIVS_TIME ) / GLfloat( fabs( scope->getMarker( 1 ) - scope->getMarker( 0 ) ) ), 1.0f, 1.0f ) );
        m.translate( -GLfloat( scope->getMarker( 0 ) + scope->getMarker( 1 ) ) / 2, 0.0f, 0.0f );
        m_program->setUniformValue( matrixLocation, pmvMatrix * m );
    }

    drawMarkers();

    unsigned historyIndex = 0;
    for ( Graph &graph : m_GraphHistory ) {
        for ( ChannelID channel = 0; channel < scope->voltage.size(); ++channel ) {
            if ( scope->horizontal.format == Dso::GraphFormat::TY ) {
                drawSpectrumChannelGraph( channel, graph, int( historyIndex ) );
                if ( scope->histogram ) {
                    drawHistogramChannelGraph( channel, graph, int( historyIndex ) );
                }
            }
            drawVoltageChannelGraph( channel, graph, int( historyIndex ) );
        }
        ++historyIndex;
    }

    if ( zoomed ) {
        m_program->setUniformValue( matrixLocation, pmvMatrix );
    }

    drawGrid();
    m_program->release();
}

void GlScope::resizeGL( int width, int height ) {
    if ( !shaderCompileSuccess )
        return;
    auto *gl = context()->functions();
    gl->glViewport( 0, 0, GLint( width ), GLint( height ) );

    // Set axes to div-scale and apply correction for exact pixelization
    float pixelizationWidthCorrection = float( width ) / ( width - 1 );
    float pixelizationHeightCorrection = float( height ) / ( height - 1 );

    pmvMatrix.setToIdentity();
    pmvMatrix.ortho( -float( DIVS_TIME ) / 2.0f * pixelizationWidthCorrection,
                     float( DIVS_TIME ) / 2.0f * pixelizationWidthCorrection,
                     -float( DIVS_VOLTAGE ) / 2.0f * pixelizationHeightCorrection,
                     float( DIVS_VOLTAGE ) / 2.0f * pixelizationHeightCorrection, -1.0f, 1.0f );

    m_program->bind();
    m_program->setUniformValue( matrixLocation, pmvMatrix );
    m_program->release();
}

// draw 4 small crosses @ (x,y), (-x,y), (x,-y) and (-x,-y)
// section 0:grid, 1:axes, 2:border
void GlScope::draw4Cross( std::vector< QVector3D > &va, int section, float x, float y ) {
    const float d = 0.05f; // cross size
    for ( int xSign : {-1, 1} ) {
        for ( int ySign : {-1, 1} ) {
            gridDrawCounts[ section ] += 4;
            va.push_back( QVector3D( xSign * ( x - d ), ySign * y, 0 ) );
            va.push_back( QVector3D( xSign * ( x + d ), ySign * y, 0 ) );
            va.push_back( QVector3D( xSign * x, ySign * ( y - d ), 0 ) );
            va.push_back( QVector3D( xSign * x, ySign * ( y + d ), 0 ) );
        }
    }
}

void GlScope::generateGrid( QOpenGLShaderProgram *program ) {
    gridDrawCounts[ 0 ] = 0;
    gridDrawCounts[ 1 ] = 0;
    gridDrawCounts[ 2 ] = 0;

    m_grid.create();
    m_grid.bind();
    m_grid.setUsagePattern( QOpenGLBuffer::StaticDraw );

    std::vector< QVector3D > vaGrid;

    { // Bind draw vertical lines
        m_vaoGrid[ 0 ].create();
        QOpenGLVertexArrayObject::Binder b( &m_vaoGrid[ 0 ] );
        m_grid.bind();
        program->enableAttributeArray( vertexLocation );
        program->setAttributeBuffer( vertexLocation, GL_FLOAT, 0, 3, 0 );
    }

    // Draw vertical dot lines
    for ( int vDiv = 1; vDiv < DIVS_TIME / 2; ++vDiv ) {
        for ( int dot = 1; dot < DIVS_VOLTAGE / 2 * DIVS_SUB; ++dot ) {
            float dotPosition = float( dot ) / DIVS_SUB;
            gridDrawCounts[ 0 ] += 4;
            vaGrid.push_back( QVector3D( -vDiv, -dotPosition, 0 ) );
            vaGrid.push_back( QVector3D( -vDiv, dotPosition, 0 ) );
            vaGrid.push_back( QVector3D( vDiv, -dotPosition, 0 ) );
            vaGrid.push_back( QVector3D( vDiv, dotPosition, 0 ) );
        }
    }
    // Draw horizontal dot lines
    for ( int hDiv = 1; hDiv < DIVS_VOLTAGE / 2; ++hDiv ) {
        for ( int dot = 1; dot < DIVS_TIME / 2 * DIVS_SUB; ++dot ) {
            if ( dot % DIVS_SUB == 0 )
                continue; // Already done by vertical lines
            float dotPosition = float( dot ) / DIVS_SUB;
            gridDrawCounts[ 0 ] += 4;
            vaGrid.push_back( QVector3D( -dotPosition, -hDiv, 0 ) );
            vaGrid.push_back( QVector3D( dotPosition, -hDiv, 0 ) );
            vaGrid.push_back( QVector3D( -dotPosition, hDiv, 0 ) );
            vaGrid.push_back( QVector3D( dotPosition, hDiv, 0 ) );
        }
    }

    { // Bind draw axes
        m_vaoGrid[ 1 ].create();
        QOpenGLVertexArrayObject::Binder b( &m_vaoGrid[ 1 ] );
        m_grid.bind();
        program->enableAttributeArray( vertexLocation );
        program->setAttributeBuffer( vertexLocation, GL_FLOAT, int( vaGrid.size() * sizeof( QVector3D ) ), 3 );
    }

    // Axes
    // Horizontal axis
    gridDrawCounts[ 1 ] += 2;
    vaGrid.push_back( QVector3D( -DIVS_TIME / 2, 0, 0 ) );
    vaGrid.push_back( QVector3D( DIVS_TIME / 2, 0, 0 ) );
    // Vertical axis
    gridDrawCounts[ 1 ] += 2;
    vaGrid.push_back( QVector3D( 0, -DIVS_VOLTAGE / 2, 0 ) );
    vaGrid.push_back( QVector3D( 0, DIVS_VOLTAGE / 2, 0 ) );
    // Subdiv lines on horizontal axis
    for ( int line = 1; line < DIVS_TIME / 2 * DIVS_SUB; ++line ) {
        float linePosition = float( line ) / DIVS_SUB;
        gridDrawCounts[ 1 ] += 4;
        vaGrid.push_back( QVector3D( linePosition, -0.05f, 0 ) );
        vaGrid.push_back( QVector3D( linePosition, 0.05f, 0 ) );
        vaGrid.push_back( QVector3D( -linePosition, -0.05f, 0 ) );
        vaGrid.push_back( QVector3D( -linePosition, 0.05f, 0 ) );
    }
    // Subdiv lines on vertical axis
    for ( int line = 1; line < DIVS_VOLTAGE / 2 * DIVS_SUB; ++line ) {
        float linePosition = float( line ) / DIVS_SUB;
        gridDrawCounts[ 1 ] += 4;
        vaGrid.push_back( QVector3D( -0.05f, linePosition, 0 ) );
        vaGrid.push_back( QVector3D( 0.05f, linePosition, 0 ) );
        vaGrid.push_back( QVector3D( -0.05f, -linePosition, 0 ) );
        vaGrid.push_back( QVector3D( 0.05f, -linePosition, 0 ) );
    }

    // Draw vertical cross lines
    for ( int vDiv = 1; vDiv < DIVS_TIME / 2; ++vDiv ) {
        for ( int hDiv = 1; hDiv < DIVS_VOLTAGE / 2; ++hDiv ) {
            draw4Cross( vaGrid, 1, vDiv, hDiv );
        }
    }
    // Draw horizontal cross lines
    for ( int hDiv = 1; hDiv < DIVS_VOLTAGE / 2; ++hDiv ) {
        for ( int vDiv = 1; vDiv < DIVS_TIME / 2; ++vDiv ) {
            if ( vDiv % DIVS_SUB == 0 )
                continue; // Already done by vertical lines
            draw4Cross( vaGrid, 1, vDiv, hDiv );
        }
    }

    {
        m_vaoGrid[ 2 ].create();
        QOpenGLVertexArrayObject::Binder b( &m_vaoGrid[ 2 ] );
        m_grid.bind();
        program->enableAttributeArray( vertexLocation );
        program->setAttributeBuffer( vertexLocation, GL_FLOAT, int( vaGrid.size() * sizeof( QVector3D ) ), 3 );
    }

    // Border
    gridDrawCounts[ 2 ] += 4;
    vaGrid.push_back( QVector3D( -DIVS_TIME / 2, -DIVS_VOLTAGE / 2, 0 ) );
    vaGrid.push_back( QVector3D( DIVS_TIME / 2, -DIVS_VOLTAGE / 2, 0 ) );
    vaGrid.push_back( QVector3D( DIVS_TIME / 2, DIVS_VOLTAGE / 2, 0 ) );
    vaGrid.push_back( QVector3D( -DIVS_TIME / 2, DIVS_VOLTAGE / 2, 0 ) );

    m_grid.allocate( &vaGrid[ 0 ], int( vaGrid.size() * sizeof( QVector3D ) ) );
    m_grid.release();
}

void GlScope::drawGrid() {
    auto *gl = context()->functions();
    gl->glLineWidth( 1 );

    // Grid
    m_vaoGrid[ 0 ].bind();
    m_program->setUniformValue( colorLocation, view->colors->grid );
    gl->glDrawArrays( GL_POINTS, 0, gridDrawCounts[ 0 ] );
    m_vaoGrid[ 0 ].release();

    // Axes and div crosses
    m_vaoGrid[ 1 ].bind();
    m_program->setUniformValue( colorLocation, view->colors->axes );
    gl->glDrawArrays( GL_LINES, 0, gridDrawCounts[ 1 ] );
    m_vaoGrid[ 1 ].release();

    // Border
    m_vaoGrid[ 2 ].bind();
    m_program->setUniformValue( colorLocation, view->colors->border );
    gl->glDrawArrays( GL_LINE_LOOP, 0, gridDrawCounts[ 2 ] );
    m_vaoGrid[ 2 ].release();
}

void GlScope::drawVertices( QOpenGLFunctions *gl, unsigned marker, QColor color ) {
    m_program->setUniformValue( colorLocation, ( marker == selectedCursor ) ? color : color.darker() );
    gl->glDrawArrays( GL_LINE_LOOP, GLint( marker * VERTICES_ARRAY_SIZE ), GLint( VERTICES_ARRAY_SIZE ) );
    if ( cursorInfo[ marker ]->shape == DsoSettingsScopeCursor::RECTANGULAR ) {
        color.setAlphaF( 0.25 );
        m_program->setUniformValue( colorLocation, color.darker() );
        gl->glDrawArrays( GL_TRIANGLE_FAN, GLint( marker * VERTICES_ARRAY_SIZE ), GLint( VERTICES_ARRAY_SIZE ) );
    }
}

void GlScope::drawMarkers() {
    auto *gl = context()->functions();

    m_vaoMarker.bind();

    unsigned marker = 0;
    drawVertices( gl, marker, view->colors->markers );
    ++marker;

    if ( view->cursorsVisible ) {
        gl->glDepthMask( GL_FALSE );
        for ( ChannelID channel = 0; channel < scope->voltage.size(); ++channel, ++marker ) {
            if ( scope->voltage[ channel ].used ) {
                drawVertices( gl, marker, view->colors->voltage[ channel ] );
            }
        }
        for ( ChannelID channel = 0; channel < scope->spectrum.size(); ++channel, ++marker ) {
            if ( scope->spectrum[ channel ].used ) {
                drawVertices( gl, marker, view->colors->spectrum[ channel ] );
            }
        }
        gl->glDepthMask( GL_TRUE );
    }

    m_vaoMarker.release();
}

void GlScope::drawVoltageChannelGraph( ChannelID channel, Graph &graph, int historyIndex ) {
    if ( !scope->voltage[ channel ].used )
        return;

    m_program->setUniformValue( colorLocation, view->colors->voltage[ channel ].darker( 100 + 10 * historyIndex ) );
    Graph::VaoCount &v = graph.vaoVoltage[ channel ];

    QOpenGLVertexArrayObject::Binder b( v.first );
    const GLenum dMode = ( view->interpolation == Dso::INTERPOLATION_OFF ) ? GL_POINTS : GL_LINE_STRIP;
    context()->functions()->glDrawArrays( dMode, 0, v.second );
}

void GlScope::drawHistogramChannelGraph( ChannelID channel, Graph &graph, int historyIndex ) {
    if ( !scope->voltage[ channel ].used )
        return;

    m_program->setUniformValue( colorLocation, view->colors->voltage[ channel ].darker( 100 + 10 * historyIndex ) );
    Graph::VaoCount &h = graph.vaoHistogram[ channel ];

    QOpenGLVertexArrayObject::Binder b( h.first );
    const GLenum dMode = ( view->interpolation == Dso::INTERPOLATION_OFF ) ? GL_POINTS : GL_LINES;
    context()->functions()->glDrawArrays( dMode, 0, h.second );
}

void GlScope::drawSpectrumChannelGraph( ChannelID channel, Graph &graph, int historyIndex ) {
    if ( !scope->spectrum[ channel ].used )
        return;

    m_program->setUniformValue( colorLocation, view->colors->spectrum[ channel ].darker( 100 + 10 * historyIndex ) );
    Graph::VaoCount &v = graph.vaoSpectrum[ channel ];

    QOpenGLVertexArrayObject::Binder b( v.first );
    const GLenum dMode = ( view->interpolation == Dso::INTERPOLATION_OFF ) ? GL_POINTS : GL_LINE_STRIP;
    context()->functions()->glDrawArrays( dMode, 0, v.second );
}