--- ray/src/common/xf.c 1989/10/05 07:54:31 1.4 +++ ray/src/common/xf.c 1989/10/18 23:15:37 1.7 @@ -12,8 +12,9 @@ static char SCCSid[] = "$SunId$ LBL"; #define PI 3.14159265358979323846 +#define d2r(a) ((PI/180.)*(a)) -#define checkarg(a,n) if (strcmp(av[i],a) || i+n >= ac) goto done +#define checkarg(a,n) if (av[i][a] || i+n >= ac) goto done int @@ -25,14 +26,12 @@ char *av[]; { double atof(), sin(), cos(); double xfmat[4][4], m4[4][4]; - double xfsca, theta; - int i, j; - int rept; + double xfsca, dtmp; + int i, icnt; setident4(retmat); *retsca = 1.0; - rept = 1; setident4(xfmat); xfsca = 1.0; @@ -43,7 +42,7 @@ char *av[]; switch (av[i][1]) { case 't': /* translate */ - checkarg("-t",3); + checkarg(2,3); m4[3][0] = atof(av[++i]); m4[3][1] = atof(av[++i]); m4[3][2] = atof(av[++i]); @@ -52,22 +51,22 @@ char *av[]; case 'r': /* rotate */ switch (av[i][2]) { case 'x': - checkarg("-rx",1); - theta = PI/180.0 * atof(av[++i]); - m4[1][1] = m4[2][2] = cos(theta); - m4[2][1] = -(m4[1][2] = sin(theta)); + checkarg(3,1); + dtmp = d2r(atof(av[++i])); + m4[1][1] = m4[2][2] = cos(dtmp); + m4[2][1] = -(m4[1][2] = sin(dtmp)); break; case 'y': - checkarg("-ry",1); - theta = PI/180.0 * atof(av[++i]); - m4[0][0] = m4[2][2] = cos(theta); - m4[0][2] = -(m4[2][0] = sin(theta)); + checkarg(3,1); + dtmp = d2r(atof(av[++i])); + m4[0][0] = m4[2][2] = cos(dtmp); + m4[0][2] = -(m4[2][0] = sin(dtmp)); break; case 'z': - checkarg("-rz",1); - theta = PI/180.0 * atof(av[++i]); - m4[0][0] = m4[1][1] = cos(theta); - m4[1][0] = -(m4[0][1] = sin(theta)); + checkarg(3,1); + dtmp = d2r(atof(av[++i])); + m4[0][0] = m4[1][1] = cos(dtmp); + m4[1][0] = -(m4[0][1] = sin(dtmp)); break; default: return(i); @@ -75,27 +74,30 @@ char *av[]; break; case 's': /* scale */ - checkarg("-s",1); + checkarg(2,1); + dtmp = atof(av[i+1]); + if (dtmp == 0.0) goto done; + i++; xfsca *= m4[0][0] = m4[1][1] = - m4[2][2] = atof(av[++i]); + m4[2][2] = dtmp; break; case 'm': /* mirror */ switch (av[i][2]) { case 'x': - checkarg("-mx",0); + checkarg(3,0); xfsca *= m4[0][0] = -1.0; break; case 'y': - checkarg("-my",0); + checkarg(3,0); xfsca *= m4[1][1] = -1.0; break; case 'z': - checkarg("-mz",0); + checkarg(3,0); xfsca *= m4[2][2] = -1.0; break; @@ -105,12 +107,12 @@ char *av[]; break; case 'i': /* iterate */ - checkarg("-i",1); - for (j = 0; j < rept; j++) { + checkarg(2,1); + icnt = atoi(av[++i]); + while (icnt-- > 0) { multmat4(retmat, retmat, xfmat); *retsca *= xfsca; } - rept = atoi(av[++i]); setident4(xfmat); xfsca = 1.0; break; @@ -122,10 +124,8 @@ char *av[]; multmat4(xfmat, xfmat, m4); } done: - for (j = 0; j < rept; j++) { - multmat4(retmat, retmat, xfmat); - *retsca *= xfsca; - } + multmat4(retmat, retmat, xfmat); + *retsca *= xfsca; return(i); } @@ -140,14 +140,12 @@ char *av[]; { double atof(), sin(), cos(); double xfmat[4][4], m4[4][4]; - double xfsca, theta; - int i, j; - int rept; + double xfsca, dtmp; + int i, icnt; setident4(retmat); *retsca = 1.0; - rept = 1; setident4(xfmat); xfsca = 1.0; @@ -158,7 +156,7 @@ char *av[]; switch (av[i][1]) { case 't': /* translate */ - checkarg("-t",3); + checkarg(2,3); m4[3][0] = -atof(av[++i]); m4[3][1] = -atof(av[++i]); m4[3][2] = -atof(av[++i]); @@ -167,22 +165,22 @@ char *av[]; case 'r': /* rotate */ switch (av[i][2]) { case 'x': - checkarg("-rx",1); - theta = -PI/180.0 * atof(av[++i]); - m4[1][1] = m4[2][2] = cos(theta); - m4[2][1] = -(m4[1][2] = sin(theta)); + checkarg(3,1); + dtmp = -d2r(atof(av[++i])); + m4[1][1] = m4[2][2] = cos(dtmp); + m4[2][1] = -(m4[1][2] = sin(dtmp)); break; case 'y': - checkarg("-ry",1); - theta = -PI/180.0 * atof(av[++i]); - m4[0][0] = m4[2][2] = cos(theta); - m4[0][2] = -(m4[2][0] = sin(theta)); + checkarg(3,1); + dtmp = -d2r(atof(av[++i])); + m4[0][0] = m4[2][2] = cos(dtmp); + m4[0][2] = -(m4[2][0] = sin(dtmp)); break; case 'z': - checkarg("-rz",1); - theta = -PI/180.0 * atof(av[++i]); - m4[0][0] = m4[1][1] = cos(theta); - m4[1][0] = -(m4[0][1] = sin(theta)); + checkarg(3,1); + dtmp = -d2r(atof(av[++i])); + m4[0][0] = m4[1][1] = cos(dtmp); + m4[1][0] = -(m4[0][1] = sin(dtmp)); break; default: return(i); @@ -190,28 +188,31 @@ char *av[]; break; case 's': /* scale */ - checkarg("-s",1); - *xfsca *= + checkarg(2,1); + dtmp = atof(av[i+1]); + if (dtmp == 0.0) goto done; + i++; + xfsca *= m4[0][0] = m4[1][1] = - m4[2][2] = 1.0 / atof(av[++i]); + m4[2][2] = 1.0 / dtmp; break; case 'm': /* mirror */ switch (av[i][2]) { case 'x': - checkarg("-mx",0); - *xfsca *= + checkarg(3,0); + xfsca *= m4[0][0] = -1.0; break; case 'y': - checkarg("-my",0); - *xfsca *= + checkarg(3,0); + xfsca *= m4[1][1] = -1.0; break; case 'z': - checkarg("-mz",0); - *xfsca *= + checkarg(3,0); + xfsca *= m4[2][2] = -1.0; break; default: @@ -220,12 +221,12 @@ char *av[]; break; case 'i': /* iterate */ - checkarg("-i",1); - for (j = 0; j < rept; j++) { + checkarg(2,1); + icnt = atoi(av[++i]); + while (icnt-- > 0) { multmat4(retmat, xfmat, retmat); *retsca *= xfsca; } - rept = atoi(av[++i]); setident4(xfmat); xfsca = 1.0; break; @@ -237,10 +238,8 @@ char *av[]; multmat4(xfmat, m4, xfmat); /* left multiply */ } done: - for (j = 0; j < rept; j++) { - multmat4(retmat, xfmat, retmat); - *retsca *= xfsca; - } + multmat4(retmat, xfmat, retmat); + *retsca *= xfsca; return(i); } #endif