#include "gluos.h"
#include <stdlib.h>
#include <stdio.h>
#include "glimports.h"
#include "zlassert.h"
#include <GL/gl.h>
#include "rectBlock.h"
rectBlock::rectBlock(gridBoundaryChain* left, gridBoundaryChain* right, Int beginVline, Int endVline)
{
Int i;
upGridLineIndex = left->getVlineIndex(beginVline);
lowGridLineIndex = left->getVlineIndex(endVline);
Int n = upGridLineIndex-lowGridLineIndex+1; leftIndices = (Int*) malloc(sizeof(Int) * n);
assert(leftIndices);
rightIndices = (Int*) malloc(sizeof(Int) * n);
assert(rightIndices);
for(i=0; i<n; i++)
{
leftIndices[i] = left->getInnerIndex(i+beginVline);
rightIndices[i] = right->getInnerIndex(i+beginVline);
}
}
rectBlock::~rectBlock()
{
free(leftIndices);
free(rightIndices);
}
void rectBlock::print()
{
Int i;
printf("block:\n");
for(i=upGridLineIndex; i >= lowGridLineIndex; i--)
{
printf("gridline %i, (%i,%i)\n", i, leftIndices[upGridLineIndex-i], rightIndices[upGridLineIndex-i]);
}
}
void rectBlock::draw(Real* u_values, Real* v_values)
{
Int i,j,k;
#ifdef DEBUG
printf("upGridLineIndex=%i, lowGridLineIndex=%i\n", upGridLineIndex, lowGridLineIndex);
#endif
for(k=0, i=upGridLineIndex; i > lowGridLineIndex; i--, k++)
{
glBegin(GL_QUAD_STRIP);
for(j=leftIndices[k+1]; j<= rightIndices[k+1]; j++)
{
glVertex2f(u_values[j], v_values[i]);
glVertex2f(u_values[j], v_values[i-1]);
}
glEnd();
}
}
Int rectBlock::num_quads()
{
Int ret=0;
Int k,i;
for(k=0, i=upGridLineIndex; i>lowGridLineIndex; i--, k++)
{
ret += (rightIndices[k+1]-leftIndices[k+1]);
}
return ret;
}
Int rectBlockArray::num_quads()
{
Int ret=0;
for(Int i=0; i<n_elements; i++)
ret += array[i]->num_quads();
return ret;
}
rectBlockArray::rectBlockArray(Int s)
{
Int i;
n_elements = 0;
size = s;
array = (rectBlock**) malloc(sizeof(rectBlock*) * s);
assert(array);
for(i=0; i<s; i++)
array[i] = NULL;
}
rectBlockArray::~rectBlockArray()
{
Int i;
for(i=0; i<size; i++)
{
if(array[i] != NULL)
delete array[i];
}
free(array);
}
void rectBlockArray::insert(rectBlock* newBlock)
{
Int i;
if(n_elements == size) {
rectBlock** temp = (rectBlock**) malloc(sizeof(rectBlock) * (2*size+1));
assert(temp);
for(i=0; i<2*size+1; i++)
temp[i] = NULL;
for(i=0; i<n_elements; i++)
temp[i] = array[i];
free(array);
array = temp;
size = 2*size + 1;
}
array[n_elements++] = newBlock;
}
void rectBlockArray::print()
{
Int i;
for(i=0; i<n_elements; i++)
array[i]->print();
}
void rectBlockArray::draw(Real* u_values, Real* v_values)
{
Int i;
for(i=0; i<n_elements; i++)
array[i]->draw(u_values, v_values);
}