#include <fcntl.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <gl/gl.h>
#include <gl/device.h>
#include <midi.h>
#include <midiio.h>
#include <malloc.h>
#include <unistd.h>
#include <sys/types.h>

/*** (^_^;) Bug Recovery of SGI Media Library --> Original MIDI Defines !! (^_^;) ***/
#define MES(x)			( (x[0]).mm.msgbuf )
#define midi_get_status(x)	( ( ( ( MES(x) & 0xf0000000 ) == 0x60000000 ) || ( ( MES(x) & 0xf0000000 ) == 0x40000000 ) ) ? ( ( MES(x) & 0x00f00000 ) >> 16 ) : 0 ) /* Port #2 */
#define midi_get_channel(x)	( ( MES(x) & 0x000f0000 ) >> 16 ) /* 1 - 16 channel Only ! */
#define midi_get_keyno(x)	( ( MES(x) & 0x00007f00 ) >> 8 ) /* device ID : #1 Only ! */
#define midi_get_velocity(x)	( ( MES(x) & 0x0000007f ) )


#define SCREEN_X_MAXSIZE 950
#define SCREEN_Y_MAXSIZE 700
/*
#define SCREEN_X_MAXSIZE 1200
#define SCREEN_Y_MAXSIZE 900
*/
#define R 500.0
#define XYMAX 50
#define XZMAX 50
#define PI 3.1415926535
#define DEL 0.01

float now[XYMAX][XZMAX][4][3],next[XYMAX][XZMAX][4][3];
int nowc[XYMAX][XZMAX][3],nextc[XYMAX][XZMAX][3],count=1,cg_param[20],cg_cnt;

open_window(){
	prefposition(100,SCREEN_X_MAXSIZE,100,SCREEN_Y_MAXSIZE); winopen("CG and Sound [Y.Nagashima] : with [Morph y.yura 1996]");
	RGBmode(); RGBcolor(0,0,0); clear(); doublebuffer(); qdevice(ESCKEY); gconfig(); blendfunction(BF_ONE, BF_MDC);
	doublebuffer(); mmode(MVIEWING); perspective(900,1.33,-100,0);
}

draw_plane(float pos[XYMAX][XZMAX][4][3], int col[XYMAX][XZMAX][3]){
	int i,j,k;
	float global_size;
	global_size = (float)cg_param[0] / 100.0; scale(global_size,global_size,global_size);
	pushmatrix();
	for(i=0; i<cg_param[3]; i++){
		for(j=0; j<cg_param[4]; j++){
			RGBcolor(col[i][j][0]+cg_param[5]*cos((double)(count-i*50)/150*j),col[i][j][1]+cg_param[6]*cos((double)(count-i*50)/150*j),col[i][j][2]+150.0*cos((double)(count-i*50)/150*j));
			bgnpolygon(); for(k=0; k<4; k++){ v3f(pos[i][j][k]); } endpolygon();
		}
	}
	popmatrix();
}

yota(){
	int xrr, yrr, zrr;
	lookat(cg_param[1],cg_param[1],cg_param[1],0,0,0,0);
	xrr = 500.0 * cos((double)count/100.0); yrr = 500.0 * sin((double)count/100.0); zrr = 500.0 * cos((double)count/100.0);
	rotate (xrr,'x'); rotate (yrr,'y'); rotate (zrr,'z');
}

cg1(float pos1[XYMAX][XZMAX][4][3], int col1[XYMAX][XZMAX][3]){
	int j,k;
	float jj,kk,xy,xz;
	xy=1000.0/(float)cg_param[3]; xz=1000.0/(float)cg_param[4];
	for (j=0; j<cg_param[3]; j++){
		jj = (float)j;
		for (k=0; k<cg_param[4]; k++){
			kk = (float)k;
			pos1[j][k][0][0] = kk*xz-500.0; pos1[j][k][0][1] = 0.0; pos1[j][k][0][2] = jj*xy-500.0;
			pos1[j][k][1][0] = kk*xz-(500.0-xz*0.8); pos1[j][k][1][1] = 0.0; pos1[j][k][1][2] = jj*xy-500.0;
			pos1[j][k][2][0] = kk*xz-(500.0-xz*0.8); pos1[j][k][2][1] = 0.0; pos1[j][k][2][2] = jj*xy-(500.0-xy*0.8);
			pos1[j][k][3][0] = kk*xz-500.0; pos1[j][k][3][1] = 0.0; pos1[j][k][3][2] = jj*xy-(500.0-xy*0.8);
			col1[j][k][0] = 100; col1[j][k][1] = 100; col1[j][k][2] = 100;
		}		
	}	
}

