Sh

Functions
void	sh3dmask (float *x, float *y, float *z, float *nrm, int *imask)
	Determines if a point is masked.
void	sharrow (float *x0, float *y0, float *z0, int *idx0, int *idy0, float *x1, float *y1, float *z1, int *idx1, int *idy1)
	Draws a curved arrow from one point to another.
void	shchar (int *ix, int *iy, float *z, char *chr, int *ired, int *igrn, int *iblu, int *iw)
	Draws a character string at a point.
void	shclippg (int *n0, float *x0, float *y0, float *z0, float *nrm0, int *ic0, int *n1, float **x1, float **y1, float **z1, float **nrm1, int **ic1)
	Clips a polygon against planes added with ?
void	shclr ()
	Clears the frame buffer to the current background color.
void	shcolor (float *x, float *y, float *z, float *n, int *rf, int *gf, int *bf, int *rb, int *gb, int *bb, int *r, int *g, int *b)
	Shades and renders a point.
void	shcs (int *ifnt, int *isx, int *isy)
	Sets the current font and shadow offset for characters.
void	shcube (float *xmin, float *xmax, float *ymin, float *ymax, float *zmin, float *zmax)
	Draw a wireframe cube.
void	shdopg (int *nt, float *xt, float *yt, float *zt, float *n)
	Shades and z-buffers a polygon.
void	shdraw (float *x, float *y, float *z, int *ip)
	A pen-up/pen-down style line drawing.
void	shfreeStrings (void)
	Clears the list of strings.
void	shend ()
	Releases any storage that has been allocated internally. shinit must be called before sny drawing is done.
void	shgetz (int *ipix, int *jpix, float *zz)
	Reads the value in the z-buffer at a pixel.
void	shinit (int *r, int *g, int *b)
	Initializes the storage needed by sh and sets the current background color.
void	shlinc (int *r, int *g, int *b)
	Sets the current line color.
void	shqlinc (int *r, int *g, int *b)
	Reads the current line color.
void	shline (float *x0, float *y0, float *z0, float *x1, float *y1, float *z1)
	shades and renders a line.
void	shline2s (int *i0, int *j0, int *i1, int *j1, int *ixoff, int *iyoff)
	Draws a shadowed line segment in screen coordinates.
void	shlineoff (float *x0, float *y0, float *z0, float *x1, float *y1, float *z1, int *ixoff, int *iyoff, int *idoff)
	Draws a line with an offset. Currently this routine just draws a line.
void	shlit (int *m, float *x, float *y, float *z, int *r, int *g, int *b)
	Assigns values to a light source.
void	shSetLineWidth (int w)
	Sets the width (in pixels) used when lines are drawn.
void	shlnonrm (float *x0, float *y0, float *z0, float *n0, float *x1, float *y1, float *z1, float *n1, int *ixoff, int *iyoff, int *idoff)
	Draws a shaded line with an offset. The offset is used to determine visability, but not applied to the line when drawn.
void	shlss (char *name, int *lname, int *pmag, int *ifnt, int *ierr, int ln)
	Loads a character set into the table of character sets.
void	shmask (int *swtch)
	Turns the 3d mask on and off.
void	shSetOutputFilename (char *name)
	Sets the root of the output file name.
void	shnlit (int *n)
	Sets the number of light sources to be used. The first n are used, the rest are not deleted, just not used.
void	shnpln (int *n)
	Sets the number of clipping planes to use. The first n are used, the rest are not deleted, just no used.
void	shpause ()
	On some systems the program can pause for the user to see the frame buffer. This just saves a version of the frame buffer with a number that is incremented by one each time pause is called.
void	shpers (float *x, float *y, float *z, float *s, float *t, float *depth)
	Performs the perspective transformation.
void	shunpers (float *ss, float *tt, float *dd, float *xx, float *yy, float *zz)
	Performs the inverse of the perspective transformation.
void	shvw (float *dist, float *alpha, float *beta, float *xmin, float *xmax, float *ymin, float *ymax, float *zmin, float *zmax)
	Sets the view point for the perspective transformation.
void	shqeye (float *xe, float *ye, float *ze)
	Returns the position of the eye.
void	shsize (float *scl, float *s0, float *t0, float *dsc, float *d0)
	Directly sets the scales used.
void	shpg (int *nv, float *xt, float *yt, float *zt)
	renders a polygon.
void	shpgc (int *r, int *g, int *b)
	Sets the color of the front face (determined by the normal) of polygons.
void	shpec (int *r, int *g, int *b)
	Sets the color of the front facing edges of polygons.
void	shbgc (int *r, int *g, int *b)
	Sets the color of the back face polygons.
void	shbec (int *r, int *g, int *b)
	Sets the color of the back facing edges of polygons.
void	shqpgc (int *r, int *g, int *b)
	Gets the color of the front face of polygons.
void	shqpec (int *r, int *g, int *b)
	Gets the color of the front facing edge of polygons.
void	shqbgc (int *r, int *g, int *b)
	Gets the color of the back face of polygons.
void	shqbec (int *r, int *g, int *b)
	Gets the color of the back facing edge of polygons.
void	shpgnrm (int *nv, float *xt, float *yt, float *zt, float *nrmt)
	Shades and z-buffers a polygon with specified normals.
void	shpl (int *n, float *x, float *y, float *z)
	Shades and draws a polyline.
void	shpln (int *i, float *ox, float *oy, float *oz, float *nx, float *ny, float *nz, int *iside)
	Sets a clipping plane.
void	shplnrm (int *nv, float *xt, float *yt, float *zt)
	Draws a polyline that is shaded according to the normals given by triplets of points (periodic). If any three consecutive points are colinear nothing is drawn.
void	shpnt (float *x, float *y, float *z)
	Draws a single point (one pixel).
void	shpntc (int *r, int *g, int *b)
	Sets the color used to draw points.
void	shqpntc (int *r, int *g, int *b)
	Reads the color used to draw points.
void	shputps (char *file, int *ln, int *m, int *n, unsigned char *imageR, unsigned char *imageG, unsigned char *imageB, int fln)
	writes an image to a postscript file.
void	shputtiff (char *file, int nrows, int ncols, unsigned char *r, unsigned char *g, unsigned char *b)
	writes an image to a tiff file (if libtiff was installed when sh was configured).
void	shqlit (int *m, float *x, float *y, float *z, int *r, int *g, int *b)
	Reads back a light source.
void	shqnpln (int *n)
	Reads back the number of clipping planes.
void	shrefsh ()
	In some graphics systems the screen can be refreshed to allow the user to see intermediate frames. In sh this routine is a noop (does nothing).
void	shSetOutputResolution (int x, int y)
	Sets the size of the screen buffer, if called BEFORE shinit. Otherwise it does nothing until the next call of shinit.
void	shsave (char *name, int *ln, int nln)
	Saves the current frame buffer. If shOutputFormat is "postscript" a postscript file is written (duh). If shOutputFormat is "tiff" and libtiff was installed before sh was configured, a tiff file is written.
void	shscal (float *xmin, float *xmax, float *ymin, float *ymax, float *zmin, float *zmax)
	Sets the visible volume.
void	shsetp (int *ipix, int *jpix, int *ired, int *igrn, int *iblu, float *zz)
	Sets a pixel in the frame buffer if zz is less than the value in the z-buffer. If the pixel is changed the the value in the z-buffer is set to zz.
void	shsphere (float *xo, float *yo, float *zo, float *r)
	Draws a tesselated sphere, 8*4**4 triangles, smoothed.
void	shsrfp (float *a, float *d, float *s, int *n)
	Set the surface properties for shading.
void	shstr (float x, float y, float z, char *str)
	Draws a horizontal character string at a point.
void	shdoStrings (FILE *fid, float x, float y)
	Writes the code in a postscript file for the strings that have been made with shstr.
void	shstrs (int xx, int yy, char *str)
	Adds a horizontal character string in screen coordinates.
void	shsync ()
	In some graphics systems this routine would wait until the physical screen buffer has been refreshed. In sh it is a noop.
void	shtri (float *xt1, float *yt1, float *zt1, float *xt2, float *yt2, float *zt2, float *xt3, float *yt3, float *zt3)
	A convenience routine to draw a triangle.
void	shtric (int *r, int *g, int *b)
	A convenience routine which sets the edge and face colors for polygons to the same color.
void	shSetOutputFormat (char *type)
	Sets the format of the file written by shsave. The default is "postscript", but if libtiff is installed when sh is configured, this routine can be used to change the format to "tiff".
void	shview (float *dist, float *alpha, float *beta, float *xmin, float *xmax, float *ymin, float *ymax, float *zmin, float *zmax)
	A convenience routine that calls shsize, shvw and shscal.

---

Function Documentation

void sh3dmask (  float *  x, float *  y, float *  z, float *  nrm, int *  imask )

Determines if a point is masked. Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. nrm  The normal to the surface the point lies on. An array of length at least 3. imask  TRUE if the point is masked, FALSE otherwise. Definition at line 31 of file sh3dmask.c. Referenced by shlnonrm(). { int cx,cy,cz; float tx,ty,tz; int n; n=5; tx=(*x)*n-floor((*x)*n); ty=(*y)*n-floor((*y)*n); tz=(*z)*n-floor((*z)*n); if(tx<0)tx=tx+1.; if(ty<0)ty=ty+1.; if(tz<0)tz=tz+1.; cx=fabs(tx-.5)<.4; cy=fabs(ty-.5)<.4; cz=fabs(tz-.5)<.4; *imask=0; if(!cx&&fabs(nrm[0])<.9)*imask=1; if(!cy&&fabs(nrm[1])<.9)*imask=1; if(!cz&&fabs(nrm[2])<.9)*imask=1; return; }

void sharrow (  float *  x0, float *  y0, float *  z0, int *  idx0, int *  idy0, float *  x1, float *  y1, float *  z1, int *  idx1, int *  idy1 )

Draws a curved arrow from one point to another. Parameters: x0  The first coordinate of the starting point. y0  The second coordinate of the starting point. z0  The third coordinate of the starting point. idx0  The first coordinate of the starting tangent to the arrow in screen coordinates. idy0  The second coordinate of the starting tangent to the arrow in screen coordinates. x1  The first coordinate of the end point. y1  The second coordinate of the end point. z1  The third coordinate of the end point. idx1  The first coordinate of the end tangent to the arrow in screen coordinates. idy1  The second coordinate of the end tangent to the arrow in screen coordinates. Definition at line 45 of file sharrow.c. References shlinc(), shline2s(), shpers(), and shqlinc(). { int rl,gl,bl; int orl=0; int ogl=0; int obl=0; float d,dx,dy; float s0=0.; float t0=0.; float d0=0.; float s1=0.; float t1=0.; float d1=0.; int ip; int i0,j0; int i1,j1; float dx0,dy0; float dx1,dy1; float s; int i,j,io,jo; int zero=0; int one=1; shpers(x0,y0,z0,&s0,&t0,&d0); i0=floor(s0*shMax+.5); j0=floor(t0*shMax+.5); shpers(x1,y1,z1,&s1,&t1,&d1); i1=floor(s1*shMax+.5); j1=floor(t1*shMax+.5); shqlinc(&orl,&ogl,&obl); d=1./sqrt((*idx0)*(*idx0)+(*idy0)*(*idy0)); dx0=(*idx0)*d; dy0=(*idy0)*d; d=1./sqrt((*idx1)*(*idx1)+(*idy1)*(*idy1)); dx1=(*idx1)*d; dy1=(*idy1)*d; ip=3; for(s=0.;s<1.;s+=.01) { i=f0(s)*i0+f1(s)*i1+f2(s)*dx0+f3(s)*dx1; j=f0(s)*j0+f1(s)*j1+f2(s)*dy0+f3(s)*dy1; if(ip==2) { shlinc(&rl,&gl,&bl); shline2s(&io,&jo,&i,&j,&one,&one); }else{ rl=orl; gl=ogl; bl=obl; shline2s(&io,&jo,&i,&j,&zero,&zero); } ip=2; io=i; jo=j; } rl=0; gl=0; bl=0; io=i1; jo=j1; i=i1-(*idx1)+.3*(*idy0); j=j1-(*idy0)-.3*(*idx0); shline2s(&io,&jo,&i,&j,&one,&one); rl=orl; gl=ogl; bl=obl; shline2s(&io,&jo,&i,&j,&zero,&zero); rl=0; gl=0; bl=0; i=i1-(*idx1)-.15*(*idy0); j=j1-(*idy0)+.15*(*idx0); shline2s(&io,&jo,&i,&j,&one,&one); rl=orl; gl=ogl; bl=obl; shline2s(&io,&jo,&i,&j,&zero,&zero); rl=0; gl=0; bl=0; io=i1-(*idx1)-.15*(*idy0); jo=j1-(*idy0)+.15*(*idx0); i=i1-(*idx1)+.15*(*idy0); j=j1-(*idy0)-.15*(*idx0); shline2s(&io,&jo,&i,&j,&one,&one); rl=orl; gl=ogl; bl=obl; shline2s(&io,&jo,&i,&j,&zero,&zero); return; } Here is the call graph for this function:

void shbec (  int *  r, int *  g, int *  b )

Sets the color of the back facing edges of polygons. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 93 of file shpgc.c. Referenced by shtric(). { sh_rd[3] =*r; sh_gd[3] =*g; sh_bd[3] =*b; sh_ram[3]=*r; sh_gam[3]=*g; sh_bam[3]=*b; return; }

void shbgc (  int *  r, int *  g, int *  b )

Sets the color of the back face polygons. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 75 of file shpgc.c. Referenced by shtric(). { sh_rd[2] =*r; sh_gd[2] =*g; sh_bd[2] =*b; sh_ram[2]=*r; sh_gam[2]=*g; sh_bam[2]=*b; return; }

void shchar (  int *  ix, int *  iy, float *  z, char *  chr, int *  ired, int *  igrn, int *  iblu, int *  iw )

Draws a character string at a point. Parameters: ix  The first coordinate of the screen coordinates of the character string. iy  The second coordinate of the screen coordinates of the character string. z  The depth of the character screen. chr  The character string. ired  The red component of the color of the string. igrn  The red component of the color of the string. iblu  The red component of the color of the string. iw  The width in pixels of the string. Definition at line 33 of file shchar.c. References shsetp(). { /* places a character string in the image */ /* Make a list, with ix,iy,iz and the string */ /* At the end of the postsript file, loop */ /* over list, if iz is visible, place postscript commands to draw the string in the output file */ /* places a character in the sh_image */ int xref=0; int yref=0; int ipwid=0; int iphgt=0; int width=0; int ibyte; int ichr; int ifnt; int idir; int xx,yy; float zz; int zero=0; int i,j; int jx,jy; int ib; int level; float conv; /* conv converts from sp to pixels */ conv=240./72.27/65536.; ifnt=sh_kfont[sh_cFont]; /* Blank is width of lower case m (?) */ ichr=chr[0]; if(chr[0]==' ')ichr='m'; idir=sh_dir[ifnt]+4+ichr*16; sscanf(sh_image+idir,"%4.4d",&width); *iw=width*(sh_fmag[ifnt])*0.001*1.2*conv; if(chr[0]==' ')return; sscanf(sh_image+idir+ 4,"%4.4d",&xref); sscanf(sh_image+idir+ 6,"%4.4d",&yref); sscanf(sh_image+idir+ 8,"%4.4d",&ipwid); sscanf(sh_image+idir+10,"%4.4d",&iphgt); sscanf(sh_image+idir+12,"%4.4d",&ibyte); ibyte=ibyte+sh_dir[ifnt]-1; jx=(*ix)-xref+.5; jy=(*iy)+yref+.5; /* put shadow in frame */ if(shadow_idx!=0||shadow_idy!=0) { ib=ibyte; for(j=0;j<iphgt;j++) { for(i=0;i<ipwid;i++) { level=sh_image[ib]; if(level!=0) { zz=*z-.01; xx=jx+i+shadow_idx; yy=jy-j+shadow_idy; shsetp(&xx,&yy,&zero,&zero,&zero,&zz); } ib++; } } } /* put character in frame */ ib=ibyte; for(j=0;j<iphgt;j++) { for(i=0;i<ipwid;i++) { level=sh_image[ib]; if(level!=0) { zz=*z-.01; xx=jx+i; yy=jy-j; shsetp(&xx,&yy,ired,igrn,iblu,&zz); } ib++; } } return; } Here is the call graph for this function:

