3d моделирование

#include <iostream>
#include <stdlib.h>
#include <allegro.h>
#include <math.h>
#define w 1024
#define h 768
#define kybe 1
#define treygl 2
#define chetireygl 3 
 
void draw();
int xx(int);
intyy(int);
intsdvig(int,int,int);
voidmashtab(double,double,double);
voidymnozhenie(double b[4][4]);
voidpovorot_x(double );
voidpovorot_y(double );
voidpovorot_z(double );
voidvrashenie_cikl();
voidknopki(int,char*);
void knopki0();
intproverka_koordinat(int,char*);
void proverka_koordinat0();
void focus(int,char*);
voidproverka_vrashenie();
void line1(int,int,int,int);
intkoeficenti(int);
intzakraska(int);
 
 
 
using namespace std;
intdx,dy,dz;
inttipe;
intflag_vrash=0;
int *xl,*yl,*zl,*gl;
int c;
int I1,I2,I3,I4,Iu,Iw,Iv;
int L1,L2,L3,L4,L0;
 
 
double n=1;
int o;
intaa=0,ab=0,ac=0;
BITMAP *bmp;
double a[4][4]={1,0,0,0,
                0,1,0,0,
                0,0,1,0,
0,0,0,1 };
 
intxk[]={0,0,100,100,0,0,100,100};
intyk[]={0,100,100,0,0,100,100,0};
intzk[]={100,100,100,100,0,0,0,0};
intgk[]={0,1,2,3,-1,
           4,7,6,5,-1,
           1,5,6,2,-1,
           0,3,7,4,-1,
           3,2,6,7,-1,
           0,4,5,1,-1,100};
 
intxt[]={0,0,100,0};
intyt[]={0,100,0,0};
intzt[]={0,0,0,100};
intgt[]={0,1,2,-1,
          0,2,3,-1,
          1,3,2,-1,
          0,3,1,-1,100};
 
int xc[]={-50,100,100,100,100,250};
intyc[]={50,0,100,100,0,50};
intzc[]={50,100,100,0,0,50};
intgc[]={1,0,2,-1,
          0,3,2,-1,
          3,0,4,-1,
          1,4,0,-1,
          1,2,5,-1,
          5,2,3,-1,
          4,5,3,-1,
          4,1,5,-1,100};         
 
 
intvektor_s[]={0,0,-1,0};
intvektor_l[]={0,0,-1,0};
intvektor_r[]={0,0,-1,0};
intnorm_vektor[]={1,1,1,1};
intsx=350;
intsy=300;
 
 
int main(intargc, char *argv[])
{
tipe=-1;
allegro_init();
install_mouse();
install_keyboard();
set_color_depth(16); 
bmp = create_bitmap(w,h);
set_gfx_mode(GFX_AUTODETECT_FULLSCREEN , w,  h, 0, 0);
show_mouse(bmp);
c=makecol(225,225,0);
L0=makecol(225,225,0);
L4=makecol(190,250,210);
knopki0();
 
do{
proverka_vrashenie();
    proverka_koordinat0();
draw();
if(n/2==int(n/2))vrashenie_cikl();
if (key[KEY_RIGHT]){sdvig(1,0,0);}
if (key[KEY_LEFT]){sdvig(-1,0,0);}
if (key[KEY_UP]){sdvig(0,-1,0);}
if (key[KEY_DOWN]){sdvig(0,1,0);}
if (key[KEY_1]){mashtab(1.01,1.01,1.01);}
if (key[KEY_2]){mashtab(0.99,0.99,0.99);}
 
 
if (key[KEY_X]){povorot_x(0.002);}  
if (key[KEY_Y]){povorot_y(0.002);} 
if (key[KEY_Z]){povorot_z(0.002);} 
 
blit(bmp,screen,0,0,0,0,w,h);
 
}
 
while(!key[KEY_ENTER]);
 
allegro_exit();
return 0;
}
END_OF_MAIN();
 
