【OpenGL】细分三角形

#include <GL/glut.h>
typedef GLfloat point2[2];
point2 v[]={{-1.0,-0.58},{1.0,-0.58},{0.0,1.15}};
int n;

void mydisplay()
{
	glClear(GL_COLOR_BUFFER_BIT);
	glColor3f(0.5,0.5,0.4);
	glBegin(GL_LINE);
		glVertex2d(-0.5,-0.5);
		glVertex2d(-0.5,0.5);
		glVertex2d(0.5,0.5);
		glVertex2d(0.5,-0.5);
	glEnd();
	glFlush();
}

void triangle(point2 a,point2 b,point2 c)
{
	glVertex2fv(a);
	glVertex2fv(b);
	glVertex2fv(c);
}

void divide_triangle(point2 a,point2 b,point2 c,int m)
{
	point2 v0,v1,v2;
	int j;
	if(m>0)
	{
		for (j=0;j<2;j++)
		{
			v0[j]=(a[j]+b[j])/2;
		}
		for (j=0;j<2;j++)
		{
			v1[j]=(a[j]+c[j])/2;
		}
		for (j=0;j<2;j++)
		{
			v2[j]=(b[j]+c[j])/2;
		}
		divide_triangle(a,v0,v1,m-1);
		divide_triangle(c,v1,v2,m-1);
		divide_triangle(b,v2,v0,m-1);
		//divide_triangle(v0,v1,v2,m-1);
	}
	else
	{
		triangle(a,b,c);
	}
}

void display(void)
{
	glClear(GL_COLOR_BUFFER_BIT);
	glBegin(GL_TRIANGLES);
		divide_triangle(v[0],v[1],v[2],n);
	glEnd();
	glFlush();
}

void myinit()
{
	glMatrixMode(GL_PROJECTION);
	glLoadIdentity();
	gluOrtho2D(-2.0,2.0,-2.0,2.0);
	glMatrixMode(GL_MODELVIEW);
	glClearColor(1.0,1.0,1.0,1.0);
	glColor3f(0.0,0.0,0.0);
}

int main(int argc,char* argv[])
{
	n=4;
	glutInit(&argc,argv);
	glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB);
	glutInitWindowSize(500,500);
	glutCreateWindow("Arthur");
	glutDisplayFunc(display);
	myinit();
	glutMainLoop();
	return 1;
}