void shclippg (  int *  n0, float *  x0, float *  y0, float *  z0, float *  nrm0, int *  ic0, int *  n1, float **  x1, float **  y1, float **  z1, float **  nrm1, int **  ic1 )

Clips a polygon against planes added with ? Parameters: n0  The number of vertices in the polygon. x0  An array of first coordinates the vertices. An array of length at least n0. y0  An array of second coordinates the vertices. An array of length at least n0. z0  An array of third coordinates the vertices. An array of length at least n0. nrm0  The normals to the surface at the the vertices lies on. An array of length at least 3*n0. ic0  An array ? n1  The number of vertices in the clipped polygon. x1  An array of first coordinates the vertices of the clipped polygon. An array of length n1 allocated in the routine. y1  An array of second coordinates the vertices of the clipped polygon. An array of length n1 allocated in the routine. z1  An array of third coordinates the vertices of the clipped polygon. An array of length n1 allocated in the routine. nrm1  The normals to the surface at the the vertices lies on. An array of length at least 3*n0 allocated in the routine. ic1  An array ? allocated in the routine. Definition at line 45 of file shclippg.c. Referenced by shpgnrm(). { /* clip a polygon */ float *x2,*y2,*z2; float *nrm2; int *ic2; int in0,in1; int m1,m2; int i,j,ip; float direc; float t; int i0; int verbose; verbose=0; m2=maxpol; x2=(float*)malloc(m2*sizeof(float)); y2=(float*)malloc(m2*sizeof(float)); z2=(float*)malloc(m2*sizeof(float)); nrm2=(float*)malloc(3*m2*sizeof(float)); ic2=(int*)malloc(m2*sizeof(int)); m1=(*n0)+1; if(*x1==(float*)NULL)*x1=(float*)malloc(m1*sizeof(float)); if(*y1==(float*)NULL)*y1=(float*)malloc(m1*sizeof(float)); if(*z1==(float*)NULL)*z1=(float*)malloc(m1*sizeof(float)); if(*nrm1==(float*)NULL)*nrm1=(float*)malloc(3*m1*sizeof(float)); if(*ic1==(int*)NULL)*ic1=(int*)malloc(m1*sizeof(int)); /* Clipping Planes */ /* NPLNS Number of clipping planes */ /* PLNO[2+3*M) point on plane M */ /* PLNN[2+3*M) normal to plane M */ /* IPLN(M) side of plane M to clip (IPLN*PLNN<0 ==> clipped) */ /* oper(M) whether clipping is ored or anded. */ /* copy polygon 0 into 1, removing near duplicates. */ *n1=0; (*x1)[*n1]=x0[0]; (*y1)[*n1]=y0[0]; (*z1)[*n1]=z0[0]; (*nrm1)[0+3*(*n1)]=nrm0[0]; (*nrm1)[1+3*(*n1)]=nrm0[1]; (*nrm1)[2+3*(*n1)]=nrm0[2]; (*ic1)[*n1]=ic0[0]; (*n1)++; for(i=1;i<*n0;i++) { if(fabs(x0[i]-x0[i-1])+fabs(y0[i]-y0[i-1])+fabs(z0[i]-z0[i-1])>1.e-6) { (*x1)[*n1]=x0[i]; (*y1)[*n1]=y0[i]; (*z1)[*n1]=z0[i]; (*nrm1)[ 3*(*n1)]=nrm0[ 3*i]; (*nrm1)[1+3*(*n1)]=nrm0[1+3*i]; (*nrm1)[2+3*(*n1)]=nrm0[2+3*i]; (*ic1)[(*n1)]=ic0[i]; (*n1)++; } } /* Duplicate first point */ (*x1)[*n1]=x0[0]; (*y1)[*n1]=y0[0]; (*z1)[*n1]=z0[0]; (*nrm1)[ 3*(*n1)]=nrm0[ 3*0]; (*nrm1)[1+3*(*n1)]=nrm0[1+3*0]; (*nrm1)[2+3*(*n1)]=nrm0[2+3*0]; (*ic1)[(*n1)]=ic0[0]; (*n1)++; /* Loop over the clipping planes */ /* in0=this point visible, in1=next point visible */ if(sh_nplns==0) { free(x2); free(y2); free(z2); free(nrm2); free(ic2); return; } for(ip=0;ip<sh_nplns;ip++) { if(verbose)printf("plane %d is ((x,y,z)-(%f,%f,%f)).(%f,%f,%f)*%d>0\n",ip,sh_plno[ 3*ip],sh_plno[1+3*ip],sh_plno[2+3*ip],sh_plnn[ 3*ip],sh_plnn[1+3*ip],sh_plnn[2+3*ip],sh_ipln[ip]); in0=0; i0=0; for(i=0;i<*n1;i++) { if(verbose){printf("vertex %d(%d) (%f,%f,%f) is ",i,*n1,(*x1)[i],(*y1)[i],(*z1)[i]);fflush(stdout);} direc=sh_plnn[ 3*ip]*((*x1)[i]-sh_plno[ 3*ip])+sh_plnn[1+3*ip]*((*y1)[i]-sh_plno[1+3*ip])+sh_plnn[2+3*ip]*((*z1)[i]-sh_plno[2+3*ip]); in1=sh_ipln[ip]*direc>0; if(verbose){if(in1)printf("in\n");else printf("out\n");} if(!in0)i0=i; in0=in0||in1; } if(verbose)fflush(stdout); if(!in0) { if(verbose){printf("Polygon completely outside\n\n");fflush(stdout);} *n1=0; free(x2); free(y2); free(z2); free(nrm2); free(ic2); return; } if(i0>0)i0--; if(verbose){printf("Starting at vertex %d\n",i0);fflush(stdout);} #define JUNK #ifdef JUNK j=0; direc=sh_plnn[ 3*ip]*((*x1)[i0]-sh_plno[ 3*ip])+sh_plnn[1+3*ip]*((*y1)[i0]-sh_plno[1+3*ip])+sh_plnn[2+3*ip]*((*z1)[i0]-sh_plno[2+3*ip]); in0=sh_ipln[ip]*direc>0; i=i0; while(i+1<*n1) { if(verbose){printf("Vertex %d(%d) (%f,%f,%f) next is %d (%f,%f,%f)\n",i,*n1,(*x1)[i],(*y1)[i],(*z1)[i],i+1,(*x1)[i+1],(*y1)[i+1],(*z1)[i+1]);fflush(stdout);} direc=sh_plnn[ 3*ip]*((*x1)[i+1]-sh_plno[ 3*ip]) +sh_plnn[1+3*ip]*((*y1)[i+1]-sh_plno[1+3*ip]) +sh_plnn[2+3*ip]*((*z1)[i+1]-sh_plno[2+3*ip]); in1=sh_ipln[ip]*direc>0; if(verbose) { if(in0){printf(" %d is in\n",i);fflush(stdout);} if(!in0){printf(" %d is out\n",i);fflush(stdout);} if(in1){printf(" %d is in\n",i+1);fflush(stdout);} if(!in1){printf(" %d is out\n",i+1);fflush(stdout);} } if(!in0&&!in1) { if(verbose){printf(" !in0&&!in1\n");fflush(stdout);} in0=in1; i++; }else if(in0&&in1) { if(verbose){printf(" in0&&in1\n");fflush(stdout);} in0=in1; if(j+1>m2) { m2+=5; x2=(float*)realloc((void*)x2,m2*sizeof(float)); y2=(float*)realloc((void*)y2,m2*sizeof(float)); z2=(float*)realloc((void*)z2,m2*sizeof(float)); nrm2=(float*)realloc((void*)nrm2,3*m2*sizeof(float)); ic2=(int*)realloc((void*)ic2,m2*sizeof(int)); } x2[j]=(*x1)[i]; y2[j]=(*y1)[i]; z2[j]=(*z1)[i]; nrm2[ 3*j]=(*nrm1)[ 3*i]; nrm2[1+3*j]=(*nrm1)[1+3*i]; nrm2[2+3*j]=(*nrm1)[2+3*i]; ic2[j]=(*ic1)[i]; j++; i++; }else if(!in0&&in1) { if(verbose){printf(" !in0&&in1\n");fflush(stdout);} t=-( sh_plnn[ 3*ip]*((*x1)[i]-sh_plno[ 3*ip]) +sh_plnn[1+3*ip]*((*y1)[i]-sh_plno[1+3*ip]) +sh_plnn[2+3*ip]*((*z1)[i]-sh_plno[2+3*ip]) ) /( sh_plnn[ 3*ip]*((*x1)[i+1]-(*x1)[i]) +sh_plnn[1+3*ip]*((*y1)[i+1]-(*y1)[i]) +sh_plnn[2+3*ip]*((*z1)[i+1]-(*z1)[i]) ); (*x1)[i]=(*x1)[i]+t*((*x1)[i+1]-(*x1)[i]); (*y1)[i]=(*y1)[i]+t*((*y1)[i+1]-(*y1)[i]); (*z1)[i]=(*z1)[i]+t*((*z1)[i+1]-(*z1)[i]); (*nrm1)[ 3*i]=(*nrm1)[ 3*i]+t*((*nrm1)[ 3*(i+1)]-(*nrm1)[ 3*i]); (*nrm1)[1+3*i]=(*nrm1)[1+3*i]+t*((*nrm1)[1+3*(i+1)]-(*nrm1)[1+3*i]); (*nrm1)[2+3*i]=(*nrm1)[2+3*i]+t*((*nrm1)[2+3*(i+1)]-(*nrm1)[2+3*i]); (*ic1)[i]=-1; if(j+1>m2) { m2+=5; x2=(float*)realloc((void*)x2,m2*sizeof(float)); y2=(float*)realloc((void*)y2,m2*sizeof(float)); z2=(float*)realloc((void*)z2,m2*sizeof(float)); nrm2=(float*)realloc((void*)nrm2,3*m2*sizeof(float)); ic2=(int*)realloc((void*)ic2,m2*sizeof(int)); } x2[j]=(*x1)[i]; y2[j]=(*y1)[i]; z2[j]=(*z1)[i]; nrm2[ 3*j]=(*nrm1)[ 3*i]; nrm2[1+3*j]=(*nrm1)[1+3*i]; nrm2[2+3*j]=(*nrm1)[2+3*i]; ic2[j]=(*ic1)[i]; j++; in0=1; i++; }else if(in0&&!in1) { if(verbose){printf(" in0&&!in1\n");fflush(stdout);} t=-( sh_plnn[ 3*ip]*((*x1)[i]-sh_plno[ 3*ip]) +sh_plnn[1+3*ip]*((*y1)[i]-sh_plno[1+3*ip]) +sh_plnn[2+3*ip]*((*z1)[i]-sh_plno[2+3*ip]) ) /( sh_plnn[ 3*ip]*((*x1)[i+1]-(*x1)[i]) +sh_plnn[1+3*ip]*((*y1)[i+1]-(*y1)[i]) +sh_plnn[2+3*ip]*((*z1)[i+1]-(*z1)[i]) ); if(j+2>m2) { m2+=5; x2=(float*)realloc((void*)x2,m2*sizeof(float)); y2=(float*)realloc((void*)y2,m2*sizeof(float)); z2=(float*)realloc((void*)z2,m2*sizeof(float)); nrm2=(float*)realloc((void*)nrm2,3*m2*sizeof(float)); ic2=(int*)realloc((void*)ic2,m2*sizeof(int)); } x2[j]=(*x1)[i]; y2[j]=(*y1)[i]; z2[j]=(*z1)[i]; nrm2[ 3*j]=(*nrm1)[ 3*i]; nrm2[1+3*j]=(*nrm1)[1+3*i]; nrm2[2+3*j]=(*nrm1)[2+3*i]; ic2[j]=(*ic1)[i]; j++; (*x1)[i]=(*x1)[i]+t*((*x1)[i+1]-(*x1)[i]); (*y1)[i]=(*y1)[i]+t*((*y1)[i+1]-(*y1)[i]); (*z1)[i]=(*z1)[i]+t*((*z1)[i+1]-(*z1)[i]); (*nrm1)[ 3*i]=(*nrm1)[ 3*i]+t*((*nrm1)[ 3*(i+1)]-(*nrm1)[ 3*i]); (*nrm1)[1+3*i]=(*nrm1)[1+3*i]+t*((*nrm1)[1+3*(i+1)]-(*nrm1)[1+3*i]); (*nrm1)[2+3*i]=(*nrm1)[2+3*i]+t*((*nrm1)[2+3*(i+1)]-(*nrm1)[2+3*i]); (*ic1)[i]=1; x2[j]=(*x1)[i]; y2[j]=(*y1)[i]; z2[j]=(*z1)[i]; nrm2[ 3*j]=(*nrm1)[ 3*i]; nrm2[1+3*j]=(*nrm1)[1+3*i]; nrm2[2+3*j]=(*nrm1)[2+3*i]; ic2[j]=(*ic1)[i]; j++; in0=0; } } /* End of while(i+1<*n1) */ if(verbose){printf("done plane %d, Caught %d vertices\n",ip,j);fflush(stdout);} /* Copy back from 2 into 1 */ if(j+1>m1) { m1=j+3; *x1=(float*)realloc((void*)(*x1),m1*sizeof(float)); *y1=(float*)realloc((void*)(*y1),m1*sizeof(float)); *z1=(float*)realloc((void*)(*z1),m1*sizeof(float)); *nrm1=(float*)realloc((void*)(*nrm1),3*m1*sizeof(float)); *ic1=(int*)realloc((void*)(*ic1),m1*sizeof(int)); } for(i=0;i<j;i++) { (*x1)[i]=x2[i]; (*y1)[i]=y2[i]; (*z1)[i]=z2[i]; (*nrm1)[ 3*i]=nrm2[ 3*i]; (*nrm1)[1+3*i]=nrm2[1+3*i]; (*nrm1)[2+3*i]=nrm2[2+3*i]; (*ic1)[i]=ic2[i]; } (*x1)[j]=x2[0]; (*y1)[j]=y2[0]; (*z1)[j]=z2[0]; (*nrm1)[ 3*j]=nrm2[ 3*0]; (*nrm1)[1+3*j]=nrm2[1+3*0]; (*nrm1)[2+3*j]=nrm2[2+3*0]; (*ic1)[j]=ic2[0]; (*n1)=j+1; #endif } /* End of for(ip=0;ip<sh_nplns;ip++) */ if(verbose){printf("Part of Polygon remains unclipped, %d vertices\n\n",j);fflush(stdout);} free(x2); free(y2); free(z2); free(nrm2); free(ic2); return; }

void shcolor (  float *  x, float *  y, float *  z, float *  n, int *  rf, int *  gf, int *  bf, int *  rb, int *  gb, int *  bb, int *  r, int *  g, int *  b )

Shades and renders a point. Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. n  The normal to the surface the point lies on. An array of length at least 3. rf  The red component of the front face color. gf  The green component of the front face color. bf  The blue component of the front face color. rb  The red component of the back face color. gb  The green component of the back face color. bb  The blue component of the back face color. r  (output) The red component of the shaded point. g  (output) The green component of the shaded point. b  (output) The blue component of the shaded point. Definition at line 39 of file shcolor.c. Referenced by shlnonrm(). { /* Shade and render a point. */ int rd[2]={0,0}; int gd[2]={0,0}; int bd[2]={0,0}; /* Geometric Information */ /* EYE[2] position of the viewer */ /* NLIT number of light sources */ /* LIT(3,M) position of light source m */ /* RS(M) color of light source m */ /* GS(M) */ /* BS(M) */ /* AM intensity of ambient illumination */ /* AD intensity of diffuse illumination */ /* AS intensity of specular illumination */ /* ND exponent of specular illumination */ /* Clipping Planes */ /* NPLNS Number of clipping planes */ /* PLNO(3,M) point on plane M */ /* PLNN(3,M) normal to plane M */ /* IPLN(M) side of plane M to clip (IPLN*PLNN<0 ==> clipped) */ /* oper(M) whether clipping is ored or anded. */ float e[3]={0.,0.,0.}; float l[3]={0.,0.,0.}; float h[3]={0.,0.,0.}; float an; float ah; float ae; float sd; float ssr; float ssg; float ssb; int k; float xxx=0.; float yyy=0.; float zzz=0.; float al; int icol; an=sqrt(n[0]*n[0]+n[1]*n[1]+n[2]*n[2]); if(an<1.e-7) { *r=*rf; *g=*gf; *b=*bf; return; } /* e is the vector from the pixel to the eye */ e[0]=shv_eye[0]-*x; e[1]=shv_eye[1]-*y; e[2]=shv_eye[2]-*z; ae=1./sqrt(e[0]*e[0]+e[1]*e[1]+e[2]*e[2]); e[0]=e[0]*ae; e[1]=e[1]*ae; e[2]=e[2]*ae; sd =0.; ssr=0.; ssg=0.; ssb=0.; for(k=0;k<sh_nlit;k++) { /* l is the vector from the pixel to light source k */ if(sh_type[k]==0) { l[0]=sh_lit[ 3*k]-xxx; l[1]=sh_lit[1+3*k]-yyy; l[2]=sh_lit[2+3*k]-zzz; }else{ l[0]=sh_lit[ 3*k]; l[1]=sh_lit[1+3*k]; l[2]=sh_lit[2+3*k]; } al=1./sqrt(l[0]*l[0]+l[1]*l[1]+l[2]*l[2]); l[0]=l[0]*al; l[1]=l[1]*al; l[2]=l[2]*al; /* printf("light %d, l=(%f,%f,%f) n=(%f,%f,%f\n",k,l[0],l[1],l[2],n[0],n[1],n[2]);*/ /* h will control the specular reflection */ h[0]=e[0]+l[0]; h[1]=e[1]+l[1]; h[2]=e[2]+l[2]; ah=1./sqrt(h[0]*h[0]+h[1]*h[1]+h[2]*h[2]); h[0]=h[0]*ah; h[1]=h[1]*ah; h[2]=h[2]*ah; /* printf("h=(%f,%f,%f)\n",h[0],h[1],h[2]);*/ /* printf("sd=%f sscol=(%f,%f,%f)\n",fabs(n[0]*l[0]+n[1]*l[1]+n[2]*l[2]),sh_rs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd),sh_gs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd),sh_bs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd));*/ sd =sd +fabs(n[0]*l[0]+n[1]*l[1]+n[2]*l[2]); ssr=ssr+sh_rs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd); ssg=ssg+sh_gs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd); ssb=ssb+sh_bs[k]*pow(fabs(n[0]*h[0]+n[1]*h[1]+n[2]*h[2]),sh_nd); } /*printf("total sd=%f sscol=(%f,%f,%f)\n",sd,ssr,ssg,ssb);*/ icol=0; if(n[0]*(shv_eye[0]-*x)+n[1]*(shv_eye[1]-*y)+n[2]*(shv_eye[2]-*z)<0.)icol=1; /*printf("side %d\n",icol);*/ rd[0]=*rf; gd[0]=*gf; bd[0]=*bf; rd[1]=*rb; gd[1]=*gb; bd[1]=*bb; *r=sh_am*rd[icol]+sh_ad*rd[icol]*sd/sh_nlit+sh_as*ssr/sh_nlit; *g=sh_am*gd[icol]+sh_ad*gd[icol]*sd/sh_nlit+sh_as*ssg/sh_nlit; *b=sh_am*bd[icol]+sh_ad*bd[icol]*sd/sh_nlit+sh_as*ssb/sh_nlit; if(*r>255)*r=255; if(*g>255)*g=255; if(*b>255)*b=255; if(*r<0)*r=0; if(*g<0)*g=0; if(*b<0)*b=0; return; }