void draw()
{ 
 
rectfill(bmp,25,50,750,h-100,makecol(0,0,0));
rect(bmp,25,50,750,h-100,makecol(255,255,255)); 
 
 
if (tipe==-1)return;
inti=0;
int h0=i;
int v=koeficenti(h0);
do
 {
 
if(gl[i+1]!=-1){
if (v>0) line1(xx(gl[i]),yy(gl[i]),xx(gl[i+1]),yy(gl[i+1]));
 
}
 
else{
 
if (v>0){line1(xx(gl[i]),yy(gl[i]),xx(gl[h0]),yy(gl[h0]));
zakraska(h0);}
h0=i+2;
v=koeficenti(h0);
i++;
if (gl[i]==100) break;
  }
i++;
   }
while (gl[i]!=100);
}
int xx(int j){return int(xl[j]*a[0][0]+yl[j]*a[1][0]+zl[j]*a[2][0]+a[3][0])+sx;}
intyy(int j){return int(xl[j]*a[0][1]+yl[j]*a[1][1]+zl[j]*a[2][1]+a[3][1])+sy;}
 
intsdvig(intdx,intdy,intdz)
{
double b[4][4]={1,0,0,0,
                    0,1,0,0,
                    0,0,1,0,     
dx,dy,dz,1 };
 
ymnozhenie(b);
}
voidymnozhenie(double b[4][4])
{
double  f[4][4]={0,0,0,0,
                  0,0,0,0,
                  0,0,0,0,
0,0,0,0 };
 
double d=0;
for(int k=0;k<4;k++)
{
for(int l=0;l<4;l++)
{d=0;
for(inti=0;i<4;i++)
  {
      d=d+a[k][i]*b[i][l];
f[k][l]=d;    
  }
 
}    
}          
for(inti=0;i<4;i++)
 {
for(int j=0;j<4;j++)
 {
a[i][j]=f[i][j];
 }
 }
}
voidmashtab(double dx,doubledy,doubledz)
{   aa=a[3][0];
ab=a[3][1];
ac=a[3][2];
sdvig(-aa,-ab,-ac);
double b[4][4]={ dx,0,0,0,
                     0,dy,0,0,
                     0,0,dz,0,
0,0,0,1 };
ymnozhenie(b);
sdvig(aa,ab,ac);
}  
voidpovorot_x(double alfa)
 
{   aa=a[3][0];
ab=a[3][1];
ac=a[3][2];
sdvig(-aa,-ab,-ac);
double b[4][4]={  1,0,0,0,
                     0,cos(alfa),sin(alfa),0,
                     0,-sin(alfa),cos(alfa),0,
0,0,0,1 };
ymnozhenie(b);
sdvig(aa,ab,ac);
}
 
voidpovorot_y(double alfa)
{ aa=a[3][0];
ab=a[3][1];
ac=a[3][2];
sdvig(-aa,-ab,-ac);
     double b[4][4]={cos(alfa),0,sin(alfa),0,
                     0,1,0,0,
                     -sin(alfa),0,cos(alfa),0,
0,0,0,1 };
ymnozhenie(b);
sdvig(aa,ab,ac);
                 }    
 
voidpovorot_z(double alfa)
{ aa=a[3][0];
ab=a[3][1];
ac=a[3][2];
sdvig(-aa,-ab,-ac);
     double b[4][4]={  cos(alfa),sin(alfa),0,0,
                     -sin(alfa),cos(alfa),0,0,
                     0,0,1,0,
0,0,0,1 };
 
ymnozhenie(b);
sdvig(aa,ab,ac);
                 }    
 