cg2(float pos2[XYMAX][XZMAX][4][3], int col2[XYMAX][XZMAX][3]){
	int i,j;
	double p1,p2,p3,p4,q1,q2,q3,q4;
	for(i=0; i<cg_param[3]; i++){
		p1 = R * cos( PI * ( 0.5 - (double)(i*2)/(double)cg_param[3] - pow( -1.0 , (double)(i*2) ) * DEL ) );
		p2 = R * sin( PI * ( 0.5 - (double)(i*2)/(double)cg_param[3] - pow( -1.0 , (double)(i*2) ) * DEL ) );
		p3 = R * cos( PI * ( 0.5 - (double)(i*2+1)/(double)cg_param[3] - pow( -1.0 , (double)(i*2+1) ) * DEL ) );
		p4 = R * sin( PI * ( 0.5 - (double)(i*2+1)/(double)cg_param[3] - pow( -1.0 , (double)(i*2+1) ) * DEL ) );
		for(j=0; j<cg_param[4]; j++){
			q1 = 2.0 * (double)j / (double)cg_param[4]; q2 = pow( -1.0 , (double)(j*2) ) * DEL;
			q3 = 2.0 * (double)(j+1) / (double)cg_param[4]; q4 = pow( -1.0 , (double)(j*2+1) ) * DEL;
			pos2[i][j][0][0] = p1 * cos( PI * ( q1 + q2 ) ); pos2[i][j][0][1] = p2; pos2[i][j][0][2] = p1 * sin( PI * ( q1 + q2 ) );
			pos2[i][j][1][0] = p1 * cos( PI * ( q3 + q4 ) ); pos2[i][j][1][1] = p2; pos2[i][j][1][2] = p1 * sin( PI * ( q3 + q4 ) );
			pos2[i][j][2][0] = p3 * cos( PI * ( q3 + q4 ) ); pos2[i][j][2][1] = p4; pos2[i][j][2][2] = p3 * sin( PI * ( q3 + q4 ) );
			pos2[i][j][3][0] = p3 * cos( PI * ( q1 + q2 ) ); pos2[i][j][3][1] = p4; pos2[i][j][3][2] = p3 * sin( PI * ( q1 + q2 ) );
			col2[i][j][0] = j*2*i; col2[i][j][1] = i*j/5; col2[i][j][2] = i*j/15;
		}
	}
}

cg3(float pos3[XYMAX][XZMAX][4][3], int col3[XYMAX][XZMAX][3]){	
	int i,j;
	float xy, p1, p2, q1, q2, q3, q4, bet = 0.8;
	xy = 800.0 / (float)cg_param[3];
	for (i=0; i<cg_param[3]; i++){
		p1 = (float)i * xy - 400.0; p2 = (float)i * xy - ( 400.0 - xy * bet );
		for (j=0; j<cg_param[4]; j++){
			q1 = 2.0 * (double)j / (double)cg_param[4]; q2 = pow( -1.0 , (double)(j*2) ) * DEL;
			q3 = 2.0 * (double)(j+1) / (double)cg_param[4]; q4 = pow( -1.0 , (double)(j*2+1) ) * DEL;
			pos3[i][j][0][0] = R * 0.8 * cos( PI * ( q1 + q2 ) ); pos3[i][j][0][1] = p1; pos3[i][j][0][2] = R * 0.8 * sin( PI * ( q1 + q2 ) );
			pos3[i][j][1][0] = R * 0.8 * cos( PI * ( q3 + q4 ) ); pos3[i][j][1][1] = p1; pos3[i][j][1][2] = R * 0.8 * sin( PI * ( q3 + q4 ) );
			pos3[i][j][2][0] = R * 0.8 * cos( PI * ( q3 + q4 ) ); pos3[i][j][2][1] = p2; pos3[i][j][2][2] = R * 0.8 * sin( PI * ( q3 + q4 ) );
			pos3[i][j][3][0] = R * 0.8 * cos( PI * ( q1 + q2 ) ); pos3[i][j][3][1] = p2; pos3[i][j][3][2] = R * 0.8 * sin( PI * ( q1 + q2 ) );
			col3[i][j][0] = 250.0*cos((double)(count-i*50)/(150+j)); col3[i][j][1] = 150.0*cos((double)(count+i*50)/(200+j));
			col3[i][j][2] = 150.0*cos((double)(count-i*50)/(150+j));
		}
	}
}

