meshrender.cc

MatlabVisualization Package

#include "meshrender.h"
#include <GL/glu.h>
#include <iostream>
#include <math.h>

MeshRender::MeshRender(int w, int h, int rank, double *data, int nrows, 
		       int ncols, DistLayout *data_layout, double angleX, 
		       double angleZ, double camera_zoom) : 
  TDRender(w, h, rank, data, nrows, ncols, data_layout, 
	   angleX, angleZ, camera_zoom) {

}

bool MeshRender::draw_help(int offset) {

  offset = offset * 6;

  lx = blocks_data[offset + 0];
  ly = blocks_data[offset + 1];
  data_nr = blocks_data[offset + 2] - lx + 1;
  data_nc = blocks_data[offset + 3] - ly + 1;

  if (data_nr <= 1 &&
      data_nc <= 1)
    return true;

  int rOffset = blocks_data[offset + 4];
  int cOffset = blocks_data[offset + 5];

  glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);

  draw_surface(rOffset, cOffset);

  return false;
 
}

void MeshRender::draw_surface(int rOffset, int cOffset) {

  double x_const = (double) axes_length / (double) (nc + 1);
  double y_const = (double) axes_length / (double) (nr + 1);
  double z_interval = maxZ - minZ;
  double z_const;
  if (z_interval != 0) 
    z_const = 0.8 * (double) axes_length / z_interval;
  else 
    z_const = 0.;

  double color_offset = (NUM_COLORS - 1)/ z_interval;
  double z1, z2, z3, z4;
  int row, col, color_index;
  
  glBegin(GL_QUAD_STRIP);

  row = rOffset;
  col = cOffset;
  
  z1 = data_buffer[col * lRows + row];
  z2 = data_buffer[(col + data_nc - 1) * lRows + row];
  z3 = data_buffer[col * lRows + row + data_nr - 1];
  z4 = data_buffer[(col + data_nc - 1) * lRows + row + data_nr - 1];
      
  color_index = FIRST_COLOR + (int) ((z1 - minZ) * color_offset);
  assert(color_index < FIRST_COLOR + NUM_COLORS);
  glColor3ub(color_index, 0, 0);
	  
  glVertex3f((ly + 1) * x_const, (lx + 1) * y_const, 
	     (z1 - minZ) * z_const);
  glVertex3f((ly + data_nc + 1) * x_const, (lx + 1) * y_const, 
	     (z2 - minZ) * z_const);
  glVertex3f((ly + 1) * x_const, (lx + data_nr + 1) * y_const,
	     (z3 - minZ) * z_const);
  glVertex3f((ly + data_nc + 1) * x_const, (lx + data_nr + 1) * y_const,
	     (z4 - minZ) * z_const);
  
  glEnd();

}