[OpenGL]课后案例02：生成Sierpinski镂垫的递归程序

//A.2  生成Sierpinski镂垫的递归程序
/* Recursive subdivision of triangle to form Sierpinski gasket */
/* number of recursive steps given on command line */
#include <GL/glut.h>
#include<stdlib.h>
/* initial triangle */
GLfloat v[3][2]={{-1.0, -0.58}, {1.0, -0.58}, {0.0, 1.15}};
int n;
void triangle( GLfloat *a, GLfloat *b, GLfloat *c)
    /* display one triangle  */
{
    glVertex2fv(a);
    glVertex2fv(b);
    glVertex2fv(c);
}
void divide_triangle(GLfloat  *a, GLfloat *b, GLfloat *c, int m)
{
    /* triangle subdivision using vertex numbers */
    GLfloat v0[2], v1[2], v2[2];
    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);
    }
    else triangle(a,b,c); /* draw triangle at end of recursion */
}
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);
}
void  main(int argc, char **argv)
{
    n=10;
    //n=atoi(argv[1]); /* or set number of subdivision steps here */
    //注意，上句若是在VS或VC++中直接使用会报错，因为需要输入一个参数。
    //可以考虑使用命令提示行，在windows下打开项目的debug文件夹，按住shift键右击
    //选择：在当前位置打开命令提示行。
    //输入：Main.exe 10
    //其中，Main.exe根据自己情况确定，后面的10是递归次数，可以自定义
    //实际检测，当数字大于15会出现卡死状态，数字大于100直接死机
    glutInit(&argc, argv);
    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB);
    glutInitWindowSize(500, 500);
    glutCreateWindow("Sierpinski Gasket");
    glutDisplayFunc(display);
    myinit();
    glutMainLoop();

}