cg4(float pos4[XYMAX][XZMAX][4][3], int col4[XYMAX][XZMAX][3]){
	int i,j,k;
	float p1, p2, q1, q2, q3, q4, xy, bet = 0.8, hole = 0.5;
	xy = 800.0 / ( (float)cg_param[3] - 2.0 );
	for (i=0; i<cg_param[3]; i++){
		p1 = (float)i * xy - 400.0; p2 = (float)i * xy - ( 400.0 - xy * bet );
		for (j=0; j<cg_param[4]; j++){
			q1 = 2.0 * (double)j / (double)cg_param[4]; q2 = pow( -1.0 , (double)(j*2) ) * DEL;
			q3 = 2.0 * (double)(j+1) / (double)cg_param[4]; q4 = pow( -1.0 , (double)(j*2+1) ) * DEL;
			if(i==0){
				pos4[i][j][0][0] = hole * R * 0.8 * cos( PI * ( q1 + q2 ) ); pos4[i][j][0][1] = p2; pos4[i][j][0][2] = hole * R * 0.8 * sin( PI * ( q1 + q2 ) );
	 			pos4[i][j][1][0] = hole * R * 0.8 * cos( PI * ( q3 + q4 ) ); pos4[i][j][1][1] = p2; pos4[i][j][1][2] = hole * R * 0.8 * sin( PI * ( q3 + q4 ) );
			}
			else{
				pos4[i][j][0][0] = R * 0.8 * cos( PI * ( q1 + q2 ) ); pos4[i][j][0][1] = p1; pos4[i][j][0][2] = R * 0.8 * sin( PI * ( q1 + q2 ) );
				pos4[i][j][1][0] = R * 0.8 * cos( PI * ( q3 + q4 ) ); pos4[i][j][1][1] = p1; pos4[i][j][1][2] = R * 0.8 * sin( PI * ( q3 + q4 ) );
			}
			if(i==cg_param[3]-1){
				pos4[i][j][2][0] = hole * R * 0.8 * cos( PI * ( q3 + q4 ) ); pos4[i][j][2][1] = p1; pos4[i][j][2][2] = hole * R * 0.8 * sin( PI * ( q3 + q4 ) );
				pos4[i][j][3][0] = hole * R * 0.8 * cos( PI * ( q1 + q2 ) ); pos4[i][j][3][1] = p1; pos4[i][j][3][2] = hole * R * 0.8 * sin( PI * ( q1 + q2 ) );
			}
			else{
				pos4[i][j][2][0] = R * 0.8 * cos( PI * ( q3 + q4 ) ); pos4[i][j][2][1] = p2; pos4[i][j][2][2] = R * 0.8 * sin( PI * ( q3 + q4 ) );
				pos4[i][j][3][0] = R * 0.8 * cos( PI * ( q1 + q2 ) ); pos4[i][j][3][1] = p2; pos4[i][j][3][2] = R * 0.8 * sin( PI * ( q1 + q2 ) );
			}
			col4[i][j][0] = 100.0*cos((double)(count-i*200)/(200+j)); col4[i][j][1] = 250.0*cos((double)(count+i*100)/(150+j));
			col4[i][j][2] = 250.0*cos((double)(count-i*150)/(250+j));
		}
	}
}

cg5(float pos5[XYMAX][XZMAX][4][3], int col5[XYMAX][XZMAX][3]){
	int j,k;
	float jj,kk,xy,xz,k1,k2,k3;
	xy = 700.0 / (float)cg_param[3]; xz = 2800.0 / (float)cg_param[4];
	k1 = -0.25 * (float)cg_param[4]; k2 = -0.5 * (float)cg_param[4]; k3 = -0.75 * (float)cg_param[4];
	for (j=0; j<cg_param[3]; j++){
		jj = (float)j;
		for (k=0; k<cg_param[4]; k++){
			kk = (float)k;
			if(k<cg_param[4]/4){
				pos5[j][k][0][0] = kk*xz-350.0; pos5[j][k][0][1] = 350.0; pos5[j][k][0][2] = jj*xy-350.0;
				pos5[j][k][1][0] = kk*xz-(350.0-xz*0.8); pos5[j][k][1][1] = 350.0; pos5[j][k][1][2] = jj*xy-350.0;
				pos5[j][k][2][0] = kk*xz-(350.0-xz*0.8); pos5[j][k][2][1] = 350.0; pos5[j][k][2][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][3][0] = kk*xz-350.0; pos5[j][k][3][1] = 350.0; pos5[j][k][3][2] = jj*xy-(350.0-xy*0.8);
			}
			else if(k<cg_param[4]/2){
				pos5[j][k][0][0] = (kk+k1)*xz-350.0; pos5[j][k][0][1] = -350.0; pos5[j][k][0][2] = jj*xy-350.0;
				pos5[j][k][1][0] = (kk+k1)*xz-(350.0-xz*0.8); pos5[j][k][1][1] = -350.0; pos5[j][k][1][2] = jj*xy-350.0;
				pos5[j][k][2][0] = (kk+k1)*xz-(350.0-xz*0.8); pos5[j][k][2][1] = -350.0; pos5[j][k][2][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][3][0] = (kk+k1)*xz-350.0; pos5[j][k][3][1] = -350.0; pos5[j][k][3][2] = jj*xy-(350.0-xy*0.8);
			}
			else if(k<3*cg_param[4]/4){
				pos5[j][k][0][0] = 350.0; pos5[j][k][0][1] = (kk+k2)*xz-350.0; pos5[j][k][0][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][1][0] = 350; pos5[j][k][1][1] = (kk+k2)*xz-(350.0-xz*0.8); pos5[j][k][1][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][2][0] = 350.0; pos5[j][k][2][1] = (kk+k2)*xz-(350.0-xz*0.8); pos5[j][k][2][2] = jj*xy-350.0;
				pos5[j][k][3][0] = 350.0; pos5[j][k][3][1] = (kk+k2)*xz-350.0; pos5[j][k][3][2] = jj*xy-350.0;
			}
			else{
				pos5[j][k][0][0] = -350.0; pos5[j][k][0][1] = (kk+k3)*xz-350.0; pos5[j][k][0][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][1][0] = -350.0; pos5[j][k][1][1] = (kk+k3)*xz-(350.0-xz*0.8); pos5[j][k][1][2] = jj*xy-(350.0-xy*0.8);
				pos5[j][k][2][0] = -350.0; pos5[j][k][2][1] = (kk+k3)*xz-(350.0-xz*0.8); pos5[j][k][2][2] = jj*xy-350.0;
				pos5[j][k][3][0] = -350.0; pos5[j][k][3][1] = (kk+k3)*xz-350.0; pos5[j][k][3][2] = jj*xy-350.0;
			}
			col5[j][k][0] = 100.0*cos((double)(count-j*200)/(200+k)); col5[j][k][1] = 250.0*cos((double)(count+j*100)/(150+k));
			col5[j][k][2] = 250.0*cos((double)(count-j*150)/(250+k));
		}
	}
}

