--- ray/src/rt/srcsupp.c 1991/10/21 12:58:16 1.14 +++ ray/src/rt/srcsupp.c 2014/06/22 18:05:28 2.20 @@ -1,13 +1,14 @@ -/* Copyright (c) 1991 Regents of the University of California */ - #ifndef lint -static char SCCSid[] = "$SunId$ LBL"; +static const char RCSid[] = "$Id: srcsupp.c,v 2.20 2014/06/22 18:05:28 greg Exp $"; #endif - /* * Support routines for source objects and materials + * + * External symbols declared in source.h */ +#include "copyright.h" + #include "ray.h" #include "otypes.h" @@ -18,7 +19,7 @@ static char SCCSid[] = "$SunId$ LBL"; #include "face.h" -#define SRCINC 4 /* realloc increment for array */ +#define SRCINC 32 /* realloc increment for array */ SRCREC *source = NULL; /* our list of sources */ int nsources = 0; /* the number of sources */ @@ -26,13 +27,10 @@ int nsources = 0; /* the number of sources */ SRCFUNC sfun[NUMOTYPE]; /* source dispatch table */ -initstypes() /* initialize source dispatch table */ +void +initstypes(void) /* initialize source dispatch table */ { extern VSMATERIAL mirror_vs, direct1_vs, direct2_vs; - extern int fsetsrc(), ssetsrc(), sphsetsrc(), cylsetsrc(), rsetsrc(); - extern int nopart(), flatpart(), cylpart(); - extern double fgetplaneq(), rgetplaneq(); - extern double fgetmaxdisk(), rgetmaxdisk(); static SOBJECT fsobj = {fsetsrc, flatpart, fgetplaneq, fgetmaxdisk}; static SOBJECT ssobj = {ssetsrc, nopart}; static SOBJECT sphsobj = {sphsetsrc, nopart}; @@ -51,27 +49,32 @@ initstypes() /* initialize source dispatch table */ int -newsource() /* allocate new source in our array */ +newsource(void) /* allocate new source in our array */ { if (nsources == 0) source = (SRCREC *)malloc(SRCINC*sizeof(SRCREC)); else if (nsources%SRCINC == 0) - source = (SRCREC *)realloc((char *)source, + source = (SRCREC *)realloc((void *)source, (unsigned)(nsources+SRCINC)*sizeof(SRCREC)); if (source == NULL) return(-1); source[nsources].sflags = 0; source[nsources].nhits = 1; - source[nsources].ntests = 2; /* initial hit probability = 1/2 */ + source[nsources].ntests = 2; /* initial hit probability = 50% */ +#if SHADCACHE + source[nsources].obscache = NULL; +#endif return(nsources++); } -setflatss(src) /* set sampling for a flat source */ -register SRCREC *src; +void +setflatss( /* set sampling for a flat source */ + SRCREC *src +) { double mult; - register int i; + int i; src->ss[SV][0] = src->ss[SV][1] = src->ss[SV][2] = 0.0; for (i = 0; i < 3; i++) @@ -86,18 +89,22 @@ register SRCREC *src; } -fsetsrc(src, so) /* set a face as a source */ -register SRCREC *src; -OBJREC *so; +void +fsetsrc( /* set a face as a source */ + SRCREC *src, + OBJREC *so +) { - register FACE *f; - register int i, j; + FACE *f; + int i, j; double d; src->sa.success = 2*AIMREQT-1; /* bitch on second failure */ src->so = so; /* get the face */ f = getface(so); + if (f->area == 0.0) + objerror(so, USER, "zero source area"); /* find the center */ for (j = 0; j < 3; j++) { src->sloc[j] = 0.0; @@ -106,7 +113,7 @@ OBJREC *so; src->sloc[j] /= (double)f->nv; } if (!inface(src->sloc, f)) - objerror(so, USER, "cannot hit center"); + objerror(so, USER, "cannot hit source center"); src->sflags |= SFLAT; VCOPY(src->snorm, f->norm); src->ss2 = f->area; @@ -119,8 +126,7 @@ OBJREC *so; } src->srad = sqrt(src->srad); /* compute size vectors */ - if (f->nv == 4 || (f->nv == 5 && /* parallelogram case */ - dist2(VERTEX(f,0),VERTEX(f,4)) <= FTINY*FTINY)) + if (f->nv == 4) /* parallelogram case */ for (j = 0; j < 3; j++) { src->ss[SU][j] = .5*(VERTEX(f,1)[j]-VERTEX(f,0)[j]); src->ss[SV][j] = .5*(VERTEX(f,3)[j]-VERTEX(f,0)[j]); @@ -130,9 +136,11 @@ OBJREC *so; } -ssetsrc(src, so) /* set a source as a source */ -register SRCREC *src; -register OBJREC *so; +void +ssetsrc( /* set a source as a source */ + SRCREC *src, + OBJREC *so +) { double theta; @@ -140,7 +148,7 @@ register OBJREC *so; src->so = so; if (so->oargs.nfargs != 4) objerror(so, USER, "bad arguments"); - src->sflags |= SDISTANT; + src->sflags |= (SDISTANT|SCIR); VCOPY(src->sloc, so->oargs.farg); if (normalize(src->sloc) == 0.0) objerror(so, USER, "zero direction"); @@ -156,42 +164,49 @@ register OBJREC *so; } -sphsetsrc(src, so) /* set a sphere as a source */ -register SRCREC *src; -register OBJREC *so; +void +sphsetsrc( /* set a sphere as a source */ + SRCREC *src, + OBJREC *so +) { - register int i; + int i; src->sa.success = 2*AIMREQT-1; /* bitch on second failure */ src->so = so; if (so->oargs.nfargs != 4) objerror(so, USER, "bad # arguments"); if (so->oargs.farg[3] <= FTINY) - objerror(so, USER, "illegal radius"); + objerror(so, USER, "illegal source radius"); + src->sflags |= SCIR; VCOPY(src->sloc, so->oargs.farg); src->srad = so->oargs.farg[3]; src->ss2 = PI * src->srad * src->srad; for (i = 0; i < 3; i++) src->ss[SU][i] = src->ss[SV][i] = src->ss[SW][i] = 0.0; for (i = 0; i < 3; i++) - src->ss[i][i] = .886227 * so->oargs.farg[3]; + src->ss[i][i] = 0.7236 * so->oargs.farg[3]; } -rsetsrc(src, so) /* set a ring (disk) as a source */ -register SRCREC *src; -OBJREC *so; +void +rsetsrc( /* set a ring (disk) as a source */ + SRCREC *src, + OBJREC *so +) { - register CONE *co; + CONE *co; src->sa.success = 2*AIMREQT-1; /* bitch on second failure */ src->so = so; /* get the ring */ co = getcone(so, 0); + if (CO_R1(co) <= FTINY) + objerror(so, USER, "illegal source radius"); VCOPY(src->sloc, CO_P0(co)); if (CO_R0(co) > 0.0) - objerror(so, USER, "cannot hit center"); - src->sflags |= SFLAT; + objerror(so, USER, "cannot hit source center"); + src->sflags |= (SFLAT|SCIR); VCOPY(src->snorm, co->ad); src->srad = CO_R1(co); src->ss2 = PI * src->srad * src->srad; @@ -199,22 +214,27 @@ OBJREC *so; } -cylsetsrc(src, so) /* set a cylinder as a source */ -register SRCREC *src; -OBJREC *so; +void +cylsetsrc( /* set a cylinder as a source */ + SRCREC *src, + OBJREC *so +) { - register CONE *co; - register int i; + CONE *co; + int i; src->sa.success = 4*AIMREQT-1; /* bitch on fourth failure */ src->so = so; /* get the cylinder */ co = getcone(so, 0); + if (CO_R0(co) <= FTINY) + objerror(so, USER, "illegal source radius"); if (CO_R0(co) > .2*co->al) /* heuristic constraint */ objerror(so, WARNING, "source aspect too small"); + src->sflags |= SCYL; for (i = 0; i < 3; i++) src->sloc[i] = .5 * (CO_P1(co)[i] + CO_P0(co)[i]); - src->srad = co->al; + src->srad = .5*co->al; src->ss2 = 2.*CO_R0(co)*co->al; /* set sampling vectors */ for (i = 0; i < 3; i++) @@ -227,36 +247,74 @@ OBJREC *so; fcross(src->ss[SW], src->ss[SV], co->ad); normalize(src->ss[SW]); for (i = 0; i < 3; i++) - src->ss[SW][i] *= .886227 * CO_R0(co); + src->ss[SW][i] *= .8559 * CO_R0(co); fcross(src->ss[SV], src->ss[SW], co->ad); } SPOT * -makespot(m) /* make a spotlight */ -register OBJREC *m; +makespot( /* make a spotlight */ + OBJREC *m +) { - register SPOT *ns; + SPOT *ns; + if ((ns = (SPOT *)m->os) != NULL) + return(ns); if ((ns = (SPOT *)malloc(sizeof(SPOT))) == NULL) return(NULL); + if (m->oargs.farg[3] <= FTINY) + objerror(m, USER, "zero angle"); ns->siz = 2.0*PI * (1.0 - cos(PI/180.0/2.0 * m->oargs.farg[3])); VCOPY(ns->aim, m->oargs.farg+4); if ((ns->flen = normalize(ns->aim)) == 0.0) objerror(m, USER, "zero focus vector"); + m->os = (char *)ns; return(ns); } +int +spotout( /* check if we're outside spot region */ + RAY *r, + SPOT *s +) +{ + double d; + FVECT vd; + + if (s == NULL) + return(0); + if (s->flen < -FTINY) { /* distant source */ + vd[0] = s->aim[0] - r->rorg[0]; + vd[1] = s->aim[1] - r->rorg[1]; + vd[2] = s->aim[2] - r->rorg[2]; + d = DOT(r->rdir,vd); + /* wrong side? + if (d <= FTINY) + return(1); */ + d = DOT(vd,vd) - d*d; + if (PI*d > s->siz) + return(1); /* out */ + return(0); /* OK */ + } + /* local source */ + if (s->siz < 2.0*PI * (1.0 + DOT(s->aim,r->rdir))) + return(1); /* out */ + return(0); /* OK */ +} + + double -fgetmaxdisk(ocent, op) /* get center and squared radius of face */ -FVECT ocent; -OBJREC *op; +fgetmaxdisk( /* get center and squared radius of face */ + FVECT ocent, + OBJREC *op +) { double maxrad2; double d; - register int i, j; - register FACE *f; + int i, j; + FACE *f; f = getface(op); if (f->area == 0.) @@ -281,11 +339,12 @@ OBJREC *op; double -rgetmaxdisk(ocent, op) /* get center and squared radius of ring */ -FVECT ocent; -OBJREC *op; +rgetmaxdisk( /* get center and squared radius of ring */ + FVECT ocent, + OBJREC *op +) { - register CONE *co; + CONE *co; co = getcone(op, 0); VCOPY(ocent, CO_P0(co)); @@ -294,11 +353,12 @@ OBJREC *op; double -fgetplaneq(nvec, op) /* get plane equation for face */ -FVECT nvec; -OBJREC *op; +fgetplaneq( /* get plane equation for face */ + FVECT nvec, + OBJREC *op +) { - register FACE *fo; + FACE *fo; fo = getface(op); VCOPY(nvec, fo->norm); @@ -307,11 +367,12 @@ OBJREC *op; double -rgetplaneq(nvec, op) /* get plane equation for ring */ -FVECT nvec; -OBJREC *op; +rgetplaneq( /* get plane equation for ring */ + FVECT nvec, + OBJREC *op +) { - register CONE *co; + CONE *co; co = getcone(op, 0); VCOPY(nvec, co->ad); @@ -319,9 +380,12 @@ OBJREC *op; } -commonspot(sp1, sp2, org) /* set sp1 to intersection of sp1 and sp2 */ -register SPOT *sp1, *sp2; -FVECT org; +int +commonspot( /* set sp1 to intersection of sp1 and sp2 */ + SPOT *sp1, + SPOT *sp2, + FVECT org +) { FVECT cent; double rad2, cos1, cos2; @@ -342,20 +406,20 @@ FVECT org; } -commonbeam(sp1, sp2, dir) /* set sp1 to intersection of sp1 and sp2 */ -register SPOT *sp1, *sp2; -FVECT dir; +int +commonbeam( /* set sp1 to intersection of sp1 and sp2 */ + SPOT *sp1, + SPOT *sp2, + FVECT dir +) { FVECT cent, c1, c2; double rad2, d; - register int i; /* move centers to common plane */ d = DOT(sp1->aim, dir); - for (i = 0; i < 3; i++) - c1[i] = sp1->aim[i] - d*dir[i]; + VSUM(c1, sp1->aim, dir, -d); d = DOT(sp2->aim, dir); - for (i = 0; i < 3; i++) - c2[i] = sp2->aim[i] - d*dir[i]; + VSUM(c2, sp2->aim, dir, -d); /* compute overlap */ rad2 = intercircle(cent, c1, c2, sp1->siz/PI, sp2->siz/PI); if (rad2 <= FTINY) @@ -366,9 +430,11 @@ FVECT dir; } -checkspot(sp, nrm) /* check spotlight for behind source */ -register SPOT *sp; /* spotlight */ -FVECT nrm; /* source surface normal */ +int +checkspot( /* check spotlight for behind source */ + SPOT *sp, /* spotlight */ + FVECT nrm /* source surface normal */ +) { double d, d1; @@ -382,15 +448,15 @@ FVECT nrm; /* source surface normal */ double -spotdisk(oc, op, sp, pos) /* intersect spot with object op */ -FVECT oc; -OBJREC *op; -register SPOT *sp; -FVECT pos; +spotdisk( /* intersect spot with object op */ + FVECT oc, + OBJREC *op, + SPOT *sp, + FVECT pos +) { FVECT onorm; double offs, d, dist; - register int i; offs = getplaneq(onorm, op); d = -DOT(onorm, sp->aim); @@ -399,46 +465,45 @@ FVECT pos; dist = (DOT(pos, onorm) - offs)/d; if (dist < 0.) return(0.); - for (i = 0; i < 3; i++) - oc[i] = pos[i] + dist*sp->aim[i]; + VSUM(oc, pos, sp->aim, dist); return(sp->siz*dist*dist/PI/(d*d)); } double -beamdisk(oc, op, sp, dir) /* intersect beam with object op */ -FVECT oc; -OBJREC *op; -register SPOT *sp; -FVECT dir; +beamdisk( /* intersect beam with object op */ + FVECT oc, + OBJREC *op, + SPOT *sp, + FVECT dir +) { FVECT onorm; double offs, d, dist; - register int i; offs = getplaneq(onorm, op); d = -DOT(onorm, dir); if (d >= -FTINY && d <= FTINY) return(0.); dist = (DOT(sp->aim, onorm) - offs)/d; - for (i = 0; i < 3; i++) - oc[i] = sp->aim[i] + dist*dir[i]; + VSUM(oc, sp->aim, dir, dist); return(sp->siz/PI/(d*d)); } double -intercircle(cc, c1, c2, r1s, r2s) /* intersect two circles */ -FVECT cc; /* midpoint (return value) */ -FVECT c1, c2; /* circle centers */ -double r1s, r2s; /* radii squared */ +intercircle( /* intersect two circles */ + FVECT cc, /* midpoint (return value) */ + FVECT c1, /* circle centers */ + FVECT c2, + double r1s, /* radii squared */ + double r2s +) { double a2, d2, l; FVECT disp; - register int i; - for (i = 0; i < 3; i++) - disp[i] = c2[i] - c1[i]; + VSUB(disp, c2, c1); d2 = DOT(disp,disp); /* circle within overlap? */ if (r1s < r2s) { @@ -458,93 +523,6 @@ double r1s, r2s; /* radii squared */ return(0.); /* overlap, compute center */ l = sqrt((r1s - a2)/d2); - for (i = 0; i < 3; i++) - cc[i] = c1[i] + l*disp[i]; + VSUM(cc, c1, disp, l); return(a2); -} - - -sourcehit(r) /* check to see if ray hit distant source */ -register RAY *r; -{ - int first, last; - register int i; - - if (r->rsrc >= 0) { /* check only one if aimed */ - first = last = r->rsrc; - } else { /* otherwise check all */ - first = 0; last = nsources-1; - } - for (i = first; i <= last; i++) - if ((source[i].sflags & (SDISTANT|SVIRTUAL)) == SDISTANT) - /* - * Check to see if ray is within - * solid angle of source. - */ - if (2.0*PI * (1.0 - DOT(source[i].sloc,r->rdir)) - <= source[i].ss2) { - r->ro = source[i].so; - if (!(source[i].sflags & SSKIP)) - break; - } - - if (r->ro != NULL) { - for (i = 0; i < 3; i++) - r->ron[i] = -r->rdir[i]; - r->rod = 1.0; - r->rox = NULL; - return(1); - } - return(0); -} - - -#define wrongsource(m, r) (r->rsrc>=0 && \ - source[r->rsrc].so!=r->ro && \ - (m->otype!=MAT_ILLUM || \ - objptr(source[r->rsrc].so->omod)->otype==MAT_ILLUM)) - -#define distglow(m, r) (m->otype==MAT_GLOW && \ - r->rot > m->oargs.farg[3]) - -#define badambient(m, r) ((r->crtype&(AMBIENT|SHADOW))==AMBIENT && \ - !distglow(m, r)) - -#define passillum(m, r) (m->otype==MAT_ILLUM && \ - !(r->rsrc>=0&&source[r->rsrc].so==r->ro)) - -#define srcignore(m, r) (directinvis && !(r->crtype&SHADOW) && \ - !distglow(m, r)) - - -m_light(m, r) /* ray hit a light source */ -register OBJREC *m; -register RAY *r; -{ - /* check for over-counting */ - if (wrongsource(m, r) || badambient(m, r)) - return; - /* check for passed illum */ - if (passillum(m, r)) { - if (m->oargs.nsargs < 1 || !strcmp(m->oargs.sarg[0], VOIDID)) - raytrans(r); - else - rayshade(r, modifier(m->oargs.sarg[0])); - return; - } - /* otherwise treat as source */ - /* check for behind */ - if (r->rod < 0.0) - return; - /* check for invisibility */ - if (srcignore(m, r)) - return; - /* get distribution pattern */ - raytexture(r, m->omod); - /* get source color */ - setcolor(r->rcol, m->oargs.farg[0], - m->oargs.farg[1], - m->oargs.farg[2]); - /* modify value */ - multcolor(r->rcol, r->pcol); }