OpenGL: создание 2D-текстуры из массива данных для отображения на квадроцикле

Я написал программу, которая выполняет вычисления конечной разности во временной области на двумерной сетке. Я хочу взять данные о величине одного из компонентов поля (скажем, z-компонент магнитного поля, Гц), преобразовать величины в цветовую шкалу, а затем отобразить окончательную цветовую карту на экране. Следующее изображение в ссылке imgur было сгенерировано с помощью Matlab, но я хотел бы получить тот же тип изображения в OpenGL.

http://imgur.com/8Rj77MK

Я использую GLFW для создания окна OpenGL (800x480) с ортогональной проекцией, используя следующий код:

//-->Construct GLFW Window<--//
int WINDOW_WIDTH = 800*(xIndex/4000);
int WINDOW_HEIGHT = 480*(yIndex/2400);

GLFWwindow* window;
window = glfwCreateWindow(WINDOW_WIDTH, WINDOW_HEIGHT, "Hz Component Magnitude", NULL, NULL);

//GLFW terminate and error function setup not shown

glfwMakeContextCurrent(window);
glfwSetKeyCallback(window, key_callback);

glewInit();

glfwGetFramebufferSize(window, &WINDOW_WIDTH, &WINDOW_HEIGHT);

glViewport(0, 0, WINDOW_WIDTH, WINDOW_HEIGHT);
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);

//Create orthographic view
glMatrixMode(GL_PROJECTION);
glOrtho(0, 0, -1.f, 1.f, 1.f, -1.f);
glMatrixMode(GL_MODELVIEW);
glEnable(GL_TEXTURE_2D);
glLoadIdentity();

glfwSwapBuffers(window);
glfwPollEvents();
//-->END Construct GLFW Window<--//

Данные о цвете для поля Hz сохраняются:

• Создание указателя 1D int длиной 800x480*(3 цвета)
• Вычислите и сохраните значения цвета поля Hz в указателе 1D

Из того, что я прочитал, типичным способом отображения 2D-текстуры в окне является создание четырехугольника и применение текстуры к четырехугольнику. Я могу успешно создать четырехугольник и успешно сгенерировать данные о цвете (проверено с помощью Matlab), но у меня возникают проблемы с созданием текстуры для применения к четырехугольнику.

Мой код для создания текстуры:

int num_colors = 3;    //Number of colors per pixel
int tex_size = num_colors*WINDOW_WIDTH*WINDOW_HEIGHT;    //Number of entries in tex_Hz

int size = tex_size*sizeof(int);

//Create pointer for storing color map
int* tex_Hz = (int*)malloc(size);

for(int i = 0; i < tex_size; i++)    //Initialize colormap to 0
    tex_Hz[i] = 0;

//OpenCL buffer for the color map
cl::Buffer texture_mem(context, CL_MEM_READ_WRITE | CL_MEM_COPY_HOST_PTR, size, tex_Hz, &error);

//Generate a texture ID
GLuint tex_id;

glGenTextures(1, &tex_id);
glBindTexture(GL_TEXTURE_2D, tex_id);

    //*A bunch of irrelevant code*//

    //-->Start of time loop<--//

    //*A bunch of irrelevant code*//

//-->Update the Hz field<--//
time_loop_queue.enqueueNDRangeKernel(timeLoopKernels[1], 0, global_size, local_size);
time_loop_queue.finish();

//-->Generate 2D image every 100 time steps<--//
if(step%100 == 0)
{
    //Read colormap values from OpenCL buffer
    time_loop_queue.enqueueReadBuffer(texture_mem, NULL, 0, size, tex_Hz, NULL, NULL);
    time_loop_queue.finish();

    //Used only as a debug runtime check to make sure data copied back
    cout << tex_Hz[3*800*240 + 3*350] << "\n";

    glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB, WINDOW_WIDTH, WINDOW_HEIGHT, 0, GL_RGB, GL_INT, tex_Hz);

    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
    glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);

    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
    glLoadIdentity();

    glBindTexture(GL_TEXTURE_2D, tex_id);

    //Draw quad and set texture coords
    glBegin(GL_QUADS);
    glTexCoord2d(0.0,0.0); glVertex2d(-1.0,-1.0);
    glTexCoord2d(1.0,0.0); glVertex2d(1.0,-1.0);
    glTexCoord2d(1.0,1.0); glVertex2d(1.0,1.0);
    glTexCoord2d(0.0,1.0); glVertex2d(-1.0,1.0);
    glEnd();

    glfwSwapBuffers(window);
}

Проблема заключается в том, что четырехугольник остается черным.

*** Я знаю, что копирование данных обратно в память хоста для последующего создания GL-текстуры не является оптимальным с точки зрения производительности, но на данный момент все, что имеет значение, — это реализация отображения поля в кратчайшие сроки, а не достижение оптимальной производительности (пока). .