cg6(float pos6[XYMAX][XZMAX][4][3], int col6[XYMAX][XZMAX][3]){
	int i,j;
	float p1, p2, p3, q1, q2, q3, q4, q5, r1, r2, r3, r4, ra, xy, bet = 0.8, hole = 0.7;
	xy = 800.0 / ( (float)cg_param[3] - 2.0 ); ra = R * 1.3 / ( (float)cg_param[3] * 2 );
	for (i=0; i<cg_param[3]; i++){
		p1 = (float)i * xy - 500.0; p2 = (float)i * xy - ( 500.0 - xy * bet );
		q1 = ra * (float)( i * 2 ); q2 = ra * hole; q3 = ra * (float)( i * 2 + 1 );
		for (j=0; j<cg_param[4]; j++){
			r1 = 2.0 * (double)j / (double)cg_param[4]; r2 = pow( -1.0 , (double)(j*2) ) * DEL;
			r3 = 2.0 * (double)(j+1) / (double)cg_param[4]; r4 = pow( -1.0 , (double)(j*2+1) ) * DEL;
			pos6[i][j][0][0] = q1 * cos(PI*(2*(double)(j)/cg_param[4]+pow(-1,(double)(j*2))*DEL)); pos6[i][j][0][1] = p1;
			pos6[i][j][0][2] = q1 * sin(PI*(2*(double)(j)/cg_param[4]+pow(-1,(double)(j*2))*DEL));
			pos6[i][j][1][0] = q1 * cos(PI*(2*(double)(j+1)/cg_param[4]+pow(-1,(double)(j*2+1))*DEL)); pos6[i][j][1][1] = p1;
			pos6[i][j][1][2] = q1 * sin(PI*(2*(double)(j+1)/cg_param[4]+pow(-1,(double)(j*2+1))*DEL));
			if(i==cg_param[3]-2){
				pos6[i][j][2][0] = q2 * cos( PI * ( r3 + r4 ) ); pos6[i][j][2][1] = p1; pos6[i][j][2][2] = q2 * sin( PI * ( r3 + r4 ) );
				pos6[i][j][3][0] = q2 * cos( PI * ( r1 + r2 ) ); pos6[i][j][3][1] = p1; pos6[i][j][3][2] = q2 * sin( PI * ( r1 + r2 ) );
			}
			else{
				pos6[i][j][2][0] = q3 * cos( PI * ( r3 + r4 ) ); pos6[i][j][2][1] = p2; pos6[i][j][2][2] = q3 * sin( PI * ( r3 + r4 ) );
				pos6[i][j][3][0] = q3 * cos( PI * ( r1 + r2 ) ); pos6[i][j][3][1] = p2; pos6[i][j][3][2] = q3 * sin( PI * ( r1 + r2 ) );
			}
			if(i==cg_param[3]-1){
				p3 = (float)(i-1) * xy - 500.0;
				q4 = ra * hole * bet * (float)((i-1)*2); q5 = ra * hole * bet * (float)((i-1)*2+1);
				pos6[i][j][0][0] = q4 * cos( PI * ( r1 + r2 ) ); pos6[i][j][0][1] = p3; pos6[i][j][0][2] = q4 * sin( PI * ( r1 + r2 ) );
				pos6[i][j][1][0] = q4 * cos( PI * ( r3 + r4 ) ); pos6[i][j][1][1] = p3; pos6[i][j][1][2] = q4 * sin( PI * ( r3 + r4 ) );
				pos6[i][j][2][0] = q5 * cos( PI * ( r3 + r4 ) ); pos6[i][j][2][1] = p3; pos6[i][j][2][2] = q5 * sin( PI * ( r3 + r4 ) );
				pos6[i][j][3][0] = q5 * cos( PI * ( r1 + r2 ) ); pos6[i][j][3][1] = p3; pos6[i][j][3][2] = q5 * sin( PI * ( r1 + r2 ) );
			}
			col6[i][j][0] = 150.0*cos((double)(count-i*200)/(100+j)); col6[i][j][1] = 100.0*cos((double)(count+i*150)/(150+j));
			col6[i][j][2] = 130.0*cos((double)(count-i*150)/(250+j));
		}
	}
}

