ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/Development/ray/src/rt/rv3.c
(Generate patch)

Comparing ray/src/rt/rv3.c (file contents):
Revision 1.12 by greg, Fri Jan 19 00:00:31 1990 UTC vs.
Revision 1.18 by greg, Tue Mar 19 13:14:03 1991 UTC

# Line 18 | Line 18 | static char SCCSid[] = "$SunId$ LBL";
18  
19   #include  "random.h"
20  
21 + #ifndef WFLUSH
22 + #define WFLUSH          30              /* flush after this many rays */
23 + #endif
24  
25 +
26   getrect(s, r)                           /* get a box */
27   char  *s;
28   register RECT  *r;
# Line 89 | Line 93 | double  *mp;
93                  (*dev->comout)("Pick view center\n");
94                  if ((*dev->getcur)(&x, &y) == ABORT)
95                          return(-1);
96 <                viewray(thisray.rorg, thisray.rdir, &ourview,
97 <                                (x+.5)/hresolu, (y+.5)/vresolu);
96 >                if (viewray(thisray.rorg, thisray.rdir, &ourview,
97 >                                (x+.5)/hresolu, (y+.5)/vresolu) < 0) {
98 >                        error(COMMAND, "not on image");
99 >                        return(-1);
100 >                }
101                  if (!direc || ourview.type == VT_PAR) {
102                          rayorigin(&thisray, NULL, PRIMARY, 1.0);
103                          if (!localhit(&thisray, &thescene)) {
# Line 130 | Line 137 | paint(p, xmin, ymin, xmax, ymax)       /* compute and paint
137   register PNODE  *p;
138   int  xmin, ymin, xmax, ymax;
139   {
140 +        extern long  nrays;
141 +        static long  lastflush = 0;
142          static RAY  thisray;
143          double  h, v;
135 –        register int  i;
144  
145          if (xmax - xmin <= 0 || ymax - ymin <= 0) {     /* empty */
146                  p->x = xmin;
# Line 141 | Line 149 | int  xmin, ymin, xmax, ymax;
149                  return;
150          }
151                                                  /* jitter ray direction */
152 <        p->x = h = xmin + (xmax-xmin)*frandom();
153 <        h /= hresolu;
146 <        p->y = v = ymin + (ymax-ymin)*frandom();
147 <        v /= vresolu;
152 >        h = xmin + (xmax-xmin)*frandom();
153 >        v = ymin + (ymax-ymin)*frandom();
154          
155 <        viewray(thisray.rorg, thisray.rdir, &ourview, h, v);
155 >        if (viewray(thisray.rorg, thisray.rdir, &ourview,
156 >                        h/hresolu, v/vresolu) < 0) {
157 >                setcolor(thisray.rcol, 0.0, 0.0, 0.0);
158 >        } else {
159 >                rayorigin(&thisray, NULL, PRIMARY, 1.0);
160 >                rayvalue(&thisray);
161 >        }
162  
163 <        rayorigin(&thisray, NULL, PRIMARY, 1.0);
164 <        
153 <        rayvalue(&thisray);
154 <
163 >        p->x = h;
164 >        p->y = v;
165          copycolor(p->v, thisray.rcol);
156 –
166          scalecolor(p->v, exposure);
167  
168          (*dev->paintr)(greyscale?greyof(p->v):p->v, xmin, ymin, xmax, ymax);
169 +
170 +        if (dev->flush != NULL && nrays - lastflush >= WFLUSH) {
171 +                lastflush = nrays;
172 +                (*dev->flush)();
173 +        }
174   }
175  
176  
# Line 349 | Line 363 | register VIEW  *vp;
363          } else if (bcmp((char *)vp, (char *)&ourview, sizeof(VIEW))) {
364                  copystruct(&oldview, &ourview);
365                  copystruct(&ourview, vp);
366 <                newimage();             /* newimage() calls with vp=&ourview! */
366 >                newimage();
367          }
368   }
369  
# Line 389 | Line 403 | FVECT  vc;
403   }
404  
405  
406 < spinvector(vres, vorig, vnorm, theta)   /* rotate vector around normal */
407 < FVECT  vres, vorig, vnorm;
408 < double  theta;
406 > zoomview(vp, zf)                        /* zoom in our out */
407 > register VIEW  *vp;
408 > double  zf;
409   {
410 <        extern double  sin(), cos();
411 <        double  sint, cost, dotp;
412 <        FVECT  vperp;
413 <        register int  i;
414 <        
401 <        if (theta == 0.0) {
402 <                VCOPY(vres, vorig);
410 >        switch (vp->type) {
411 >        case VT_PAR:                            /* parallel view */
412 >        case VT_ANG:                            /* angular fisheye */
413 >                vp->horiz /= zf;
414 >                vp->vert /= zf;
415                  return;
416 +        case VT_PER:                            /* perspective view */
417 +                vp->horiz = atan(tan(vp->horiz*(PI/180./2.))/zf) /
418 +                                (PI/180./2.);
419 +                vp->vert = atan(tan(vp->vert*(PI/180./2.))/zf) /
420 +                                (PI/180./2.);
421 +                return;
422 +        case VT_HEM:                            /* hemispherical fisheye */
423 +                vp->horiz = sin(vp->horiz*(PI/180./2.))/zf;
424 +                if (vp->horiz >= 1.0-FTINY)
425 +                        vp->horiz = 180.;
426 +                else
427 +                        vp->horiz = asin(vp->horiz) / (PI/180./2.);
428 +                vp->vert = sin(vp->vert*(PI/180./2.))/zf;
429 +                if (vp->vert >= 1.0-FTINY)
430 +                        vp->vert = 180.;
431 +                else
432 +                        vp->vert = asin(vp->vert) / (PI/180./2.);
433 +                return;
434          }
405 –        sint = sin(theta);
406 –        cost = cos(theta);
407 –        dotp = DOT(vorig, vnorm);
408 –        fcross(vperp, vnorm, vorig);
409 –        for (i = 0; i < 3; i++)
410 –                vres[i] = vnorm[i]*dotp*(1.-cost) +
411 –                                vorig[i]*cost + vperp[i]*sint;
435   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines (old)
> Changed lines (new)