home *** CD-ROM | disk | FTP | other *** search
- /* TC-13.C */
- /* This program demonstrates a mirror reflection of
- a computer generated 3d solid cube */
-
- /* ----------------------------------------------------------------------- */
- /* INCLUDE FILES */
- #include <process.h>
- #include <bios.h>
- #include <stdio.h>
- #include <graphics.h>
- #include <math.h>
-
- /* ----------------------------------------------------------------------- */
- /* DECLARATIONS */
-
- float x=0.0,y=0.0,z=0.0; /* world coordinates */
- float x1=0.0,x2=0.0,x3=0.0; /* polygon vertices */
- float y01=0.0,y2=0.0,y3=0.0; /* polygon vertices */
- float z1=0.0,z2=0.0,z3=0.0; /* polygon vertices */
- float sx1=0.0,sx2=0.0,sx3=0.0,sx4=0.0,sx5=0.0; /* display coordinaes */
- float sy1=0.0,sy2=0.0,sy3=0.0,sy4=0.0,sy5=0.0; /* display coordinates */
- float sx=0.0,sy=0.0; /* output of 3d perspective formulas */
- float xa=0.0,ya=0.0,za=0.0; /* temporary values in 3d formulas */
- float d=1200.0; /* angular perspective value */
- double r1=5.68319; /* yaw angle in radians */
- double r2=6.28319; /* roll angle in radians */
- double r3=5.79778; /* pitch angle in radians */
- double sr1=0.0,sr2=0.0,sr3=0.0; /* sine rotation factors */
- double cr1=0.0,cr2=0.0,cr3=0.0; /* cosine rotation factors */
- float mx=0.0,my=0.0,mz=-450.0; /* viewpoint position */
- int maxx=639,minx=0,maxy=199,miny=0; /* scaling viewport */
- float screen_x=639,screen_y=199; /* dimensions of screen mode */
- int C0=0,C1=1,C2=2,C3=3,C4=4,C5=5,C6=6,C7=7,C8=8, /* color varariables */
- C9=9,C10=10,C11=11,C12=12,C13=13,C14=14,C15=15,
- mode_flag=0;
- float rx=0.0,ry=0.0; /* scaling values used in mapping routine */
- int t=0,t1=0,t2=0; /* loop counters */
- int p1=0; /* array indexer */
- int key_matte_clr=6; /* exclusive key matte color */
- int edge_clr=7; /* used to draw edges on models */
- int solid_clr=0; /* used to fill surfaces of solid models */
- float x_res,y_res; /* used for 2d mappinf from 640*480 template */
- float sp=0.0; /* visibility factor in hidden surface routine */
- float sp1=0.0,sp2=0.0,sp3=0.0; /* temporary values of sp */
-
- int NORM_HUE=0; /* halftoning color */
- int BASE_CLR=0; /* prep cooor for dithering */
- int CC4=0; /* underlay color for dithering */
- int CC5=0; /* overlay color for dithering */
- unsigned short CC6=0xffff; /* pattern for dithering */
- float xx1=.5,yy1=.70711,zz1=.5; /* location of light source - elevation 45 degrees, compas 135 degrees {{REDIFINED !!}}*/
- float v4=0.0; /* illumination factor 0 to 1 range */
- int v5=0; /* illumination factor base 1 range */
- int v6=11; /* VGA and EGA illumination range */
- int v7=0; /* temporary varaible used for illumination range */
- float xu=0.0,yu=0.0,zu=0.0; /* vector from vertex one to vertex two */
- float xv=0.0,yv=0.0,zv=0.0; /* vector from vertex 1 to vertex 3 */
- float xn=0.0,y0n=0.0,zn=0.0; /* surface perpendicular factor */
- float v1=0.0,v2=0.0; /* length of surface perpendicular vector */
- float v3=0.0; /* ratio of length to unit vector magnitude */
- float xw=0.0,yw=0.0,zw=0.0; /* surface perpendicual unit vector */
- int polary[4][2]; /* used in my 3d solid polygon draw routine */
-
- /* define halftone codes */
- char fill_0[]={0,0,0,0,0,0,0,0}; /* 0% fill */
- char fill_3[]={0,32,0,0,0,2,0,0}; /* 3% fill */
- char fill_6[]={32,0,2,0,128,0,8,0}; /* 6% fill */
- char fill_12[]={32,2,128,8,32,2,128,8}; /* 12% fill */
- char fill_25[]={68,17,68,17,68,17,68,17}; /* 25% fill */
- char fill_37[]={170,68,170,17,170,68,170,17}; /* 37% fill */
- char fill_50[]={85,170,85,170,85,170,85,170}; /* 50% fill */
- char fill_62[]={85,187,85,238,85,187,85,238}; /* 62% fill */
- char fill_75[]={187,238,187,238,187,238,187,238}; /* 75% fill */
- char fill_87[]={223,253,127,247,223,253,127,247}; /* 87% fill */
- char fill_93[]={255,223,255,223,255,223,255,223}; /* 93% fill */
- char fill_100[]={255,255,255,255,255,255,255,255}; /* 100% fill */
-
- /* database of vertices for 3d model */
- int array1[][3]={
- 30,-30,30, 30,30,30, -30,30,30, -30,-30,30, 30,30,-30,
- -30,30,-30, -30,-30,-30, 30,-30,-30};
-
- /* database of fill origins for 3d model */
- int array2[][3]={
- 0,-30,0, 0,0,-30, -30,0,0, 0,0,30, 30,0,0, 0,30,0};
-
- /* database of mirror image coordinates */
- int array3[][3]={
- -30,-90,30, 30,-90,30, 30,-90,-30};
-
- /* database of mirror fill coordinates */
- int array4[][3]={
- 0,-60,30, 30,-60,0};
-
- float B1[8][3]; /* array of 8 sets of xyz coordinates */
- float B2[8][2]; /* array of 8 sets of sx,sy diaplay coordinates */
- float B3[6][2]; /* array of 6 sets of sx,sy fill coordinates */
- float B5[3][2]; /* array of 3 sets of sx,sy mirror image coords */
- float B6[2][2]; /* array of 2 sets of sx,sy mirrir fill coordinates */
-
- void keyboard(void);void quit_pgm(void);void calc_3d(void);
- void rotation(void);void window(void);void graphics_setup(void);
- void coords(void);void draw_poly(void);void notice(int x,int y);
- void visibility_test(void);void store_coords(void);
- void illumination(void);void dither(void);void shade(void);
- void values(void);
-
- /* ----------------------------------------------------------------------- */
- /* MAIN ROUTINES */
-
- main(){
- graphics_setup();
- setviewport(0,0,maxx,maxy,1);
- key_matte_clr=C11;
- edge_clr=C7;solid_clr=C0;
- NORM_HUE=C11;
- store_coords();
-
- setcolor(C7);setfillpattern(fill_50,C7); /* draw the mirror */
- x=-80;y=-30;z=80;calc_3d();window();moveto(sx,sy);putpixel(sx,sy,C7);
- x=80;y=-30;z=80;calc_3d();window();lineto(sx,sy);
- x=80;y=-30;z=-80;calc_3d();window();lineto(sx,sy);
- x=-80;y=-30;z=-80;calc_3d();window();lineto(sx,sy);
- x=-80;y=-30;z=80;calc_3d();window();lineto(sx,sy);
- x=0;y=-30;z=0;calc_3d();window();floodfill(sx,sy,C7);
-
- surface3:
- x1=B1[0][0];y01=B1[0][1];z1=B1[0][2];x2=B1[1][0];y2=B1[1][1];
- z2=B1[1][2];x3=B1[2][0];y3=B1[2][1];z3=B1[2][2];visibility_test();
- if (sp>0) goto surface4;
- sx1=B2[0][0];sy1=B2[0][1];sx2=B2[1][0];sy2=B2[1][1];sx3=B2[2][0];
- sy3=B2[2][1];sx4=B2[3][0];sy4=B2[3][1];sx5=B3[3][0];sy5=B3[3][1];
- draw_poly();
- illumination();v5=4;
- shade();dither();
- mirror_surface3:
- sx1=B2[0][0];sy1=B2[0][1];sx2=B2[3][0];sy2=B2[3][1];sx3=B5[0][0];
- sy3=B5[0][1];sx4=B5[1][0];sy4=B5[1][1];sx5=B6[0][0];sy5=B6[0][1];
- draw_poly();
- v5=v5-1; /* reduce light to account for mirror */
- shade();dither();
-
- surface4:
- x1=B1[7][0];y01=B1[7][1];z1=B1[7][2];x2=B1[4][0];y2=B1[4][1];
- z2=B1[4][2];x3=B1[1][0];y3=B1[1][1];z3=B1[1][2];visibility_test();
- if (sp>0) goto surface5;
- sx1=B2[7][0];sy1=B2[7][1];sx2=B2[4][0];sy2=B2[4][1];sx3=B2[1][0];
- sy3=B2[1][1];sx4=B2[0][0];sy4=B2[0][1];sx5=B3[4][0];sy5=B3[4][1];
- draw_poly();illumination();shade();dither();
- mirror_surface4:
- sx1=B2[7][0];sy1=B2[7][1];sx2=B2[0][0];sy2=B2[0][1];sx3=B5[1][0];
- sy3=B5[1][1];sx4=B5[2][0];sy4=B5[2][1];sx5=B6[1][0];sy5=B6[1][1];
- draw_poly();
- v5=v5-1; /* reduce light to account for mirror */
- shade();dither();
-
- surface5:
- x1=B1[1][0];y01=B1[1][1];z1=B1[1][2];x2=B1[4][0];y2=B1[4][1];
- z2=B1[4][2];x3=B1[5][0];y3=B1[5][1];z3=B1[5][2];visibility_test();
- if (sp>0) goto surfaces_done;
- sx1=B2[1][0];sy1=B2[1][1];sx2=B2[4][0];sy2=B2[4][1];sx3=B2[5][0];
- sy3=B2[5][1];sx4=B2[2][0];sy4=B2[2][1];sx5=B3[5][0];sy5=B3[5][1];
- draw_poly();illumination();shade();dither();
-
- /* tidy up the background */
- setcolor(C6);setlinestyle(USERBIT_LINE,0xffff,NORM_WIDTH);
- x=-80;y=-30;z=80;calc_3d();window();moveto(sx,sy);putpixel(sx,sy,C6);
- x=80;y=-30;z=80;calc_3d();window();lineto(sx,sy);sx2=sx;sy2=sy;
- x=80;y=-30;z=180;calc_3d();window();lineto(sx,sy);sx3=sx;sy3=sy;
- x=-80;y=-30;z=180;calc_3d();window();lineto(sx,sy);sx4=sx;sy4=sy;
- x=-80;y=-30;z=80;calc_3d();window();lineto(sx,sy);sx1=sx;sy1=sy;
- x=0;y=-30;z=90;calc_3d();window();
- setfillstyle(SOLID_FILL,C0);floodfill(sx,sy,C6);
- setfillstyle(SOLID_FILL,C0);floodfill(sx,sy,C6);
- setcolor(C0);
- moveto(sx1,sy1);lineto(sx2,sy2);lineto(sx3,sy3);
- lineto(sx4,sy4);lineto(sx1,sy1);
- setcolor(C7);moveto(sx1,sy1);lineto(sx2,sy2);
-
- surfaces_done:setcolor(C7);notice(0,0);
- for (t1=1;t1!=2;) keyboard();
- quit_pgm();
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: CALCULATE SIN,COS FACTORS */
-
- void rotation(void){
- sr1=sin(r1);sr2=sin(r2);sr3=sin(r3);cr1=cos(r1);cr2=cos(r2);
- cr3=cos(r3);return;}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: STANDARD 3D FORMULAS */
- /* Pass: x,y,z cartesian world coordinates.
- Returns: sx,sy cartesian display coordinates.
- x,y,z catesian view coordinates */
-
- void calc_3d(void){
- x=(-1)*x;xa=cr1*x-sr1*z;za=sr1*x+cr1*z;x=cr2*xa+sr2*y;
- ya=cr2*y-sr2*xa;z=cr3*za-sr3*ya;y=sr3*za+cr3*ya;x=x+mx;y=y+my;
- z=z+mz;sx=d*x/z;sy=d*y/z;return;}
-
- void store_coords(void){
- rotation(); /* store roll,pitch yaw rotation factors */
- for (t=0;t<=7;t++){
- x=array1[t][0];y=array1[t][1];z=array1[t][2];
- calc_3d();window();
- B1[t][0]=x;B1[t][1]=y;B1[t][2]=z;
- B2[t][0]=sx;B2[t][1]=sy;
- };
-
- for (t=0;t<=5;t++){
- x=array2[t][0];y=array2[t][1];z=array2[t][2];
- calc_3d();window();
- B3[t][0]=sx;B3[t][1]=sy; /* store fill display coordinates */
- };
-
- for (t=0;t<=2;t++){
- x=array3[t][0];y=array3[t][1];z=array3[t][2];
- calc_3d();window();
- B5[t][0]=sx;B5[t][1]=sy; /* store mirrored display coordinates */
- };
-
- for (t=0;t<=1;t++){
- x=array4[t][0];y=array4[t][1];z=array4[t][2];
- calc_3d();window();
- B6[t][0]=sx;B6[t][1]=sy; /* store mirrored fill display coordinates */
- };
-
- return;
- }
-
- /* ----------------------------------------------------------------------- */
- /* HIDDEN SURFACE VISIBILTY TEST */
- void visibility_test(void){
- sp1=x1*(y2*z3-y3*z2);sp1=(-1)*sp1;sp2=x2*(y3*z1-y01*z3);
- sp3=x3*(y01*z2-y2*z1);sp=sp1-sp2-sp3;return;}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: DRAW 4-SIDED SOLID POLYGON IN 3D SPACE */
-
- void draw_poly(void){
- setlinestyle(USERBIT_LINE,0xffff,NORM_WIDTH);
- setfillstyle(SOLID_FILL,key_matte_clr);
- setcolor(edge_clr);
-
- polary[0][0]=sx1;polary[0][1]=sy1;polary[1][0]=sx2;polary[1][1]=sy2;
- polary[2][0]=sx3;polary[2][1]=sy3;polary[3][0]=sx4;polary[3][1]=sy4;
-
- fillpoly (4,(int far*) polary);
- return;
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINES: CALCULATE ILLUMINATION LEVEL */
-
- void illumination(void){
- xu=x2-x1;yu=y2-y01;zu=z2-z1; /* vector from vertex 1 to vertex 2 */
- xv=x3-x1;yv=y3-y01;zv=z3-z1; /* vector from vertex 1 to vertex 3 */
- xn=(yu*zv)-(zu*yv);y0n=(zu*xv)-(xu*zv);
- zn=(xu*yv)-(yu*xv); /* surface perpendicular vector */
- y0n=y0n*(-1);zn=zn*(-1); /* convert to cartesian system */
- v1=(xn*xn)+(y0n*y0n)+(zn*zn);
- v2=sqrt(v1); /* magnitude of surface perpendicular vector */
- v3=1/v2; /* ratio of magnitude to length of unit vector */
- xw=v3*xn;yw=v3*y0n;zw=v3*zn; /* surface perpendicular to unit vector */
- v4=(xw*xx1)+(yw*yy1)+(zw*zz1); /* illumination factor 0 to 1 */
- v4=v4*v6; /* expand illumination range from base 0 */
- v7=v4; /* convert illumination range to integer */
- v5=v7+1; /* illumination range from base 1 */
- return;
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: ILLUMINATION MATRIX */
-
- void shade(void){
- switch (v5){
- case 1:setfillpattern(fill_6,NORM_HUE);CC6=0x1010;values();return;
- case 2:setfillpattern(fill_6,NORM_HUE);CC6=0x1010;values();return;
- case 3:setfillpattern(fill_6,NORM_HUE);CC6=0x1010;values();return;
- case 4:setfillpattern(fill_12,NORM_HUE);CC6=0x2020;values();return;
- case 5:setfillpattern(fill_25,NORM_HUE);CC6=0x2222;values();return;
- case 6:setfillpattern(fill_37,NORM_HUE);CC6=0xaaaa;values();return;
- case 7:setfillpattern(fill_50,NORM_HUE);CC6=0xaaaa;values();return;
- case 8:setfillpattern(fill_62,NORM_HUE);CC6=0xaaaa;values();return;
- case 9:setfillpattern(fill_75,NORM_HUE);CC6=0xbbbb;values();return;
- case 10:setfillpattern(fill_87,NORM_HUE);CC6=0xdddd;values();return;
- case 11:setfillpattern(fill_93,NORM_HUE);CC6=0xefef;values();return;
- case 12:setfillpattern(fill_100,NORM_HUE);CC6=0xffff;values();return;
- defualt:setfillpattern(fill_6,NORM_HUE);CC6=0x1010;values();}
- return;}
-
- /* ----------------------------------------------------------------------- */
- /* LOCAL SUBROUTINE: prep fill and set dithering variables */
- void values(void){
- floodfill(sx5,sy5,edge_clr);CC4=BASE_CLR;CC5=NORM_HUE;return;}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: APPLY DITHERING */
-
- void dither(void){
- setlinestyle(USERBIT_LINE,0xffff,NORM_WIDTH);setcolor(CC4);
- moveto (sx1,sy1);lineto(sx2,sy2);lineto(sx3,sy3);
- lineto(sx4,sy4);lineto(sx1,sy1);
- setlinestyle(USERBIT_LINE,CC6,NORM_WIDTH);setcolor(CC5);
- moveto(sx1,sy1);lineto(sx2,sy2);lineto(sx3,sy3);
- lineto(sx4,sy4);lineto(sx1,sy1);
- return;
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: MAP CARTESIAN COORDS TO PHYSICAL SCREEN COORDS */
-
- void window(void){
- sx=sx+399;sy=sy+299;rx=screen_x/799;ry=screen_y/599;sx=sx*rx;
- sy=sy*ry;return;}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: CHACK THE KEYBOARD BUFFER */
- void keyboard(void){
- if (bioskey(1)==0) return; else quit_pgm();}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: GRACEFUL EXIT FROM PROGRAM */
-
- void quit_pgm(void){
- cleardevice();restorecrtmode();exit(0);}
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: VGA/EGA/MCGA/CGA COMPATIBILITY MODULE */
-
- void graphics_setup(void){
- int graphics_adapter,graphics_mode;
- detectgraph(&graphics_adapter,&graphics_mode);
- if (graphics_adapter==VGA) goto VGA_mode;
- if (graphics_mode==EGAHI) goto EGA_ECD_mode;
- if (graphics_mode==EGALO) goto EGA_SCD_mode;
- if (graphics_adapter==CGA) goto CGA_mode;
- if (graphics_adapter==MCGA) goto CGA_mode;
- goto abort_message;
-
- VGA_mode:
- graphics_adapter=VGA;graphics_mode=VGAHI;
- initgraph(&graphics_adapter,&graphics_mode,"");
- x_res=640;y_res=480;mode_flag=1;
- maxx=639;minx=0;maxy=479;miny=0;screen_x=639;screen_y=479;
- setcolor(7);moveto(0,30);
- outtext("Revisions by A. Helder");
- moveto(472,472);
- outtext("Press any key to quit");
- moveto(160,0);
- outtext("USING C TO GENERATE MIRROR REFLECTIONS");
- return;
-
- EGA_ECD_mode:
- graphics_adapter=EGA;graphics_mode=EGAHI;
- initgraph(&graphics_adapter,&graphics_mode,"");
- x_res=640;y_res=350;mode_flag=2;
- maxx=639;minx=0;maxy=349;miny=0;screen_x=639;screen_y=349;
- setcolor(7);moveto(0,20);
- outtext("Revisions by A. Helder");
- moveto(472,342);
- outtext ("Press any key to quit");
- moveto(160,0);
- outtext("USING C TO GENERATE MIRROR REFLECTIONS");
- return;
-
- EGA_SCD_mode:
- graphics_adapter=EGA;graphics_mode=EGALO;
- initgraph(&graphics_adapter,&graphics_mode,"");
- x_res=640;y_res=200;mode_flag=3;
- maxx=639;minx=0;maxy=199;miny=0;screen_x=639;screen_y=199;
- setcolor(7);moveto(0,19);
- outtext("Revisions by A. Helder");
- moveto(472,192);
- outtext("PRESS ANY KEY TO QUIT");
- moveto(160,0);
- outtext("USING C TO GENERATE MIRROR REFLECTIONS");
- return;
-
- CGA_mode:
- graphics_adapter=CGA;graphics_mode=CGAC3;
- initgraph(&graphics_adapter,&graphics_mode,"");
- x_res=320;y_res=200;mode_flag=1;C7=3;
- maxx=319;minx=0;maxy=199;miny=0;screen_x=319;screen_y=199;
- setcolor(3);moveto(48,10);
- outtext("Revisions by A. Helder");
- moveto(88,0);
- outtext ("MIRROR REFLECTIONS");
- return;
-
- abort_message:
- printf("\n\nUnable to proceed - Requires VGA,EGA,CGA or MCGA adapter");
- printf("\nWith appropriate monitor");
- exit(0);
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: MAP 640*480 TEMPLATE TO 2D SCREEN */
- void coords(void)
- {
- sx=sx*(x_res/640);sy=sy*(y_res/480);return;
- }
-
- /* ----------------------------------------------------------------------- */
- /* SUBROUTINE: COPYRIGHT NOTICE */
-
- int copyright[][3]={0x7c00,0x0000,0x0000,0x8231,
- 0x819c,0x645e,0xba4a,0x4252,0x96d0,0xa231,0x8252,0x955e,0xba4a,
- 0x43d2,0xf442,0x8231,0x825c,0x945e,0x7c00,0x0000,0x0000};
-
- void notice(int x, int y){
- int a,b,c; int t1=0;
-
- for (t1=0;t1<=6;t1++)
- {
- a=copyright[t1][0];b=copyright[t1][1];
- c=copyright[t1][2];
- setlinestyle(USERBIT_LINE,a,NORM_WIDTH);
- moveto(x,y);lineto(x+15,y);
- setlinestyle(USERBIT_LINE,b,NORM_WIDTH);
- moveto(x+16,y);lineto(x+31,y);
- setlinestyle(USERBIT_LINE,c,NORM_WIDTH);
- moveto(x+32,y);lineto(x+47,y);y++;
- };
- setlinestyle(USERBIT_LINE,0xFFFF,NORM_WIDTH);
- return;}
-
- /* end of source code */