cg7(float pos7[XYMAX][XZMAX][4][3], int col7[XYMAX][XZMAX][3]){
	int i,j;
	for(i=0; i<cg_param[3]; i++){
		for (j=0; j<cg_param[4]; j++){
			pos7[i][j][0][0] = (random()%1001)/1000.0 * 2000.0 - 1000.0; pos7[i][j][0][1] = (random()%1001)/1000.0 * 2000.0 - 1000.0;
			pos7[i][j][0][2] = (random()%1001)/1000.0 * 2000.0 - 1000.0;
			pos7[i][j][1][0] = pos7[i][j][0][0]+50.0; pos7[i][j][1][1] = pos7[i][j][0][1]; pos7[i][j][1][2] = pos7[i][j][0][2];
			pos7[i][j][2][0] = pos7[i][j][0][0]+50.0; pos7[i][j][2][1] = pos7[i][j][0][1]+50.0; pos7[i][j][2][2] = pos7[i][j][0][2];
			pos7[i][j][3][0] = pos7[i][j][0][0]; pos7[i][j][3][1] = pos7[i][j][0][1]+50.0; pos7[i][j][3][2] = pos7[i][j][0][2];
			col7[i][j][0] = 150.0*cos((double)(count-i*200)/(100+j)); col7[i][j][1] = 100.0*cos((double)(count+i*150)/(150+j));
			col7[i][j][2] = 130.0*cos((double)(count-i*150)/(250+j));
		}
	}
}

cg8(float pos8[XYMAX][XZMAX][4][3], int col8[XYMAX][XZMAX][3]){
	int i,j;
	float cn, buf;
	cn = 30.0 * sin((double)count/5000.0);
	for (i=0; i<cg_param[3]; i++){
		buf = 100.0 * (float)(i+1);
		for (j=0; j<cg_param[4]-1; j++){
			col8[i][j][0] = 250.0*cos((double)(i+1)*cn+(double)count/800.0); col8[i][j][1] = 200.0*sin((double)(i+1)*cn*(double)count/900.0);
			col8[i][j][2] = 200.0*cos((double)(i+1)*cn*(double)count/1000.0);
			pos8[i][j][0][0] = buf*cos((double)j*cn) * cos((double)i*cn/10.0); pos8[i][j][0][1] = buf*sin((double)j*cn);
			pos8[i][j][0][2] = buf*cos((double)j*cn) * sin((double)i*cn/10.0);
			pos8[i][j][1][0] = buf*cos((double)(j+1)*cn) * cos((double)i*cn/10.0); pos8[i][j][1][1] = buf*sin((double)(j+1)*cn);
			pos8[i][j][1][2] = buf*cos((double)(j+1)*cn) * sin((double)i*cn/10.0);
			pos8[i][j][2][0] = pos8[i][j][1][0]; pos8[i][j][2][1] = pos8[i][j][1][1]+6.0; pos8[i][j][2][2] = pos8[i][j][1][2];
			pos8[i][j][3][0] = pos8[i][j][0][0]; pos8[i][j][3][1] = pos8[i][j][0][1]+6.0; pos8[i][j][3][2] = pos8[i][j][0][2];
		}
		col8[i][cg_param[4]-1][0] = col8[i][0][0]; col8[i][cg_param[4]-1][1] = col8[i][0][1]; col8[i][cg_param[4]-1][2] = col8[i][0][2];
		pos8[i][cg_param[4]-1][0][0] = buf*cos((double)(cg_param[4]-1)*cn) * cos((double)i*cn/10.0);
		pos8[i][cg_param[4]-1][0][1] = buf*sin((double)(cg_param[4]-1)*cn);
		pos8[i][cg_param[4]-1][0][2] = buf*cos((double)(cg_param[4]-1)*cn) * sin((double)i*cn/10.0);
		pos8[i][cg_param[4]-1][1][0] = pos8[i][0][0][0]; pos8[i][cg_param[4]-1][1][1] = pos8[i][0][0][1]; pos8[i][cg_param[4]-1][1][2] = pos8[i][0][0][2];
		pos8[i][cg_param[4]-1][2][0] = pos8[i][cg_param[4]-1][1][0]; pos8[i][cg_param[4]-1][2][1] = pos8[i][cg_param[4]-1][1][1]+6.0;
		pos8[i][cg_param[4]-1][2][2] = pos8[i][cg_param[4]-1][1][2];
		pos8[i][cg_param[4]-1][3][0] = pos8[i][cg_param[4]-1][0][0]; pos8[i][cg_param[4]-1][3][1] = pos8[i][cg_param[4]-1][0][1]+6.0;
		pos8[i][cg_param[4]-1][3][2] = pos8[i][cg_param[4]-1][0][2];
	}
}