void shcs (  int *  ifnt, int *  isx, int *  isy )

Sets the current font and shadow offset for characters. Parameters: ifnt  The font number. isx  The horizontal shadow offset in pixels. isy  The vertical shadow offset in pixels. Definition at line 28 of file shcs.c. { sh_cFont=*ifnt; shadow_idx=*isx; shadow_idy=*isy; return; }

void shcube (  float *  xmin, float *  xmax, float *  ymin, float *  ymax, float *  zmin, float *  zmax )

Draw a wireframe cube. Parameters: xmax  The first coordinate of the furthest extreme of the cube. ymax  The second coordinate of the furthest extreme of the cube. zmax  The third coordinate of the furthest extreme of the cube. xmin  The first coordinate of the least extreme of the cube. ymin  The second coordinate of the least extreme of the cube. zmin  The third coordinate of the least extreme of the cube. Definition at line 31 of file shcube.c. References shline(). { /* Draw a wireframe cube. */ shline(xmin,ymin,zmin,xmax,ymin,zmin); shline(xmax,ymin,zmin,xmax,ymin,zmax); shline(xmax,ymin,zmax,xmin,ymin,zmax); shline(xmin,ymin,zmax,xmin,ymin,zmin); shline(xmin,ymax,zmin,xmax,ymax,zmin); shline(xmax,ymax,zmin,xmax,ymax,zmax); shline(xmax,ymax,zmax,xmin,ymax,zmax); shline(xmin,ymax,zmax,xmin,ymax,zmin); shline(xmin,ymin,zmin,xmin,ymax,zmin); shline(xmin,ymin,zmax,xmin,ymax,zmax); shline(xmax,ymin,zmin,xmax,ymax,zmin); shline(xmax,ymin,zmax,xmax,ymax,zmax); return; } Here is the call graph for this function:

void shdopg (  int *  nt, float *  xt, float *  yt, float *  zt, float *  n )