voidvrashenie_cikl()
{
for(inti=0;i<1;i++)
       {
povorot_x(0.001);
povorot_y(0.001);             
povorot_z(0.001);
 
        }
}
voidknopki(int y ,char *str)
{
rectfill(bmp,800,y,915,y-50,0);
rect(bmp,800,y,915,y-50,makecol(255,255,255));
 
textout(bmp,font,str,830,y-25,255);
textout(screen, font, "VRASHENIE-X,Y,Z", 160, 700, 255); 
textout(screen, font, "MASHTABIROVANIE-1,2", 160, 730, 255);
textout(screen, font, "PEREMESHENIE-RIGHT,LEFT,UP,DOWN", 430, 700, 255); 
}
void knopki0()
{
knopki(100,"kyb");
knopki(200,"tetraedr");
knopki(300,"oktaedr");
knopki(400,"vrashenie");
knopki(500,"exit");
}
void focus(int y ,char *str)
{
rectfill(bmp,800,y,915,y-50,makecol(0,0,255));
textout(bmp,font,str,1050,y-25,255);
}
intproverka_koordinat(inty,char *str)
{
if((mouse_x>800)&& (mouse_x<900) && (mouse_y>y-50) && (mouse_y<y) && (mouse_b& 1)) 
 {   focus(y,str); 
return 1;
  }
else
  {
knopki(y,str);
return 0;    
 
 
}
 
void proverka_koordinat0()
{ 
 
if (proverka_koordinat(100,"kyb"))       {xl=xk;yl=yk;zl=zk;gl=gk;tipe=kybe;draw();}
if (proverka_koordinat(200,"tetraedr"))  {xl=xt;yl=yt;zl=zt;gl=gt;tipe=treygl;draw();}
if (proverka_koordinat(300,"oktaedr"))   {xl=xc;yl=yc;zl=zc;gl=gc;tipe=chetireygl;draw();}
if (proverka_koordinat(500,"exit")) exit(0);
 
}
voidproverka_vrashenie()
{ if(mouse_x>800&&mouse_x<900 &&mouse_y>350 &&mouse_y<400)
 
{ if(mouse_b& 1)flag_vrash=1;
else if(!(mouse_b& 1)&& (flag_vrash))
        {
        n=n+1;
flag_vrash=0;   
      }    
        }    
}
void line1(int x1,int y1,int x2,int y2)
{
 
if(x1>750)x1=750;
if(x2>750)x2=750; 
if(x1<50)x1=50;
if(x2<50)x2=50;
if(y1<50)y1=50; 
if(y2<50)y2=50; 
if(y1>h-100)y1=h-100; 
if(y2>h-100)y2=h-100;  
 line(bmp,x1,y1,x2,y2,makecol(0,0,255));
}    
intkoeficenti(int t)
{
int a_k=(xx(gl[t+1])-xx(gl[t]))*(yy(gl[t+2])-yy(gl[t]))-(xx(gl[t+2])-xx(gl[t]))*(yy(gl[t+1])-yy(gl[t]));
if(a_k<0)return 1;
if(a_k>=0)return 0;
}
intzakraska(int t)
{ 
double beta=M_PI/1.5;
double R=255*L1+0.3*L0*cos(beta);
double G=155*L2+0.3*L0*cos(beta);
double B=115*L4+0.3*L0*cos(beta);
 
doublex_x=(xx(gl[t])+xx(gl[t+1])+xx(gl[t+2]))/3;
doubley_y=(yy(gl[t])+yy(gl[t+1])+yy(gl[t+2]))/3;
 
doublex_g=xx(gl[t]);double x_g1=xx(gl[t+2]);
doubley_g=yy(gl[t]);double y_g1=yy(gl[t+2]);
 
  //  line(bmp,x_g,y_g,x_g1,y_g1,c);
   // line(bmp,xx(gl[t+1]),yy(gl[t+1]),xx(gl[t+3]),yy(gl[t+3]),c);
 
 
floodfill(bmp,x_x,y_y,Iw);
 
 
    I1=0.7*L1+0.3*L0*cos(beta);
    I2=0.3*L2+0.3*L0*cos(beta);
    I3=0.7*L3+0.3*L0*cos(beta);
    I4=0.8*L4+0.3*L0*cos(beta);
 
    Iv=0.00001*L4+0.00001*L1+0.3*L0*cos(beta);
Iu=0.00001*L3+0.00001*L2+0.3*L0*cos(beta);
 
Iw=(((1-x_x)*Iu+x_x*Iv)+((1-y_y)*Iu+y_y*Iv))/32;
 
}