cg9(float pos9[XYMAX][XZMAX][4][3], int col9[XYMAX][XZMAX][3]){
	int i,j;
	float cn, buf;
	cn = 30.0 * sin((double)count/5000.0);
	for (i=0; i<cg_param[3]; i++){
		buf = 100.0 * (float)(i+1);
		for (j=0; j<cg_param[4]-1; j++){
			col9[i][j][0] =250.0*cos((double)(i+1)+cn*(double)count/1800.0); col9[i][j][1] =200.0*sin((double)(i+1)+cn*(double)count/1900.0);
			col9[i][j][2] =200.0*cos((double)(i+1)+cn*(double)count/1500.0);
			pos9[i][j][0][0] = (buf*cos(PI/10.0*(double)j) + 100.0*cos((double)j*cn)) * cos((double)i*cn/10.0);
			pos9[i][j][0][1] = (buf*sin(PI/10.0*(double)j) + 100.0*sin((double)j*cn)) * sin((double)i*cn/10.0); pos9[i][j][0][2] = 200.0*sin((double)j*cn);
			pos9[i][j][1][0] = (buf*cos(PI/10.0*(double)(j+1)) + 100.0*cos((double)(j+1)*cn)) * cos((double)i*cn/10.0);
			pos9[i][j][1][1] = (buf*sin(PI/10.0*(double)(j+1)) + 100.0*sin((double)(j+1)*cn)) * sin((double)i*cn/10.0); pos9[i][j][1][2] = 200.0*sin((double)(j+1)*cn);
			pos9[i][j][2][0] = pos9[i][j][1][0]; pos9[i][j][2][1] = pos9[i][j][1][1]+16.0; pos9[i][j][2][2] = pos9[i][j][1][2];
			pos9[i][j][3][0] = pos9[i][j][0][0]; pos9[i][j][3][1] = pos9[i][j][0][1]+16.0; pos9[i][j][3][2] = pos9[i][j][0][2]; 
		}
		col9[i][cg_param[4]-1][0] = col9[i][0][0]; col9[i][cg_param[4]-1][1] = col9[i][0][1]; col9[i][cg_param[4]-1][2] = col9[i][0][2];
		pos9[i][cg_param[4]-1][0][0] = (buf*cos(PI/10.0*(double)(cg_param[4]-1)) + 100.0*cos((double)(cg_param[4]-1)*cn)) * cos((double)i*cn/10.0);
		pos9[i][cg_param[4]-1][0][1] = (buf*sin(PI/10.0*(double)(cg_param[4]-1)) + 100.0*sin((double)(cg_param[4]-1)*cn)) * sin((double)i*cn/10.0);
		pos9[i][cg_param[4]-1][0][2] = 200.0*sin((double)(cg_param[4]-1)*cn);
		pos9[i][cg_param[4]-1][1][0] = pos9[i][0][0][0]; pos9[i][cg_param[4]-1][1][1] = pos9[i][0][0][1]; pos9[i][cg_param[4]-1][1][2] = pos9[i][0][0][2];
		pos9[i][cg_param[4]-1][2][0] = pos9[i][cg_param[4]-1][1][0]; pos9[i][cg_param[4]-1][2][1] = pos9[i][cg_param[4]-1][1][1]+16.0;
		pos9[i][cg_param[4]-1][2][2] = pos9[i][cg_param[4]-1][1][2];
		pos9[i][cg_param[4]-1][3][0] = pos9[i][cg_param[4]-1][0][0]; pos9[i][cg_param[4]-1][3][1] = pos9[i][cg_param[4]-1][0][1]+16.0;
		pos9[i][cg_param[4]-1][3][2] = pos9[i][cg_param[4]-1][0][2];
	}
}