Shades and z-buffers a polygon. Parameters: nt  The number of vertices in the polygon. xt  The first coordinate of the vertices of the polygon. An array of length at least nt. yt  The second coordinate of the vertices of the polygon. An array of length at least nt. zt  The third coordinate of the vertices of the polygon. An array of length at least nt. n  The normal to the surface the vertices lie on. An array of length at least 3*nt. Definition at line 63 of file shdopg.c. Referenced by shpgnrm(). { /* shade and z-buffer a polygon defined by the vertices xt,yt,zt. */ /* with specified normal */ float nn[3]={1.,0.,0.}; float t; float de0,de1; float xx,yy,zz; float xxx=0.; float yyy=0.; float zzz=0.; float zbuf=0.; int r=0; int g=0; int b=0; int nv,imax,jmax,imin,jmin; int m,mi,mj; int mleft,mright; int kt; int nedges; float an; int imask=0; int ipix,jpix; /* Project the vertices into Perspective coordinates */ /* Transform into raster coordinates */ /* SHk0[m] points to one point of each edge, */ /* SHk1[m] points to the other */ #ifdef VERBOSE printf(" shdopg n=%d\n",*nt); for(m=0;m<*nt;m++) { printf(" %d (%f,%f,%f) (%f,%f,%f)\n",m,xt[m],yt[m],zt[m],n[3*m],n[1+3*m],n[2+3*m]); fflush(stdout); } #endif if(*nt>SHmaxpol) { SHmaxpol+=1000; SHx=(float*)realloc((void*)SHx,SHmaxpol*sizeof(float)); if(SHx==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHy=(float*)realloc((void*)SHy,SHmaxpol*sizeof(float)); if(SHy==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHz=(float*)realloc((void*)SHz,SHmaxpol*sizeof(float)); if(SHz==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHi=(int*)realloc((void*)SHi,SHmaxpol*sizeof(int)); if(SHi==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHj=(int*)realloc((void*)SHj,SHmaxpol*sizeof(int)); if(SHj==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHk0=(int*)realloc((void*)SHk0,SHmaxpol*sizeof(int)); if(SHk0==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHk1=(int*)realloc((void*)SHk1,SHmaxpol*sizeof(int)); if(SHk1==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHke=(int*)realloc((void*)SHke,SHmaxpol*sizeof(int)); if(SHke==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHo=(float*)realloc((void*)SHo,SHmaxpol*sizeof(float)); if(SHo==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHs=(float*)realloc((void*)SHs,SHmaxpol*sizeof(float)); if(SHs==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHxo=(float*)realloc((void*)SHxo,SHmaxpol*sizeof(float)); if(SHxo==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHyo=(float*)realloc((void*)SHyo,SHmaxpol*sizeof(float)); if(SHyo==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHzo=(float*)realloc((void*)SHzo,SHmaxpol*sizeof(float)); if(SHzo==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHno=(float*)realloc((void*)SHno,(2+3*SHmaxpol)*sizeof(float)); if(SHno==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHxs=(float*)realloc((void*)SHxs,SHmaxpol*sizeof(float)); if(SHxs==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHys=(float*)realloc((void*)SHys,SHmaxpol*sizeof(float)); if(SHys==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHzs=(float*)realloc((void*)SHzs,SHmaxpol*sizeof(float)); if(SHzs==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHns=(float*)realloc((void*)SHns,(2+3*SHmaxpol)*sizeof(float)); if(SHns==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHi0=(int*)realloc((void*)SHi0,SHmaxpol*sizeof(int)); if(SHi0==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHx0=(float*)realloc((void*)SHx0,SHmaxpol*sizeof(float)); if(SHx0==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHy0=(float*)realloc((void*)SHy0,SHmaxpol*sizeof(float)); if(SHy0==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHz0=(float*)realloc((void*)SHz0,SHmaxpol*sizeof(float)); if(SHz0==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHn0=(float*)realloc((void*)SHn0,(2+3*SHmaxpol)*sizeof(float)); if(SHn0==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHs0=(float*)realloc((void*)SHs0,SHmaxpol*sizeof(float)); if(SHs0==(float*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} SHiflg=(int*)realloc((void*)SHiflg,SHmaxpol*sizeof(int)); if(SHiflg==(int*)NULL){fprintf(stderr,"shdopg: out of memory!");fflush(stderr);return;} } nv=*nt; if(fabs(xt[0]-xt[nv-1])+fabs(yt[0]-yt[nv-1])+fabs(zt[0]-zt[nv-1])<3.e-7)nv--; imax=0; jmax=0; imin=shIMax; jmin=shJMax; for(m=0;m<nv;m++) { shpers(xt+m,yt+m,zt+m,SHx+m,SHy+m,SHz+m); SHi[m]=SHx[m]*shMax+.5; SHj[m]=SHy[m]*shMax+.5; SHk0[m]=m; SHk1[m]=m+1; if(m==nv-1)SHk1[m]=0; if(SHi[m]>imax)imax=SHi[m]; if(SHj[m]>jmax)jmax=SHj[m]; if(SHi[m]<imin)imin=SHi[m]; if(SHj[m]<jmin)jmin=SHj[m]; } #ifdef VERBOSE printf(" shdopg n=%d\n",nv); for(m=0;m<nv;m++) { printf(" %d (%d,%d)\n",m,SHi[m],SHj[m]); fflush(stdout); } #endif /* Quick clip if polygon is not in image */ if((imax<0||imin>=shIMax||jmax<0||jmin>=shJMax)) { #ifdef VERBOSE printf(" done shdopg -- polygon out of range (%d,%d)->(%d,%d), (0,0)->(%d,%d)\n",imin,jmin,imax,jmax,shIMax,shJMax); #endif return; } if(jmax==jmin||nv<3) { shpl(&nv,xt,yt,zt); #ifdef VERBOSE printf(" done shdopg -- too few vertices, or horiSHzontal line, nv=%d, j in (%d,%d)\n",nv,jmin,jmax); #endif return; } /* flip edges so that SHy[SHk0[m]]<SHy[SHk1[m]] */ for(m=0;m<nv;m++) { if(SHj[SHk0[m]]>SHj[SHk1[m]]) { kt=SHk0[m]; SHk0[m]=SHk1[m]; SHk1[m]=kt; } } /* sort edges by increasing SHy[SHk0[m]] */ for(mi=0;mi<nv-1;mi++) { for(mj=mi+1;mj<nv;mj++) { if(SHj[SHk0[mj]]<SHj[SHk0[mi]]) { kt=SHk0[mi]; SHk0[mi]=SHk0[mj]; SHk0[mj]=kt; kt=SHk1[mi]; SHk1[mi]=SHk1[mj]; SHk1[mj]=kt; } } } #ifdef VERBOSE printf(" shdopg edges\n"); for(m=0;m<nv;m++) { printf(" %d (%d,%d) (%f,%f,%f)->(%f,%f,%f)\n",m,SHk0[m],SHk1[m],SHx[SHk0[m]],SHy[SHk0[m]],SHz[SHk0[m]],SHx[SHk1[m]],SHy[SHk1[m]],SHz[SHk1[m]]); fflush(stdout); } #endif /* Normalize the normal vector at each vertex */ for(m=0;m<nv;m++) { an=sqrt(n[3*m]*n[3*m]+n[1+3*m]*n[1+3*m]+n[2+3*m]*n[2+3*m]); if(fabs(an)<1.e-6) { shpl(&nv,xt,yt,zt); #ifdef VERBOSE printf(" done shdopg -- flat line m=%d, n=(%f,%f,%f), |n|=%f\n",m,n[3*m],n[1+3*m],n[2+3*m],an); #endif return; }else{ an=1./an; n[3*m]=n[3*m]*an; n[1+3*m]=n[1+3*m]*an; n[2+3*m]=n[2+3*m]*an; } /* Compute the equations of the edges of the polygon */ /* (x,y,z)=(SHxo[m],SHyo[m],SHzo[m])+j*(SHxs[m],SHys[m],SHzs[m]) */ if(SHj[SHk1[m]]==SHj[SHk0[m]]) { SHs[m]=0.; SHxs[m]=0.; SHys[m]=0.; SHzs[m]=0.; SHns[3*m]=0; SHns[1+3*m]=0; SHns[2+3*m]=0; }else{ t=1./(float)(SHj[SHk1[m]]-SHj[SHk0[m]]); SHs[m]= (SHi[SHk1[m]]-SHi[SHk0[m]])*t; SHxs[m]= (SHx[SHk1[m]]-SHx[SHk0[m]])*t; SHys[m]= (SHy[SHk1[m]]-SHy[SHk0[m]])*t; SHzs[m]= (SHz[SHk1[m]]-SHz[SHk0[m]])*t; SHns[3*m]= (n[ 3*SHk1[m]]-n[ 3*SHk0[m]])*t; SHns[1+3*m]=(n[1+3*SHk1[m]]-n[1+3*SHk0[m]])*t; SHns[2+3*m]=(n[2+3*SHk1[m]]-n[2+3*SHk0[m]])*t; } SHo[m]= SHi[SHk0[m]]-SHj[SHk0[m]]*SHs[m]; SHxo[m]= SHx[SHk0[m]]-SHj[SHk0[m]]*SHxs[m]; SHyo[m]= SHy[SHk0[m]]-SHj[SHk0[m]]*SHys[m]; SHzo[m]= SHz[SHk0[m]]-SHj[SHk0[m]]*SHzs[m]; SHno[3*m]= n[ 3*SHk0[m]]-SHj[SHk0[m]]*SHns[3*m]; SHno[1+3*m]=n[1+3*SHk0[m]]-SHj[SHk0[m]]*SHns[1+3*m]; SHno[2+3*m]=n[2+3*SHk0[m]]-SHj[SHk0[m]]*SHns[2+3*m]; } #ifdef VERBOSE printf(" shdopg equations\n"); for(m=0;m<nv;m++) { printf(" %4.4d s %f, SHxs (%f,%f,%f), SHns (%f,%f,%f)\n",m,SHs[m],SHxs[m],SHys[m],SHzs[m],SHns[3*m],SHns[1+3*m],SHns[2+3*m]); printf(" o %f, SHxo (%f,%f,%f), SHno (%f,%f,%f)\n",SHo[m],SHxo[m],SHyo[m],SHzo[m],SHno[3*m],SHno[1+3*m],SHno[2+3*m]); fflush(stdout); } #endif /* Now fill the polygon, bottom to top */ for(jpix=jmin;jpix<=jmax;jpix++) { /* Find the range of edges in k which cross this value of jpix */ mleft=0; mright=nv-1; for(m=0;m<nv;m++) { if(SHj[SHk1[mleft]]<jpix)mleft++; if(SHj[SHk0[mright]]>jpix)mright--; } /* Compute the intersection points on each edge */ #ifdef VERBOSE printf("scan line %d, left %d, right %d\n",jpix,mleft,mright); fflush(stdout); #endif nedges=0; for(m=mleft;m<=mright;m++) { if(SHj[SHk1[m]]>=jpix) { if(SHj[SHk1[m]]==jpix&&SHj[SHk0[m]]==jpix) { SHi0[nedges]= SHi[SHk0[m]]; SHiflg[nedges]=3; SHs0[nedges]=fabs( SHi[SHk1[m]] - SHi[SHk0[m]] )+.5; SHx0[nedges]=SHx[SHk0[m]]; SHy0[nedges]=SHy[SHk0[m]]; SHz0[nedges]=SHz[SHk0[m]]; SHn0[3*nedges]=n[3*SHk0[m]]; SHn0[1+3*nedges]=n[1+3*SHk0[m]]; SHn0[2+3*nedges]=n[2+3*SHk0[m]]; SHke[nedges]=nedges; nedges++; SHi0[nedges]= SHi[SHk1[m]]; SHiflg[nedges]=3; SHs0[nedges]=fabs( SHi[SHk1[m]] - SHi[SHk0[m]] )+.5; SHx0[nedges]=SHx[SHk1[m]]; SHy0[nedges]=SHy[SHk1[m]]; SHz0[nedges]=SHz[SHk1[m]]; SHn0[3*nedges]=n[3*SHk1[m]]; SHn0[1+3*nedges]=n[1+3*SHk1[m]]; SHn0[2+3*nedges]=n[2+3*SHk1[m]]; SHke[nedges]=nedges; nedges++; }else{ SHiflg[nedges]=0; if(SHj[SHk1[m]]==jpix)SHiflg[nedges]=1; if(SHj[SHk0[m]]==jpix)SHiflg[nedges]=2; if(SHo[m]+jpix*SHs[m]>=0.) SHi0[nedges]= SHo[m]+jpix* SHs[m]+.5; else SHi0[nedges]= SHo[m]+jpix* SHs[m]-.5; if(SHj[SHk1[m]]==SHj[SHk0[m]]) { fprintf(stderr," divide by zero in shpgnrm, SHj[%d]=SHj[%d]=%d, jpix=%d\n",SHk1[m],SHk0[m],SHj[SHk0[m]],jpix); abort(); } SHs0[nedges]=fabs((SHi[SHk1[m]]-SHi[SHk0[m]])/(float)(SHj[SHk1[m]]-SHj[SHk0[m]])); SHx0[nedges]=SHxo[m]+jpix*SHxs[m]; SHy0[nedges]=SHyo[m]+jpix*SHys[m]; SHz0[nedges]=SHzo[m]+jpix*SHzs[m]; SHn0[3*nedges]=SHno[3*m]+jpix*SHns[3*m]; SHn0[1+3*nedges]=SHno[1+3*m]+jpix*SHns[1+3*m]; SHn0[2+3*nedges]=SHno[2+3*m]+jpix*SHns[2+3*m]; SHke[nedges]=nedges; nedges++; } } } /* Sort SHke by increasing SHi0 (intersection points left to right.) */ for(mi=0;mi<nedges;mi++) { for(mj=mi+1;mj<nedges;mj++) { if(SHi0[SHke[mj]]<SHi0[SHke[mi]]) { kt=SHke[mi]; SHke[mi]=SHke[mj]; SHke[mj]=kt; } if(SHi0[SHke[mj]]==SHi0[SHke[mi]]&&SHiflg[SHke[mj]]<SHiflg[SHke[mi]]) { kt=SHke[mi]; SHke[mi]=SHke[mj]; SHke[mj]=kt; } } } /* Remove duplicates */ mi=0; while(mi<nedges-1) { if(SHi0[SHke[mi]]==SHi0[SHke[mi+1]]&&SHiflg[SHke[mi]]!=0&&SHiflg[SHke[mi+1]]!=0&&SHiflg[SHke[mi]]!=SHiflg[SHke[mi+1]]) { if(SHs0[SHke[mi+1]]>SHs0[SHke[mi]])SHs0[SHke[mi]]=SHs0[SHke[mi+1]]; SHs0[SHke[mi+1]]=SHs0[SHke[mi]]; for(mj=mi+1;mj<nedges;mj++)SHke[mj]=SHke[mj+1]; nedges--; }else{ mi++; } } /* Fill the Slabs */ for(m=0;m<nedges;m+=2) { de0=.5+SHs0[SHke[m]]; de1=.5+SHs0[SHke[m+1]]; for(ipix=SHi0[SHke[m]];ipix<=SHi0[SHke[m+1]];ipix++) { if(ipix>-1&&ipix<shIMax&&jpix>-1&&jpix<shJMax) { /* Compute perspective coordinates of the pixel being filled (xx,yy,zz) */ if(SHi0[SHke[m]]==SHi0[SHke[m+1]]) { xx=SHx0[SHke[m]]; yy=SHy0[SHke[m]]; zz=SHz0[SHke[m]]; nn[0]=SHn0[3*SHke[m]]; nn[1]=SHn0[1+3*SHke[m]]; nn[2]=SHn0[2+3*SHke[m]]; }else{ t=(float)(ipix-SHi0[SHke[m]])/(float)(SHi0[SHke[m+1]]-SHi0[SHke[m]]); xx=SHx0[SHke[m]]+t*(SHx0[SHke[m+1]]-SHx0[SHke[m]]); yy=SHy0[SHke[m]]+t*(SHy0[SHke[m+1]]-SHy0[SHke[m]]); zz=SHz0[SHke[m]]+t*(SHz0[SHke[m+1]]-SHz0[SHke[m]]); nn[0]=SHn0[3*SHke[m]]+t*(SHn0[3*SHke[m+1]]-SHn0[3*SHke[m]]); nn[1]=SHn0[1+3*SHke[m]]+t*(SHn0[1+3*SHke[m+1]]-SHn0[1+3*SHke[m]]); nn[2]=SHn0[2+3*SHke[m]]+t*(SHn0[2+3*SHke[m+1]]-SHn0[2+3*SHke[m]]); } an=sqrt(nn[0]*nn[0]+nn[1]*nn[1]+nn[2]*nn[2]); if(fabs(an)<1.e-7) { nn[0]=SHn0[3*SHke[m]]; nn[1]=SHn0[1+3*SHke[m]]; nn[2]=SHn0[2+3*SHke[m]]; an=sqrt(nn[0]*nn[0]+nn[1]*nn[1]+nn[2]*nn[2]); nn[0]=nn[0]*an; nn[1]=nn[1]*an; nn[2]=nn[2]*an; }else{ an=1./an; nn[0]=nn[0]*an; nn[1]=nn[1]*an; nn[2]=nn[2]*an; } shunpers(&xx,&yy,&zz,&xxx,&yyy,&zzz); /* Check against the Z-buffer */ shgetz(&ipix,&jpix,&zbuf); #ifdef VERBOSE printf("shgetz %d,%d %f %f\n",ipix,jpix,zbuf,zz);fflush(stdout); #endif if(zbuf>zz) { /* If pixel is visible, shade it and put it in the image */ shcolor(&xxx,&yyy,&zzz,nn,&(sh_rd[0]),&(sh_gd[0]),&(sh_bd[0]),&(sh_rd[2]),&(sh_gd[2]),&(sh_bd[2]),&r,&g,&b); #ifdef VERBOSE printf("shcolor, x-( %f %f %f ) n-( %f %f %f ) fc-(%d %d %d ) bc-(%d %d %d) c %d %d %d\n",xxx,yyy,zzz,nn[0],nn[1],nn[2],sh_rd[0],sh_gd[0],sh_bd[0],sh_rd[2],sh_gd[2],sh_bd[2],r,g,b); #endif if(shMask) { sh3dmask(&xxx,&yyy,&zzz,nn,&imask); if(imask==1)shsetp(&ipix,&jpix,&r,&g,&b,&zz); }else{ shsetp(&ipix,&jpix,&r,&g,&b,&zz); } } } } } } #ifdef VERBOSE printf(" done shdopg -- normal end\n"); #endif return; }

void shdoStrings (  FILE *  fid, float  x, float  y )

Writes the code in a postscript file for the strings that have been made with shstr. Parameters: fid  The File descriptor of the postscript file for output. x  The width of the screen. y  The height of the screen. Definition at line 69 of file shstr.c. { int i; int ix,iy; fprintf(fid,"initclip\n"); fprintf(fid,"/Helvetica findfont %f scalefont setfont\n",20.*shIMax/1024); fprintf(fid,"0. setgray\n"); for(i=0;i<SHnstr;i++) { ix=SHxstr[i]/x; iy=SHystr[i]/y; fflush(stdout); if((SHstr[i])[0]!='/') fprintf(fid,"%d %d moveto (%s) show\n",ix,iy,SHstr[i]); else fprintf(fid,"%s\n",SHstr[i]); } return; }

void shdraw (  float *  x, float *  y, float *  z, int *  ip )

A pen-up/pen-down style line drawing. Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. ip  If ip=3, move to the point with the pen up. If ip=2 move to the point with the pen down. Definition at line 29 of file shdraw.c. References shline(). { /* Calcomp type interface for line drawing routine. */ if(*ip==2) shline(&sh_x0,&sh_y0,&sh_z0,x,y,z); sh_x0=*x; sh_y0=*y; sh_z0=*z; return; } Here is the call graph for this function:

void shgetz (  int *  ipix, int *  jpix, float *  zz )

Reads the value in the z-buffer at a pixel. Parameters: ipix  The first coordinate of the point in pixels. jpix  The second coordinate of the point in pixels. zz  A place to store the z-buffer value.l Definition at line 28 of file shgetz.c. Referenced by shlnonrm(). { int index; /* Return the z-buffer. */ *zz=0.; if(*ipix<0||*jpix<0)return; if(*ipix>=shIMax||*jpix>=shJMax)return; index=(*ipix)+shIMax*(*jpix); *zz=shZBuffer[index]; return; }

void shinit (  int *  r, int *  g, int *  b )

Initializes the storage needed by sh and sets the current background color. Parameters: r  The red component of the background color. g  The green component of the background color. b  The blue component of the background color. Definition at line 28 of file shinit.c. References shlinc(), shlit(), shmask(), shnlit(), shnpln(), shpntc(), shsetp(), shsrfp(), shtric(), and shview(). { int i,j; int full=255; int zero=0; float am; float ad; float as; int ae; float dist,alpha,beta,xmin,xmax; float ftwo=2.; int one=1; int two=2; float xp,yp,zp; /* Initialization Routine. r,g,b Background Color, integer 0-&full. This routine fills the image buffer iax with the supplied color. Then sets the z-buffer to be infinitly far. (z is scaled from zero to 1., so it is set to 2.) Two light sources are defined. No clipping planes. Default color is white (&full,&full,&full) Default view is looking down the x-axis, distance=10. box=(0.,1.)**3 */ dist=10.; alpha=0.; beta=0.; xmin=0.; xmax=1.; shview(&dist,&alpha,&beta,&xmin,&xmax,&xmin,&xmax,&xmin,&xmax); shnlit(&two); xp=.3; yp=.3; zp=.5; shlit(&one,&xp,&yp,&zp,&full,&full,&full); xp=.7; yp=.3; zp=-1.; shlit(&two,&xp,&yp,&zp,&full,&full,&full); shnpln(&zero); shmask(&zero); shRedBackground=*r; shGreenBackground=*g; shBlueBackground=*b; shtric(&full,&full,&full); shlinc(&full,&full,&full); shpntc(&full,&full,&full); /* shsrfp(.3,.7,.5,5)*/ am=.3; ad=.9; as=.3; ae=5; shsrfp(&am,&ad,&as,&ae); shRedBuffer=(unsigned char*)malloc(shIMax*shJMax*sizeof(unsigned char)); if(shRedBuffer==(unsigned char*)NULL) { printf("Out of memory in shinit, allocating %d, line %d in file %s\n",shIMax*shJMax*sizeof(unsigned char),__LINE__,__FILE__);fflush(stdout); abort(); } shGreenBuffer=(unsigned char*)malloc(shIMax*shJMax*sizeof(unsigned char)); if(shGreenBuffer==(unsigned char*)NULL) { printf("Out of memory in shinit, allocating %d, line %d in file %s\n",shIMax*shJMax*sizeof(unsigned char),__LINE__,__FILE__);fflush(stdout); abort(); } shBlueBuffer=(unsigned char*)malloc(shIMax*shJMax*sizeof(unsigned char)); if(shBlueBuffer==(unsigned char*)NULL) { printf("Out of memory in shinit, allocating %d, line %d in file %s\n",shIMax*shJMax*sizeof(unsigned char),__LINE__,__FILE__);fflush(stdout); abort(); } shZBuffer=(float*)malloc(shIMax*shJMax*sizeof(float)); if(shZBuffer==(float*)NULL) { printf("Out of memory in shinit, allocating %d, line %d in file %s\n",shIMax*shJMax*sizeof(float),__LINE__,__FILE__);fflush(stdout); abort(); } /*printf("shinit, shZBuffer=0x%8.8x\n",shZBuffer);fflush(stdout);*/ for(j=0;j<shJMax;j++) for(i=0;i<shIMax;i++)shsetp(&i,&j,r,g,b,&ftwo); /* shlss('CMR10',5,1000,0,ierr) */ /* shcs(0,-3,3) */ sh_frame=0; return; } Here is the call graph for this function:

void shlinc (  int *  r, int *  g, int *  b )

Sets the current line color. Parameters: r  The red component of the line color. g  The green component of the line color. b  The blue component of the line color. Definition at line 36 of file shlinc.c. Referenced by sharrow(), and shinit(). { sh_rl=*r; sh_gl=*g; sh_bl=*b; }

void shline (  float *  x0, float *  y0, float *  z0, float *  x1, float *  y1, float *  z1 )

shades and renders a line. Parameters: x0  The first coordinate of the beginning of the line. y0  The second coordinate of the beginning of the line. z0  The third coordinate of the beginning of the line. x1  The first coordinate of the end of the line. y1  The second coordinate of the end of the line. z1  The third coordinate of the end of the line. Definition at line 31 of file shline.c. References shlnonrm(). Referenced by shcube(), shdraw(), shlineoff(), and shpl(). { /* Shade and render a line segment. */ float nn[3]={0.,0.,0.}; int zero=0; shlnonrm(x0,y0,z0,nn,x1,y1,z1,nn,&zero,&zero,&zero); return; } Here is the call graph for this function:

void shline2s (  int *  i0, int *  j0, int *  i1, int *  j1, int *  ixoff, int *  iyoff )

Draws a shadowed line segment in screen coordinates. Parameters: i0  The first coordinate of the beginning of the line in pixels. j0  The second coordinate of the beginning of the line in pixels. i1  The first coordinate of the end of the line in pixels. j1  The second coordinate of the end of the line in pixels. ixoff  The first coordinate of the offset of the shadow in pixels. iyoff  The second coordinate of the offest of the shadow in pixels. Definition at line 32 of file shline2s.c. References shsetp(). Referenced by sharrow(). { float zz; float xc,yc; float dx,dy; int ipix,jpix; /* draw a line segment */ if((*i0)==(*i1)&&(*j0)==(*j1)) { if((*i0)>0&&(*i0)<shIMax&&(*j0)>0&&(*j0)<shJMax) { zz=0.; shsetp(i0,j0,&sh_rl,&sh_gl,&sh_bl,&zz); } return; } xc=(*i0); yc=(*j0); dx=(*i1)-xc; dy=(*j1)-yc; zz=0.; while(fabs(dx)>.5||fabs(dy)>.5) { ipix=xc+.5+(*ixoff); jpix=yc+.5+(*iyoff); if(ipix>0&&ipix<shIMax&&jpix>0&&jpix<shJMax)shsetp(&ipix,&jpix,&sh_rl,&sh_gl,&sh_bl,&zz); dx=(*i1)-xc; dy=(*j1)-yc; if(abs(dx)>.5||abs(dy)>.5) { if(abs(dx)>abs(dy)) { if(dx>0) { xc=xc+1; yc=yc+dy/dx; }else{ xc=xc-1; yc=yc-dy/dx; } }else{ if(dy>0) { yc=yc+1; xc=xc+dx/dy; }else{ yc=yc-1; xc=xc-dx/dy; } } } dx=(*i1)-xc; dy=(*j1)-yc; } return; } Here is the call graph for this function:

void shlineoff (  float *  x0, float *  y0, float *  z0, float *  x1, float *  y1, float *  z1, int *  ixoff, int *  iyoff, int *  idoff )

Draws a line with an offset. Currently this routine just draws a line. Parameters: x0  The first coordinate of the beginning of the line. y0  The second coordinate of the beginning of the line. z0  The third coordinate of the beginning of the line. x1  The first coordinate of the end of the line. y1  The second coordinate of the end of the line. z1  The third coordinate of the end of the line. ixoff  The first coordinate of the offset in pixels. iyoff  The second coordinate of the offset in pixels. idoff  The depth coordinate of the offset in pixels. Definition at line 34 of file shlineoff.c. References shline(). { /* Line routine with offset */ shline(x0,y0,z0,x1,y1,z1); return; } Here is the call graph for this function:

void shlit (  int *  m, float *  x, float *  y, float *  z, int *  r, int *  g, int *  b )

Assigns values to a light source. Parameters: m  The number of the light source being changed. x  The first coordinate of the position of the light source. y  The second coordinate of the position of the light source. z  The third coordinate of the position of the light source. r  The red component of the color of the light. g  The green component of the color of the light. b  The blue component of the color of the light. Definition at line 44 of file shlit.c. References shpers(). Referenced by shinit(). { float xt=0.; float yt=0.; float zt=0.; /* This Routine Sets the parameters associated with light source m. */ /* x,y,z (real) position. */ /* r,g,b (integer 0-255) color. */ /* type (integer 0-1) 0=point,1=direction. */ if((*m)>=sh_mlit) { if(sh_mlit==0) { sh_mlit=(*m)+10; sh_rs=(int*)malloc(sh_mlit*sizeof(int)); sh_gs=(int*)malloc(sh_mlit*sizeof(int)); sh_bs=(int*)malloc(sh_mlit*sizeof(int)); sh_type=(int*)malloc(sh_mlit*sizeof(int)); sh_lit=(double*)malloc(3*sh_mlit*sizeof(double)); }else{ sh_mlit=(*m)+10; sh_rs=(int*)realloc((void*)sh_rs,sh_mlit*sizeof(int)); sh_gs=(int*)realloc((void*)sh_gs,sh_mlit*sizeof(int)); sh_bs=(int*)realloc((void*)sh_bs,sh_mlit*sizeof(int)); sh_type=(int*)realloc((void*)sh_type,sh_mlit*sizeof(int)); sh_lit=(double*)realloc((void*)sh_lit,sh_mlit*sizeof(double)); } } sh_rs[(*m)-1]=*r; sh_gs[(*m)-1]=*g; sh_bs[(*m)-1]=*b; shpers(x,y,z,&xt,&yt,&zt); sh_lit[0+3*((*m)-1)]=xt; sh_lit[1+3*((*m)-1)]=yt; sh_lit[2+3*((*m)-1)]=zt; sh_type[(*m)-1]=1; return; } Here is the call graph for this function:

void shlnonrm (  float *  x0, float *  y0, float *  z0, float *  n0, float *  x1, float *  y1, float *  z1, float *  n1, int *  ixoff, int *  iyoff, int *  idoff )

Draws a shaded line with an offset. The offset is used to determine visability, but not applied to the line when drawn. Parameters: x0  The first coordinate of the beginning of the line. y0  The second coordinate of the beginning of the line. z0  The third coordinate of the beginning of the line. n0  The normal to the surface at the beginning of the line. An array of length at least 3. x1  The first coordinate of the end of the line. y1  The second coordinate of the end of the line. z1  The third coordinate of the end of the line. n1  The normal to the surface at the end of the line. An array of length at least 3. ixoff  The horizontal coordinate (in pixels) of the offest. iyoff  The vertical coordinate (in pixels) of the offest. idoff  The depth coordinate (in pixels) of the offest. Definition at line 52 of file shlnonrm.c. References sh3dmask(), shcolor(), shgetz(), shpers(), shsetp(), and shunpers(). Referenced by shline(), and shpgnrm(). { int rr=0; int gg=0; int bb=0; int i0,j0,i1,j1; int imask=0; float nn[3]={0.,0.,0.}; float zbuf=0.; float dx,dy; int ipix,jpix; float s0=0.; float t0=0.; float d0; float s1=0.; float t1=0.; float d1; float xx,yy,zz; float xxx=0.; float yyy=0.; float zzz=0.; float xc,yc; int clipped; int ip; float direc; /* Project into perspective coordinates. */ shpers(x0,y0,z0,&s0,&t0,&d0); shpers(x1,y1,z1,&s1,&t1,&d1); if((*idoff)<0) { d0=1.005*d0; d1=1.005*d1; } if((*idoff)>0) { d0=.995*d0; d1=.995*d1; } /* Transform to pixel coordinates. */ i0=s0*shMax+.5; j0=t0*shMax+.5; i1=s1*shMax+.5; j1=t1*shMax+.5; if(i0==i1&&j0==j1) { if(i0>0&&i0<shIMax&&j0>0&&j0<shJMax) { shgetz(&i0,&j0,&zbuf); if(zbuf>d0) { if(shMask) { sh3dmask(x0,y0,z0,n0,&imask); if(imask==1) { shcolor(x0,y0,z0,n0,&sh_rl,&sh_gl,&sh_bl,&sh_rl,&sh_gl,&sh_bl,&rr,&gg,&bb); if(shlw==0)shsetp(&i0,&j0,&rr,&gg,&bb,&d0); else shtmpsetp(i0,j0,rr,gg,bb,d0,shlw); } }else{ shcolor(x0,y0,z0,n0,&sh_rl,&sh_gl,&sh_bl,&sh_rl,&sh_gl,&sh_bl,&rr,&gg,&bb); if(shlw==0)shsetp(&i0,&j0,&rr,&gg,&bb,&d0); else shtmpsetp(i0,j0,rr,gg,bb,d0,shlw); } } } return; } /* Trace and shade the line. */ xc=i0; yc=j0; dx=i1-xc; dy=j1-yc; while(fabs(dx)>.5||fabs(dy)>.5) { if(fabs(i1-i0)>fabs(j1-j0)) { xx=s0+(float)(xc-i0)*(s1-s0)/(float)(i1-i0); yy=t0+(float)(xc-i0)*(t1-t0)/(float)(i1-i0); zz=d0+(float)(xc-i0)*(d1-d0)/(float)(i1-i0); nn[0]=n0[0]+(float)(xc-i0)*(n1[0]-n0[0])/(float)(i1-i0); nn[1]=n0[1]+(float)(xc-i0)*(n1[1]-n0[1])/(float)(i1-i0); nn[2]=n0[2]+(float)(xc-i0)*(n1[2]-n0[2])/(float)(i1-i0); }else{ xx=s0+(float)(yc-j0)*(s1-s0)/(float)(j1-j0); yy=t0+(float)(yc-j0)*(t1-t0)/(float)(j1-j0); zz=d0+(float)(yc-j0)*(d1-d0)/(float)(j1-j0); nn[0]=n0[0]+(float)(yc-j0)*(n1[0]-n0[0])/(float)(j1-j0); nn[1]=n0[1]+(float)(yc-j0)*(n1[1]-n0[1])/(float)(j1-j0); nn[2]=n0[2]+(float)(yc-j0)*(n1[2]-n0[2])/(float)(j1-j0); } /* clip */ clipped=0; shunpers(&xx,&yy,&zz,&xxx,&yyy,&zzz); if(sh_nplns>0) { for(ip=0;ip<sh_nplns;ip++) { direc=sh_plnn[ 3*ip]*(xxx-sh_plno[ 3*ip])+sh_plnn[1+3*ip]*(yyy-sh_plno[1+3*ip])+sh_plnn[2+3*ip]*(zzz-sh_plno[2+3*ip]); if(sh_oper[ip]==0) clipped=clipped||(sh_ipln[ip]*direc<0); else clipped=clipped&&(sh_ipln[ip]*direc<0); } } if(!clipped) { ipix=xc+.5; jpix=yc+.5; if(ipix>0&&ipix<shIMax&&jpix>0&&jpix<shJMax) { shgetz(&ipix,&jpix,&zbuf); if(zbuf>zz-.005) { if(shMask) { sh3dmask(&xxx,&yyy,&zzz,nn,&imask); if(imask==1) { shcolor(&xxx,&yyy,&zzz,nn,&sh_rl,&sh_gl,&sh_bl,&sh_rl,&sh_gl,&sh_bl,&rr,&gg,&bb); if(shlw==0)shsetp(&ipix,&jpix,&rr,&gg,&bb,&zz); else shtmpsetp(ipix,jpix,rr,gg,bb,zz,shlw); } }else{ shcolor(&xxx,&yyy,&zzz,nn,&sh_rl,&sh_gl,&sh_bl,&sh_rl,&sh_gl,&sh_bl,&rr,&gg,&bb); if(shlw==0)shsetp(&ipix,&jpix,&rr,&gg,&bb,&zz); else shtmpsetp(ipix,jpix,rr,gg,bb,zz,shlw); } } } } dx=i1-xc; dy=j1-yc; if(fabs(dx)>.5||fabs(dy)>.5) { if(fabs(dx)>fabs(dy)) { if(dx>0) { xc=xc+1; yc=yc+dy/dx; }else{ xc=xc-1; yc=yc-dy/dx; } }else{ if(dy>0) { yc=yc+1; xc=xc+dx/dy; }else{ yc=yc-1; xc=xc-dx/dy; } } } } return; } Here is the call graph for this function:

void shlss (  char *  name, int *  lname, int *  pmag, int *  ifnt, int *  ierr, int  ln )

Loads a character set into the table of character sets. Parameters: name  The root of the name of the character set (.iax is appended). lname  The lenght of the character string containing the name (needed for fortran). pmag  The size of the character set (a la TeX, 1000 is nominal size). ifnt  The table entry for this font. ierr  TRUE if the font loaded. FALSE otherwise. ln  The length of the name. (not used). Definition at line 40 of file shlss.c. { /* loads a symbol set */ char *ssname; char *cname; char *file; int i,k; int len=0; int irem; int irec,nrec; FILE *fid; int xref,yref; char fontnm[20]; /* loads a symbol set */ *ierr=0; cname=(char*)malloc(((*lname)+1)*sizeof(char)); ssname=(char*)malloc(((*lname)+5)*sizeof(char)); strncpy(cname,name,(*lname)); cname[(*lname)]=0x0; sprintf(ssname,"%s.%4.4d",cname,*pmag); file=(char*)malloc((strlen(QUOTE( FONTPATH ))+(*lname)+10)*sizeof(char)); strcpy(file,QUOTE( FONTPATH ) ); strcat(file,ssname); strcat(file,"IAX"); free(cname); free(ssname); /* Is it loaded already? */ if(*ifnt>=sh_kFonts) { if(sh_kfont==(int*)NULL) { sh_kFonts=10; sh_kfont=(int*)malloc(sh_kFonts*sizeof(int)); }else{ sh_kFonts=*ifnt+10; sh_kfont=(int*)realloc((void*)sh_kfont,sh_kFonts*sizeof(int)); } } if(sh_nFonts>0) { k=0; for(i=0;i<sh_nFonts;i++) { if(!strcmp(sh_fontnm[i],ssname))k=i; } } if(k!=0) { sh_kfont[*ifnt]=k; *ierr=1; return; } if((fid=fopen(file,"r"))==(FILE*)NULL) { fprintf(stderr,"shSoft Font file -->%s<-- not found\n",file); free(file); return; } free(file); /* read in character set */ if(sh_nFonts>=sh_mFonts) { if(sh_mFonts==0) { sh_mFonts=10; sh_dir=(int*)malloc(sh_mFonts*sizeof(int)); sh_fmag=(int*)malloc(sh_mFonts*sizeof(int)); sh_fontnm=(char**)malloc(sh_mFonts*sizeof(char*)); }else{ sh_mFonts=*ifnt+10; sh_dir=(int*)realloc((void*)sh_dir,sh_mFonts*sizeof(int)); sh_fmag=(int*)realloc((void*)sh_fmag,sh_mFonts*sizeof(int)); sh_fontnm=(char**)realloc((void*)sh_fontnm,sh_mFonts*sizeof(char*)); } } sh_kfont[*ifnt]=sh_nFonts; fscanf(fid,"%d",&len); if(sh_nimg+len>=sh_lenimg) { if(sh_lenimg==0) { sh_lenimg=sh_nimg+len; sh_image=(char*)malloc(sh_lenimg*sizeof(char)); }else{ sh_lenimg+=sh_nimg+len; sh_image=(char*)realloc((void*)sh_image,sh_lenimg*sizeof(char)); } } nrec=floor(len/512); for(irec=0;irec<nrec;irec++); { for(i=0;i<512;i++)sh_image[sh_nimg+512*irec+i]=fgetc(fid); fgetc(fid); } irem=len%512; if(irem>0) for(i=0;i<irem;i++)sh_image[sh_nimg+512*nrec+i]=fgetc(fid); fclose(fid); sh_dir[sh_nFonts]=sh_nimg; sh_fmag[sh_nFonts]=*pmag; sh_fontnm[sh_nFonts]=(char*)malloc((strlen(ssname)+1)*sizeof(char)); strcpy(sh_fontnm[sh_nFonts],ssname); sh_nimg+=len; sh_nFonts++; return; }

void shmask (  int *  swtch  )

Turns the 3d mask on and off. Parameters: swtch  If TRUE (>0) the mask is used, o.w. it is not used. Definition at line 32 of file shmask.c. Referenced by shinit(). { shMask=*swtch; return; }

void shnlit (  int *  n  )

Sets the number of light sources to be used. The first n are used, the rest are not deleted, just not used. Parameters: n  The number of light sources to use. Definition at line 26 of file shnlit.c. Referenced by shinit(). { /* This Routine Sets the number of light sources to be used. */ sh_nlit=*n; return; }

void shnpln (  int *  n  )

Sets the number of clipping planes to use. The first n are used, the rest are not deleted, just no used. Parameters: n  The number of clipping planes to use. Definition at line 26 of file shnpln.c. Referenced by shinit(), and shpgnrm(). { /* Set the number of clipping planes. */ sh_nplns=*n; return; }

void shpec (  int *  r, int *  g, int *  b )

Sets the color of the front facing edges of polygons. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 57 of file shpgc.c. Referenced by shtric(). { sh_rd[1] =*r; sh_gd[1] =*g; sh_bd[1] =*b; sh_ram[1]=*r; sh_gam[1]=*g; sh_bam[1]=*b; return; }

void shpers (  float *  x, float *  y, float *  z, float *  s, float *  t, float *  depth )

Performs the perspective transformation. Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. s  (output) The first coordinate of the transformed point. t  (output) The second coordinate of the transformed point. depth  (output) The depth coordinate of the transformed point. Definition at line 32 of file shpers.c. Referenced by sharrow(), shlit(), shlnonrm(), shscal(), and shstr(). { float c1,c2,c3; float dtemp; /* Project x,y,z to perspective coordinates. */ c1= (*x-shv_x0)*shv_n11+(*y-shv_y0)*shv_n12+(*z-shv_z0)*shv_n13; c2= (*x-shv_x0)*shv_n21+(*y-shv_y0)*shv_n22+(*z-shv_z0)*shv_n23; c3= (*x-shv_x0)*shv_n31+(*y-shv_y0)*shv_n32+(*z-shv_z0)*shv_n33; dtemp=(shv_l-shv_d)/(c1-shv_d); *s=dtemp*c2*shv_pcale+shv_soff; *t=dtemp*c3*shv_pcale+shv_toff; *depth=-dtemp*c1*shv_dscl+shv_doff; return; }

void shpg (  int *  nv, float *  xt, float *  yt, float *  zt )

renders a polygon. Parameters: nv  The number of vretices. xt  An array giving the first coordinates of the vertices. yt  An array giving the first coordinates of the vertices. zt  An array giving the first coordinates of the vertices. Definition at line 30 of file shpg.c. Referenced by shtri(). { int r,g,b; int r0,g0,b0; int i; float an; /* shade and z-buffer a polygon defined by the vertices xt,yt,zt. */ float nx,ny,nz; float *nrm; nrm=(float*)malloc(3*(*nv)*sizeof(float)); /* find normal */ nx=(yt[1]-yt[0])*(zt[2]-zt[0])-(yt[2]-yt[0])*(zt[1]-zt[0]); ny=(xt[2]-xt[0])*(zt[1]-zt[0])-(xt[1]-xt[0])*(zt[2]-zt[0]); nz=(xt[1]-xt[0])*(yt[2]-yt[0])-(xt[2]-xt[0])*(yt[1]-yt[0]); an=sqrt(nx*nx+ny*ny+nz*nz); if(an==0.)return; an=1./an; nx=nx*an; ny=ny*an; nz=nz*an; for(i=0;i<*nv;i++) { nrm[ 3*i]=nx; nrm[1+3*i]=ny; nrm[2+3*i]=nz; } shpgnrm(nv,xt,yt,zt,nrm); free(nrm); return; }

void shpgc (  int *  r, int *  g, int *  b )

Sets the color of the front face (determined by the normal) of polygons. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 39 of file shpgc.c. Referenced by shtric(). { sh_rd[0] =*r; sh_gd[0] =*g; sh_bd[0] =*b; sh_ram[0]=*r; sh_gam[0]=*g; sh_bam[0]=*b; return; }

void shpgnrm (  int *  nv, float *  xt, float *  yt, float *  zt, float *  nrmt )

Shades and z-buffers a polygon with specified normals. Parameters: nv  The number of vertices. xt  The first coordinates of the vertices. yt  The second coordinates of the vertices. zt  The third coordinates of the vertices. nrmt  The normals at the vertices. The normal for vertex i is (nrm[0+3*i],nrm[1+3*i],nrm[2+3*i]) Definition at line 32 of file shpgnrm.c. References shclippg(), shdopg(), shlnonrm(), shnpln(), and shqnpln(). { int r=0; int g=0; int b=0; int r0=0; int g0=0; int b0=0; int zero=0; int one=1; int full=255; int i; /* shade and z-buffer a polygon defined by the vertices xt,yt,zt. */ float *x,*y,*z; float nx,ny,nz; int *ict,*ic; float *nrm; int nn=0; int np=0; x=(float*)NULL; y=(float*)NULL; z=(float*)NULL; nrm=(float*)NULL; ic=(int*)NULL; ict=(int*)malloc((*nv)*sizeof(int)); for(i=0;i<*nv;i++)ict[i]=1; shclippg(nv,xt,yt,zt,nrmt,ict,&nn,&x,&y,&z,&nrm,&ic); shdopg(&nn,x,y,z,nrm); #ifdef FGHJKL x=(float*)realloc((void*)x,(nn+1)*sizeof(float)); y=(float*)realloc((void*)y,(nn+1)*sizeof(float)); z=(float*)realloc((void*)z,(nn+1)*sizeof(float)); ic=(int*)realloc((void*)ic,(nn+1)*sizeof(int)); x[nn]=x[0]; y[nn]=y[0]; z[nn]=z[0]; ic[nn]=ic[0]; shqnpln(&np); shnpln(&zero); for(i=0;i<nn;i++) { if(0||ic[i+1]!=1) { /* shqlinc(&r,&g,&b); shlinc(&full,&full,&full);*/ if(fabs(x[i]-x[i+1])>1.e-6||fabs(y[i]-y[i+1])>1.e-6||fabs(z[i]-z[i+1])>1.e-6) { shlnonrm(&(x[i]),&(y[i]),&(z[i]),&(nrm[3*i]),&(x[i+1]),&(y[i+1]),&(z[i+1]),&(nrm[3*i+3]),&zero,&zero,&one); /* shlinc(&r,&g,&b);*/ } }else{ if(fabs(x[i]-x[i+1])>1.e-6||fabs(y[i]-y[i+1])>1.e-6||fabs(z[i]-z[i+1])>1.e-6) { shqlinc(&r0,&g0,&b0); shqpec(&r,&g,&b); if(r0!=r&&g0!=g&&b0!=b) { shlinc(&r,&g,&b); shlnonrm(&(x[i]),&(y[i]),&(z[i]),&(nrm[3*i]),&(x[i+1]),&(y[i+1]),&(z[i+1]),&(nrm[3*i+3]),&zero,&zero,&one); shlinc(&r0,&g0,&b0); } } } } shnpln(&np); #endif free(x); free(y); free(z); free(nrm); free(ict); free(ic); return; } Here is the call graph for this function:

void shpl (  int *  n, float *  x, float *  y, float *  z )

Shades and draws a polyline. Parameters: n  The number of points in the polyline. x  The first coordinates of the points on the polyline. y  The second coordinates of the points on the polyline. z  The third coordinates of the points on the polyline. Definition at line 29 of file shpl.c. References shline(). { int i; /* Shade and render a polyline */ for(i=0;i<*n-1;i++) shline(x+i,y+i,z+i,x+i+1,y+i+1,z+i+1); return; } Here is the call graph for this function:

void shpln (  int *  i, float *  ox, float *  oy, float *  oz, float *  nx, float *  ny, float *  nz, int *  iside )

Sets a clipping plane. Parameters: i  The number of the plane in the table of clipping planes. ox  The first coordinate of a point on the plane. oy  The second coordinate of a point on the plane. oz  The third coordinate of a point on the plane. nx  The first coordinate of the normal to the plane. ny  The second coordinate of the normal to the plane. nz  The third coordinate of the normal to the plane. iside  If positive the normal points out of the half space that is kept. If negative the opposite half space is kept. Definition at line 33 of file shpln.c. { int ip; ip=(*i)-1; /* Set the clipping plane info. */ if(*i>=sh_mplns) { if(sh_mplns==0) { sh_mplns=10; sh_plno=(float*)malloc(3*sh_mplns*sizeof(float)); sh_plnn=(float*)malloc(3*sh_mplns*sizeof(float)); sh_ipln=(int*)malloc(sh_mplns*sizeof(int)); sh_oper=(int*)malloc(sh_mplns*sizeof(int)); }else{ sh_mplns+=10; sh_plno=(float*)realloc((void*)sh_plno,3*sh_mplns*sizeof(float)); sh_plnn=(float*)realloc((void*)sh_plnn,3*sh_mplns*sizeof(float)); sh_ipln=(int*)realloc((void*)sh_ipln,sh_mplns*sizeof(int)); sh_oper=(int*)realloc((void*)sh_oper,sh_mplns*sizeof(int)); } } sh_plno[ 3*ip]=*ox; sh_plno[1+3*ip]=*oy; sh_plno[2+3*ip]=*oz; sh_plnn[ 3*ip]=*nx; sh_plnn[1+3*ip]=*ny; sh_plnn[2+3*ip]=*nz; sh_ipln[ip]=*iside; sh_oper[ip]=0; /*printf("plane %d is ((x,y,z)-(%f,%f,%f)).(%f,%f,%f)*%d>0\n",ip,sh_plno[ 3*ip],sh_plno[1+3*ip],sh_plno[2+3*ip],sh_plnn[ 3*ip],sh_plnn[1+3*ip],sh_plnn[2+3*ip],sh_ipln[ip]);fflush(stdout);*/ return; }

void shplnrm (  int *  nv, float *  xt, float *  yt, float *  zt )

Draws a polyline that is shaded according to the normals given by triplets of points (periodic). If any three consecutive points are colinear nothing is drawn. Parameters: nv  The number of points in the polyline xt  An array of length at leas tnv with the first coordinates of the points on the polyline. yt  An array of length at leas tnv with the first coordinates of the points on the polyline. zt  An array of length at leas tnv with the first coordinates of the points on the polyline. Definition at line 31 of file shplnrm.c. { float *n; float an; int m,m1,m2,m3; int i; int zero=0; /* Shade and render a polyline with normals */ n=(float*)malloc(3*(*nv)*sizeof(float)); for(m=0;m<*nv;m++) { m1=m; m3=m-1; if(m3<1)m3=*nv; m2=m+1; if(m2>*nv)m2=1; m1=1; m2=2; m3=3; n[3*m]=(yt[m2]-yt[m1])*(zt[m3]-zt[m1])-(yt[m3]-yt[m1])*(zt[m2]-zt[m1]); n[1+3*m]=(xt[m3]-xt[m1])*(zt[m2]-zt[m1])-(xt[m2]-xt[m1])*(zt[m3]-zt[m1]); n[2+3*m]=(xt[m2]-xt[m1])*(yt[m3]-yt[m1])-(xt[m3]-xt[m1])*(yt[m2]-yt[m1]); an=sqrt(n[3*m]*n[3*m]+n[1+3*m]*n[1+3*m]+n[2+3*m]*n[2+3*m]); if(an==0.)return; an=1./an; n[3*m]=n[3*m]*an; n[1+3*m]=n[1+3*m]*an; n[2+3*m]=n[2+3*m]*an; } for(i=0;i<*nv-1;i++) shlnonrm(xt+i,yt+i,zt+i,&(n[3*i]),xt+i+1,yt+i+1,zt+i+1,&(n[3*i+3]),&zero,&zero,&zero); return; }

void shpnt (  float *  x, float *  y, float *  z )

Draws a single point (one pixel). Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. Definition at line 28 of file shpnt.c. { /* Shade and Render a point. */ float s=0.; float t=0.; float d=0.; int clipped; int ip; float direc; int ipix,jpix; float zbuf=0.; /* clip */ clipped=0; if(sh_nplns>0) { for(ip==0;ip<sh_nplns;ip++) { direc=sh_plnn[ 3*ip]*((*x)-sh_plno[ 3*ip])+sh_plnn[1+3*ip]*((*y)-sh_plno[1+3*ip])+sh_plnn[2+3*ip]*((*z)-sh_plno[2+3*ip]); if(sh_oper[ip]==0) clipped=clipped||(sh_ipln[ip]*direc==0); else clipped=clipped&&(sh_ipln[ip]*direc==0); } if(clipped)return; } shpers(x,y,z,&s,&t,&d); ipix=s*shMax+.5; jpix=t*shMax+.5; if(ipix>0&&ipix<shIMax&&jpix>0&&jpix<shJMax) { shgetz(&ipix,&jpix,&zbuf); if(zbuf>d)shsetp(&ipix,&jpix,&sh_rp,&sh_gp,&sh_bp,&d); } return; }

void shpntc (  int *  r, int *  g, int *  b )

Sets the color used to draw points. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 36 of file shpntc.c. Referenced by shinit(). { sh_rp=*r; sh_gp=*g; sh_bp=*b; return; }

void shputps (  char *  file, int *  ln, int *  m, int *  n, unsigned char *  imageR, unsigned char *  imageG, unsigned char *  imageB, int  fln )

writes an image to a postscript file. Parameters: file  The root of the name of the postscript file. (.ps is added). ln  The length of the character string with the name. (Needed for FORTRAN). m  The width of the image in pixels. n  The height of the image in pixels. imageR  The red component of the image. Stored by horizontal scan lines, starting at the top. imageG  The green component of the image. imageB  The blue component of the image. fln  The length of the file name. (not used). Definition at line 34 of file shputps.c. Referenced by shsave(). { char title[255]=""; FILE *fid; int i,j; int i0,j0,i1,j1; int xll,yll,xur,yur; float w; i0=*m; j0=*n; i1=-1; j1=-1; for(j=0;j<*n;j++) { for(i=0;i<*m;i++) { if(shRedBuffer[i+j*shIMax]!=255||shGreenBuffer[i+j*shIMax]!=255||shBlueBuffer[i+j*shIMax]!=255) { if(j<=j0)j0=j; if(i<=i0)i0=i; if(j>=j1)j1=j; if(i>=i1)i1=i; } } } w=(*n)*504./(*m); xll=i0*w/(*m); yll=j0*w/(*n); xur=i1*w/(*m); yur=j1*w/(*n); strncpy(title,file,*ln); title[*ln]=0x0; fid=fopen(file,"w"); if(fid==(FILE*)NULL) { fprintf(stderr,"shputps: Problem opening file -->%s<-- for writing.\n",file); fflush(stderr); return; } fprintf(fid,"%s\n","%!PS-Adobe-3.0"); fprintf(fid,"%s\n","%%Creator: SH by Michael E. Henderson, 03/17/93"); fprintf(fid,"%s%s\n","%%Title: ",title); fprintf(fid,"%s %d %d %d %d\n","%%BoundingBox:",xll+50,yll+50,xur+52,yur+52); fprintf(fid,"%s\n","%%EndComments"); fprintf(fid,"%s\n","%%BeginProlog"); fprintf(fid,"%s\n","%%BeginResource: "); fprintf(fid,"%s\n"," "); fprintf(fid,"%s\n","% see if we have the ""colorimage"" operator."); fprintf(fid,"%s\n","% define one if we don""t"); fprintf(fid,"%s\n","/colorimage where % do we know about ""colorimage""?"); fprintf(fid,"%s\n"," { pop } % yes: pop off the ""dict"" returned"); fprintf(fid,"%s\n"," { % no: define one"); fprintf(fid,"%s\n"," /str1 1 string def"); fprintf(fid,"%s\n"," /str3 3 string def"); fprintf(fid,"%s\n"," /colorimage"); fprintf(fid,"%s\n"," { pop pop % pop off ""false"", ""3"" operands"); fprintf(fid,"%s\n"," pop % pop off old ""readhexstring"" proc"); fprintf(fid,"%s\n"," % and define a new one for ""image"""); fprintf(fid,"%s\n"," { currentfile str3 readhexstring pop pop"); fprintf(fid,"%s\n"," str1 0 % for the ""put"" below"); fprintf(fid,"%s\n"," str3 0 get 20 mul % Red"); fprintf(fid,"%s\n"," str3 1 get 32 mul % Green"); fprintf(fid,"%s\n"," str3 2 get 12 mul % Blue"); fprintf(fid,"%s\n"," add add 64 idiv % I = .5G + .31R + .18B"); fprintf(fid,"%s\n"," put str1 % str1 = intensity of R,G,B triplet"); fprintf(fid,"%s\n"," } image"); fprintf(fid,"%s\n"," } def % end of colorimage def"); fprintf(fid,"%s\n"," } ifelse % end of ""false"" case"); fprintf(fid,"%s\n"," "); fprintf(fid,"%s\n","%%EndResource"); fprintf(fid,"%s\n","%%EndProlog"); fprintf(fid,"%s\n"," "); fprintf(fid,"%s\n","%%Page: 1 1"); fprintf(fid,"%s\n"," "); fprintf(fid,"%s\n","%%BeginPageSetup"); fprintf(fid,"%s\n"," "); fprintf(fid," newpath\n"); fprintf(fid," %d %d moveto\n",xll+50,yll+50); fprintf(fid," %d %d lineto\n",xur+52,yll+50); fprintf(fid," %d %d lineto\n",xur+52,yur+52); fprintf(fid," %d %d lineto\n",xll+50,yur+52); fprintf(fid," closepath clip\n"); fprintf(fid,"%s\n","% Image size (1/72 inch coords,"); fprintf(fid,"%d %f %s\n",504,w," scale"); fprintf(fid,"%s\n"," "); fprintf(fid,"%s\n","%%EndPageSetup"); fprintf(fid,"%s\n"," "); fprintf(fid,"%f %f %s\n",.1,.1," translate"); fprintf(fid,"%s %d %s\n","/readstr ",3*(*m)," string def"); fprintf(fid,"%d %d %d %s %d %d %d %d %d %d %s\n",*m,*n,8,"[ ", *m, 0,0,*n,0,0, " ]"); fprintf(fid,"%s\n","{currentfile readstr readhexstring pop}"); fprintf(fid,"%s\n","false 3 colorimage"); for(j=0;j<*n;j++) { for(i=0;i<*m;i++) { fprintf(fid,"%2.2x%2.2x%2.2x",shRedBuffer[i+j*shIMax],shGreenBuffer[i+j*shIMax],shBlueBuffer[i+j*shIMax]); if(i%20==19)fprintf(fid,"\n"); } if((*m-1)%20!=19)fprintf(fid,"\n"); } fprintf(fid,"%s\n"," "); fprintf(fid,"%f %f %s\n",1./504.,1./w," scale"); shdoStrings(fid,shIMax/512.,shJMax/w); fprintf(fid,"%s\n"," showpage"); fclose(fid); return; }

void shputtiff (  char *  file, int  nrows, int  ncols, unsigned char *  r, unsigned char *  g, unsigned char *  b )

writes an image to a tiff file (if libtiff was installed when sh was configured). Parameters: file  The root of the name of the postscript file. (.ps is added). nrows  The height of the image in pixels. ncols  The width of the image in pixels. r  The red component of the image. Stored by horizontal scan lines, starting at the top. g  The green component of the image. b  The blue component of the image. Definition at line 40 of file shputtiff.c. Referenced by shsave(). { #ifdef HAVE_LIBTIFF TIFF *out = NULL; char *fn; char *buf; unsigned int i,j; unsigned red=0; unsigned green=1; unsigned blue=2; int rc; int ln; int separate; separate=0; ln=nrows*ncols; fn=(char *)malloc((strlen(file)+1)*sizeof(char)); strcpy(fn,file); /*printf(" puttiff for ==>%s<==\n",fn);*/ /*printf(" Image is %dx%d\n",nrows,ncols);*/ out = TIFFOpen(fn, "w"); free(fn); if (out == NULL)return; TIFFSetField(out, TIFFTAG_IMAGEWIDTH, (unsigned long) ncols); TIFFSetField(out, TIFFTAG_IMAGELENGTH, (unsigned long) nrows); TIFFSetField(out, TIFFTAG_PHOTOMETRIC, PHOTOMETRIC_RGB); TIFFSetField(out, TIFFTAG_SAMPLESPERPIXEL, 3); TIFFSetField(out, TIFFTAG_BITSPERSAMPLE, 8,8,8); TIFFSetField(out, TIFFTAG_COMPRESSION, COMPRESSION_NONE); TIFFSetField(out, TIFFTAG_ROWSPERSTRIP, 1); /*TIFFSetField(out, TIFFTAG_STRIPBYTECOUNTS, ncols);*/ TIFFSetField(out, TIFFTAG_ORIENTATION, ORIENTATION_TOPLEFT); if(separate) { TIFFSetField(out, TIFFTAG_PLANARCONFIG, PLANARCONFIG_SEPARATE); buf=(char *)malloc(ncols*sizeof(char)); }else{ TIFFSetField(out, TIFFTAG_PLANARCONFIG, PLANARCONFIG_CONTIG); buf=(char *)malloc(3*ncols*sizeof(char)); } TIFFSetField(out, TIFFTAG_SOFTWARE, "puttiff"); for(i=0;i<nrows;i++) { if(separate) { for(j=0;j<ncols;j++)buf[j]=r[j+i*ncols]; /* rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)i, red);*/ rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)(nrows-i-1), red); for(j=0;j<ncols;j++)buf[j]=g[j+i*ncols]; /* rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)i, green);*/ rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)(nrows-i-1), green); for(j=0;j<ncols;j++)buf[j]=b[j+i*ncols]; /* rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)i, blue);*/ rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)(nrows-i-1), blue); }else{ for(j=0;j<ncols;j++)buf[3*j ]=r[j+i*ncols]; for(j=0;j<ncols;j++)buf[3*j+1]=g[j+i*ncols]; for(j=0;j<ncols;j++)buf[3*j+2]=b[j+i*ncols]; /* rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)i, 0);*/ rc=TIFFWriteScanline(out, (unsigned char*)buf, (unsigned)(nrows-i-1), 0); } } free(buf); TIFFWriteDirectory(out); TIFFClose(out); #endif return; }

void shqbec (  int *  r, int *  g, int *  b )

Gets the color of the back facing edge of polygons. Parameters: r  (output) The red component of the color. g  (output) The green component of the color. b  (output) The blue component of the color. Definition at line 156 of file shpgc.c. { *r=sh_rd[3]; *g=sh_gd[3]; *b=sh_bd[3]; return; }

void shqbgc (  int *  r, int *  g, int *  b )

Gets the color of the back face of polygons. Parameters: r  (output) The red component of the color. g  (output) The green component of the color. b  (output) The blue component of the color. Definition at line 141 of file shpgc.c. { *r=sh_rd[2]; *g=sh_gd[2]; *b=sh_bd[2]; return; }

void shqeye (  float *  xe, float *  ye, float *  ze )

Returns the position of the eye. Parameters: xe  (output) The first coordinate of the eye. ye  (output) The second coordinate of the eye. ze  (output) The third coordinate of the eye. Definition at line 132 of file shpers.c. { /* Entry to return viewpoint. */ *xe=shv_d*shv_n11; *ye=shv_d*shv_n12; *ze=shv_d*shv_n13; return; }

void shqlinc (  int *  r, int *  g, int *  b )

Reads the current line color. Parameters: r  A place to sotre red component of the line color. g  A place to sotre green component of the line color. b  A place to sotre blue component of the line color. Definition at line 50 of file shlinc.c. Referenced by sharrow(). { *r=sh_rl; *g=sh_gl; *b=sh_bl; }

void shqlit (  int *  m, float *  x, float *  y, float *  z, int *  r, int *  g, int *  b )

Reads back a light source. Parameters: m  The number of the light source. x  (output) The first coordinate of the position of the light. y  (output) The second coordinate of the position of the light. z  (output) The third coordinate of the position of the light. r  (output) The red component of the color of the light. g  (output) The green component of the color of the light. b  (output) The blue component of the color of the light. Definition at line 32 of file shqlit.c. { /* This Routine Returns the parameters associated with light source */ /* x,y,z (real) position. */ /* r,g,b (integer 0-255) color. */ /* type (integer 0-1) 0=point,1=direction. */ *r=sh_rs[*m]; *g=sh_gs[*m]; *b=sh_bs[*m]; *x=sh_lit[ +3*(*m)]; *y=sh_lit[1+3*(*m)]; *z=sh_lit[2+3*(*m)]; return; }

void shqnpln (  int *  n  )

Reads back the number of clipping planes. Parameters: n  (output) The number of clipping planes. Definition at line 26 of file shqnpln.c. Referenced by shpgnrm(). { /* Return the number of clipping planes. */ *n=sh_nplns; return; }

void shqpec (  int *  r, int *  g, int *  b )

Gets the color of the front facing edge of polygons. Parameters: r  (output) The red component of the color. g  (output) The green component of the color. b  (output) The blue component of the color. Definition at line 126 of file shpgc.c. { *r=sh_rd[1]; *g=sh_gd[1]; *b=sh_bd[1]; return; }

void shqpgc (  int *  r, int *  g, int *  b )

Gets the color of the front face of polygons. Parameters: r  (output) The red component of the color. g  (output) The green component of the color. b  (output) The blue component of the color. Definition at line 111 of file shpgc.c. { *r=sh_rd[0]; *g=sh_gd[0]; *b=sh_bd[0]; return; }

void shqpntc (  int *  r, int *  g, int *  b )

Reads the color used to draw points. Parameters: r  (output) The red component of the color. g  (output) The green component of the color. b  (output) The blue component of the color. Definition at line 51 of file shpntc.c. { *r=sh_rp; *g=sh_gp; *b=sh_bp; return; }

void shsave (  char *  name, int *  ln, int  nln )

Saves the current frame buffer. If shOutputFormat is "postscript" a postscript file is written (duh). If shOutputFormat is "tiff" and libtiff was installed before sh was configured, a tiff file is written. Parameters: name  The root of the file name. ln  The length of the root of the file name (needed for FORTRAN). nln  The lenth of the root of the file name (not used). Definition at line 36 of file shsave.c. References shputps(), and shputtiff(). Referenced by shpause(). { /* Write image to disk.*/ if(!strcmp(shOutputFormat,"postscript")) shputps(name,ln,&shIMax,&shJMax,shRedBuffer,shGreenBuffer,shBlueBuffer,nln); #ifdef HAVE_LIBTIFF else if(!strcmp(shOutputFormat,"tiff")) shputtiff(name, shIMax,shJMax, shRedBuffer,shGreenBuffer,shBlueBuffer); #endif else shputps(name,ln,&shIMax,&shJMax,shRedBuffer,shGreenBuffer,shBlueBuffer,nln); return; } Here is the call graph for this function:

void shscal (  float *  xmin, float *  xmax, float *  ymin, float *  ymax, float *  zmin, float *  zmax )

Sets the visible volume. Parameters: xmin  The lower limit of the the first coordinate of the region. ymin  The lower limit of the the second coordinate of the region. zmin  The lower limit of the the third coordinate of the region. xmax  The upper limit of the the first coordinate of the region. ymax  The upper limit of the the second coordinate of the region. zmax  The upper limit of the the third coordinate of the region. Definition at line 31 of file shscal.c. References shpers(), and shsize(). Referenced by shview(). { int i,j; float one=1.; float zero=0.; float smax,smin; float tmax,tmin; float dmax,dmin; float s=0.; float t=0.; float d=0.; float x,y,z; float as,at,pcale; float soff; float toff; float dcale; float doff; float I[8*6]={0.,0.,0.,0.,1.,1.,1.,1., 1.,1.,1.,1.,0.,0.,0.,0., 0.,0.,1.,1.,0.,0.,1.,1., 1.,1.,0.,0.,1.,1.,0.,0., 0.,1.,0.,1.,0.,1.,0.,1., 1.,0.,1.,0.,1.,0.,1.,0.}; /* Perform Scaling for the perspective transform. */ shsize(&one,&zero,&zero,&one,&zero); smax=-1.e5; smin=1.e5; tmax=-1.e5; tmin=1.e5; dmax=-1.e5; dmin=1.e5; for(j=0;j<8;j++) { x=I[j+8*0]*(*xmax)+I[j+8*1]*(*xmin); y=I[j+8*2]*(*ymax)+I[j+8*3]*(*ymin); z=I[j+8*4]*(*zmax)+I[j+8*5]*(*zmin); shpers(&x,&y,&z,&s,&t,&d); if(s>smax)smax=s; if(s<smin)smin=s; if(t>tmax)tmax=t; if(t<tmin)tmin=t; if(d>dmax)dmax=d; if(d<dmin)dmin=d; } as=1.*shIMax/shMax; at=1.*shJMax/shMax; pcale=.9*as/(smax-smin); if(pcale*(tmax-tmin)>.9*at)pcale=.9*at/(tmax-tmin); soff=.5*(as-pcale*(smax+smin)); toff=.5*(at-pcale*(tmax+tmin)); dcale=.9/(dmax-dmin); doff=.5*(1.-dcale*(dmax+dmin)); shsize(&pcale,&soff,&toff,&dcale,&doff); return; } Here is the call graph for this function:

void shSetLineWidth (  int  w  )

Sets the width (in pixels) used when lines are drawn. Parameters: w  The line width in pixels Definition at line 29 of file shlnonrm.c. { shlw=w; return; }

void shSetOutputFilename (  char *  name  )

Sets the root of the output file name. Parameters: name  The root of the output file name (.ps or .tiff is appended). Definition at line 35 of file shname.c. { shOutputName=(char*)realloc(shOutputName,(strlen(name)+10)*sizeof(char)); strcpy(shOutputName,name); if(!strcmp(shOutputFormat,"tiff")) strcat(shOutputName,".tiff"); else strcat(shOutputName,".ps"); return; }

void shSetOutputFormat (  char *  type  )

Sets the format of the file written by shsave. The default is "postscript", but if libtiff is installed when sh is configured, this routine can be used to change the format to "tiff". Parameters: type  Either "postscript" or "tiff". Definition at line 35 of file shtype.c. { shOutputFormat=(char*)realloc(shOutputFormat,(strlen(type)+1)*sizeof(char)); strcpy(shOutputFormat,type); return; }

void shSetOutputResolution (  int  x, int  y )

Sets the size of the screen buffer, if called BEFORE shinit. Otherwise it does nothing until the next call of shinit. Parameters: x  The number of pixels in the horizontal direction. y  The number of pixels in the vertical direction. Definition at line 31 of file shres.c. { shIMax=x; shJMax=y; return; }

void shsetp (  int *  ipix, int *  jpix, int *  ired, int *  igrn, int *  iblu, float *  zz )

Sets a pixel in the frame buffer if zz is less than the value in the z-buffer. If the pixel is changed the the value in the z-buffer is set to zz. Parameters: ipix  The horizontal coordinate of the pixel. jpix  The vertical coordinate of the pixel. ired  The red commponent of the color. igrn  The green commponent of the color. iblu  The blue commponent of the color. zz  The z value (depth) of the pixel. Definition at line 33 of file shsetp.c. Referenced by shchar(), shclr(), shinit(), shline2s(), and shlnonrm(). { /* Set the pixel color and z-buffer. */ int index; /*printf("shsetp %d %d %d %d %d %f\n",*ipix,*jpix,*ired,*igrn,*iblu,*zz);fflush(stdout); printf("shsetp, shZBuffer=0x%8.8x\n",shZBuffer);fflush(stdout);*/ if(*ipix<0||*jpix<0)return; if(*ipix>=shIMax||*jpix>=shJMax)return; index=*ipix+shIMax*(*jpix); shRedBuffer[index]=*ired; shGreenBuffer[index]=*igrn; shBlueBuffer[index]=*iblu; shZBuffer[index]=*zz; return; }

void shsize (  float *  scl, float *  s0, float *  t0, float *  dsc, float *  d0 )

Directly sets the scales used. Parameters: scl  The scale for the screen coordinates (which will be in [0,1] before being transformed to pixels). s0  The origin for the horizontal screen coordinate. t0  The origin for the vertical screen coordinate. dsc  The scale for the depth coordinate (which will be in [0,1]) d0  The origin for the depth coordinate. Definition at line 152 of file shpers.c. Referenced by shscal(), and shview(). { /* Entry to set scale factors. */ shv_pcale=*scl; shv_soff=*s0; shv_toff=*t0; shv_doff=*d0; shv_dscl=*dsc; return; }

void shsphere (  float *  xo, float *  yo, float *  zo, float *  r )

Draws a tesselated sphere, 8*4**4 triangles, smoothed. Parameters: xo  The first coordinate of the center of the sphere. yo  The second coordinate of the center of the sphere. zo  The third coordinate of the center of the sphere. r  The radius of the sphere. Definition at line 31 of file shsphere.c. { float *xt,*yt,*zt; float x[3]={0.,0.,0.}; float y[3]={0.,0.,0.}; float z[3]={0.,0.,0.}; float nrm[9]={0.,0.,0.,0.,0.,0.,0.,0.,0.}; float xa,ya,za,xb,yb,zb,xc,yc,zc; float x0,y0,z0,x1,y1,z1,x2,y2,z2; float d; int maxtri; int i; int itri,ntri,mtri; int n,ismooth; int nlevel; /* Draws a tesselated sphere, 8*4**n triangles, if ismooth^=0 smoothed. */ ismooth=1; nlevel=4; maxtri=8; for(i=0;i<nlevel;i++)maxtri=maxtri*4; xt=(float*)malloc(3*maxtri*sizeof(float)); yt=(float*)malloc(3*maxtri*sizeof(float)); zt=(float*)malloc(3*maxtri*sizeof(float)); /* 1 */ ntri=0; xt[ 3*ntri]= 0.; yt[ 3*ntri]= 1.; zt[ 3*ntri]= 0.; xt[1+3*ntri]= 0.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 1.; xt[2+3*ntri]= 1.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 0.; ntri++; /* 2 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]= 1.; zt[ 3*ntri]= 0.; xt[1+3*ntri]=-1.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 0.; xt[2+3*ntri]= 0.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 1.; ntri++; /* 3 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]= 1.; zt[ 3*ntri]= 0.; xt[1+3*ntri]= 0.; yt[1+3*ntri]= 0.; zt[1+3*ntri]=-1.; xt[2+3*ntri]=-1.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 0.; ntri++; /* 4 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]= 1.; zt[ 3*ntri]= 0.; xt[1+3*ntri]= 1.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 0.; xt[2+3*ntri]= 0.; yt[2+3*ntri]= 0.; zt[2+3*ntri]=-1.; ntri++; /* 5 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]=-1.; zt[ 3*ntri]= 0.; xt[2+3*ntri]= 0.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 1.; xt[1+3*ntri]= 1.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 0.; ntri++; /* 6 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]=-1.; zt[ 3*ntri]= 0.; xt[2+3*ntri]=-1.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 0.; xt[1+3*ntri]= 0.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 1.; ntri++; /* 7 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]=-1.; zt[ 3*ntri]= 0.; xt[2+3*ntri]= 0.; yt[2+3*ntri]= 0.; zt[2+3*ntri]=-1.; xt[1+3*ntri]=-1.; yt[1+3*ntri]= 0.; zt[1+3*ntri]= 0.; ntri++; /* 8 */ xt[ 3*ntri]= 0.; yt[ 3*ntri]=-1.; zt[ 3*ntri]= 0.; xt[2+3*ntri]= 1.; yt[2+3*ntri]= 0.; zt[2+3*ntri]= 0.; xt[1+3*ntri]= 0.; yt[1+3*ntri]= 0.; zt[1+3*ntri]=-1.; ntri++; if(nlevel>0) { for(i=0;i<nlevel;i++) { mtri=ntri; for(itri=0;itri<ntri;itri++) { x0=xt[3*itri]; y0=yt[3*itri]; z0=zt[3*itri]; x1=xt[1+3*itri]; y1=yt[1+3*itri]; z1=zt[1+3*itri]; x2=xt[2+3*itri]; y2=yt[2+3*itri]; z2=zt[2+3*itri]; xa=.5*(x0+x1); ya=.5*(y0+y1); za=.5*(z0+z1); d=1./sqrt(xa*xa+ya*ya+za*za); xa=xa*d; ya=ya*d; za=za*d; xb=.5*(x1+x2); yb=.5*(y1+y2); zb=.5*(z1+z2); d=1./sqrt(xb*xb+yb*yb+zb*zb); xb=xb*d; yb=yb*d; zb=zb*d; xc=.5*(x2+x0); yc=.5*(y2+y0); zc=.5*(z2+z0); d=1./sqrt(xc*xc+yc*yc+zc*zc); xc=xc*d; yc=yc*d; zc=zc*d; xt[ 3*itri]=x0; yt[ 3*itri]=y0; zt[ 3*itri]=z0; xt[1+3*itri]=xa; yt[1+3*itri]=ya; zt[1+3*itri]=za; xt[2+3*itri]=xc; yt[2+3*itri]=yc; zt[2+3*itri]=zc; xt[ 3*mtri]=xa; yt[ 3*mtri]=ya; zt[ 3*mtri]=za; xt[1+3*mtri]=x1; yt[1+3*mtri]=y1; zt[1+3*mtri]=z1; xt[2+3*mtri]=xb; yt[2+3*mtri]=yb; zt[2+3*mtri]=zb; mtri++; xt[ 3*mtri]=xb; yt[ 3*mtri]=yb; zt[ 3*mtri]=zb; xt[1+3*mtri]=x2; yt[1+3*mtri]=y2; zt[1+3*mtri]=z2; xt[2+3*mtri]=xc; yt[2+3*mtri]=yc; zt[2+3*mtri]=zc; mtri++; xt[ 3*mtri]=xa; yt[ 3*mtri]=ya; zt[ 3*mtri]=za; xt[1+3*mtri]=xb; yt[1+3*mtri]=yb; zt[1+3*mtri]=zb; xt[2+3*mtri]=xc; yt[2+3*mtri]=yc; zt[2+3*mtri]=zc; mtri++; } ntri=mtri; } } for(itri=0;itri<ntri;itri++) { for(i=0;i<3;i++) { nrm[0+3*i]= xt[i+3*itri]; nrm[1+3*i]= yt[i+3*itri]; nrm[2+3*i]= zt[i+3*itri]; x[i]=(*xo)+(*r)*xt[i+3*itri]; y[i]=(*yo)+(*r)*yt[i+3*itri]; z[i]=(*zo)+(*r)*zt[i+3*itri]; } n=3; if(ismooth!=0)shpgnrm(&n,x,y,z,nrm); else if(ismooth<0)shpl(&n,x,y,z); else shpg(&n,x,y,z); } free(xt); free(yt); free(zt); return; }

void shsrfp (  float *  a, float *  d, float *  s, int *  n )

Set the surface properties for shading. Parameters: a  The ambient intensity. d  The diffusive intensity. s  The specular intensity. n  The specular exponent (larger means a smaller highlight) Definition at line 29 of file shsrfp.c. Referenced by shinit(). { /* This Routine Sets the surface properties used for the Shading */ /* a ambient intensity. */ /* d difuse intensity. */ /* s specular intensity. */ /* n specular exponent. */ sh_am=*a; sh_ad=*d; sh_as=*s; sh_nd=*n; return; }

void shstr (  float  x, float  y, float  z, char *  str )

Draws a horizontal character string at a point. Parameters: x  The first coordinate of the point. y  The second coordinate of the point. z  The third coordinate of the point. str  The character string. Definition at line 39 of file shstr.c. References shpers(). { float xx=0.; float yy=0.; float zz=0.; SHstr=(char**)realloc(SHstr,(SHnstr+1)*sizeof(char*)); SHxstr=(int*)realloc(SHxstr,(SHnstr+1)*sizeof(int)); SHystr=(int*)realloc(SHystr,(SHnstr+1)*sizeof(int)); SHzstr=(double*)realloc(SHzstr,(SHnstr+1)*sizeof(double)); shpers(&x,&y,&z,&xx,&yy,&zz); SHxstr[SHnstr]=xx*shMax+.5; SHystr[SHnstr]=yy*shMax+.5; SHzstr[SHnstr]=zz; SHstr[SHnstr]=(char*)malloc((strlen(str)+1)*sizeof(char)); strcpy(SHstr[SHnstr],str); SHnstr++; return; } Here is the call graph for this function:

void shstrs (  int  xx, int  yy, char *  str )

Adds a horizontal character string in screen coordinates. Parameters: xx  The first coordinate of the point in pixels. yy  The second coordinate of the point in pixels. str  The character string. Definition at line 34 of file shstrs.c. { SHstr=(char**)realloc(SHstr,(SHnstr+1)*sizeof(char*)); SHxstr=(int*)realloc(SHxstr,(SHnstr+1)*sizeof(int)); SHystr=(int*)realloc(SHystr,(SHnstr+1)*sizeof(int)); SHzstr=(double*)realloc(SHzstr,(SHnstr+1)*sizeof(double)); SHxstr[SHnstr]=xx; SHystr[SHnstr]=yy; SHzstr[SHnstr]=-1.; SHstr[SHnstr]=(char*)malloc((strlen(str)+1)*sizeof(char)); strcpy(SHstr[SHnstr],str); SHnstr++; return; }

void shtri (  float *  xt1, float *  yt1, float *  zt1, float *  xt2, float *  yt2, float *  zt2, float *  xt3, float *  yt3, float *  zt3 )

A convenience routine to draw a triangle. Parameters: xt1  The first coordinate of the first corner. yt1  The second coordinate of the first corner. zt1  The third coordinate of the first corner. xt2  The first coordinate of the second corner. yt2  The second coordinate of the second corner. zt2  The third coordinate of the second corner. xt3  The first coordinate of the third corner. yt3  The second coordinate of the third corner. zt3  The third coordinate of the third corner. Definition at line 34 of file shtri.c. References shpg(). { /* Shade a triangle. Note this just calls shpg.*/ float x[3]={0.,0.,0.}; float y[3]={0.,0.,0.}; float z[3]={0.,0.,0.}; int n; x[0]=*xt1; x[1]=*xt2; x[2]=*xt3; y[0]=*yt1; y[1]=*yt2; y[2]=*yt3; z[0]=*zt1; z[1]=*zt2; z[2]=*zt3; n=3; shpg(&n,x,y,z); return; } Here is the call graph for this function:

void shtric (  int *  r, int *  g, int *  b )

A convenience routine which sets the edge and face colors for polygons to the same color. Parameters: r  The red component of the color. g  The green component of the color. b  The blue component of the color. Definition at line 28 of file shtric.c. References shbec(), shbgc(), shpec(), and shpgc(). Referenced by shinit(). { shpgc(r,g,b); shpec(r,g,b); shbgc(r,g,b); shbec(r,g,b); return; } Here is the call graph for this function:

void shunpers (  float *  ss, float *  tt, float *  dd, float *  xx, float *  yy, float *  zz )

Performs the inverse of the perspective transformation. Parameters: ss  The first coordinate of the point. tt  The second coordinate of the point. dd  The depth coordinate of the point. xx  The first coordinate of the transformed point. yy  The second coordinate of the transformed point. zz  The third coordinate of the transformed point. Definition at line 60 of file shpers.c. Referenced by shlnonrm(). { float c1,c2,c3; /* Entry to invert perspective transform. */ c1=shv_d*(*dd-shv_doff)/(*dd-shv_doff+shv_dscl*(shv_l-shv_d)); c2=(*ss-shv_soff)*(c1-shv_d)/(shv_l-shv_d)/shv_pcale; c3=(*tt-shv_toff)*(c1-shv_d)/(shv_l-shv_d)/shv_pcale; *xx=shv_x0+c1*shv_n11+c2*shv_n21+c3*shv_n31; *yy=shv_y0+c1*shv_n12+c2*shv_n22+c3*shv_n32; *zz=shv_z0+c1*shv_n13+c2*shv_n23+c3*shv_n33; return; }

void shview (  float *  dist, float *  alpha, float *  beta, float *  xmin, float *  xmax, float *  ymin, float *  ymax, float *  zmin, float *  zmax )

A convenience routine that calls shsize, shvw and shscal. Parameters: dist  The distance of the eye from the center of the visible cube. alpha  The angle around the equator (in degrees) of the line from the cube center to the eye. alpha=0 is the y=0 plane. beta  The angle up from the equator (in degrees) of the line from the cube center to the eye. beta=90 is the north pole. xmin  The lower limit of the first coordinate of the visible cube. ymin  The lower limit of the second coordinate of the visible cube. zmin  The lower limit of the third coordinate of the visible cube. xmax  The upper limit of the first coordinate of the visible cube. ymax  The upper limit of the second coordinate of the visible cube. zmax  The upper limit of the third coordinate of the visible cube. Definition at line 34 of file shview.c. References shscal(), shsize(), and shvw(). Referenced by shinit(). { float zero=0.; float one=0.; /* Set view point and scale. */ shsize(&one,&zero,&zero,&one,&zero); shvw(dist,alpha,beta,xmin,xmax,ymin,ymax,zmin,zmax); shscal(xmin,xmax,ymin,ymax,zmin,zmax); return; } Here is the call graph for this function:

void shvw (  float *  dist, float *  alpha, float *  beta, float *  xmin, float *  xmax, float *  ymin, float *  ymax, float *  zmin, float *  zmax )

Sets the view point for the perspective transformation. Parameters: dist  The distance of the view point from the center of the cube of limits. alpha  The angle in degrees around the equator of the line from the view point to the center of the cube. alpha=0 is the y=0 plane. beta  The angle in degrees up from the equator of the line from the view point to the center of the cube. beta=0 is the z=0. plane, beta=90 is looking down on the north pole. xmin  The lower limit of the first coordinate of points in the cube. ymin  The lower limit of the second coordinate of points in the cube. zmin  The lower limit of the third coordinate of points in the cube. xmax  The upper limit of the first coordinate of points in the cube. ymax  The upper limit of the second coordinate of points in the cube. zmax  The upper limit of the third coordinate of points in the cube. Definition at line 90 of file shpers.c. Referenced by shview(). { float a,b,pi; /* Entry to set view point. */ shv_x0=((*xmin)+(*xmax))/2.; shv_y0=((*ymin)+(*ymax))/2.; shv_z0=((*zmin)+(*zmax))/2.; pi=3.1415926; a=(*alpha-90.)*pi/180.; b=*beta*pi/180.; shv_n11=cos(b)*cos(a); shv_n12=cos(b)*sin(a); shv_n13=sin(b); shv_eye[0]=*dist*shv_n11; shv_eye[1]=*dist*shv_n12; shv_eye[2]=*dist*shv_n13; shv_n21=-sin(a); shv_n22=cos(a); shv_n23=0.0; shv_n31=-sin(b)*cos(a); shv_n32=-sin(b)*sin(a); shv_n33=cos(b); shv_d=*dist; shv_l=0.75*shv_d; return; }

---