cg10(float pos10[XYMAX][XZMAX][4][3], int col10[XYMAX][XZMAX][3]){
	int i,j;
	float cn, buf;	
	cn = 30.0 * sin((double)count/5000.0);
	for (i=1; i<cg_param[3]+1; i++){
		buf = 100.0 * (float)i;
		for (j=1; j<cg_param[4]; j++){
			col10[i-1][j-1][0] =250.0*cos((double)i+cn*count/2000.0); col10[i-1][j-1][1] =200.0*sin((double)i+cn*count/1500.0);
			col10[i-1][j-1][2] =200.0*cos((double)i+cn*count/2300.0);
			pos10[i-1][j-1][0][0] = (buf*cos(PI/10.0*(double)j) - 50.0*sin((double)j*cn/(double)i*3.0)) * cos((double)i*cn/50.0);
			pos10[i-1][j-1][0][1] = (buf*sin(PI/10.0*(double)j) - 60.0*cos((double)j*cn/(double)i*3.0)) * sin((double)i*cn/50.0);
			pos10[i-1][j-1][0][2] = 200.0*sin((double)(j*i)*cn/5.0);
			pos10[i-1][j-1][1][0] = (buf*cos(PI/10.0*(double)(j+1)) - 50.0*sin((double)(j+1)*cn/(double)i*3.0)) * cos((double)i*cn/50.0);
			pos10[i-1][j-1][1][1] = (buf*sin(PI/10.0*(double)(j+1)) - 60.0*cos((double)(j+1)*cn/(double)i*3.0)) * sin((double)i*cn/50.0);
			pos10[i-1][j-1][1][2] = 200.0*sin((double)((j+1)*i)*cn/10.0);
			pos10[i-1][j-1][2][0] = pos10[i-1][j-1][1][0]; pos10[i-1][j-1][2][1] = pos10[i-1][j-1][1][1]+8.0;
			pos10[i-1][j-1][2][2] = pos10[i-1][j-1][1][2];
			pos10[i-1][j-1][3][0] = pos10[i-1][j-1][0][0]; pos10[i-1][j-1][3][1] = pos10[i-1][j-1][0][1]+8.0;
			pos10[i-1][j-1][3][2] = pos10[i-1][j-1][0][2];
		}
		col10[i-1][cg_param[4]-1][0] = 250.0*cos((double)(i+1)+cn*count/2000.0); col10[i-1][cg_param[4]-1][1] = 200.0*sin((double)(i+1)+cn*count/1500.0);
		col10[i-1][cg_param[4]-1][2] = 200.0*cos((double)(i+1)+cn*count/2300.0);
		pos10[i-1][cg_param[4]-1][0][0] = (buf*cos(PI/10.0*(double)cg_param[4]) - 50.0*sin((double)cg_param[4]*cn/(double)i*3.0)) * cos((double)i*cn/50.0);
		pos10[i-1][cg_param[4]-1][0][1] = (buf*sin(PI/10.0*(double)cg_param[4]) - 60.0*cos((double)cg_param[4]*cn/(double)i*3.0)) * sin((double)i*cn/50.0);
		pos10[i-1][cg_param[4]-1][0][2] = 200.0*sin((double)cg_param[4]*cn*(double)i/5.0); pos10[i-1][cg_param[4]-1][1][0] = pos10[i-1][0][0][0];
		pos10[i-1][cg_param[4]-1][1][1] = pos10[i-1][0][0][1]; pos10[i-1][cg_param[4]-1][1][2] = 200.0*sin((double)cg_param[4]*cn*(double)i/10.0);
		pos10[i-1][cg_param[4]-1][2][0] = pos10[i-1][cg_param[4]-1][1][0]; pos10[i-1][cg_param[4]-1][2][1] = pos10[i-1][cg_param[4]-1][1][1]+8.0;
		pos10[i-1][cg_param[4]-1][2][2] = pos10[i-1][cg_param[4]-1][1][2]; pos10[i-1][cg_param[4]-1][3][0] = pos10[i-1][cg_param[4]-1][0][0];
		pos10[i-1][cg_param[4]-1][3][1] = pos10[i-1][cg_param[4]-1][0][1]+8.0; pos10[i-1][cg_param[4]-1][3][2] = pos10[i-1][cg_param[4]-1][0][2];
	}
}

cgend(float pose[XYMAX][XZMAX][4][3], int cole[XYMAX][XZMAX][3]){
	int i,j,k;
	for (i=0; i<cg_param[3]; i++){
		for (j=0; j<cg_param[4]; j++){
			for (k=0; k<4; k++){ pose[i][j][k][0] = 0.0; pose[i][j][k][1] = 0.0; pose[i][j][k][2] = 0.0; }
			cole[i][j][0] = 0; cole[i][j][1] = 0; cole[i][j][2] = 0;
		}
	}
}

cg_change_now(){
	switch(cg_cnt){
		case 1:	cg1(now, nowc);	break;
		case 2:	cg2(now, nowc); break;
		case 3:	cg3(now, nowc); break;
		case 4:	cg4(now, nowc); break;
		case 5:	cg5(now, nowc); break;
		case 6:	cg6(now, nowc); break;
		case 7:	cg7(now, nowc); break;
		case 8:	cg8(now, nowc); break;
		case 9:	cg9(now, nowc); break;
		case 10:	cg10(now, nowc); break;
	}
}

cg_change_next(){
	switch(cg_cnt){
		case 0:	cgend(next,nextc); break;
		case 1:	cg1(next, nextc); break;
		case 2:	cg2(next, nextc); break;
		case 3:	cg3(next, nextc); break;
		case 4:	cg4(next, nextc); break;
		case 5:	cg5(next, nextc); break;
		case 6:	cg6(next, nextc); break;
		case 7:	cg7(next, nextc); break;
		case 8:	cg8(next, nextc); break;
		case 9:	cg9(next, nextc); break;
		case 10:	cg10(next, nextc); break;
	}
}

sound_play(int no){
	switch(no){
		case 0:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s00.aiff", NULL ); exit(7); }	break;
		case 1:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s01.aiff", NULL ); exit(7); }	break;
		case 2:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s02.aiff", NULL ); exit(7); }	break;
		case 3:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s03.aiff", NULL ); exit(7); }	break;
		case 4:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s04.aiff", NULL ); exit(7); }	break;
		case 5:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s05.aiff", NULL ); exit(7); }	break;
		case 6:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s06.aiff", NULL ); exit(7); }	break;
		case 7:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s07.aiff", NULL ); exit(7); }	break;
		case 8:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s08.aiff", NULL ); exit(7); }	break;
		case 9:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s09.aiff", NULL ); exit(7); }	break;
		case 10:	if(fork()==0){ execl( "sounds/mini_play", "mini_play", "sounds/s10.aiff", NULL ); exit(7); }	break;
	}
}

morph(int cg_cnt){
	int i, j, k, l, cnt=0, changec[XYMAX][XZMAX][3], vc[XYMAX][XZMAX][4];
	float change[XYMAX][XZMAX][4][3], vec[XYMAX][XZMAX][4][3], VELOCITY;
	VELOCITY = (float)cg_param[2];
	for(i=0; i<cg_param[3]; i++){
		for(j=0; j<cg_param[4]; j++){
			for(k=0; k<4; k++){ for(l=0; l<3; l++) vec[i][j][k][l] = next[i][j][k][l] - now[i][j][k][l]; }	
			for(l=0; l<3; l++) vc[i][j][l]= nextc[i][j][l] - nowc[i][j][l];
		}
	}
	while((cnt++) < (int)VELOCITY){
		RGBcolor(0,0,0);  clear();	
		for(i=0; i<cg_param[3]; i++){
			for(j=0; j<cg_param[4]; j++){
				for(k=0; k<4; k++){ for(l=0; l<3; l++) change[i][j][k][l] = vec[i][j][k][l]/2.0 * (cos(PI-(double)cnt*PI/VELOCITY) + 1.0) + now[i][j][k][l]; }
				for(l=0; l<3; l++) changec[i][j][l] = vc[i][j][l]/2.0 * (cos(PI-(double)cnt*PI/VELOCITY) + 1.0) + nowc[i][j][l];	
			}
		}
		pushmatrix(); yota(); draw_plane(change, changec); popmatrix();
		if( cg_cnt < 8 ) count++;
		swapbuffers();
	}
	for(i=0; i<cg_param[3]; i++){
		for(j=0; j<cg_param[4]; j++){
			for(k=0; k<4; k++){ for(l=0; l<3; l++) now[i][j][k][l] = next[i][j][k][l]; }
			for(l=0; l<3; l++) nowc[i][j][l] = nextc[i][j][l];	
		}
	}
}

main()
{
	unsigned char status, channel, keyno;
	short val;
	MIport *midi_port;
	MIevent midi[10000];
	MIconfig *c;
	c = MInewconfig(); midi_port = MInewport();
	if( MIopen( midi_port, "r", &c ) < 0 ) exit(-1);
	cg_param[0] = 0; cg_param[1] = 0; cg_param[2] = 15; cg_param[3] = 2; cg_param[4] = 2;
	cg_param[5] = 100; cg_param[6] = 100;
	cgend(now,nowc); open_window();
	while (!(qtest()&&qread(&val)==ESCKEY&&val==0)) {
		if( MIreceive( midi_port, midi, 1 ) > 0 ){
			status = midi_get_status( midi ); channel = midi_get_channel( midi ); keyno = midi_get_keyno( midi );
			if( status == 0xd0 ){
				switch(channel) {
					case 0x00:	cg_param[0] = keyno + 5;		break;
					case 0x01:	cg_param[1] = (128-keyno) * 9;		break;
					case 0x02:	cg_param[2] = keyno + 1;		break;
					case 0x03:	cg_param[3] = keyno/6 + 2; cg_change_now();	break;
					case 0x04:	cg_param[4] = keyno/6 + 2; cg_change_now();	break;
					case 0x05:	cg_param[5] = keyno * 2;		break;
					case 0x06:	cg_param[6] = keyno * 2;		break;
					case 0x0e:	sound_play(keyno);		break;
					case 0x0f:	if( keyno != cg_cnt ){
							cg_cnt = keyno;
							cg_change_next();
							morph( cg_cnt );
						}
						break;
				}
			}
		}
		else if( cg_cnt > 7 ){
			switch(cg_cnt){
				case 8:	cg8(now, nowc);		break;
				case 9:	cg9(now, nowc);		break;
				case 10:	cg10(now, nowc);	break;
			}
		}
		RGBcolor(0,0,0); clear(); pushmatrix(); yota(); draw_plane(now, nowc); popmatrix(); count ++; swapbuffers();
	}
	MIclose( midi_port ); MIfreeport( midi_port ); gexit(); return 0;
}